Holomorphic Flow Equations

• Holomorphic flow equations are evolution equations defined by holomorphic vector fields that enforce analytic continuation and enable rigid classifications in complex domains.
• They provide a unified framework for analyzing planar dynamics, singularity structures, and global linearization in ℂⁿ through precise rectification and phase portrait techniques.
• Applications extend to complex geometry, spectral Hamiltonian reformulations, and power-flow modeling, demonstrating their versatility in both pure and applied mathematical contexts.

Searching arXiv for papers on holomorphic flow equations and related holomorphic flows. arxiv_search(query="holomorphic flow equations holomorphic flows complex dynamics", max_results=10, sort_by="submittedDate") Holomorphic flow equations are evolution equations whose defining vector fields, symmetry groups, or recursive amplitudes are constrained by holomorphy. In current research they occur in several technically distinct settings: planar systems

$\dot z=f(z),$

with $f$ holomorphic and real-time parameterization; holomorphic conformal flows on locally conformally Kähler manifolds; perturbed holomorphic cylinders converging to configurations joined by a flow line; Hamiltonian systems built from the holomorphic $\xi$-flow; positive-time flow invariant domains and global linearization problems on $\mathbb C^n$; refined holomorphic anomaly equations for Wilson-loop BPS sectors; and holomorphic embedding formulations of AC power flow (Rondón et al., 6 Jan 2026, Kainz et al., 2024, Ornea et al., 2010, Cant, 2021, Lebiedz, 2020, Chatterjee et al., 29 Jun 2026, Wang, 2023, Trias, 2015).

1. Analytic framework and model equations

A basic holomorphic flow equation in one complex variable is

$\dot z=f(z), \qquad f \text{ holomorphic}.$

A closely related classical construction starts from a holomorphic complex potential

$\Omega(z)=\phi(x,y)+i\psi(x,y),$

with associated velocity field

