Modified Ricci-Harmonic Flow

• Modified Ricci-Harmonic flow is a family of coupled geometric evolution equations that blend Ricci flow with harmonic map heat flow to analyze curvature behavior and singularity formation.
• The approach integrates auxiliary fields such as scalar potentials and weighted measures to define effective curvature tensors and corresponding entropy functionals.
• It leverages soliton analysis, gauge fixing, and singularity classification to elucidate the stability, convergence, and structure of diverse geometric settings.

Modified Ricci-Harmonic flow denotes a family of coupled geometric evolution equations rather than a single canonical PDE. In the standard usage, it is the Ricci flow coupled to harmonic map heat flow, usually written for an evolving metric $g(t)$, a map $\phi(t)$, and often a time-dependent coupling $\alpha(t)$. In parallel, the literature also uses closely related “modified” terminology for Perelman/List-type systems with an auxiliary scalar potential or diffusivity $f$, and for weighted-measure formulations on manifolds with boundary. The subject therefore combines a metric evolution, a heat-type evolution of an auxiliary field, and a spectrum of entropy, compactness, and soliton structures whose precise form depends on which modification is being used (Matteo, 2018, Chatterjee et al., 2013, Gomes et al., 27 Oct 2025).

1. Terminology and basic equations

The standard coupled system is the harmonic Ricci flow, also called Ricci-harmonic flow, on a domain $(M^n,g(t))$ with target $(N^k,\gamma)$, coupling function $\alpha(t)\ge 0$, and map $\phi(t):M\to N$: $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$ A standard derived quantity is the modified curvature tensor

$\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$

which plays the role of an effective Ricci tensor and effective scalar curvature in the coupled theory. When $\phi(t)$0 and $\phi(t)$1 is real-valued, one obtains List’s flow; in one scalar normalization this appears as

$\phi(t)$2

with

$\phi(t)$3

These formulations are the baseline meaning of Ricci-harmonic flow in the modern geometric-analysis literature (Matteo, 2018, Li, 2010).

A second, genuinely modified, formulation introduces a potential $\phi(t)$4 and preserves the weighted measure $\phi(t)$5. On a compact manifold with boundary, the modified $\phi(t)$6 flow is

$\phi(t)$7

$\phi(t)$8

together with boundary conditions

$\phi(t)$9

This system is the negative gradient flow of a weighted Gibbons-Hawking-York type functional and is explicitly designed to encode Perelman–Lott–List-type weighted geometry in the Ricci-harmonic setting (Gomes et al., 27 Oct 2025).

A recurrent source of ambiguity is the Perelman-inspired scalar-diffusivity flow studied on asymptotically non-flat spaces: $\alpha(t)$0 That system has a metric and a scalar obeying a heat equation, but no separate evolution equation for a map $\alpha(t)$1. It is therefore not the standard Ricci-harmonic flow in the metric-map sense, even though it is structurally analogous to scalar-coupled Ricci-harmonic-type systems (Chatterjee et al., 2013).

2. Solitons, fixed points, and entropy functionals

Self-similar solutions organize the subject. For the metric-map formulation, a gradient harmonic Ricci soliton with potential $\alpha(t)$2 satisfies

$\alpha(t)$3

with $\alpha(t)$4, $\alpha(t)$5, and $\alpha(t)$6 corresponding respectively to shrinking, steady, and expanding solitons. The trace identities

$\alpha(t)$7

supply the elliptic structure behind soliton rigidity and normalization. In canonical shrinking form one has

$\alpha(t)$8

$\alpha(t)$9

The same pattern reappears in weighted variants, where critical points of the weighted functional are gradient steady $f$0 solitons subject to the boundary conditions $f$1 and $f$2 (Matteo, 2018, Gomes et al., 27 Oct 2025).

Perelman-type functionals have several coupled analogues. In the scalar case,

$f$3

and, along the appropriate coupled system,

$f$4

Equality is equivalent to the gradient steady harmonic-Ricci soliton equations. The same paper develops $f$5, $f$6, shrinker and expander $f$7-type functionals, and the associated $f$8-entropies; the critical point equations are exactly the expanding or shrinking harmonic-Ricci soliton equations (Li, 2010).

For the boundary-weighted theory, the central object is

$f$9

where

$(M^n,g(t))$0

Its gradient is the modified $(M^n,g(t))$1 system above, and its first variation separates into bulk Euler-Lagrange terms

$(M^n,g(t))$2

plus boundary terms involving the second fundamental form and normal derivative of $(M^n,g(t))$3 (Gomes et al., 27 Oct 2025).

3. Singularity formation, compactness, and reduced volume

The modern singularity theory of harmonic Ricci flow parallels the Ricci-flow theory of Type I blow-up. A solution on $(M^n,g(t))$4, $(M^n,g(t))$5, is Type I if

$(M^n,g(t))$6

A decisive refinement is that a Type I curvature bound automatically controls the map sector: if

$(M^n,g(t))$7

