The idea that black holes possess genuine thermodynamic properties dates back to the foundational work of Jacob Bekenstein, Stephen Hawking and their collaborators, who established that black holes carry a well-defined temperature and entropy proportional to the area of their event horizons. That analogy, however, raised a deeper question: do black holes exhibit the full machinery of thermodynamics, including phase structure and critical phenomena? Research beginning in the 1970s by P.C.W. Davies and P. Hut suggested they might, and the discovery gained fresh momentum after Juan Maldacena’s 1997 AdS/CFT correspondence made asymptotically AdS black holes central to modern theoretical physics.

A pivotal conceptual advance came with the reinterpretation of the cosmological constant as a thermodynamic pressure. In this extended framework, the phase behaviour of charged and rotating AdS black holes closely mirrors the van der Waals liquid-gas transition of ordinary fluids: a small black hole phase corresponds roughly to the gas, a large black hole phase to the liquid, and a first-order transition connects them, complete with a critical point where the distinction between the phases dissolves. The new study asks what happens to the dynamics of such transitions as the critical point is approached, treating the black hole not as a static equilibrium object but as a stochastic system buffeted by thermal fluctuations.

The researchers build on a framework known as free energy landscape dynamics, which has proven powerful in physics, chemistry and biology, from protein folding to chemical reactions. In this picture, the thermodynamic states of a system are valleys on a landscape defined by its free energy, and thermal noise jiggles the system between them. Previously, Rong-Gen Li and Jin Wang and their collaborators applied this framework to the Hawking-Page transition and to the small-large black hole transition of Reissner-Nordström AdS (RN-AdS) black holes, showing that stochastic switching between phases can occur in both directions. More recently, it was shown that near the critical and spinodal points of the RN-AdS system, the relaxation dynamics slows dramatically.

The Dibrugarh team extends this analysis to rotating Kerr-AdS black holes, and in doing so makes a deliberate technical choice: rather than using the horizon radius as the fluctuating order parameter, they treat the black hole’s entropy as the dynamical variable evolving on the free energy landscape. For Kerr-AdS black holes, the generalized Gibbs free energy in the canonical ensemble is most naturally written as a function of the entropy, and the entropy uniquely labels each equilibrium macrostate. The authors also show that this choice does not affect the universal long-time behaviour, because entropy and horizon radius are related by a smooth transformation that merely relabels coordinates on the same thermodynamic manifold.

The mathematical core of the analysis is a Langevin equation: a stochastic differential equation in which the entropy evolves under a deterministic driving force generated by the slope of the free energy, a friction term describing dissipation into the thermal environment, and a Gaussian white noise term constrained by the fluctuation-dissipation relation. Far from criticality, the free energy well surrounding a stable state is approximately parabolic, and perturbations decay exponentially with a characteristic time given by the damping coefficient divided by the curvature of the free energy at equilibrium. At the critical point, however, the first three derivatives of the free energy with respect to entropy vanish, and the quadratic approximation collapses entirely.

That collapse is the origin of critical slowing down. As the landscape flattens, the restoring force that pulls fluctuations back toward equilibrium weakens and eventually disappears, so fluctuations persist for ever longer times. The researchers demonstrate this both analytically and numerically. Simulating the Langevin equation with a Heun predictor-corrector scheme, they extract the autocorrelation time and the variance of the entropy trajectories, finding that both rise sharply near the critical point and near the spinodal lines where one of the black hole phases ceases to exist. Independently, they solve the associated Fokker-Planck equation, which describes the probability distribution of the entropy, and find that its smallest nonzero eigenvalue, which sets the slowest relaxation rate, is strongly suppressed near criticality, confirming the same physics from the spectral side.

The study’s most striking result concerns universality. Fitting the numerically obtained relaxation time to a power law of the form tau proportional to the reduced distance from criticality raised to a negative exponent, the team recovers a dynamical critical exponent of approximately two-thirds along every path they examined, whether varying the temperature at fixed pressure or angular momentum, varying the pressure at fixed temperature, or varying the angular momentum at fixed temperature, and regardless of whether the critical point is approached from above or below. Analytically, this exponent follows from a mean-field argument: near a critical inflection point the order parameter scales as the cube root of the distance from criticality, so the free energy curvature scales as the two-thirds power, and the relaxation time, its inverse, diverges as the minus two-thirds power.

Remarkably, the same exponent emerges for three physically distinct black hole families: charged RN-AdS black holes, rotating Kerr-AdS black holes, and Bardeen black holes, an early example of regular black holes whose cores are nonsingular and which satisfy the weak energy condition. Despite radically different spacetime geometries and thermodynamic variables, all three systems realize identical dynamical scaling, placing them in the same mean-field dynamical universality class. The conclusion is that critical slowing down is governed not by the microscopic details of the black hole solution but by the generic structure of the free energy landscape, specifically its quartic form at a critical inflection point, echoing the logic of universality that underpins conventional critical phenomena.

