---
title: "Beyond the branching exponent"
date: 2026-10-05
version: "1.1.0-candidate"
doi: 10.5281/zenodo.23160583
pdf: https://github.com/ipitchford/horton-tokunaga-kingman-asymptotics/releases/download/v1.1.0-candidate/paper.pdf
repository: https://github.com/ipitchford/horton-tokunaga-kingman-asymptotics
archive: https://zenodo.org/records/23160583
license: CC0-1.0
status: unrefereed (internally replayed; not peer reviewed, not independently reproduced, not formally verified)
---

# Beyond the branching exponent

## Summary

A tree’s number of leaves does not tell us how deeply its branches are nested. Horton–Strahler order measures that nesting: leaves start at one; joining equal orders raises the order by one, while joining unequal orders retains the larger value. A balanced four-leaf tree has order three; a four-leaf comb has order two.

This candidate sharpens two mathematical descriptions of random branching. For critical Tokunaga trees, it describes the order when the leaf count is fixed **exactly**, including fluctuations around logarithmic growth. For Kingman’s coalescent, it shows that limiting branch frequencies have a positive geometric prefactor—not merely a known exponential rate.

These are two model-specific theorems, not a claim that natural river networks obey either model.

## Summary for specialists

For every $c>1$, put $R=2c$. In the critical Tokunaga model, with leaves of order one, the candidate proves

$$
\begin{aligned}
\mathbb E[K\mid L=n]&=\log_R n\\
&\quad+D_c(\log_R n)+o(1),
\end{aligned}
$$

where $D_c$ is continuous and period one. The centered conditional law converges in total variation along its actual logarithmic phases, with every fixed absolute moment. No nonconstancy assertion is made for $D_c$.

For the deterministic limiting Kingman branch densities $\mathcal N_j$, it establishes

$$
\mathcal N_j\rho^j\longrightarrow H\in(0,\infty),\qquad 2\le\rho\le4.
$$

The finite-tree size limit is taken first, at each fixed order; then $j$ tends to infinity. Arbitrary joint limits are not included.

## Technical account

The Tokunaga reduction matches leaf-count marginals to a Galton–Watson process with immigration. It does not identify the original tree filtration with that auxiliary Markov chain. A uniform lattice local limit theorem, supplemented by weighted point-probability bounds, makes conditioning on an exact leaf count legitimate. The logarithmic phase survives this conditioning.

For Kingman, endpoint-normalized profiles produce a monotone integral limit. A separate fractional-linear comparison gives a supermultiplicative bound and controls the normalization. That second step is essential: convergent ratios alone do not rule out a vanishing geometric prefactor.

## Evidence, assurance and limitations

The universal conclusions rest on the written arguments. Symbolic identities and rational probabilities through 64 leaves provide finite checks. Numerical ODE calculations are diagnostics, not proof certificates or interval-certified decimals.

The supplied review exposed a roughly 0.245% drift in one original 60-order estimate despite passing finite inequalities. The revised code uses a stable residual and controlled integration window, reports quadrature warnings and errors, and separately checks agreement between 40- and 60-order runs. The analytic tail bound does not bound all ODE or roundoff errors.

The review has been actioned, but its external identity is unauthenticated. Internal checking, archived availability, formal verification and independent review remain separate assurance dimensions.

## Relationship to earlier work

The $c=2$ binary/Catalan case, including periodic register-function behavior, is classical. The manuscript translates leaf-versus-vertex and order-zero-versus-order-one conventions explicitly. Modern conditioned Galton–Watson results provide important context but use different model assumptions or observables.

The Kingman exponent and endpoint comparison also predate this work. The claimed advance is the positive geometric prefactor. The source audit is targeted, not an exhaustive guarantee of priority.

## Who should care, and why

| Audience | Potential use | Required caution |
| --- | --- | --- |
| Probability researchers | Exact-conditioning and normalization arguments | Audit the analytic proofs, not just the scripts |
| Branching and algorithm researchers | Compare order conventions and fluctuation laws | The classical binary case is not a new result |
| Mathematical modelers | More precise predictions within specified random-tree models | No empirical calibration or model validation is supplied |

## Why the problem matters

A leading exponent leaves substantial information unresolved. It may not describe fluctuations at a fixed size, and it need not determine whether a rescaled quantity has a nonzero limit. This package addresses those two gaps in distinct settings.

## How to inspect or reproduce the recorded checks

Start with `AI_INDEX.md` and the manuscript. In an isolated Python environment, install `requirements.txt`, then run the exact checker normally and with `-O`. Run the numerical checker at 40 orders and at 60 orders with step 0.04, then compare their JSON records using `verify_stability.py`. Run `negative_controls.py` in both modes: all six deliberate errors must be rejected. Full commands are in the package README.

Successful replay checks the shipped diagnostics. It does not mechanically verify the local-limit theorem or the infinite ODE hierarchy.

## The most valuable next projects

Independent scrutiny of the Fourier bounds and positive-normalization argument is the first priority. Validated ODE enclosures could certify numerical constants. Joint finite-size/order limits and empirical model assessment are separate research questions, not missing conditions silently assumed here.

## What is in the evidence package

The archive includes manuscript source and PDF, an agent-readable index, review dispositions, source comparison, pinned dependencies, exact and numerical checkers, retained revised diagnostics, negative controls and a complete checksum manifest. The original supplied review audit is retained privately rather than relicensed as part of the research output.




## Verification status

Unrefereed analytical candidate with supplied review actioned and internal finite replay. No authenticated external specialist review, formal proof or independent reproduction.

## References

1. Kovchegov and Zaliapin (2020). Random self-similar trees: A mathematical theory of Horton laws. Model definitions, inverse pruning and open questions. <https://doi.org/10.1214/19-PS331>
2. Kovchegov and Zaliapin (2017). Horton self-similarity of Kingman’s coalescent tree. Existing root-Horton and endpoint comparison framework. <https://doi.org/10.1214/16-AIHP748>
3. Flajolet, Raoult and Vuillemin (1979). The number of registers required for evaluating arithmetic expressions. Classical binary case. <https://doi.org/10.1016/0304-3975(79)90009-4>
4. Brandenberger, Devroye and Reddad (2021). The Horton–Strahler number of conditioned Galton–Watson trees. <https://doi.org/10.1214/21-EJP678>
5. Khanfir (2025). The Horton–Strahler number of Galton–Watson trees with possibly infinite variance. <https://doi.org/10.1214/25-AAP2204>
6. Khanfir (2025, v2). Fluctuations of the Horton–Strahler number of stable Galton–Watson trees. <https://arxiv.org/abs/2401.13771v2>
