A Combinatorial Framework for the Pons-Batle Identity: Young Tableaux, Lattice Paths, and Limit Laws

TL;DR

This paper introduces a combinatorial and analytic framework using Young tableaux and lattice paths to confirm a conjecture about counting tree-child networks, revealing their asymptotic probabilistic structures.

Contribution

It provides explicit bijections, generating functions, recurrence relations, and limit laws for the enumeration of tree-child networks with bounded reticulation nodes.

Findings

Confirmed the Pons-Batle conjecture for bounded reticulation nodes.

Derived explicit formulas and recurrences for counting networks.

Identified Beta and Uniform distributions as limit laws for structural parameters.

Abstract

Tree-child networks are an important class of phylogenetic network used to model reticulate evolutionary processes. These networks have attracted increasing attention from researchers with interests in both combinatorics and algorithms. A fundamental open problem posed by Pons and Batle asks whether the number $TC_{n,k}$ of bicombining tree-child networks with $n$ leaves and $k$ reticulation nodes equals the number of certain constrained words, now called Pons-Batle words. In this paper, we confirm the conjecture for tree-child networks with a bounded number of reticulation nodes. Our approach is combinatorial and analytic. We introduce families of Young tableaux with walls and holes and construct explicit bijections with Pons-Batle words, yielding a direct combinatorial explanation of the identities. These tableaux encode structural features of the underlying networks, including the placement of reticulation nodes. By projecting them to decorated Dyck paths, we obtain algebraic generating functions with differential operators encoding step weights, leading to explicit recurrence relations and closed-form formulas for $TC_{n,k}$. Beyond finite verification for moderate $k$, the framework reveals an underlying probabilistic structure. For $k=1$, natural structural parameters, such as the position and value of distinguished cells, converge, after rescaling, to $\mathrm{Beta}(2,1)$, $\mathrm{Beta}(1,2)$, and Uniform (i.e., $\mathrm{Beta}(1,1)$) distributions. These limit laws arise from a coalescence of singularities at the dominant square-root singularity, producing a non-analytic transition in the local expansion. Overall, our results provide both combinatorial insight and a unified analytic perspective on the asymptotic behavior of tree-child networks, showing how algebraic generating functions with interacting singularities systematically produce Beta limit laws.