---
title: Combinatorial Ricci Flow
url: https://www.emergentmind.com/topics/combinatorial-ricci-flow
type: topic
---

# Combinatorial Ricci Flow

Combinatorial Ricci flow denotes a family of discrete curvature flows on triangulated or cell-decomposed spaces in which discrete conformal parameters evolve by curvature defects. On surfaces, the basic parameters are typically circle-packing radii or logarithmic conformal factors, and the curvature is usually an angle defect of the form $K_i=2\pi-\sum\theta_i^{jk}$; for surfaces with boundary or in spherical background geometry, the prescribed quantity may instead be a boundary length or a total geodesic curvature; in three dimensions, the curvature becomes the cone-angle deficit along edges of ideal or hyper-ideal tetrahedral decompositions. Across these settings, the theory is organized by convex or proper energy functionals, symmetry and definiteness properties of curvature Jacobians, extensions through degeneracies, and convergence theorems that make the flow an effective algorithm for constructing discrete uniformizations and hyperbolic structures [1806.05812] [2208.03643] [2303.09191] [2009.02496].

## 1. Discrete curvature models on triangulated surfaces

On a weighted triangulation $(T,\Theta)$ of a connected, oriented, closed surface, a circle-packing metric is a vector of radii $r=(r_v)_{v\in V}\in \mathbb R_{>0}^V$. The edge length attached to an edge $e=uv$ is
\[
\ell(e;r)=
\begin{cases}
\sqrt{r_u^2+r_v^2+2\,r_u r_v\cos\Theta(e)},&\chi(S)=0,\\[4pt]
\operatorname{arccosh}\!\bigl(\cosh r_u\cosh r_v+\sinh r_u\sinh r_v\cos\Theta(e)\bigr),&\chi(S)<0,
\end{cases}
\]
and each face becomes a geodesic triangle with inner angles $\phi_{f,v}(r)$. The combinatorial curvature at $v$ is then
\[
K_v(r)=2\pi-\sum_{f\ni v}\phi_{f,v}(r).
\]
A circle-pattern metric is characterized by $K_v(r)=0$ for all vertices [1806.05812].

A broader discrete conformal framework assigns metrics through vertex functions $f:V\to\mathbb R$ together with coefficients $\varepsilon:V\to\{-1,0,1\}$ and $\eta:E\to\mathbb R$. In this setting the classical curvature remains $K_i=2\pi-\sum \theta_i^{jk}$, while the parameterized combinatorial $\alpha$-curvature is
\[
R_i^{(\alpha)}=u_i^{-\alpha}K_i,
\]
with $u_i=e^{2f_i}$ in Euclidean background and the corresponding hyperbolic $u$-coordinates determined by the chosen discrete conformal structure. The case $\alpha=0$ recovers the classical curvature [2105.14714].

For ideally triangulated surfaces with boundary, the basic object is a type $(-1,-1,-1)$ generalized circle-packing metric. If $\partial\Sigma$ has components labeled by $i=1,\dots,N$, with positive radii $r_i$ and edge weights $\Phi_{ij}>0$, then the edge lengths satisfy
\[
\cosh l_{ij}=\cosh\Phi_{ij}\,\sinh r_i\,\sinh r_j-\cosh r_i\,\cosh r_j.
\]
Gluing the corresponding right-angled hyperbolic hexagons produces a hyperbolic structure on $\Sigma$ with totally geodesic boundary. In this case the combinatorial curvature at boundary component $i$ is
\[
K_i=\sum_{\{i,j,k\}\in F}\theta_i^{jk},
\]
where $\theta_i^{jk}$ is the boundary-arc length at $i$ in the relevant hexagon [2208.03643].

## 2. Canonical flow equations and normalization conventions

The literature uses several equivalent or analogous formulations, depending on the discrete metric variable, the background geometry, and whether a target curvature is prescribed.