The work connects black hole thermodynamics to a web of ideas spanning the Kibble-Zurek mechanism of defect formation in cosmological phase transitions, early-warning signals of critical transitions in complex systems, and the time-dependent Ginzburg-Landau theory of near-critical dynamics. The authors suggest several future directions, including extensions to higher-dimensional black holes, modified gravity theories and multicritical systems, and possible links to quasinormal modes, thermodynamic geometry, Lyapunov exponents and holographic nonequilibrium phenomena. For now, the message is conceptually simple and profound: a black hole poised at a thermodynamic critical point forgets its past ever more slowly, and the way it forgets obeys a law that charged, rotating and even regular black holes all share. In the slow drift toward a phase transition, gravity and statistical mechanics appear to speak the same universal language.

Subject of Research: Universal dynamical scaling and critical slowing down in black hole phase transitions in anti-de Sitter spacetime

Article Title: Critical slowing down of black hole phase transition and universal dynamic scaling in AdS black holes

Article References: Awal, M. B., & Phukon, P. (2026). Critical slowing down of black hole phase transition and universal dynamic scaling in AdS black holes. The European Physical Journal C, 86(9), Article 1058. https://doi.org/10.1140/epjc/s10052-026-16329-0

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16329-0

Keywords: black hole thermodynamics, phase transitions, critical slowing down, AdS black holes, free energy landscape, Langevin equation, Fokker-Planck equation, Kerr-AdS black holes, RN-AdS black holes, Bardeen black holes, universality, dynamical critical exponent

The black hole in the study is the Bañados–Teitelboim–Zanelli, or BTZ, solution, a mathematically exact black hole in a universe with three spacetime dimensions and a negative cosmological constant. Its geometry is described by an anti-de Sitter, or AdS, background, in which space curves inward at large distances. In pure Einstein gravity, the local geometry of a vacuum three-dimensional spacetime is completely fixed by the cosmological constant. There are no ordinary propagating graviton degrees of freedom: unlike four-dimensional black holes, a disturbance does not travel through empty space as a local gravitational wave. The BTZ black hole can nevertheless have a horizon and nontrivial global properties, making it a powerful laboratory for testing ideas about black holes, quantum gravity and holography.

Almeida and Lima studied what happens when Einstein’s theory is replaced by quadratic f(R) gravity, whose simplest form adds an (R^2) term to the familiar Ricci-scalar action. Here, (R) measures the curvature of spacetime and the coefficient (\alpha) determines the strength of the correction. This modification does more than slightly alter Einstein’s equations. Through a mathematical reformulation known as the scalar–tensor correspondence, the extra curvature term behaves like a massive scalar field coupled to gravity. The resulting excitation is commonly called the scalaron. In three dimensions, it is especially important because it becomes the unique local propagating gravitational degree of freedom introduced by the theory. The researchers asked whether a localized topological defect could act as a source for this otherwise hidden mode.

Their source is a Maxwell–Higgs vortex, a field configuration related to the Nielsen–Olesen vortices used in particle physics and cosmology. A vortex forms when a complex field breaks a continuous symmetry while retaining a quantized winding around a narrow core. The winding number cannot change continuously, so the defect is topologically protected. Around its center, the Higgs-like field and gauge field vary rapidly, concentrating energy and stress in a finite region. That localized energy-momentum tensor generally has a nonzero trace, represented by (T). In standard Einstein gravity, the trace of the field equations fixes the Ricci scalar algebraically. In quadratic f(R) gravity, however, the trace becomes a differential equation, turning the vortex into a source that can launch a curvature disturbance.

In the regime where the higher-curvature correction is small compared with the AdS curvature scale, the researchers linearized the theory. The trace equation then reduces to a massive Klein–Gordon equation for (R), schematically ((\Box-m^2)R=-(2\pi/\alpha)T), with an effective scalaron mass (m^2=1/(4\alpha)). This relation produces an unusual but clear physical picture: reducing (\alpha) makes the scalaron heavier, causing the curvature excitation to become increasingly localized. The approximation requires (\alpha R\ll1), ensuring that nonlinear terms such as (\alpha R^2) remain subdominant. The authors interpret this small-(\alpha) limit as a controlled effective-field-theory expansion rather than an arbitrary mathematical simplification.

To calculate the response, the team treated the BTZ geometry as a fixed background and assumed a static, circularly symmetric vortex positioned outside the event horizon. The resulting radial equation has a self-adjoint Sturm–Liouville form, a structure that makes it possible to construct an exact radial Green function. One homogeneous solution is chosen to remain regular at the horizon, while the other decays at spatial infinity. The Green function combines these two solutions across the source region, allowing the curvature profile to be written as an integral over the vortex’s energy-momentum trace. In essence, every part of the vortex contributes to the scalaron field, but the black-hole geometry determines how those contributions propagate and how they are filtered at large distances.

