Limiting Curvature Condition in Gravity

• Limiting curvature condition is a theoretical framework that enforces a universal bound on curvature invariants, keeping them below a critical scale set by fundamental lengths.
• It is implemented by modifying the gravitational action with auxiliary fields or control functions that preserve Einstein gravity in low-curvature regimes and switch dynamics in high-curvature phases.
• Applications include regular black-hole interiors and cosmological bounces, where bounded invariants replace singularities with smooth or de Sitter-like cores.

Searching arXiv for papers on limiting curvature and related gravity models. arXiv search query: "limiting curvature gravity black holes cosmology mimetic cuscuton" The limiting curvature condition is a family of ideas in gravitational theory according to which curvature growth is prevented from exceeding prescribed critical values set by a fundamental scale. In one major formulation, the condition is implemented dynamically by modifying the gravitational action so that general relativity is recovered exactly in a low-curvature branch, while saturation of a chosen curvature bound activates auxiliary “control functions” and forces evolution onto a constrained high-curvature branch. In another formulation, especially in phenomenological regular-black-hole model building, it is a metric-level requirement that relevant curvature invariants remain everywhere finite and, in some constructions, remain sub-Planckian for every black-hole mass (Frolov, 2021, Aşmanoğlu et al., 24 Jul 2025).

1. Definitions and principal formulations

A common starting point is the Markov-type statement that curvature should be universally bounded by a scale set by a fundamental length $\ell$, schematically

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$

with $B_{\cal R}$ depending on the chosen invariant but not on solution parameters such as black-hole mass or the entropy of the universe (Frolov, 2021). In this sense, the limiting curvature condition is stronger than simple regularity at one point of a metric: it demands a parameter-independent curvature cap set by the theory.

The literature does not use a single universal implementation. Instead, several technically distinct realizations coexist.

Formulation	Limited quantity	Mechanism
Constrained modified gravity	Chosen curvature functions $\Phi_i$	Auxiliary fields $\chi_i,\zeta_i$ enforce $\Phi_i\le 0$
Mimetic gravity	$\chi=\Box\phi$; in FLRW, $\chi=3H$	Multi-valued potential $f(\chi)$ reconstructs bounded cosmological evolution
Cuscuton gravity	$K=\nabla_\mu u^\mu$; in FLRW, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$0	Nondynamical scalar imposes a constraint equation for the mean curvature
Phenomenological regular metrics	$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$1, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$2, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$3, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$4 or analogous invariants	Metric ansatz chosen so curvature maxima scale as $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$5 or $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$6

These formulations are related by purpose rather than by identical field content. A recurrent misconception is that the limiting curvature condition is merely the statement “replace the singularity by a smooth core.” The phenomenological black-hole literature explicitly distinguishes regularity at the center from the stronger requirement that the relevant curvature scalars remain uniformly bounded for every black-hole mass, and in some models remain sub-Planckian everywhere (Aşmanoğlu et al., 24 Jul 2025). A second misconception is that bounding one convenient scalar, such as $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$7, is automatically sufficient. In the cosmological models of Frolov and Zelnikov, bounding only the Ricci scalar is explicitly stated to be insufficient, because in a contracting universe different curvature components can compensate in $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$8 while the Kretschmann invariant still diverges (Frolov et al., 2021).

2. Constraint-based gravity and branch structure

The most explicit action-level realization in four-dimensional gravity is based on the Einstein–Hilbert action supplemented by a constraint sector,

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$9

with

$B_{\cal R}$0

Here the $B_{\cal R}$1 are scalar functions of curvature and limiting-curvature parameters $B_{\cal R}$2, while $B_{\cal R}$3 and $B_{\cal R}$4 are auxiliary or Lagrange-multiplier fields. Variation with respect to $B_{\cal R}$5 and $B_{\cal R}$6 yields

$B_{\cal R}$7

Since $B_{\cal R}$8, the first equation is solvable only if

$B_{\cal R}$9

This is the mechanism that enforces the inequality dynamically (Frolov, 2022).

The branch structure is central. In the subcritical regime, $\Phi_i$0, one has $\Phi_i$1, hence $\Phi_i$2, so the extra sector does not affect the metric equations and the theory is exactly on the Einstein branch. When the critical value is reached, $\Phi_i$3, then $\Phi_i$4 and $\Phi_i$5 may become nonzero. The solution then enters a supercritical branch, in which the metric equations are modified so that the saturated constraint remains dynamically consistent. In this language, $\Phi_i$6 acts as a control function (Frolov, 2022).

This structure was also used in closed FRW cosmology, where the reduced constraint action is written with control fields and auxiliary fields,

$\Phi_i$7

and the low-curvature Einstein regime is recovered when the control functions vanish (Frolov et al., 2021). In the two-dimensional dilaton-gravity construction, the same logic appears in

$\Phi_i$8

which gives

$\Phi_i$9

and therefore

$\chi_i,\zeta_i$0

inside the model (Frolov et al., 2021).

A major conceptual feature of these constructions is that the limiting curvature condition is not imposed as an external truncation of allowed solutions. It is encoded in the variational principle. This separates the low-curvature Einstein branch from the saturated high-curvature branch and makes the transition itself part of the dynamics.

