Papers

We investigate null geodesics and their stability in four-dimensional charged Kalb-Ramond black holes with a cosmological constant. The non-vanishing vacuum expectation value of the background antisymmetric tensor field induces spontaneous Lorentz symmetry breaking, controlled by a dimensionless parameter $l$, and deforms the Reissner-Nordström-(A)dS geometry. We derive the effective potential for massless particles and obtain analytic expressions for the photon-sphere radius, critical impact parameter, capture cross section, and deflection trajectories. The limiting cases $l\rightarrow0$ recover the RN, RN-dS and RN-AdS geodesic structures, while finite values of $l$ shift the turning points, modify the bending of light and deform both first-kind and second-kind photon trajectories. In the AdS branch, we find a special limiting limaçon-type null geodesic whose angular structure is modified by the Kalb-Ramond parameter. We also compute the Lyapunov exponent of the unstable circular null orbit and show that the cosmological constant changes the instability timescale, whereas the Lorentz-violating parameter produces a non-trivial deformation of the photon-sphere instability. Our results identify null geodesics as sensitive probes of Lorentz-violating effects in charged black hole spacetimes.

In 2025, the 4th Gravitational Wave Transient Catalog reported the discovery of over 200 gravitational wave events since 2015. These gravitational waves (GWs) are emitted by compact binary coalescenses (CBCs) between neutron stars and black holes, and detected with a global network of ground-based laser interferometers. GWs represent a new `messenger' – ripples in spacetime that propagate across cosmological distances, and carry information about the physical properties of compact objects and their location in the Universe. In this review, we consider what information we can learn about black holes and neutron stars using GWs as a messenger. We explore how the GW observation process allows us to transform the outputs of laser interferometers like LIGO, Virgo and KAGRA into knowledge about black holes and neutron stars. We show that the GW source population acts as a messenger that can lead us to breakthroughs in our understanding of stellar evolution, how the Universe has changed across cosmic time, and fundamental physics. Finally, we discuss three key opportunities to learn more about our Universe using GWs as a messenger in the coming decades.

Whether a rotating black hole (BH) immersed in an external magnetic field charges up to the Wald value $Q_{\rm W}=2aMB$, where $M$ and $a$ are the BH mass and spin parameter, and $B$ is the strength of the external field, is a long-standing open question in BH electrodynamics, with consequences for charge separation, particle acceleration and the structure of BH magnetospheres. We address it with axisymmetric general-relativistic particle-in-cell simulations performed with GRZeltron, including self-consistent pair creation. For representative values of the BH spin, we evolve the same asymptotically uniform Wald field from two opposite initial horizon charges, $Q_{0}=0$ and $Q_{0}=Q_{\rm W}$, and track the accumulated charge through the BH horizon. We find that both branches relax within a few tens of gravitational times to the same saturated charge value. Thus, the equilibrium is a dynamical attractor of the kinetic magnetosphere rather than a memory of the initial data. The attractor lies well below $Q_{\rm W}$, at $\xi_{\rm eq}\equiv Q_{\rm eq}/Q_{\rm W}\approx 0.3$ for $a\lesssim0.7$, then falling almost to $\xi_{\rm eq}\approx0$ for large spins $a\approx1$. We derive an analytic expression for $\xi_{\rm eq}$ by requiring equal magnetic fluxes through the horizon associated with positive and negative charges, and find its spin dependence to be in agreement with the simulations. The charge state of an astrophysical BH is therefore driven by kinetic plasma processes, and the Wald charge is an upper bound, not a general equilibrium value. Since $\xi_{\rm eq}<1$, charging does not quench the horizon-infinity potential drop that powers energy-extraction processes such as the Blandford-Znajek mechanism.

Null-congruences encode the causal structure of spacetime, describing a phase space of radiation channels, whose expansion scalar $\theta$ satisfies the Raychaudhuri equation. In spherical symmetry, the expansion scalar of wave fronts of area $A=4\pi r^2$ satisfies $\theta = A^\prime/A = 2/r$. About a mass $M$, this carries an encoding area $A_E = \lambda \varphi \ell_p^2 = \lambda R_g r \le A$ in a UV-IR consistent coupling $\ell_p^2= G\hbar/c^3$ to its Compton phase $\varphi = (mc/\hbar)r \propto 1/\hbar$. Here, $\lambda = 8\pi $ is fixed by the Newtonian limit at large distances and the saturation limit $A_E=A$ at horizon about $M$ based on the non-local extension of the equivalence principle in Rindler spacetime, without using the Einstein equations. In spherical symmetry, this recovers the Schwarzschild metric from a geometric embedding of $p=A_E/A$, which is directly observable in gravitational lensing. The result shows a non-thermal origin of the scaling dimension of two of gravitational phase space, previously derived based on the second law of thermodynamics by Bekenstein (1973), with possible implications for cosmological spacetime.

We employ the statefinder hierarchy along with the fractional growth parameter to investigate the recently proposed noninteracting Viaggiu holographic dark energy (VHDE) model with the future event horizon as the infrared cut-off. By adjusting different values of the VHDE model parameter $\delta$, we illustrate the evolution trajectories of several important parameters, prticularly, the statefinder hierarchy $S_3^{(1)}$, $S_3^{(2)}$, $S_4^{(1)}$, $S_4^{(2)}$ and the fractional growth parameter $\epsilon$ against the redshift $z$ and the fractional energy density of matter, $\Omega_m$. Together with the growth rate of matter perturbations, we present the statefinder hierarchy that defines a composite null diagnostic (CND), which can differentiate dynamical dark energy from $\Lambda$CDM. The joint use of the statefinder hierarchy and the fractional growth parameter offers a powerful diagnostic tool for probing VHDE, effectively resolving the degeneracy between different parameter choices within the model. We find that employing a CND pair, rather than a single diagnostic, significantly enhances the diagnostic power for the VHDE model. Finally, we perform an $\omega_d-\omega_d'$ analysis and, in addition, assess the squared sound speed $v_{s}^{2}$ to determine the model's dynamical stability under small perturbations.

Nonextensive Bose and Fermi state is well studied with counting near black-hole horizons in a one-function, nonextremal, static, spherically symmetric class. The calculation is organized around a $q$-deformed grand ensemble whose total partition function is a $q$-product of single-level partitions. This identifies the bosonic occupation-number sum as the primary level partition and allows us to separate two operations that are often conflated: replacing that sum by a geometric-series form and applying the ordinary reciprocal-log identity to a $q$-logarithm. A finite-energy tunneling prescription yields a leading near-horizon state-counting term controlled by the surface gravity and a horizon-response factor. Entropy matching then gives an energy scale proportional to the Hawking temperature; the geometry factor cancels from normalized nonextensive corrections and from the Bose–Fermi ratio. We derive the first-order coefficients analytically and evaluate the unexpanded $q$-integrals numerically. The linear approximation is accurate to better than one percent when $q$ lies within roughly two percent of unity. All prescriptions recover the extensive limit, but the exact bosonic level sum and the geometric approximations give opposite signs for the leading deformation of the conventional $7/8$ Bose–Fermi relation. The result isolates a robust geometry-factorized structure and quantifies a distinct $q$-algebra systematic uncertainty.

We revisit the Nakamura conjecture, which relates the Tomimatsu-Sato solutions of stationary axisymmetric gravity to finite Toda molecules. While the conjecture has been established partially, its general rotating sector remains an open problem. We show that the Toda determinants underlying the conjecture possess a natural weight grading. In particular, the two functions entering the Ernst potential have weights n^2 and n^2-1, and this grading extends systematically to shifted determinants labelled by partitions. In coordinates adapted to the Toda generators, each differentiation corresponds to adding one box to the associated Young diagram and increases the weight by one. The same integer n^2 also appears in the zero-order term of the Nakamura bilinear operator, revealing a compatibility between the differential equation and the determinant grading. The partition structure further explains the previously unresolved behavior of second derivatives. Repeated differentiation in one direction produces an internal sector and an external sector requiring only a one-step extension of the Wronskian hierarchy; the latter is reduced by a local three-term Pluecker relation. Thus weight grading, Young-diagram growth, Wronskian enlargement, and Pluecker reduction emerge as parts of a single determinant structure. The unit weight relation n^2 = (n^2-1) + 1 also singles out the elementary Toda seed as a natural third object, suggesting a possible route toward a genuine trilinear formulation. Although no trilinear closure is assumed here, the present construction reduces the remaining general-n Nakamura problem to definite determinant-minor identities and provides a structural framework in which such a formulation can be investigated.

In this work we establish general constitutive equations for a relativistic plasma composed of an arbitrary number of classical, charged, and chemically non-reacting species in the presence of a background electromagnetic field. The intensity of the electromagnetic force acting on either species is assumed to be comparable with the gradients of the state variables and is thus considered as a first-order driving force for dissipation within the gradient expansion. Assuming that the species share the same temperature and velocity flow when in equilibrium, we determine the entropy production and establish constitutive equations in terms of first-order frame invariant quantities. Conditions on the corresponding transport coefficients are given in order to comply with Onsager's reciprocal relations and the second law of thermodynamics. In particular, we discuss cross effects for the binary mixture which may be relevant for astrophysical and cosmological applications. Finally, we derive a complete set of evolution equations governing the dynamics of a binary mixture propagating on an arbitrary spacetime background.

Stars embedded in the disks of active galactic nuclei (AGN) can grow to hundreds of solar masses, and the same disks are expected to host a population of stellar-mass black holes. A star may therefore capture passing black holes, turning its interior into a potential factory for forming compact binary black hole systems and, ultimately, for the gravitational-wave events seen by our detectors. We test how readily this channel operates using three-dimensional hydrodynamic simulations of black hole–star encounters. As a proof of principle, we adopt a star-to-black-hole mass ratio of $34\!:\!1$, and find that the capture does not disrupt the structure of the star: it heats the star by $10$–$30\%$, and the star loses only $\lesssim1\%$ of its mass throughout. The captured black hole loses its orbital angular momentum rapidly to dynamical friction, sinking to the stellar center within a stellar dynamical time $\lesssim10^{4}\,$s. If the star already harbors a black hole at its center, the two form a bound binary whose gravitational-wave coalescence time falls below $10^{4}\,$yr. Depending on the geometry and speed of the encounter, the resulting orbit may be nearly circular or retain substantial eccentricity. Our calculations confirm that black hole–star encounters in AGN disks are indeed a viable channel for assembling stellar-mass black hole binaries and supplying sources for gravitational-wave detectors.

We study the relaxation of a complex massive scalar field around a supermassive black hole, modeling the near-horizon dynamics of Fuzzy Dark Matter (FDM) with spatial granularity. The field is initialized as a broadband, anisotropic multi-mode random field in momentum space with various spectral widths $\sigma$, decomposed into spherical harmonics up to $\ell_{\max}=100$, and numerically evolved on a finite domain using a specialized radial basis adapted to the background spatial geometry. We track the evolution of each multipole and calculate its Noether flux across both the event horizon and the outer boundary. Since angular momentum barriers suppress the horizon absorption of high-$\ell$ modes, while high-frequency radial components and fine-scale angular structures dissipate via outward radiation and horizon accretion, the Schwarzschild spacetime acts as an effective low-pass filter. Regardless of the initial spectral width, all configurations converge toward a universal late-time relaxation state dominated almost exclusively by low multipoles ($\ell \le 2$) and small radial wavenumbers, with the initial bandwidth $\sigma$ determining the timescale of global charge depletion. For concrete physical estimates, we set the boson mass to $m_b=10^{-22}\,\mathrm{eV}/c^2$ and the black hole mass to $M_{\mathrm{BH}}=6.5\times10^9 M_\odot$.

Gravitational-wave searches for subsolar-mass compact binaries have left a gap between the planetary-mass regime accessible to long-duration searches and the higher masses probed by conventional compact-binary analyses. We close this gap by searching for the inspiral of ultra-compact binaries with chirp masses $10^{-2}$–$10^{-1}\,M_\odot$ in LIGO O4ab data. Using the statistically robust BinaryGFH-v2 method, a time-frequency technique to search for rapidly evolving inspirals lasting minutes to hours, we cover this previously challenging parameter space with a computational cost of only $\sim 350$ CPU-hours, and find no significant candidates. We thus set 90% confidence upper limits on the binary merger-rate density, providing the first gravitational-wave constraints across this chirp-mass range. Assuming that primordial-black-hole binaries formed in the early Universe with no binary suppression ($f_{\rm sup}=1$), these results provide the first gravitational-wave constraints on the primordial-black-hole \dmh fraction in this mass range, yielding upper limits with $f_{\rm PBH}<1$ for equal-mass binaries with component masses $10^{-2}$–$10^{-1}\,M_\odot$ and asymmetric binaries with secondary masses $m_2=3\times10^{-4}$–$10^{-2}\,M_\odot$ for a primary mass $m_1=2.5\,M_\odot$.

We present a three-dimensional visco-magnetoacoustic framework for analog gravity in magnetohydrodynamic (MHD) flows. While standard fluid dissipation typically breaks the Lorentzian signature of acoustic metrics, we demonstrate that evaluating wave perturbations in the eikonal limit alongside the Shakura-Sunyaev $\alpha$-viscosity prescription preserves a well-defined effective spacetime geometry for the fast magnetoacoustic mode. The slow-magnetoacoustic mode and Alfvén mode do not admit a non-degenerate metric. To investigate analog horizon thermodynamics, we utilize astrophysical accretion disks as background media, specifically modeling numerical magnetized advective accretion flows around rotating black holes and analytical advection-dominated inflow-outflow solutions (ADIOS). Standard self-similar ADIOS models strictly enforce a constant Mach number, precluding horizon formation. We therefore introduce a magnetic field perturbation that breaks self-similarity, generates a dynamic Mach number, and enables the formation of a visco-magnetoacoustic horizon. By evaluating the spontaneous phonon emission at these horizons, we reveal that the analog Hawking temperature is highly sensitive to the magnetic field topology; the spatial orientation of the background magnetic gradients dictates whether the Hawking radiation is amplified or suppressed. Furthermore, we find that increasing the viscosity parameter leads to a monotonic increase in the analog Hawking temperature.

Kehle and Unger have recently shown that, under suitable conditions, a massless charged classical scalar field can produce an extremal black hole in finite time. This result appears to violate the third law of thermodynamics (TLT), since extremal black holes have zero temperature. However, because the TLT is intrinsically a quantum principle, a self-consistent analysis of this problem must account for quantum effects. In particular, we point out that the electric field of the system may spontaneously polarize the vacuum, thereby creating massless charged pairs that will discharge the black hole. This mechanism guarantees the quantum validity of the TLT in the regime $|eQ|\geq\hbar/2$ (here ${Q,e}$ denote, respectively, the electric charges of the black hole and the scalar field). This leaves open the question of what happens for lower values of $|eQ|/\hbar$. Recently, Schneider derived a classical bound showing that extremal black holes cannot be formed à la Kehle-Unger in the small-charge regime $|eQ|\leq\hbar/3$. His analysis leaves open the possibility that a stronger bound may exist. We conjecture here that this bound can indeed be strengthened such that classical considerations would prevent extremal (zero-temperature) charged Reissner-Nordström black holes from forming dynamically in the regime $|eQ|<\hbar/2$, where quantum polarization effects are too weak to preserve the validity of the TLT.

A sufficiently strong hadron-to-quark first-order phase transition can give rise to a third family of stable compact stars that are commonly referred to as hybrid hadron-quark stars. Stable hybrid stars that have the same gravitational mass as neutron stars are referred to as twin stars. Although equilibrium twin stars may exist, whether they can be dynamically formed remains an open question. We investigate this problem by examining the gravitational binding energy of competing equilibrium configurations at fixed baryonic rest mass and by performing general relativistic hydrodynamical simulations of several possible transition channels. While twin stars are more gravitationally bound than neutron stars with the same rest mass, this energetic preference alone does not determine the dynamical outcome. Compression and shocks during the evolution generate thermal pressure that can prevent the system from settling on the cold twin star branch. Our simulations show that sufficiently rapid cooling can remove this thermal support and enable twin star formation, whereas slower or no cooling generally favors a neutron star remnant. Accessing the twin star branch through the formation channels considered here requires cooling on a timescale comparable to or shorter than the stellar dynamical timescale. Since realistic cooling mechanisms operate on much longer timescales, our results suggest that in channels that conserve the total rest-mass neutron stars may be dynamically favored even when a more gravitationally bound twin star configuration exists with the same rest mass. Our results demonstrate a point of principle, at least for equations of state where the quark deconfinement density does not change appreciably for temperatures up to $10-20\,\textrm{MeV}.$

Gravitational sums over spacetime topology suggest a baby-universe field theory whose quanta are entire universes. We develop such a description from group-averaged amplitudes that impose gravitational constraints and define physical inner products in the single-universe theory. To sharpen the contrast with ordinary QFT, we start with gravity on one-dimensional worldlines, which can arise from reducing higher-dimensional cylinders. For multiple noninteracting worldlines that pair boundary components, the baby-universe Hilbert space is the symmetric Fock space over the physical one-universe Hilbert space, and the universe field expands in the on-shell modes spanning the latter. Each mode multiplies the sum of creation and annihilation operators, so the mode coefficients commute. The mode expansion thus makes explicit the field's commutativity and its Wheeler-DeWitt linear field equation. Diagonalizing the coefficients recovers a statistical description by a Gaussian path integral over on-shell field configurations labelling superselection ($\alpha$-)sectors. We take exploratory steps toward more general topologies. Disk-like topologies coherently shift the no-boundary state, turning on the Hartle-Hawking wavefunction as a one-point function. Even with nontrivial topologies represented by interacting Feynman diagrams, we argue that group averaging near each boundary still enforces the same linear field equation, while interactions instead enter through a non-Gaussian measure over $\alpha$-sectors. In the topological Marolf-Maxfield model, we promote a disk-and-cylinder truncation to an exact Fock-space representation of the full theory incorporating all topologies. We also relate the truncated and exact inner products by an insertion of interactions, whose perturbative expansion produces a sum over time evolutions generated by pairs of pants, evoking an analogy to group averaging.

We find the necessary and sufficient conditions for the gravitational path integral (GPI) to admit an interpretation in terms of (i) an inner product in the Hilbert space of open/closed universes or (ii) a statistical average of boundary observables. The conditions needed for (i) are necessary but not sufficient for (ii). Requiring the stronger condition (ii) places additional constraints on all wormhole amplitudes, giving a concrete diagnostic of whether the GPI computes an average or merely a pseudo-average. To exemplify the difference between the two constraints, we present a class of gravitational toy models that satisfy (i) but fail the stronger condition (ii); to emphasize the power of the stronger condition, we initiate a bootstrap study of statistical ensembles based on condition (ii), finding that the resulting constraints substantially sharpen the bounds on the moments of the ensemble obtained by imposing (i). More broadly, both of these positivity conditions impose nontrivial constraints on wormhole amplitudes and thereby restrict which gravitational effective field theories can be embedded in consistent theories of quantum gravity.

We confront a class of dark-energy equations of state emerging from group field theory (GFT) quantum gravity with DESI Data Release 2 baryon acoustic oscillations and Pantheon+ type-Ia supernovae. We introduce sampling parametrisations that replace microscopic initial-condition parameters by combinations more directly measured by background probes. The GFT solutions separate into logarithmic, power-law and oscillatory branches, determined by the microscopic interaction parameter $m$. The logarithmic branch is constrained to lie extremely close to a cosmological constant, while the power-law branch permits a small phantom deviation. Without perturbative-theory priors, oscillatory solutions can reproduce the mild preference of the distance data for a dip in $w(z)$ near $z\simeq0.5$–$1$. Profile-likelihood constraints favour $m\sim-2$ and $m\sim-5$ from BAO and supernovae, shifting towards $m\sim-3.5$ when CMB information is included. Conservative perturbative priors strongly suppress these deviations from $\Lambda$CDM. The quantum-gravity scale $z_q$, related to the average number of quantum gravity atoms, remains unconstrained, although its role in the time evolution makes higher-redshift observations a promising route to probing it. We further find that strong projection effects highlight the importance of performing likelihood profiling alongside our marginal posterior constraints. Our results provide a first direct test of GFT-motivated dynamical dark energy and demonstrate the potential for cosmological observations to inform quantum-gravity model building.

Next-generation gravitational-wave observatories like Cosmic Explorer will provide an unprecedented opportunity to probe the cosmological stochastic gravitational-wave background, providing insight into the early Universe that is inaccessible with other means. We investigate the prospects for measuring a Gaussian cosmological gravitational-wave background in Cosmic Explorer in the presence of a distinguishable, non-Gaussian foreground of binary black hole mergers, using a Bayesian search that simultaneously fits binary black hole and background parameters. We successfully recover a Gaussian cosmological background with a signal-to-noise ratio of 5 in one day of simulated data in Cosmic Explorer without contamination from unresolved signals or data exclusion. However, when introducing waveform systematic errors to our analysis, our inference on the stochastic background parameters becomes biased and can falsely infer the presence of a background when there is none. The bias is driven by signal mismatches at low frequencies (<20 Hz), especially when the recovery waveform is missing physical effects included in the injected signals. While we propose a metric for diagnosing if measurements will be strongly biased by waveform systematics, our results suggest improvements in waveform modeling or other mitigation techniques will be required to accurately recover a primordial stochastic gravitational-wave background in next-generation detectors.

We identify the spatial projection of a conserved string charge density with a bit-thread flow. This correspondence gives bit threads a string worldsheet description and allows us to compute black hole entropy from the maximal entropy flux of a worldsheet-derived current. We first reproduce the BTZ entropy and then extend the construction to the ten-dimensional F1-NS5-P and D1-D5-P systems. Under Type IIB S-duality, the F1 current maps to the D1 current. Although the local currents, metrics, and flow norms change, the maximal physical flux remains unchanged and gives the same entropy. This provides a nontrivial check of the correspondence. Finally, we propose a D1-D5 boundary-matching conjecture. In the coarse-grained picture, each string-charge line defines an entropy-carrying tube. The boundary CFT entanglement entropy and the bulk black hole entropy both fix the product of the number of tubes and the maximum entropy carried by each tube. For a fixed tube discretization, they assign the same capacity to each tube. They therefore describe the same conserved maximal entropy flux, providing an information-flow interpretation of the equality between entanglement entropy and black hole entropy.

Whether classical spacetime remains welldefined at arbitrarily high curvature is a central question. In general relativity, singularities signal the breakdown of the classical description, motivating modifications of spacetime's short-distance structure. Pseudocomplex general relativity (pcGR) extends spacetime geometry to pseudo-complex coordinates, naturally introducing two metric sectors, the physical metric g{\mu}{\nu} and an auxiliary field f{\mu}{\nu}. The magnitude of the pseudo imaginary component defines the invariant acceleration scale a0, a direct consequence of the pseudo-complex structure. The simultaneity condition, requiring both idempotent sectors to satisfy the pseudo-complex Einstein equations independently, fixes the auxiliary field algebraically, with no new propagating degrees of freedom, and sets a lower bound on the lapse that regularizes curvature. The regularization is achieved through the combined effect of the lapse-gap condition e{\nu}(r) gt. a0, which ensures the radial metric component remains nondegenerate, together with the regularity conditions at the areal-radius origin, B(0) = 1 and B'(0) = 0, which emerge naturally from the pseudo-complex geometry. These conditions jointly eliminate the Schwarzschildtype curvature divergence, rather than the gap condition acting alone.

We study a curvature-coupled dark energy model that can modify cosmological evolution both before recombination and during the late-time accelerated era. The model belongs to the class of scalar-tensor theories, in which a quintessence field is non-minimally coupled to the Ricci scalar. We specify the model through a shifted quartic coupling, $f(\varphi)=\alpha(\varphi^2-\varphi_{\rm today}^2)^2$, and an inverse power-law potential, $V(\varphi)=\Lambda\varphi^{-\sigma}$, where $\varphi_{\rm today}$ is a constant fixed by requiring the effective Planck mass to recover its present-day normalization. We implement the model in a modified version of $\mathtt{hi\_class}$ and compute its background and linear cosmological predictions. This specific form of non-minimal coupling allows an effective crossing of the phantom divide at late times while naturally suppressing deviations from standard gravity today and satisfying local gravity constraints. At the same time, the scalar field can modify the expansion history before recombination, shifting the acoustic scale in the direction required to alleviate the $H_0$ tension, while remaining dynamically relevant as dark energy at low redshift. For the parameter choices studied here, we find tens-of-percent deviations from $\Lambda$CDM in the expansion history, matter clustering, and metric-potential spectra. The modified evolution of the gravitational potentials leaves characteristic signatures in relativistic observables, with weak-lensing power suppressed by $\mathcal{O}(20$-$40\%)$ at low multipoles and order-unity changes in the late Integrated Sachs-Wolfe signal. Together, these results reveal a broad, scale-dependent phenomenology that motivates both a full parameter-space analysis and an extension of the present linear treatment to a dedicated non-linear $N$-body implementation.

The interpretation of gravitational wave sources depends on the specific choice of model used to interpret signals. Previous analyses of GWTC-3 event candidates using SEOBNRv4PHM, IMRPhenomXPHM, and NRSur7dq4 have sometimes arrived at notably different conclusions about event properties. Subtle analysis settings and code differences can also produce notable differences. We revisit 60 publicly available events from the first three observing runs, released across GWTC-1, GWTC-2, GWTC-2.1, and GWTC-3, using a consistent analysis framework and three waveform models: two state-of-the-art time-domain models with higher-order modes (SEOBNRv5PHM and IMRPhenomTPHM) and one older frequency-domain model without higher-order modes (IMRPhenomPv2). The two state-of-the-art models agree well overall, although about 20% of events have a visible difference in at least one one-dimensional marginal posterior. By contrast, the older model often arrives at qualitatively different conclusions for key events across the mass spectrum. We also identify a bug in published LVK GWTC-2.1 SEOBNRv4PHM results and in a later reprocessing of GW200105. At the adopted sampling rates, the model's Nyquist-frequency restriction excluded part of the allowed parameter space and artificially truncated the mass-ratio posterior away from equal mass for seven low-mass events. These results emphasize the need to use multiple state-of-the-art waveform models and carefully validated analysis settings to characterize uncertainty in gravitational-wave source properties.

We explore the thermodynamic properties of static black holes in conformal gravities in arbitrary, even dimensions. In the spherically symmetric case, up to local Weyl rescalings, these spacetimes are static and unique, due to a Birkhoff's theorem that applies for all conformal gravities whose Lagrangians are complete contractions of $k$ Weyl tensors in dimension $D=2k$. The solutions can be chosen to be asymptotically locally (A)dS or flat, and the horizon can be extended to be a surface of arbitrary constant curvature, $\gamma$. We focus on the asymptotically AdS case, for which the cosmological constant at infinity is an integration constant, and compute the mass of the spacetime in terms of the diverse hair parameters which have a purely gravitational origin. We show that the first law of thermodynamics is satisfied, and study phase transitions in the canonical ensemble, showing the existence of zeroth-order black-hole/black-hole phase transitions for arbitrary dimension.

In this paper, we derive the complete transformations of a generic first order perturbative metric in Schwarzschild spacetime to the Newman-Unti (NU) gauge in series expansion near null infinity. This allows us to determine the asymptotic shear, the mass aspect and the angular momentum aspect of the perturbative fields. As a direct application, we derive the corresponding asymptotic NU data for quasinormal modes (QNMs) of a Schwarzschild black hole. The total energy and angular momentum of QNMs obtained from the asymptotic NU data vanish but classical supertranslation charges are non-trivial, which can be applied to fix the BMS frame at first perturbative order.

In this paper we investigate accelerated gluodynamics for a broad intervals of temperature and acceleration. Our study is carried out within lattice simulation in the co-moving reference frame parameterized by the Rindler coordinates. We developed the renormalization prescription that allowed us to calculate renormalized Polyakov loop as a function of coordinate in the Rindler spacetime. Using the data for the renormalized local Polyakov loop, we calculated spatial dependence of the static quark free energy and effective mass of static quark. Besides the Rindler coordinates, it is believed that accelerated gluodynamics can be approximated utilizing non-accelerated gluodynamics with a properly adjusted temperature gradient in accordance with the Tolman-Ehrenfest law. We compared these approaches for the observables under study. It was found that they agree quite well close to the critical temperature and demonstrate disagreement at higher temperatures. We believe that this disagreement might be attributed to the Tolman-Ehrenfest law corrections which appear in the Rindler gluodynamics.

We study slowly rotating, nonlinearly charged black holes in arbitrary spacetime dimension $d \geq 4$, using a generalized Lense–Thirring ansatz. Working to linear order in the rotation parameters, we solve the Einstein–nonlinear electrodynamics (NLE) field equations for restricted theories ${\cal L}({\cal S})$ depending on single invariant. We provide an extension of the four-dimensional "No Go theorem" proven in ArXiv:2203.01919 and demonstrate the uniqueness of Maxwell theory in $d \geq 4$. Similar to four dimensions, it is proven that the standard rotation generating trick fails to produce full rotating solutions for any nontrivial NLE in higher dimensions. We identify a family of NLE theories admitting slowly rotating solutions with vector potentials encoding rotation effect in Maxwell-like manner; in four dimensions this reproduces the previously found RegMax theory, while in $d\geq4$ we obtain the corresponding electric potential in closed form and the associated Lagrangian implicitly as a function of the radial coordinate.

The magnetic reconnection process near a black hole can efficiently convert magnetic energy into the kinetic and thermal energy of plasma outflows, enabling the extraction of rotational energy from a spinning black hole. However, its impact on the surrounding accretion disk remains insufficiently explored. In this work, we develop a framework to investigate the 'heating' effect of magnetic reconnection outflow plasma on thin accretion disks around Kerr black holes. We incorporate the energy and angular momentum fluxes of reconnection outflow into the disk conservation laws as additional source terms, yielding a modified radiative flux. We investigate the dependence of the modified radiative flux on the black hole spin, reconnection radius, plasma magnetization, and outflow orientation. For a fixed magnetic field prescription, the contribution of the reconnection outflow plasma to the modified radiative flux decreases as the black hole spin increases. Its dependence on the reconnection radius differs inside and outside the ergosphere: inside, a smaller radius enhances disk heating through more efficient extraction of black hole rotational energy, whereas outside, a larger radius produces stronger heating due to reduced radial dilution of the outflow. The reconnection outflow plasma contribution also increases with plasma magnetization and the radial component of the outflow velocity. Our results demonstrate that magnetic reconnection can provide an additional source of energy and angular momentum for relativistic accretion disks, establishing a connection between plasma processes near rotating black holes and the observable accretion disk properties, providing a new perspective on the strong gravitational field regime.

We derive an extension of Chandrasekhar's pulsational stability criterion for relativistic stars with a pressure boundary condition at the surface. We then apply the new criterion to hylotropes, a model describing rapidly accreting supermassive stars with a core-envelope structure. We find that confinement typically stabilises stars, delaying the onset of the general relativistic instability when the external pressure at the surface is roughly $\sim10^{-3}$ of the stars' hydrostatic pressure scale. Beyond this value, the critical mass at marginal stability increases rapidly, exceeding the value of $\sim 10^5$ M$_{\odot}$ for helium burning, isolated $n=3$ polytropes by orders of magnitude. We discuss how the required pressure boundary conditions may be produced by the combined ram and magnetic pressure supplied by the accretion flow as it assembles the supermassive star, and provide a fitting formula for the final mass before collapse. As a rough estimate assuming equipartition between magnetic and thermal pressures in a standard $\alpha$-disk, we find that the critical mass doubles at accretion rates of order $\dot M\sim100\,{\rm M_\odot}\,{\rm yr}^{-1}$, again increasing rapidly after this value. Our results show under what conditions, potentially realized in a subset of high redshift halos, non-rotating SMS may collapse into black holes beyond the typical scale of heavy seeds.

We investigate the loss of coherence of a two-level Unruh–DeWitt detector coupled to massless real scalar fields in $(1+3)$-dimensional de Sitter spacetime. Treating the detector–field interaction perturbatively to second order, we derive analytic expressions for the real asymptotic coherence-decay coefficients for conformally and minimally coupled scalar fields. For the conformally coupled field, the coefficient remains finite in the small-gap limit and grows linearly with the detector energy gap in the large-gap regime. For the massless minimally coupled field, the logarithmic sector of the Wightman function produces an additional positive contribution. This provides an explicit analytic extension of the previously known infrared enhancement of detector transition responses to the asymptotic coefficient governing detector coherence. The explicitly average-time-dependent sector does not contribute to the real nonzero-frequency asymptotic coefficient under the adiabatically regulated half-line prescription adopted here, when the regulator is removed at fixed average time and fixed positive detector gap. This asymptotic result does not describe the complete finite-time dynamics, which depends on the switching function and observation interval. Within the stated prescription, the coefficients satisfy $\Gamma_{\rm MMC}>\Gamma_{\rm CC}$ for $\omega>0$, with $\Gamma_{\rm MMC}/\Gamma_{\rm CC} =1+(H/\omega)^2$. The distinction is most pronounced in the small-gap regime, whereas the two coefficients approach one another for large detector gaps.

We explore the cosmological dynamics and formation of structures within the Tsallis-Cirto modified gravity framework with a model parameter $\beta$ that is inspired by nonextensive statistics. By applying the gravity-thermodynamics conjecture, we modify the Friedmann equations and the evolution of the Hubble parameter. The impact of the Tsallis-Cirto entropy on different cosmographic parameters and the development of the linear matter overdensities have been analyzed by formulating perturbed field equations using the spherical collapse approach in a flat Friedmann-Lemaître-Robertson-Walker background. We present a novel and established diagnostic approach to distinguish among distinct cosmological models compared to flat and non-flat $\Lambda$CDM scenarios, discovering that the Tsallis-Cirto entropy-inspired modified cosmology ($\beta\ne 1$) successfully withstands all tests, contradicting both the flat and non-flat $\Lambda$CDM models. This model also meets the conditions for the Universe to reach thermodynamic equilibrium in the far future. Additionally, we investigate the halo mass function and cluster number counts within this modified gravity framework. All findings are compared to the fiducial $\Lambda$CDM profile, indicating that the added entropic correction affects the history of expansion, the growth rate of structures, and the dark matter halo abundance. We find that the more massive structures are less abundant and develop at later times, which aligns with the hierarchical model of formation of large-scale structures.