- The paper proves that the largest Ricci-matrix eigenvalue converges under repeated leaf attachment to the maximum of zero and the largest Dirichlet eigenvalue of the tree’s branches.
- It reduces the problem to a finite-dimensional perturbation, showing that the first-order coefficient α determines whether eigenvalues eventually increase toward or decrease toward the limit.
- A numerical example demonstrates that convergence need not be globally monotone, with eigenvalues crossing from negative to positive before eventually declining toward zero.
Background and motivation
For a finite tree T=(V,E), the Ricci matrix RT∈RE×E introduced by Bai, Cheng, and Hua encodes the Lin–Lu–Yau optimal-transport curvature of weighted trees: a discrete Einstein metric (edge weights with constant curvature) exists on T if and only if the Perron eigenvector of RT is positive, and the Einstein curvature equals κ=−λmax(RT) (Bai et al., 24 Apr 2026). The sign of λmax(RT) therefore determines the sign of the curvature. Prior work established that λmax(RT)≤0 forces T to be a caterpillar, that attaching a leaf at a vertex of degree at most two does not decrease λmax, and that when λmax(RT)<0, leaf attachment at any vertex strictly increases it; Cheng subsequently classified all trees with nonpositive Perron eigenvalue and confirmed that the only infinite zero-curvature family is RT∈RE×E0 (Cheng, 20 May 2026). These results motivate the question addressed here: under repeated attachment of pendant edges at a fixed vertex RT∈RE×E1, does RT∈RE×E2 converge, and is it eventually monotone?
The paper answers both questions affirmatively in a precise sense. The main structural observation is that the Perron eigenvector of RT∈RE×E3 must be constant on the RT∈RE×E4 new pendant edges (they form an orbit under automorphisms commuting with RT∈RE×E5), so RT∈RE×E6 reduces to the largest eigenvalue of a fixed-dimensional matrix RT∈RE×E7 admitting the affine decomposition
RT∈RE×E8
where RT∈RE×E9 and T0 is independent of T1. All subsequent results flow from this representation.
The paper first develops a variational analysis of a single leaf addition using the vertex decomposition of the quadratic form,
T2
with T3 and T4. Extending T5 by a value T6 on the new edge yields the exact difference formula
T7
Maximizing over T8 produces the T9-sharp criterion: for a unit vector RT0 with RT1 and RT2, if
RT3
then some extension attains Rayleigh quotient at least RT4, whence RT5. Applying this to the Perron vector gives a sufficient condition for one-step monotonicity. Since RT6, the worst case RT7 yields a degree-only threshold RT8 for RT9 (and κ=−λmax(RT)0 for κ=−λmax(RT)1): whenever κ=−λmax(RT)2, attaching a leaf at κ=−λmax(RT)3 cannot decrease the Perron eigenvalue. This recovers coarse degree-dependent information while retaining the sharper κ=−λmax(RT)4-based test.
Removing κ=−λmax(RT)5 splits κ=−λmax(RT)6 into branches κ=−λmax(RT)7. On the orbit-reduced subspace, each diagonal block κ=−λmax(RT)8 of κ=−λmax(RT)9 is obtained from λmax(RT)0 by setting λmax(RT)1 to zero—equivalently, imposing a Dirichlet boundary condition at λmax(RT)2 or taking λmax(RT)3. In particular, the root-edge diagonal entry becomes λmax(RT)4 rather than λmax(RT)5. The interface vectors satisfy λmax(RT)6 (supported only on the root edge), λmax(RT)7, and λmax(RT)8.
The resulting limit matrix λmax(RT)9 is block upper-triangular with diagonal blocks λmax(RT)≤00 and a scalar zero block, so its spectrum is λmax(RT)≤01, giving the limit formula
λmax(RT)≤02
Two consequences deserve emphasis. First, the growing leaf cluster contributes a floor value of exactly λmax(RT)≤03: since the symmetric mode of the new leaves has Rayleigh quotient λmax(RT)≤04, repeated leaf addition can never drive the limit negative even if every branch is subcritical (λmax(RT)≤05). Second, λmax(RT)≤06 occurs precisely when some branch, viewed as a Schrödinger operator λmax(RT)≤07 on its line graph with absorbing boundary at λmax(RT)≤08, has a potential well deep enough to bind a positive-energy state. This dichotomy is a feature not shared by adjacency matrices (whose spectral radius strictly increases under any pendant attachment) or Laplacians, where no analogous convergence-plus-sign-dependent-monotonicity phenomenon arises.
First-order expansion and tail monotonicity
Writing λmax(RT)≤09 and T0, analytic perturbation theory applies when T1 is simple. With right/left eigenvectors T2 normalized by T3, the first-order coefficient is the spectral projection
T4
and
T5
The tail monotonicity theorem then states: if T6, the sequence T7 is eventually strictly monotonic—strictly decreasing toward T8 from above when T9, strictly increasing toward it from below when λmax0. The proof is elementary once the expansion is available: λmax1 retains the sign of λmax2 near λmax3, and λmax4 decreases monotonically to λmax5.
The simplicity assumption is handled explicitly: if λmax6 is multiple (e.g., two branches attain equal maxima, or coincide with the zero block), degenerate perturbation theory replaces λmax7 by λmax8, where λmax9 is the compression of λmax(RT)<00 to the eigenspace of λmax(RT)<01; the tail monotonicity conclusion is unchanged. The derivation for this case is deferred to Kato's classical treatment rather than carried out in full, which is the paper's principal expository gap.
Numerical verification
The worked example uses a base tree with λmax(RT)<02, one branch a path of length two and the other a binary fork. Both branch blocks are subcritical, λmax(RT)<03 and λmax(RT)<04, so λmax(RT)<05. Computing left and right eigenvectors for the zero eigenvalue gives λmax(RT)<06, predicting
λmax(RT)<07
with eventual strict decrease toward λmax(RT)<08 from above. The numerics confirm this and, notably, show that early behavior is genuinely non-monotonic: λmax(RT)<09, RT∈RE×E00, RT∈RE×E01, peaking near RT∈RE×E02 before decreasing through RT∈RE×E03. Thus the sequence crosses from negative to positive before settling into its predicted tail—a concrete demonstration that eventual monotonicity need not be global monotonicity. The complementary case RT∈RE×E04 is realized by the star RT∈RE×E05, where RT∈RE×E06 exactly, so RT∈RE×E07 and the sequence increases toward RT∈RE×E08 from below, matching the theorem.
Limitations and open questions
The theory rests on the affine structure RT∈RE×E09, which holds because only the degree of RT∈RE×E10 varies; it does not address simultaneous leaf addition at multiple vertices, where the perturbation parameter would couple several degrees. The degenerate-eigenvalue case is asserted by appeal to standard perturbation theory rather than proven in detail, and the paper does not characterize when RT∈RE×E11—in that case the first-order term vanishes and eventual monotonicity is undetermined by the present methods, requiring second-order analysis. Finally, the connection between the sign of RT∈RE×E12 and combinatorial features of the branch data (beyond the computed examples) remains unexplored.
Conclusion
This paper reduces the dynamics of RT∈RE×E13 under iterated leaf attachment to a finite-dimensional rank-one spectral perturbation problem. The limit RT∈RE×E14 depends only on Dirichlet-decoupled branch data and a universal zero floor contributed by the growing leaf cluster, while the sign of the projection coefficient RT∈RE×E15 governs both the approach direction and eventual strict monotonicity. Together with the earlier classification results, this yields a complete asymptotic picture of how local leaf addition shapes the spectrum governing discrete Einstein curvature on trees.