| Setting | Variable | Flow law |
|---|---|---|
| Chow–Luo surface flow | $r_i$ | $\displaystyle \frac{dr_i}{dt}=-2\,K(v_i)\,r_i$ |
| Closed surface, log variable | $x_v$ or $u_i$ | $\displaystyle \dot x_v=-K_v,\qquad \dot u_i=-K_i$ |
| Surface with boundary, prescribed $\bar K$ | $u_i=\log\tanh(r_i/2)$ | $\displaystyle \frac{du_i}{dt}=K_i(u)-\bar K_i$ |
| Prescribed $\alpha$-curvature | $u_i$ | $\displaystyle \frac{du_i}{dt}=R_i^*-R_i^{(\alpha)}$ |
| Spherical total geodesic curvature | $K_v=\ln\cot r_v$ | $\displaystyle \frac{dK_v}{dt}=-(L_v-\hat L_v)$ |
| Three-dimensional edge flow | $l_e$ | $\displaystyle \frac{dl_e}{dt}=K_e(l)\quad\text{or}\quad \frac{dl_e}{dt}=K_e(l)\,l_e$ |

These formulas exhibit a common pattern: the state variable is driven by a discrete curvature defect. The multiplicative factors $r_i$, $\sinh r_i$, $\sin(2r_i)$, or $l_e$ reflect the chosen coordinates rather than a contradiction in the underlying idea. In particular, prescribed-curvature versions replace $K$ by $K-\bar K$, $K-\hat K$, $L-\hat L$, or their higher-dimensional analogues [1104.2033] [1806.05812] [2208.03643] [2105.14714] [2303.09191] [2009.02496] [2009.03731].

## 3. Variational structure, Jacobians, and discrete parabolicity

A central feature of combinatorial Ricci flow is its variational formulation. For generalized circle packings on surfaces with boundary, the discrete conformal coordinates
\[
u_i=\log\tanh\frac{r_i}{2}
\]
make the curvature map $K:u\mapsto (K_1(u),\dots,K_N(u))$ smooth, and
\[
\Delta_{ij}=\frac{\partial K_i}{\partial u_j}
\]
defines a symmetric negative-definite $N\times N$ matrix. The associated Ricci energy
\[
\mathcal E(u)=-\int^u\sum_i (K_i-\bar K_i)\,du_i
\]
is strictly convex, and the linearization has strictly negative spectrum at the target metric [2208.03643].

In spherical background geometry with prescribed total geodesic curvatures, the local 1-form
\[
\omega_e=L_u\,dK_u+L_w\,dK_w
\]
is closed on each edge, hence admits a local convex potential $E_e(K_u,K_w)$. These edge potentials assemble into a global strictly convex functional
\[
E(K)=\sum_{e\in E}E_e(K_u,K_w)-\sum_{v\in V}\hat L_v K_v,
\]
and the flow in $K$-variables is
\[
\frac{dK}{dt}=-\nabla E(K).
\]
Properness of this energy is characterized exactly by the coherence inequalities
\[
\sum_{v\in X}\hat L_v<2\sum_{e\in E(X)}\Theta(e)
\]
for every nonempty subset $X\subset V$ [2303.09191].

For parameterized combinatorial curvature flows, the Ricci energy
\[
F(u)= - \sum_{\{i,j,k\}\in F} F_{ijk}(u_i,u_j,u_k)+\sum_i \int^{u_i}(2\pi-R_i^*\,r)\,dr
\]
satisfies
\[
\nabla F(u)_i = K_i-R_i^*u_i^\alpha,
\]
so the modified $\alpha$-Ricci flow is the negative-gradient flow of $F$. A key technical device is the continuous extension of triangle angles across degenerate configurations, which extends $F$ to a $C^1$ function on $\mathbb R^N$ [2105.14714].

In three dimensions the analogous energy is the co-volume functional. On an ideally triangulated compact $3$-manifold with boundary one defines
\[
H(l)=\sum_{\sigma\in T}F_\sigma(l)-2\pi\sum_{ij\in E}l_{ij},
\qquad
\frac{\partial H}{\partial l_{ij}}=-K_{ij},
\]
while in the cusped case
\[
C(\ell)=\sum_{\sigma\in T}\operatorname{cov}(\sigma)-2\pi\sum_{e\in E}\ell_e,
\qquad
\nabla C(\ell)=-K(\ell).
\]
Both identities come from Schläfli’s formula, and both realize the flow as a negative-gradient evolution of a convex functional [2009.02496] [2009.05477].

A plausible implication is that the various versions of combinatorial Ricci flow are best understood not as a single discretization, but as a variational paradigm: choose a discrete metric space, identify a curvature map with symmetric derivative structure, and run the negative-gradient flow of the corresponding energy.

## 4. Long-time existence and convergence on finite surfaces