then there is $(M^n,g(t))$8 such that

$(M^n,g(t))$9

and

$(N^k,\gamma)$0

After parabolic rescaling about a Type I singularity, every subsequential pointed limit is a normalized gradient shrinking harmonic Ricci soliton in canonical form, and at a Type I singular point the limit is nontrivial, i.e. not locally Gaussian (Matteo, 2018).

The analytic infrastructure consists of refined compactness theorems, a pseudolocality theorem, and reduced length and reduced volume based at the singular time. The singular-time reduced volume is monotone, and constant reduced volume characterizes gradient shrinking harmonic Ricci solitons through the vanishing of the coupled analogue of Perelman’s $(N^k,\gamma)$1-identity. This framework yields the identification

$(N^k,\gamma)$2

so the scalar Type I set, the essential blow-up set, and the ordinary singular set coincide for Type I harmonic Ricci flows (Matteo, 2018).

A complementary formulation uses the augmented curvature

$(N^k,\gamma)$3

and classifies singularities by the growth of $(N^k,\gamma)$4 or $(N^k,\gamma)$5. In this language one has Type I, IIa, IIb, and III singularities exactly as in Hamilton’s Ricci-flow classification. Compact blow-up limits of finite-time singularities are shrinking Ricci harmonic solitons via the monotonicity of the coupled entropy $(N^k,\gamma)$6, and ancient solutions satisfy

$(N^k,\gamma)$7

both on compact manifolds and, with localized maximum-principle arguments, on complete noncompact manifolds (Shi, 2013).

4. Harmonic map heat flow as gauge and the shrinker asymptotics problem

A distinct but closely related line of work uses harmonic map heat flow as a gauge for Ricci flow near a compact gradient shrinker $(N^k,\gamma)$8. In rescaled variables, the metric satisfies the modified rescaled Ricci flow

$(N^k,\gamma)$9

while the corresponding modified harmonic map heat flow is

$\alpha(t)\ge 0$0

The gauge-fixed metric perturbation $\alpha(t)\ge 0$1 then solves a strictly parabolic equation

$\alpha(t)\ge 0$2

and the map perturbation satisfies a compatible heat-type equation driven by

$\alpha(t)\ge 0$3

on vector fields. The paper explicitly interprets this as coupling Ricci flow to a modified harmonic map heat flow and uses it to construct a global parabolic gauge (Choi et al., 3 Apr 2025).

For an ancient Ricci flow asymptotic to a compact integrable shrinker, or a Ricci flow developing a finite-time singularity modeled on such a shrinker, there exists an ancient or immortal harmonic map heat flow between the evolving metric and the shrinker for all times. This produces a global gauge in which both the map perturbation and the gauge-fixed metric perturbation decay exponentially. Two consequences are central: all ancient Ricci flows asymptotic to a compact integrable shrinker are constructed and classified, and tangent flows are unique without passing to equivalence modulo diffeomorphisms. The optimal convergence rate at a singularity is governed by the first negative eigenvalue of the stability operator for the entropy; in particular, a Ricci flow developing a round $\alpha(t)\ge 0$4 singularity converges at least at the rate

$\alpha(t)\ge 0$5

This suggests a gauge-theoretic version of modified Ricci-harmonic flow in which the harmonic map sector is not additional geometry but the mechanism that removes diffeomorphism degeneracy (Choi et al., 3 Apr 2025).

5. Asymptotically non-flat and asymptotically hyperbolic modifications

In the Perelman-inspired scalar-diffusivity variant, the fixed-point equations

$\alpha(t)\ge 0$6

were analyzed on the spatial metric induced by Marder’s cylindrically symmetric asymptotically non-flat vacuum solution. With the radial ansatz $\alpha(t)\ge 0$7, the system admits the explicit solution

$\alpha(t)\ge 0$8

and, in the full spacetime metric, this satisfies

$\alpha(t)\ge 0$9

When the full modified flow is imposed, consistency of the $\phi(t):M\to N$0, $\phi(t):M\to N$1, and $\phi(t):M\to N$2 component equations forces

$\phi(t):M\to N$3

For that special case, the area $\phi(t):M\to N$4 of a cylindrical surface obeys

$\phi(t):M\to N$5

where $\phi(t):M\to N$6 is the compactness; since $\phi(t):M\to N$7, the area decreases along the flow. The paper is explicit that it does not establish global existence, singularity formation, or uniqueness, and that the result is an existence demonstration in a particular case rather than a general theory (Chatterjee et al., 2013).

A different asymptotic regime appears in the normalized Ricci-harmonic flow adapted to $\phi(t):M\to N$8: $\phi(t):M\to N$9 With the substitution $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$0, the static vacuum equations

$\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$1

become

$\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$2