$\dot z=\overline{\Omega'(z)}.$

In that setting, $\overline{\Omega'}=\nabla\phi$, while $\psi$ is a first integral; the real and imaginary parts of $\Omega$ therefore generate orthogonal foliations, representing equipotential lines and streamlines (Rondón et al., 6 Jan 2026).

The modern generalization to arbitrary holomorphic $f$0 is to define a potential from the reciprocal field,

$f$1

If $f$2, then

$f$3

so the rectifying coordinate $f$4 transforms the system into the constant vector field $f$5. The trajectories are precisely the level curves

$f$6

Because primitives of $f$7 may require cutting the domain, the construction is often carried out on a star-shaped domain obtained by removing a ray from a punctured disk (Rondón et al., 6 Jan 2026).

This analytic picture already separates several meanings of “holomorphic flow equation.” In some papers the phrase refers to an actual real-time holomorphic ODE; in others it refers to holomorphic symmetry flows on complex manifolds or to recursive equations whose non-holomorphic dependence is encoded by propagators. This suggests that the unifying feature is not a single canonical equation but the use of holomorphy to rigidify continuation, classification, or recursion (Rondón et al., 6 Jan 2026, Wang, 2023).

2. Planar dynamics, singularities, and local classification

For planar systems

$f$8

with $f$9, holomorphy imposes strong restrictions on equilibria and orbit geometry. If $\xi$0 is a simple equilibrium with

$\xi$1

then the Jacobian has the Cauchy–Riemann form

$\xi$2

so a simple equilibrium can be a node, a focus, or a center/focus when $\xi$3, but it cannot be a saddle. If $\xi$4 is a zero of order $\xi$5, then the definite directions are exactly

$\xi$6

with $\xi$7. Every orbit tending to $\xi$8 does so in one of these definite directions, and the local phase portrait decomposes into a finite elliptic decomposition of order $\xi$9. On simply connected domains, a holomorphic Poincaré–Bendixson-type theorem shows that a bounded non-periodic orbit has a limit set consisting of exactly one equilibrium; for entire holomorphic vector fields, bounded non-periodic orbits are homoclinic or heteroclinic (Kainz et al., 2024).

The rectifying potential framework makes these restrictions explicit. For cubic monic centered polynomial systems

$\mathbb C^n$0

the global phase portrait on the Poincaré disk decomposes into $\mathbb C^n$1, $\mathbb C^n$2, or $\mathbb C^n$3 canonical regions of center-type, sepal-type, or $\mathbb C^n$4–$\mathbb C^n$5-type. In piecewise holomorphic or anti-holomorphic systems, the same complex-potential machinery converts the existence of periodic orbits into explicit matching equations on the switching line. The resulting bounds are sharp in several cases: for a mixed anti-holomorphic–holomorphic linear system there is at most one limit cycle; in the shifted equilibrium case there are at most three; and for piecewise anti-holomorphic polynomial systems the bounds are degree-dependent—linear: no limit cycles, quadratic: at most one, cubic: at most three (Rondón et al., 6 Jan 2026).

Holomorphic flows also act on singularity theory. For a singular holomorphic vector field

$\mathbb C^n$6

with flow $\mathbb C^n$7, the time-$\mathbb C^n$8 map is $\mathbb C^n$9, and the induced deformation of a plane branch $\dot z=f(z), \qquad f \text{ holomorphic}.$0 is

$\dot z=f(z), \qquad f \text{ holomorphic}.$1

The first Puiseux coefficient changed by the deformation is the contact exponent, while the tangency order is

$\dot z=f(z), \qquad f \text{ holomorphic}.$2

These notions govern which coefficients can be eliminated by holomorphic flows in Zariski’s moduli problem. The paper on analytic moduli of plane branches proves

$\dot z=f(z), \qquad f \text{ holomorphic}.$3

but also exhibits a non-complete analytic class, with explicit example

$\dot z=f(z), \qquad f \text{ holomorphic}.$4

showing that analytic equivalence need not be realizable by a single embedded holomorphic flow (Ayuso et al., 2017).

3. Holomorphic conformal flows in complex and Hermitian geometry

A complex manifold $\dot z=f(z), \qquad f \text{ holomorphic}.$5 of complex dimension $\dot z=f(z), \qquad f \text{ holomorphic}.$6 is locally conformally Kähler (LCK) if it admits a Kähler covering

$\dot z=f(z), \qquad f \text{ holomorphic}.$7

such that the deck transformation group $\dot z=f(z), \qquad f \text{ holomorphic}.$8 acts by holomorphic homotheties of the Kähler metric. Equivalently, $\dot z=f(z), \qquad f \text{ holomorphic}.$9 carries a Hermitian form $\Omega(z)=\phi(x,y)+i\psi(x,y),$0 and a closed $\Omega(z)=\phi(x,y)+i\psi(x,y),$1-form $\Omega(z)=\phi(x,y)+i\psi(x,y),$2 with

$\Omega(z)=\phi(x,y)+i\psi(x,y),$3

On a Kähler covering $\Omega(z)=\phi(x,y)+i\psi(x,y),$4, the monodromy character is defined by

$\Omega(z)=\phi(x,y)+i\psi(x,y),$5

and an automorphic potential is a function $\Omega(z)=\phi(x,y)+i\psi(x,y),$6 such that

$\Omega(z)=\phi(x,y)+i\psi(x,y),$7

Vaisman manifolds provide an explicit example, with

$\Omega(z)=\phi(x,y)+i\psi(x,y),$8

for the lifted Lee form (Ornea et al., 2010).

The central result for holomorphic conformal flows is a characterization of LCK manifolds with automorphic potential. If a compact LCK manifold admits a holomorphic conformal flow

$\Omega(z)=\phi(x,y)+i\psi(x,y),$9

which lifts to a flow of non-isometric homotheties on the minimal Kähler covering, then the closure $\dot z=\overline{\Omega'(z)}.$0 of the generated group is connected and contains the monodromy group. A crucial step is that on a Kähler manifold of complex dimension $\dot z=\overline{\Omega'(z)}.$1, every holomorphic conformal map is automatically a homothety. From this, together with a Gauduchon metric

$\dot z=\overline{\Omega'(z)}.$2

the existence of an automorphic potential follows. Conversely, if $\dot z=\overline{\Omega'(z)}.$3 admits an automorphic potential, then one can choose an LCK metric $\dot z=\overline{\Omega'(z)}.$4 with the same monodromy and construct a conformal flow of holomorphic diffeomorphisms whose lift acts by non-trivial homotheties on the covering. The proof uses an embedding into a Hopf manifold

$\dot z=\overline{\Omega'(z)}.$5

the formal logarithm $\dot z=\overline{\Omega'(z)}.$6, and averaging over an induced $\dot z=\overline{\Omega'(z)}.$7-action. The paper explicitly weakens the symmetry assumption in the earlier Kamishima–Ornea theorem: instead of a holomorphic conformal $\dot z=\overline{\Omega'(z)}.$8-action yielding a Vaisman conclusion, a holomorphic conformal flow already characterizes the broader class of LCK manifolds with automorphic potential (Ornea et al., 2010).

The same paper also records an important correction: an earlier stronger claim that any LCK manifold with potential has monodromy $\dot z=\overline{\Omega'(z)}.$9 is false. This correction matters because it separates the existence of an automorphic potential from an overly restrictive monodromy description (Ornea et al., 2010).

A different geometric use of “flow” appears in the Chern-Ricci flow,

$\overline{\Omega'}=\nabla\phi$0

which evolves Hermitian metrics rather than holomorphic diffeomorphisms. On the Hopf manifold $\overline{\Omega'}=\nabla\phi$1, an explicit solution shows that non-negativity of holomorphic bisectional curvature is not preserved along this flow, in contrast with the Kähler-Ricci case (Yang, 2015).

4. Perturbed holomorphic cylinders and flow-line breaking

In symplectic and Floer-type analysis, holomorphic flow equations arise through perturbed Cauchy–Riemann equations on long cylinders,

$\overline{\Omega'}=\nabla\phi$2

for maps

$\overline{\Omega'}=\nabla\phi$3

with $\overline{\Omega'}=\nabla\phi$4 and $\overline{\Omega'}=\nabla\phi$5 smoothly. If $\overline{\Omega'}=\nabla\phi$6 is a gradient vector field, this is the finite-cylinder version of Floer’s equation. The compactness theorem proved in this setting states that, after passing to a subsequence and choosing $\overline{\Omega'}=\nabla\phi$7 with $\overline{\Omega'}=\nabla\phi$8, the two translated ends converge to holomorphic half-cylinders with removable singularities $\overline{\Omega'}=\nabla\phi$9 and $\psi$0, while the rescaled middle converges uniformly to a flow line $\psi$1 of $\psi$2 of length $\psi$3 joining $\psi$4 to $\psi$5 (Cant, 2021).

The analytical core is the center of mass

$\psi$6

and the oscillation functional

$\psi$7

An a priori exponential estimate yields a differential inequality of the form

$\psi$8

and hence exponential decay of $\psi$9. Elliptic bootstrapping upgrades this to $\Omega$0-estimates. After rescaling by $\Omega$1, the center-of-mass equation becomes an asymptotic ODE for $\Omega$2, and the middle region collapses to a genuine flow line. The resulting limiting object is therefore a “holomorphic-flow-line-holomorphic” configuration rather than a single holomorphic cylinder (Cant, 2021).

5. Hamiltonian and spectral realizations of holomorphic flows

The holomorphic $\Omega$3-flow

$\Omega$4

where $\Omega$5 is the Riemann $\Omega$6-function, has been recast as a complex Hamiltonian system with

$\Omega$7

Hamilton’s equations give

$\Omega$8

so the $\Omega$9-equation is the original holomorphic flow and the $f$00-equation is its variational evolution. Eliminating time yields

$f$01

Using the product representation

$f$02

one obtains the logarithmic derivative

$f$03

and the implicit phase portrait

$f$04

The paper interprets this phase portrait $f$05 as a Riemann surface whose branching is governed by the zero set $f$06 (Lebiedz, 2020).

The variational differential of the flow map contains the same spectral data. One finds

$f$07

and the full differential matrix includes the sum

$f$08

This is described as a “spectral sum” structure reminiscent of trace formulas. A nonlinear reparameterization of time,

$f$09

transforms the system into the Newton-flow form

$f$10

For a simple zero $f$11, the closed-orbit period is

$f$12

Quantization on a closed orbit circle with line element $f$13 produces the Dirac-type operator

$f$14

whose spectrum is determined by the classical periods,

$f$15

The construction is explicitly presented as an analogy with trace formulas and not as a proof of the Riemann hypothesis (Lebiedz, 2020).

6. Flow-invariant Runge domains and global linearization on $f$16

For a holomorphic vector field $f$17 on a domain $f$18, the flow $f$19 solves

$f$20

If $f$21 is complete and $f$22 is a globally attracting fixed point, flow invariance can be converted into approximation-theoretic and conjugacy statements. In the linear case $f$23, with decomposition

$f$24

a positive-time invariant domain $f$25 is Runge if it contains $f$26 and the orbitwise thickness condition

$f$27

holds for large $f$28 on every compact set $f$29. A simpler corollary states that if $f$30 is positive-time invariant, contains $f$31 and $f$32, and

$f$33

then $f$34 is Runge. The paper also gives examples of non-plurisubharmonic Runge domains and a counterexample showing that without the distance hypothesis the conclusion can fail: $f$35 is positive-time invariant for a linear flow but is not Runge (Chatterjee et al., 29 Jun 2026).

The same paper proves a global linearization theorem. Write

$f$36

If, for every compact $f$37, the improper integral

$f$38

converges uniformly on $f$39, then

$f$40

exists locally uniformly, defines an automorphism $f$41, and satisfies

$f$42

A corollary gives a concrete spectral-gap criterion: if

$f$43

near $f$44 and

$f$45

then the flow is globally linearizable by an automorphism of $f$46. This produces an explicit conjugating map and shows that global holomorphic dynamics can sometimes be reduced to its linear model by a limit-normalization procedure (Chatterjee et al., 29 Jun 2026).

7. Recursive holomorphic flow equations in physics and network theory

In refined topological string theory and five-dimensional $f$47 theories on

$f$48

the phrase “holomorphic flow equation” appears in the form of refined holomorphic anomaly equations. For Wilson-loop expectation values,

$f$49

the full free energies $f$50 satisfy the same refined anomaly equation as the ordinary topological-string amplitudes. The paper’s central new result is that the BPS sectors themselves satisfy a refined holomorphic anomaly equation,

$f$51

with a rank-one specialization in one modulus $f$52. These equations are solved by direct integration, using regularity at the conifold point, regularity at the orbifold point, and large-volume asymptotics. Expanding Wilson-loop expectation values around the conifold point yields quantum spectra of the associated quantum Hamiltonians, and the same structure leads to a generalized blowup equation in which the factor $f$53 becomes a linear combination of Wilson-loop expectation values (Wang, 2023).

In power-systems analysis, the Holomorphic Embedding Load-Flow Method reframes the AC power-flow equations as a holomorphic family on an algebraic curve. Starting from

$f$54

one introduces an embedding parameter $f$55 so that $f$56 is a trivial no-load state and $f$57 is the physical problem. Because complex conjugation is not holomorphic, the embedded system uses independent holomorphic unknowns $f$58 and $f$59, with the reflection condition

$f$60

After clearing denominators, the equations define an algebraic curve; the operational solution is the analytic continuation of the “white germ” from $f$61 to $f$62. Power-series coefficients are computed recursively, and near-diagonal Padé approximants provide analytic continuation. Stahl’s theorem gives the completeness statement: if the white germ can be continued to $f$63, Padé approximants recover the correct branch, while nonconvergence signals infeasibility. The method also distinguishes physical branches from “ghost” branches, which solve the algebraic system but violate the reflection condition (Trias, 2015).

Later work sharpened the handling of PV buses and multidimensional loading. A general parametrized PV/PQ formulation was introduced to avoid the triple-product term

$f$64

which caused double convolutions and accuracy problems in an earlier PV model. The newer bilinear formulation proves, via the Complex Implicit Function Theorem, that the reflecting condition is redundant rather than an extra assumption, and it provides several practical model variants, with Model 4 reported as especially strong numerically on standard IEEE cases. The multidimensional holomorphic embedding method then assigns separate scales $f$65 to loads, powers, or load groups and represents each bus voltage by a multivariate power series

$f$66

where the number of monomials of total degree $f$67 is

$f$68

This turns the operating-condition space into an explicitly parameterized analytic object that can be prepared offline and evaluated online by substitution (Wallace et al., 2016, Liu et al., 2017).

Across these settings, holomorphic flow equations serve different immediate purposes—rectification of complex ODEs, classification of equilibria, characterization of LCK manifolds with potential, compactness and degeneration of perturbed holomorphic curves, Hamiltonian reformulations with spectral data, Runge approximation and linearization on $f$69, recursive anomaly equations in B-model physics, and constructive solution of nonlinear network equations. The common structural theme is that holomorphy constrains continuation, excludes generic planar pathologies, and replaces many local existence arguments by exact analytic or algebraic mechanisms.