For compact polyhedral surfaces without boundary and bounded combinatorial or Wald curvature, one can prove forward and backward short-time existence and uniqueness of the discrete or metric Ricci flow by approximating the polyhedral surface with smoothings, applying classical Ricci-flow results on the smooth approximants, and then passing to the limit. In the normalized flow on a compact closed surface, convergence to constant curvature occurs exponentially fast in the negative- and positive-curvature cases and at rate $O\!\bigl(\frac1{1+t}\bigr)$ in the flat case. The same source also states that on genus $g\ge 2$ polyhedral surfaces the unnormalized combinatorial flow never develops combinatorial singularities, whereas on the torus or sphere with Euclidean background triangle inequalities may fail late in the flow and local retriangulations by edge-flips may be used to continue it [1104.2033].

For hyperbolic background geometry, the modified curvature
\[
\widetilde K_i=\frac{K_i^{\rm cl}}{2\pi(\cosh r_i-1)}
\]
yields the flow $\dot u_i=-\widetilde K_i$. The convergence theorem is an if and only if statement: the flow converges exactly when a circle-packing metric with zero curvature exists. A stronger sufficient condition is also available: if the initial modified curvatures are all non-positive, then the flow exists for all time, the curvatures remain non-positive, the trajectory stays in a compact subset of $(0,\infty)^N$, and the solution converges exponentially fast to the unique zero-curvature metric. The paper also states the same equivalence for arbitrary smooth positive area elements $A_i(r)$ [1505.05076].

For ideally triangulated surfaces with boundary and prescribed boundary lengths $\bar K\in(0,\infty)^N$, the combinatorial Ricci flow
\[
\frac{du_i}{dt}=K_i(u)-\bar K_i
\]
has a particularly strong global theory. For any initial $u(0)\in\Omega(\Phi)$, the solution exists for all $t\ge 0$, remains in a compact subset of $\Omega(\Phi)$, and converges exponentially fast to the unique $\bar u\in\Omega(\Phi)$ with $K(\bar u)=\bar K$. The same paper proves the analogous result for the combinatorial Calabi flow. It also records an explicit Euler-type algorithm in the $u$-variables, with complexity $O(|F|)$ per step and $O(\log(1/\varepsilon))$ steps to reach accuracy $\varepsilon$, up to CFL-type restrictions on the step size [2208.03643].

For the parameterized $\alpha$-Ricci flow, singularities are handled by extending the angle functions by constants across degenerate triangles. Under the structure conditions
\[
\varepsilon_s\varepsilon_t+\eta_{st}>0,\qquad
\varepsilon_i\eta_{jk}+\varepsilon_j\eta_{ki}+\varepsilon_k\eta_{ij}\ge 0,
\]
and assuming the prescribed curvature is realized by some discrete conformal structure, the extended modified flow has a unique global solution. If additionally $\alpha\chi(M)\le 0$ in Euclidean background or $\alpha R^*\le 0$ in hyperbolic background, then the solution converges exponentially fast to the unique target metric. In particular, the $\alpha=0$ case recovers long-time convergence of Chow–Luo’s flow and Luo’s Yamabe flow under the usual sign conditions [2105.14714].

## 5. Degeneration, prescribed data, and spherical background geometry

A recurring issue is whether failure of the strict existence criterion destroys the asymptotic structure of the flow. For closed surfaces with $\chi(S)\le 0$, Takatsu analyzed the marginal case
\[
\Phi(U)\le 0\quad\forall\,U\subset V,
\qquad
Z_T=\{z\in V:\exists\,U\ni z,\ \Phi(U)=0\}\neq\varnothing.
\]
In this regime the combinatorial Ricci flow still satisfies $\lim_{t\to\infty}K_v(r(t))=0$ for every vertex, but some radii collapse: in the Euclidean case the full radius vector does not converge in $\mathbb R_{>0}^V$, while in the hyperbolic case the flow does converge in $\mathbb R_{>0}^V$ and the limiting metric vanishes exactly on the collapsing set. The limiting metric on the complement is the unique circle-pattern metric on the reduced complex [1806.05812].

In hyperbolic background geometry, degenerated generalized circle packings can be treated directly by prescribing the total geodesic curvature data
\[
T_i(k)=k_i\,L_i(k).
\]
For the interior target region $\Omega_T$, the map $\Phi:\Omega_k\to\Omega_T$, $k\mapsto T(k)$ is a homeomorphism, so every $T^*\in\Omega_T$ determines a unique nondegenerate metric and the Ricci flow converges exponentially fast to it. The map extends continuously to the boundary, and for $T^*$ on the partial boundary $\partial_0\Omega_T$ there is a unique degenerated metric solving $T(k)=T^*$. The corresponding flow exists for all time on the relevant boundary stratum and converges to that unique degenerated solution [2407.08496].