3. Black-hole interiors and singularity replacement

In the four-dimensional spherically symmetric black-hole model of Frolov and Zelnikov, the interior metric is written as

$\chi_i,\zeta_i$1

For this symmetry class, the orthonormal-frame Riemann tensor has four independent eigenvalues,

$\chi_i,\zeta_i$2

$\chi_i,\zeta_i$3

with

$\chi_i,\zeta_i$4

$\chi_i,\zeta_i$5

The reason these variables are decisive is that every scalar polynomial curvature invariant is a polynomial in $\chi_i,\zeta_i$6 for the ansatz under study, so bounding them bounds all scalar polynomial invariants (Frolov, 2022).

The specific limiting-curvature conditions chosen in the black-hole model are

$\chi_i,\zeta_i$7

In Schwarzschild interior one has

$\chi_i,\zeta_i$8

so singularity formation corresponds to unbounded growth of all four eigenvalues. The first constraint directly controls the curvature of the two-dimensional $\chi_i,\zeta_i$9 sector, because the corresponding two-dimensional scalar curvature is $\Phi_i\le 0$0. The second controls the dangerous growth in the $\Phi_i\le 0$1 sector, where $\Phi_i\le 0$2 and $\Phi_i\le 0$3 govern the intrinsic curvature of the $\Phi_i\le 0$4-dimensional closed homogeneous isotropic slice (Frolov, 2022).

The 2021 analysis shows that a regular first transition into the supercritical regime requires the first active constraint to involve $\Phi_i\le 0$5, and that the primary condition $\Phi_i\le 0$6 by itself is not sufficient: it caps the curvature of the $\Phi_i\le 0$7 sector, but $\Phi_i\le 0$8 still runs to zero in finite proper time and $\Phi_i\le 0$9 diverge. This motivates the second constraint $\chi=\Box\phi$0 with $\chi=\Box\phi$1 (Frolov et al., 2021). When both constraints are active, the solution no longer approaches $\chi=\Box\phi$2. Instead, in synchronous gauge,

$\chi=\Box\phi$3

the $\chi=\Box\phi$4 factor obeys a de Sitter-type equation and takes the explicit form

$\chi=\Box\phi$5

with

$\chi=\Box\phi$6

Thus the $\chi=\Box\phi$7 sector becomes de Sitter-like, while the areal radius $\chi=\Box\phi$8 decreases only to a minimum of order $\chi=\Box\phi$9, then increases, reaches a maximum which for $\chi=3H$0 is slightly larger than $\chi=3H$1, and decreases again. The angular two-sphere oscillates rather than collapsing to zero size. The resulting spacetime is described as exponentially expanding in one direction and oscillating in the other two, with all polynomial curvature invariants uniformly bounded; within the symmetry-reduced argument, the spacetime is complete (Frolov, 2022).

A two-dimensional analogue appears in modified dilaton gravity. There the classical black-hole solution

$\chi=3H$2

has curvature

$\chi=3H$3

which diverges in the interior. The limiting-curvature modification enforces

$\chi=3H$4

and once the bound is saturated the interior enters a constant-curvature branch

$\chi=3H$5

The singularity is then replaced by a de Sitter-like core. The same paper studies Vaidya-type collapse and evaporation in this two-dimensional model and uses the bounded-curvature core to analyze late-time causal structure (Frolov et al., 2021).

4. Cosmological bounces and bounded expansion

In closed homogeneous isotropic cosmology, Frolov and Zelnikov write the metric as

$\chi=3H$6

with

$\chi=3H$7

For radiation, Einstein evolution gives

$\chi=3H$8

and the Kretschmann invariant

$\chi=3H$9

diverges in the collapse. The regular linear family is

$f(\chi)$0

Once the critical point is reached, the trajectory leaves the Einstein line $f(\chi)$1 and evolves on the constraint curve. The high-curvature regime becomes “de Sitter-type,” because the trajectory approaches the region near $f(\chi)$2, which is the de Sitter condition in these variables. The contracting universe therefore reaches a finite minimum scale factor and bounces into an expanding phase. The general discussion extends this result from the linear family to more general constraint curves $f(\chi)$3 provided the curve intersects both $f(\chi)$4 and $f(\chi)$5 and satisfies

$f(\chi)$6

along the relevant segment (Frolov et al., 2021).

A different cosmological realization appears in limiting-curvature mimetic gravity. The covariant action is

$f(\chi)$7

with

$f(\chi)$8

The mimetic constraint

$f(\chi)$9

defines a preferred timelike congruence, and in the $K=\nabla_\mu u^\mu$0-time gauge one has

$K=\nabla_\mu u^\mu$1

For spatially flat FLRW,

$K=\nabla_\mu u^\mu$2

The limiting quantity is therefore the expansion scalar rather than a direct polynomial in the Riemann tensor. In the group-field-theory reconstruction, the cosmological dynamics yields a bounded function $K=\nabla_\mu u^\mu$3 with extrema

$K=\nabla_\mu u^\mu$4