so asymptotically hyperbolic static vacuum triples are fixed points, or self-similar expanders, of the normalized flow. The paper constructs a renormalized expander entropy $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$3, identifies static pairs as its critical points, proves that its gradient flow is the normalized Ricci-harmonic system modulo gauge, and shows that dynamical stability is equivalent to a local positive mass property. In this sense, static vacuum spacetimes with $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$4 are attractors of a modified Ricci-harmonic flow exactly when the relevant positive-mass statement holds nearby (Jouttijärvi et al., 24 Apr 2026).

6. Boundary-weighted theories and coupling to mean curvature flow

The boundary theory begins from the weighted scalar and mean curvatures

$\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$5

and the functional

$\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$6

Along the modified $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$7 flow, its derivative has bulk terms

$\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$8

plus a boundary contribution involving $\frac{\partial g(t)}{\partial t} = -2 \Ric_{g(t)} + 2\alpha(t)\,\nabla\phi(t)\otimes\nabla\phi(t),\qquad \frac{\partial \phi(t)}{\partial t} = \tau_{g(t)}\phi(t).$9, tangential derivatives of $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$0 and $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$1, ambient curvature terms, and the additional coupled term $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$2. This extends the weighted Gibbons-Hawking-York framework previously developed for Perelman’s modified Ricci flow to the genuinely coupled Ricci-harmonic setting (Gomes et al., 27 Oct 2025, Lott, 2011).

When the boundary or a hypersurface evolves by mean curvature flow in a gradient steady $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$3 soliton background, the resulting boundary integrand becomes an extension of Hamilton’s differential Harnack expression. In the general steady soliton background, the identity is

$\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$4

so the Harnack quantity vanishes identically on the steady soliton background. In the flat Euclidean case with constant $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$5, this reduces to Hamilton’s original condition for a translating soliton (Gomes et al., 27 Oct 2025).

The same background supports a Huisken monotonicity-type formula. For a mean curvature flow $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$6 in a gradient $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$7 soliton background with potential $\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$8, the weighted area

$\mathcal S = \Ric - \alpha \nabla\phi\otimes\nabla\phi,\qquad \Sh = \Sc - \alpha |\nabla\phi|^2,$9

in the steady case, and

$\phi(t)$00

in the shrinking or expanding cases, is non-increasing. Equality holds exactly for mean curvature solitons satisfying

$\phi(t)$01

Moreover, these mean curvature solitons are characterized by $\phi(t)$02-minimal hypersurfaces in the initial metric, so the modified Ricci-harmonic structure controls both bulk entropy and interface self-similarity (Gomes et al., 27 Oct 2025).

7. Exceptional-geometry analogues and explicit singularity models

The phrase “Ricci-harmonic flow” has also been extended to $\phi(t)$03- and Spin(7)-structures, where the evolving variable is a defining differential form rather than a metric-map pair. For a $\phi(t)$04-structure $\phi(t)$05 with torsion tensor $\phi(t)$06 and dual 4-form $\phi(t)$07, the Ricci-harmonic flow of $\phi(t)$08-structures is

$\phi(t)$09

The paper derives this from the Taylor expansion of $\phi(t)$10 in adapted normal coordinates and interprets it as the heat equation for $\phi(t)$11-structures. On compact manifolds it proves short-time existence and uniqueness, global and local Shi-type estimates, a compactness theorem, a long-time existence criterion in terms of bounded velocity, and the rigidity statement that stationary points are exactly torsion-free $\phi(t)$12-structures. Compact expanding solitons do not exist, and the only compact steady solitons are torsion-free. A parallel Spin(7) flow is derived by the same method (Dwivedi, 8 Jan 2026).

Explicit singularity models were then obtained on 7-dimensional contact Calabi-Yau manifolds and on the 7-dimensional Heisenberg group. For the natural co-closed $\phi(t)$13-structure on a contact Calabi-Yau manifold, the Ricci-harmonic flow admits an ancient solution with a finite-time Type I singularity. The same ansatz yields immortal solutions for the Ricci-like flows associated with the $\phi(t)$14-Einstein-Hilbert functional, with infinite-time singularities that are Type III when the transversal Calabi-Yau distribution is flat and Type IIb otherwise; analogous solutions are also obtained for the negative gradient flow of the torsion energy and on the Heisenberg group. These are presented as the first examples of Type I singularities for the Ricci-harmonic flow and of Type IIb and Type III singularities for the Ricci-like flows in this exceptional-holonomy context (Dwivedi et al., 23 Jan 2026).

The accumulated literature therefore supports a broad but coherent view of modified Ricci-Harmonic flow. At one pole lie the standard metric-map systems of Müller/List type, with their soliton equations, entropy functionals, reduced-volume monotonicity, and Type I singularity theory. At another lie weighted, boundary, scalar-diffusivity, gauge-theoretic, and exceptional-geometry variants, each retaining the central pattern of a Ricci-type metric evolution coupled to a diffusion equation for an auxiliary field. This suggests that “modified Ricci-Harmonic flow” is best understood as a family of Perelman-style coupled parabolic systems whose unifying themes are weighted curvature, self-similar solutions, and monotonicity-driven control of singularities.