The spherical theory replaces angle deficit by total geodesic curvature. For a circle of radius $r_v$ one has $k_v=\cot r_v$, total length $\ell_v=\alpha_v\sin r_v$, and total geodesic curvature
\[
L_v=\alpha_v\cos r_v.
\]
The spherical prescribed-curvature flow converges for every initial data if and only if the target data satisfy the coherence inequalities
\[
\sum_{v\in X}\hat L_v<2\sum_{e\in E(X)}\Theta(e)
\]
for every nonempty subset $X\subset V$; under these inequalities the flow converges exponentially fast to the unique circle pattern realizing $L_v=\hat L_v$. The same source gives an explicit or semi-implicit Euler implementation with $O(|E|)$ work per step and overall complexity $O(|E|\cdot \log 1/\varepsilon)$ [2303.09191].

A further spherical refinement treats degenerated circle pattern metrics and spherical conical metrics with prescribed total geodesic curvatures. The existence theorem identifies the admissible sets $\mathcal L_1$, $\bar\partial\mathcal L$, and $\mathcal L=\mathcal L_1\cup\bar\partial\mathcal L$ for nondegenerate, degenerate, and general spherical conical metrics, respectively. The same work proves convergence of the prescribed combinatorial Ricci flows in three forms: the nondegenerate flow for $\hat L\in\mathcal L_1$, the restricted flow on an index subset $I$, and a mixed flow on the full space in which some geodesic curvatures remain zero. It also gives two distinct constructions of degenerated solutions, one on the reduced coordinate subspace and one on the full space with vanishing components preserved [2404.05982].

## 6. Infinite triangulations and noncompact analogues

Combinatorial Ricci flow extends beyond finite complexes to infinite, locally finite triangulations. For infinite disk triangulations in Euclidean and hyperbolic background geometries, the flow
\[
\frac{du_i}{dt}=-K_i(u),\qquad i\in V,
\]
is defined on the entire infinite vertex set. Long-time existence is obtained by exhausting the triangulation by finite simply connected subcomplexes, solving Dirichlet-boundary flows on the exhaustion pieces, and then using a maximum principle together with an Arzelà–Ascoli diagonal argument. Uniqueness holds under uniformly bounded curvature. In hyperbolic background, if the initial metric satisfies $K_i(u(0))\le 0$ for all vertices, then the flow converges to a zero-curvature packing. In the Euclidean case, the regular infinite hexagonal triangulation admits a small-data convergence theorem in which the Dirichlet energy
\[
E(u(t))=\sum_{\{i,j\}\in E_H}(u_i(t)-u_j(t))^2
\]
tends to zero and $\|u(t)\|_{\ell^\infty}\to 0$ [2504.05817].

The same noncompact philosophy appears in the theory of infinite ideal circle patterns. For a locally finite infinite cell decomposition $\mathcal D=(V,E,F)$ of $S^2$ or $\mathbb R^2$, one defines vertex curvature by angle deficits and runs
\[
\frac{du_v}{dt}=-K_v(u,\Theta).
\]
Under the ideal-pattern condition
\[
\sum_{\ell=1}^m \Theta_{e_\ell}=(m-2)\pi
\]
on each face, the flow exists for all $t\ge 0$. Uniqueness holds under either uniformly bounded curvature or uniformly bounded degrees. In the hyperbolic case, if the initial metric has $K_v(0)\le 0$ for all vertices, then the flow converges to a good ideal circle pattern with $K_v(\infty)=0$. The same work uses Thurston’s correspondence between ideal circle patterns on $S^2$ and convex ideal polyhedra in the Poincaré ball to characterize infinite ideal polyhedra in $\mathbb H^3$ and thereby solve Rivin’s infinite version under the stated combinatorial conditions [2506.05036].