where

$K=\nabla_\mu u^\mu$5

The bounce occurs at

$K=\nabla_\mu u^\mu$6

and in the loop-quantum-cosmology limit $K=\nabla_\mu u^\mu$7 the effective Friedmann equation becomes

$K=\nabla_\mu u^\mu$8

The reconstructed $K=\nabla_\mu u^\mu$9 is necessarily multi-valued, and the paper emphasizes that continuity of $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$00 and $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$01 at branch points is required for an unambiguous smooth bounce (Cesare, 2018).

5. Stability, perturbations, and non-equivalence results

The limiting curvature condition at the background level does not by itself guarantee a viable theory of perturbations. This point is made explicitly in the perturbative study of nonsingular cosmologies based on Lagrange multipliers and curvature-invariant functions,

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$02

In the simplest model,

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$03

on FLRW, but the covariant completion built from $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$04, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$05, and $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$06 yields generic Ostrogradski instabilities unless the Weyl-squared term is included with the specific coefficient $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$07. Even with that tuning, the theory is equivalent to an $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$08 model and inherits instability problems in anisotropic backgrounds. The same paper therefore constructs a different invariant built from $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$09 and $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$10 that preserves the same background implementation of limiting curvature while avoiding the Ostrogradski instability and the additional anisotropic degree of freedom (Yoshida et al., 2017).

A related negative result concerns the relation between limiting-curvature mimetic gravity and polymer quantum gravity. In Bianchi I, the action

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$11

does lead to a bounded scalar

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$12

and in the isotropic sector the Hamiltonian reduces exactly to the standard effective LQC Hamiltonian with

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$13

However, after transforming to canonical LQC variables in Bianchi I, the Hamiltonian contains unavoidable terms polynomial in the connection variables $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$14, whereas the polymer Hilbert space supports only almost periodic, holonomy-type functions. The paper’s conclusion is therefore that limiting curvature mimetic gravity is not, in general, an effective polymer quantum gravity or effective LQC theory beyond the isotropic sector (Bodendorfer et al., 2018).

Cuscuton gravity provides a contrasting example in which the limiting-curvature idea is realized without higher-curvature terms. With

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$15

the scalar equation becomes

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$16

Defining the unit normal

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$17

one obtains a constraint on the mean curvature,

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$18

In FLRW this reduces to

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$19

The nonsingular bounce solutions of the model are classically stable at linear order. The paper shows that the spatially-flat gauge is ill-defined when $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$20, but this is a gauge artifact rather than a physical singularity; in appropriate gauges the stability results remain robust, the ultraviolet sound speed is generally close to unity, and curvature perturbations remain essentially constant in the infrared throughout the bounce phase (Quintin et al., 2019).

6. Phenomenological black-hole metrics, astrophysical applications, and open issues

In the recent phenomenological regular-black-hole literature, the limiting curvature condition is often adopted as a selection principle for static spherically symmetric metrics. One 2025 study uses the Hayward-type ansatz

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$21

with generalized mass-dependent regulator

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$22

and benchmark choice

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$23

The authors track

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$24

and interpret the limiting curvature condition as requiring that the maxima of these invariants remain sub-Planckian for every mass. In the de Sitter-core family this reduces to

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$25

and for the benchmark function it is enough to require

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$26

The paper also develops a Minkowski-core family, shows that the maxima of $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$27, $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$28, and $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$29 occur away from the center, and finds disconnected black-hole mass bands, vanishing Hawking temperature at band edges, stable endpoints of evaporation, and weaker extragalactic gamma-ray bounds for primordial-black-hole dark matter than in the Schwarzschild case (Aşmanoğlu et al., 24 Jul 2025).

Another 2024 construction begins with two analytic charged regular black-hole solutions from nonlinear electrodynamics and then generates two limiting-curvature-condition versions by replacing the charge parameter $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$30 with a fundamental length $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$31. The resulting metrics are

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$32

and

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$33

For the first LCC solution, the large-mass limits are

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$34

For the second, the dominant maxima are reported as

$|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$35

The same paper studies null, weak, dominant, and strong energy conditions, linear stability criteria of Moreno–Sarbach type, circular photon orbits, shadows, and eikonal quasinormal modes for these geometries (Balart et al., 2024).

Several limitations recur across the subject. The four-dimensional action-based proof of bounded polynomial curvature invariants is explicitly symmetry-reduced to the spherically symmetric interior ansatz and does not address perturbative stability, rotating or charged black holes, collapse matter, or quantum particle creation; the geometry is framed as that of an eternal black hole rather than a full collapse spacetime (Frolov, 2022). The phenomenological metric models are likewise not derived from a unique underlying ultraviolet completion, and the regulator function $|{\cal R}|<\Lambda=\frac{B_{\cal R}}{\ell^2},$36 is not uniquely selected (Aşmanoğlu et al., 24 Jul 2025). A plausible implication is that the limiting curvature condition is best viewed not as a single theory but as a program: select or dynamically constrain a finite set of curvature quantities so that high-curvature evolution is diverted away from singular geometry and toward a bounded-curvature phase.