The most striking result is that the far-field profile forgets almost everything about the vortex’s microscopic structure. Outside the vortex core, the curvature perturbation falls as (R(r)\sim r^{-(1+\nu)}), where (\nu=\sqrt{1+m^2\ell^2}) and (\ell) is the AdS radius. The exponent is fixed by the scalaron mass and the geometry, not by the detailed shape, width or internal field structure of the defect. The core affects mainly the amplitude through an integrated effective scalar charge, (Q_{\mathrm{eff}}), obtained by weighting the source profile with the radial Green-function response. This universality was checked numerically using a resolved Nielsen–Olesen vortex rather than a simplified idealized source. The calculated curves reproduced the predicted slope on a log–log plot for different vortex parameters, while their overall strengths varied as expected with the source details.

That long-distance behavior is also where the study connects to holography. In the AdS/ CFT correspondence, a massive scalar field in the bulk is associated with an operator in the lower-dimensional boundary theory. The exponent gives the operator’s conformal dimension, (\Delta=1+\nu). Because the researchers impose the normalizable boundary condition, the non-normalizable mode—which would represent an externally imposed source at the boundary—is absent. The vortex instead creates a response in the bulk that can be interpreted as a source-free expectation value of the dual operator. The boundary theory would therefore register the defect through a one-point function whose amplitude depends on (Q_{\mathrm{eff}}), while its scaling dimension depends only on (\alpha) and (\ell). The result offers a clean example of how a localized object deep in an AdS geometry could leave a universal imprint on boundary physics.

The calculations also indicate that the new black-hole hair is stable and energetically well behaved, within the assumptions of the model. For positive (\alpha), the scalaron has positive mass squared, so it is not tachyonic in the flat-space sense. In AdS, the relevant condition is the Breitenlohner–Freedman bound, which permits certain negative mass-squared fields but requires (m^2\geq -1/\ell^2). The scalaron lies safely above that threshold. Its energy, calculated from the effective massive-scalar action, contains gradient and mass contributions and remains finite both near the horizon and at infinity. Because the asymptotic field decays rapidly, the energy integral converges strongly rather than accumulating an infrared divergence.

Perhaps most importantly, the excitation does not significantly reshape the black hole in the perturbative regime. When (\alpha\ll\ell^2), the parameter (\nu) becomes large, approximately (\ell/(2\sqrt{\alpha})). For a vortex centered at radius (r_v) outside a horizon of radius (r_h), the scalar energy is estimated to contain a suppression factor of roughly ((r_h/r_v)^{2\nu}). Since (r_h/r_v<1), this becomes exponentially small as (\nu) grows. The scalaron can therefore form a mathematically extended curvature profile while carrying far less energy than the BTZ black hole itself. As (\alpha) approaches zero, its mass diverges and the mode decouples, smoothly recovering ordinary three-dimensional Einstein gravity, where no local gravitational degree of freedom remains.

The findings do not yet describe a fully backreacting black hole with a vortex, nor do they predict an immediately observable signature in an astrophysical system. The analysis assumes a static vortex outside the horizon, a fixed BTZ background and a linearized curvature equation. The authors identify rotating BTZ black holes as a more difficult next step. Rotation introduces inner and outer horizons, frame dragging and couplings between radial and angular modes; a vortex that is static in the nonrotating case may also need angular momentum inside an ergosphere. Future work could examine scalaron quasinormal modes, solve the coupled gravitational and vortex equations beyond leading order, and determine how the holographic response changes in a rotating geometry. For now, the study provides an analytically controlled demonstration of a remarkable possibility: in a minimal universe where Einstein gravity has no local gravitational waves, a topological vortex can activate a new curvature field and give a black hole a faint, stable and highly structured form of gravitational hair.

Subject of Research: Scalaron excitations induced by Maxwell–Higgs topological vortices in quadratic f(R) gravity on a BTZ black hole background

Article Title: Scalaron hair induced by topological vortices in quadratic f(R) gravity on a BTZ black hole background

Article References: Almeida, C. A. S., & Lima, F. C. E. “Scalaron hair induced by topological vortices in quadratic f(R) gravity on a BTZ black hole background.” General Relativity and Gravitation 58, 63 (2026). Original research article

Image Credits: AI Generated

DOI: 10.1007/s10714-026-03567-6

Keywords: quadratic f(R) gravity, scalaron, BTZ black hole, topological vortices, three-dimensional gravity, curvature excitation, anti-de Sitter spacetime, black-hole hair