In spherical background geometry on infinite cellular decompositions, the basic variables are radii $r_i\in(0,\pi/2)$ or logarithmic coordinates $u_i=\ln\cot r_i$, and the target quantity is the total geodesic curvature
\[
T_i=\alpha_i\cos r_i.
\]
Under the conditions
\[
\text{(S1)}\ \hat T_i>0,\qquad
\text{(S2)}\ \sum_{v\in U}\hat T_v<2\sum_{e\in E(U)}\Theta(e)
\]
for every finite subset $U\subset V$, the flow exists uniquely for all time. If in addition the initial data satisfy
\[
\text{(S3)}\ T_i(r(0))\ge \hat T_i\quad\forall i,
\]
then $r_i(t)$ converges and $\lim_{t\to\infty}T_i(r(t))=\hat T_i$ for every vertex. The proof uses finite exhaustions, uniform $C^0$–$C^2$ estimates, a maximum principle on infinite graphs, and a Lyapunov functional built from convex edge potentials [2505.05925].

## 7. Three-dimensional extensions and hyperbolization

In dimension three, combinatorial Ricci flow is formulated on ideal or hyper-ideal triangulations of pseudo-manifolds, compact manifolds with boundary, or cusped manifolds. The discrete metric variable is an edge-length vector, and the curvature at an edge is the cone-angle deficit
\[
K_e(l)=2\pi-\sum_{\sigma\ni e}\alpha_e^\sigma(l),
\]
or $\pi-\theta_e$ for certain boundary edges in the partially truncated setting. Zero curvature means that the tetrahedra glue to a complete hyperbolic metric with totally geodesic boundary or cusps, depending on the model [2009.02496] [2009.03731] [2502.06497].

For ideally triangulated compact $3$-manifolds with boundary, the basic flow is Luo’s additive equation
\[
\frac{dl_{ij}}{dt}=K_{ij}(l),
\]
which Xu shows can be extended through tetrahedral degenerations by Luo–Yang’s continuous extension of dihedral angles to all of $\mathbb R^6$. The extended flow exists uniquely for all $t\ge 0$, and if there exists an admissible zero-curvature hyper-ideal metric $l^*$, then every initial datum converges exponentially fast to $l^*$. This verifies Luo’s conjecture in that setting [2009.02496].

A multiplicative variant on pseudo-$3$-manifolds evolves
\[
\frac{d}{dt}l_e(t)=K_e(l(t))\,l_e(t),
\]
or, in the extended form, $\dot l_e=\widetilde K_e(l)\,l_e$. Under the hypothesis that every edge has valence at least $10$, there exists a unique zero-curvature hyper-ideal metric and the extended flow converges exponentially fast from arbitrary positive initial data. The paper also gives the estimate
\[
l_e^*\in \bigl[(3\max_e d_e)^{-1},\ \arccosh(3)\bigr].
\]
A later refinement lowers the valence threshold to $9$ and proves that for any compact $3$-manifold with boundary admitting an ideal triangulation with all edge valences at least $9$, there exists a unique complete hyperbolic metric with totally geodesic boundary realizing that triangulation, together with explicit upper and lower bounds contained in $(0,\arccosh 2]$ [2009.03731] [2601.15174].

For cusped $3$-manifolds, one works with decorated ideal tetrahedra and signed edge lengths. The co-volume functional on the affine cusp slice has gradient $-K$, so the extended flow
\[
\frac{d\ell}{dt}=-K(\ell)
\]
is again a negative-gradient flow of a convex energy. The main theorem states that an ideal triangulated cusped $3$-manifold admits a complete hyperbolic metric if and only if the extended combinatorial Ricci flow converges to a zero-curvature point, and when convergence occurs the rate is exponential. This provides a curvature-flow formulation dual to the Casson–Rivin volume-maximization program [2009.05477] [2009.05842].

The relation to geometric triangulations is made explicit in the partially truncated and decorated setting. There the extended combinatorial Ricci flow exists uniquely for all initial data in $\mathbb R^E$, and it converges if and only if the triangulation is geometric. Equivalently, there exists a complete hyperbolic structure on the manifold with the given ideal triangulation realized as a geometric ideal triangulation if and only if the flow converges to a zero-curvature metric. The same framework is presented as a systematic approach to Thurston’s triangulation conjecture [2502.06497].

Taken together, these results show that combinatorial Ricci flow has evolved from the original Chow–Luo surface equation into a broad discrete uniformization and hyperbolization program. The common structure is curvature-defect evolution, convex energy descent, and rigidity of the target metric; the major technical differences arise from background geometry, boundary behavior, degenerations, and the passage from finite to infinite or from two-dimensional to three-dimensional complexes.

Source: https://www.emergentmind.com/topics/combinatorial-ricci-flow