REVIEW 4 major objections 3 minor 32 references
Non-symmetric vector dyson equations
T0 review · 4 major / 3 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read Non-symmetric vector Dyson equations with symmetric or non-backtracking S have fully controlled stability: inverse size (κ+η)^(-1/2) at edges, (κ+η)^(-2/3) at cusps, bounded in the bulk except at a bipartite point.
desk verdict A genuinely new stability theory for non-symmetric vector Dyson equations, but the unproved Puiseux lemma and an implicit positivity assumption need fixing before I'd trust the qualitative framework. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the stability operator I - m(z)^2 S and its equivalent phase-modulus form R = U - F, where U is a unitary diagonal of phase ratios and F = |m| S |m| is a nonnegative matrix with the same support as S. In the symmetric theory F is self-adjoint and its norm is bounded by one; here that fails. The authors replace it with a Perron-Frobenius normalization of F into a Markov transition matrix Q_F, using a quantitative mixing-space estimate for I - λQ on the complement of the Perron direction. For bulk stability they introduce the resonance set R_Q of phase matrices diag(v_i) satisfying a cycle product condition, and prove that a small singular value of U - λQ forces U near R_
What would settle it
Compute the boundary behavior of m for a small irreducible nonnegative non-symmetric S at a point where the density vanishes: if some component diverges like a power with irrational exponent, or like log η, rather than admitting a rational-power Puiseux expansion, then Lemma 3.1 fails and the decomposition theorem collapses. A more targeted check: at a regular cusp of a non-backtracking matrix on a 3-regular graph, numerically estimate the density exponent and the inverse stability norm; observing growth faster than |ω|^(1/3) or η^(-2/3) would falsify Theorems 2.7 and 2.8.
Extended reading notes
Core claim
The paper establishes that the stability operator I - m(z)^2 S, the linearization of the vector Dyson equation -1/m(z) = z1 + a + S m(z), is fully controlled for two classes of nonnegative irreducible matrices: symmetric S and non-backtracking matrices (NB). Theorem 2.8 gives two-sided bounds: at a regular edge, ||(I - m^2 S)^(-1)|| ~ max{1,(κ+η)^(-1/2)}; at a regular cusp, under (Sym) or (NB), it is at most (κ+η)^(-2/3); in the strict bulk it is bounded unless the relevant graph is bipartite and a is a constant vector, in which case the inverse grows like |z+a|^(-1). Theorem 2.7 converts these bounds into density singularities: v_i(τ±ω)=c_i ω^(1/2)+O(ω) at edges and v_i(τ+ω)=c_i |ω|^(1/3)+O
Load-bearing premise
The framework rests on the unproved premise that, near every real point, each component of the solution expands in a convergent fractional-power series (a Puiseux expansion), and that the nonzero entries of S stay uniformly bounded away from zero; if either fails, the finite singular set, measure decomposition, interval support, and edge/cusp classification collapse.
Editorial extensions
If this is right
- At regular edges the limiting density has square-root growth, and at regular cusps cubic-root growth, in both the symmetric and the non-backtracking non-symmetric settings.
- The inverse stability operator has two-sided bound (κ+η)^(-1/2) at regular edges and upper bound (κ+η)^(-2/3) at cusps, meaning perturbations are amplified only at the scale of the density singularity.
- In the strict bulk the stability operator is bounded away from singularities, except in the bipartite case with a=a1, where z=-a is the only instability point and the inverse grows like |z+a|^(-1).
- Every irreducible nonnegative S yields a density whose singular set is finite and whose support is a finite union of intervals; the measure decomposition holds without symmetry.
- In the symmetric setting the primitivity assumption is removed, giving a complete characterization of bulk stability including the bipartite exceptional case.
Reading between the lines
- The sharp stability estimates would likely provide the deterministic core needed to extend local-law and universality proofs from symmetric to non-backtracking variance profiles, a step the paper itself does not take.
- The resonance set R_S is a computable diagnostic: for a proposed non-symmetric S, checking R_S predicts where bulk instability can occur, as the cyclic permutation example already illustrates.
- The minimum-degree-three condition in (NB) is used to force phase alignment on the non-backtracking graph; relaxing it to degree two may create new resonance modes and change the cusp picture, which the paper leaves open.
- Because the Puiseux-expansion lemma is cited rather than proved, a direct algebraic proof for irreducible nonnegative S would make the qualitative framework self-contained.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the vector Dyson equation -1/m(z) = z1 - a + S m(z) with entrywise nonnegative, not necessarily symmetric S. It proves existence/uniqueness and Stieltjes representation for all z in C+ (Thm 2.1), a measure decomposition with finite singular set and finite-union-of-intervals support (Thm 2.3), and Hölder regularity under local boundedness (Thm 2.6). The main results concern the stability operator I - m(z)^2 S: for symmetric matrices (Sym) and for weighted non-backtracking matrices (NB) arising from connected base graphs of minimum degree at least 3, the paper proves square-root growth at regular edges and cubic-root growth at regular cusps (Thm 2.7), and optimal bounds on the inverse of the stability operator in edge, cusp, bulk, and off-support regimes (Thm 2.8), including a |z+a|^(-1) bulk instability in a bipartite exceptional case. The proof strategy combines a Perron-Frobenius/Markov-chain normalization of the non-symmetric matrix F, a graph-theoretic resonance set for the phase matrix U, and a reduction to a scalar cubic equation. An appendix gives a cyclic-permutation counterexample showing that arbitrary nonnegative irreducible S can lose bulk stability.
Significance. If the results are correct, this is a meaningful extension of the Ajanki–Erdős–Krüger stability theory: it covers the non-symmetric equations arising from Green's functions on universal covers of finite graphs, and it identifies possible degeneracies through graph-theoretic resonance. The explicit Perron-Frobenius normalization (Lem. 4.2, Prop. 4.4), the resonance-set criterion (Lem. 5.4), and the concrete counterexample (App. A) are valuable concrete contributions. The paper is free of fitted parameters and does not appear circular: prior results [2,5,8,11,12,29] are used in a standard way. However, several load-bearing ingredients are delegated or implicit, and they need to be supplied before the main claims can be considered fully verified.
major comments (4)
- [§3, Lemma 3.1] Lemma 3.1 asserts, without proof, a convergent Puiseux expansion for every component m_i(z) near every real tau for arbitrary nonnegative S, citing 'Proof of Proposition 4.2 in [8]'. This lemma is load-bearing: Theorem 2.3 Step 1 uses it to make A_j finite, Step 4 uses it to make the support S a finite union of intervals, and Theorem 2.6 uses it for Hölder regularity. The cited [8] is a paper on quantum ergodicity on graphs, not a source for vector Dyson equations. The manuscript should either prove Lemma 3.1 (e.g. from algebraicity of the solution) or give a precise theorem in the literature that covers exactly this setting. As written, the qualitative framework of the paper rests on an unverified assertion.
- [§3, Thm 2.5; §4, Prop. 4.4] The proof of Theorem 2.5 says 'Since positive entries of S are uniformly bounded from below', and the proof of Proposition 4.4 uses this to conclude kappa_Q_F ~ 1. This uniform lower bound is not stated in the hypotheses of Theorem 2.5 or in Assumptions (Sym)/(NB). Without it, f_ij(z) ~ s_ij/(1+|z|)^2 does not imply that positive entries of Q_F are comparable to 1, and the constants in (2.8), (2.9), (4.14), and in Theorem 2.8 need not be uniform. Please add the lower-bound assumption explicitly, or define 'model parameters' to include min_{s_ij>0} s_ij and verify it for (NB).
- [§4, Prop. 4.8] Proposition 4.8, which yields the approximate cubic equation (4.54) and error bounds (4.55), is the bridge from the stability expansion to Theorems 2.7–2.8. Its proof consists of a statement that the argument is 'essentially identical' to [2, Prop. 6.2] after replacing the orthogonal projection Q by the oblique projection P = I - r l^T. The non-symmetric setting changes the linear algebra in a non-trivial way, and the estimates (4.55a)–(4.55b) are not derived. Please include a self-contained proof, or a detailed translation that verifies every estimate in the oblique-coordinate case.
- [§4, Prop. 4.7; §5, Prop. 5.8] In the (Sym) bipartite case of Proposition 4.7, the exclusion of the alternating mode at tau = 0 uses a reflection argument with an invalid formula for the map z -> -z and then evaluates at z = i eta with '-z = z'. This step is needed for psi(tau)>0 when sigma(tau)=0, hence for the cusp part of Theorem 2.7; Proposition 5.8 also refers back to this argument to conclude m(0) is purely imaginary. Please correct the reflection map and provide a valid proof, or the bipartite exceptional cases are not supported.
minor comments (3)
- [§4, Prop. 4.7] The displayed reflection formula contains typographical errors; it should be corrected and a proper proof of the symmetry property should be given.
- [§1.3] The set of 'model parameters' is never defined. The paper repeatedly uses constants depending on them, and the implicit lower bound on positive entries of S is part of this dependence; a brief definition would remove ambiguity.
- [Thm 2.8] In the (NB) case, the phrase 'a in R1' should be explained: a is indexed by directed edges, so the condition means a_(x,y) = a_y = a constant independent of y.
Circularity Check
No significant circularity: the derivation is self-contained, uses no self-citations, and no prediction reduces to a fitted input.
full rationale
The paper contains no fitted parameters renamed as predictions and no self-citations: all cited prior work is by other authors ([2,5] Ajanki–Erdős–Krüger, [8] Anantharaman–Sabri, [11] Avni–Breuer–Simon, [12] Banks–Garza-Vargas–Mukherjee, [29] Nagnibeda–Woess). Theorem 2.1 is quoted from [5, Section 4], an external existence theorem. Lemma 3.1, the Puiseux-expansion step, is cited to [8, Prop. 4.2]; that reference is not authored or co-authored by the present authors, so this is an ordinary external-citation dependency, not a self-citation chain. The Perron–Frobenius normalization of F to a Markov matrix Q_F (Eq. 4.15) is a change of variables, not an assumption of the conclusion. The resonance set R_Q (Definition 5.2) is defined from the support graph of S and then used to analyze degeneracies—a structural analysis rather than an input equivalent to the target result. The edge/cusp singularities are derived from a scalar cubic equation (Proposition 4.8) whose coefficients μ, σ, ψ come from the solution and Perron eigenvectors, with nondegeneracy ψ+σ²∼1 proved under (Sym)/(NB); no input equivalent to the asserted square-root or cubic-root behavior is used. The bipartite exceptional instability is proved by explicit spectral computation in Proposition 5.8 and Appendix A. The concerns raised in the skeptical summary—Lemma 3.1 is not proved in the text and may not be covered by the cited [8] result, and Theorem 2.5/Proposition 4.4 silently use a uniform lower bound on positive entries of S—are omitted-proof/assumption gaps and correctness risks, not circularity: they do not make any asserted conclusion identical to an input by definition. Therefore the appropriate circularity score is 0.
Assumptions & free parameters
assumptions (5)
- domain assumption Each component m_i(z) has a convergent Puiseux expansion near every real point (Lemma 3.1).
- standard math Perron-Frobenius theorem for irreducible nonnegative matrices and the cyclic decomposition (Lemmas 1.1 and 1.2).
- domain assumption S is nonnegative and irreducible (stated before Theorem 2.3); reducible case is claimed to be triangularizable.
- domain assumption Assumptions (Sym) and (NB): S symmetric, or S is a non-backtracking matrix of a connected base graph with minimum degree at least 3 and a terminal-dependent vector a.
- domain assumption Positive entries of S are uniformly bounded below; all constants 'depending only on model parameters' hide a positivity parameter.
Cite this review
Pith. "Pith review of Non-symmetric vector dyson equations." pith.science (2026). https://pith.science/paper/KSQSTUC3
@misc{pith2026260716333,
author = {Pith},
title = {Pith review of: Non-symmetric vector dyson equations},
year = {2026},
howpublished = {\url{https://pith.science/paper/KSQSTUC3}},
note = {Machine review of arXiv:2607.16333}
}
abstract
We study the vector Dyson equation $$-\frac{1}{m(z)}=z\mathbf{1}+\mathbf{a}+Sm(z),$$ with parameter $z$ in the complex upper half-plane $\mathbb{C}_+$, where $\mathbf{a}\in\mathbb R^d$ and $S$ is a nonnegative matrix, not necessarily symmetric. This equation has a unique vector solution $m(z)\in\mathbb{C}_+^d$, for which we establish a complete measure decomposition and prove regularity. We then develop a graph-theoretic approach to the singularity and stability problem for non-symmetric matrices $S$. The graph structure of $S$ identifies the possible degeneracies of the stability operator as $z$ approaches the real axis. In particular, for non-backtracking matrices, we prove square-root growth at regular edges, cubic-root growth at regular cusps, and complete stability estimates. We also obtain the corresponding estimates for symmetric matrices in the periodic setting.
Reference graph
Works this paper leans on
-
[8]
Recent results of quantum ergodicity on graphs and further investigation.Annales de la Faculté des Sciences de Toulouse
Nalini Anantharaman and Mostafa Sabri. Recent results of quantum ergodicity on graphs and further investigation.Annales de la Faculté des Sciences de Toulouse. Mathématiques. Série 6, 28(3):559–592, 2019. 49
2019
-
[1]
Ajanki, László Erdős, and Torben Krüger
Oskari H. Ajanki, László Erdős, and Torben Krüger. Local spectral statistics of Gaussian matrices with correlated entries.Journal of Statistical Physics, 163(2):280–302, 2016
2016
-
[2]
Ajanki, László Erdős, and Torben Krüger
Oskari H. Ajanki, László Erdős, and Torben Krüger. Singularities of solutions to quadratic vector equations on the complex upper half-plane.Communications on Pure and Applied Mathematics, 70(9):1672–1705, 2017
2017
-
[3]
Ajanki, László Erdős, and Torben Krüger
Oskari H. Ajanki, László Erdős, and Torben Krüger. Universality for general Wigner-type matrices.Probability Theory and Related Fields, 169(3-4):667–727, 2017
2017
-
[4]
Stability of the matrix dyson equation and random matrices with correlations.Probability Theory and Related Fields, 173(1):293– 373, 2019
Oskari H Ajanki, László Erdős, and Torben Krüger. Stability of the matrix dyson equation and random matrices with correlations.Probability Theory and Related Fields, 173(1):293– 373, 2019
2019
-
[5]
Quadratic vector equations on com- plex upper half-plane.Memoirs of the American Mathematical Society, 261(1261):v+133, 2019
Oskari Heikki Ajanki, László Erdős, and Torben Krüger. Quadratic vector equations on com- plex upper half-plane.Memoirs of the American Mathematical Society, 261(1261):v+133, 2019
2019
-
[6]
Local law for random Gram matrices
Johannes Alt, László Erdős, and Torben Krüger. Local law for random Gram matrices. Electronic Journal of Probability, 22(none):1 – 41, 2017
2017
-
[7]
Correlated random matrices: Band rigidity and edge universality.The Annals of Probability, 48(2):963 – 1001, 2020
Johannes Alt, László Erdős, Torben Krüger, and Dominik Schröder. Correlated random matrices: Band rigidity and edge universality.The Annals of Probability, 48(2):963 – 1001, 2020
2020
Show all 32 references
-
[9]
Anderson and Ofer Zeitouni
Greg W. Anderson and Ofer Zeitouni. A CLT for a band matrix model.Probability Theory and Related Fields, 134(2):283–338, 2006
2006
-
[10]
Anderson and Ofer Zeitouni
Greg W. Anderson and Ofer Zeitouni. A law of large numbers for finite-range dependent random matrices.Communications on Pure and Applied Mathematics, 61(8):1118–1154, 2008
2008
-
[11]
Periodic jacobi matrices on trees.Advances in Mathematics, 370:107241, 2020
Nir Avni, Jonathan Breuer, and Barry Simon. Periodic jacobi matrices on trees.Advances in Mathematics, 370:107241, 2020
2020
-
[12]
Point spectrum of peri- odic operators on universal covering trees.International Mathematics Research Notices, 2022(22):17713–17744, 2022
Jess Banks, Jorge Garza-Vargas, and Satyaki Mukherjee. Point spectrum of peri- odic operators on universal covering trees.International Mathematics Research Notices, 2022(22):17713–17744, 2022
2022
-
[13]
F. A. Berezin. Some remarks on Wigner distribution.Theoretical and Mathematical Physics, 17(3):1163–1171, 1973
1973
-
[14]
Plemmons.Nonnegative Matrices in the Mathematical Sciences, volume 9 ofClassics in Applied Mathematics
Abraham Berman and Robert J. Plemmons.Nonnegative Matrices in the Mathematical Sciences, volume 9 ofClassics in Applied Mathematics. Society for Industrial and Applied Mathematics, Philadelphia, PA, 1994
1994
-
[15]
Brualdi and Herbert J
Richard A. Brualdi and Herbert J. Ryser.Combinatorial Matrix Theory, volume 39 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 1991
1991
-
[16]
Cusp universality for random matrices, ii: The real symmetric case.Pure and Applied Analysis, 1(4):615–707, 2019
Giorgio Cipolloni, László Erdős, Torben Krüger, and Dominik Schröder. Cusp universality for random matrices, ii: The real symmetric case.Pure and Applied Analysis, 1(4):615–707, 2019
2019
-
[17]
The matrix Dyson equation and its applications for random matrices
László Erdős. The matrix Dyson equation and its applications for random matrices. In Alexei Borodin, Ivan Corwin, and Alice Guionnet, editors,Random Matrices, volume 26 of IAS/Park City Mathematics Series, pages 75–158. American Mathematical Society, Provi- dence, RI, 2019
2019
-
[18]
Random matrices with slow corre- lation decay
László Erdős, Torben Krüger, and Dominik Schröder. Random matrices with slow corre- lation decay. InForum of Mathematics, Sigma, volume 7, page e8. Cambridge University Press, 2019
2019
-
[19]
Cusp universality for random matrices i: local law and the complex hermitian case.Communications in Mathematical Physics, 378(2):1203–1278, 2020
László Erdős, Torben Krüger, and Dominik Schröder. Cusp universality for random matrices i: local law and the complex hermitian case.Communications in Mathematical Physics, 378(2):1203–1278, 2020
2020
-
[20]
Girko.Theory of Stochastic Canonical Equations
Vyacheslav L. Girko.Theory of Stochastic Canonical Equations. Vol. I, volume 535 of Mathematics and its Applications. Kluwer Academic Publishers, Dordrecht, 2001
2001
-
[21]
Large deviations upper bounds and central limit theorems for non- commutative functionals of Gaussian large random matrices.Annales de l’Institut Henri Poincaré
Alice Guionnet. Large deviations upper bounds and central limit theorems for non- commutative functionals of Gaussian large random matrices.Annales de l’Institut Henri Poincaré. Probabilités et Statistiques, 38(3):341–384, 2002
2002
-
[22]
Hachem, P
W. Hachem, P. Loubaton, and J. Najim. The empirical distribution of the eigenvalues of a gram matrix with a given variance profile.Annales de l’Institut Henri Poincare (B) Probability and Statistics, 42(6):649–670, 2006
2006
-
[23]
William Helton, Reza Rashidi Far, and Roland Speicher
J. William Helton, Reza Rashidi Far, and Roland Speicher. Operator-valued semicircular elements: solving a quadratic matrix equation with positivity constraints.International Mathematics Research Notices, 2007:rnm086, 2007. 50
2007
-
[24]
Horn and Charles R
Roger A. Horn and Charles R. Johnson.Matrix Analysis. Cambridge University Press, Cambridge, 2 edition, 2013
2013
-
[25]
On the spectral theory of trees with finite cone type.Israel Journal of Mathematics, 194(1):107–135, 2013
Matthias Keller, Daniel Lenz, and Simone Warzel. On the spectral theory of trees with finite cone type.Israel Journal of Mathematics, 194(1):107–135, 2013
2013
-
[26]
An invitation to trees of finite cone type: random and deterministic operators.Markov Processes and Related Fields, 21(3):557–574, 2015
Matthias Keller, Daniel Lenz, and Simone Warzel. An invitation to trees of finite cone type: random and deterministic operators.Markov Processes and Related Fields, 21(3):557–574, 2015
2015
-
[27]
A. M. Khorunzhy and L. A. Pastur. On the eigenvalue distribution of the deformed Wigner ensemble of random matrices. InSpectral Operator Theory and Related Topics, volume 19 ofAdvances in Soviet Mathematics, pages 97–127. American Mathematical Society, Provi- dence, RI, 1994
1994
-
[28]
The Dyson equation for2-positive maps and Hölder bounds for the Lévy distance of densities of states, 2022
Tobias Mai. The Dyson equation for2-positive maps and Hölder bounds for the Lévy distance of densities of states, 2022
2022
-
[29]
Random walks on trees with finitely many cone types.J
Tatiana Nagnibeda and Wolfgang Woess. Random walks on trees with finitely many cone types.J. Theoret. Probab., 15(2):383–422, 2002
2002
-
[30]
Absolutely continuous spectrum for random Schrödinger operators on tree- strips of finite cone type.Annales Henri Poincaré, 14(4):737–773, 2013
Christian Sadel. Absolutely continuous spectrum for random Schrödinger operators on tree- strips of finite cone type.Annales Henri Poincaré, 14(4):737–773, 2013
2013
-
[31]
Random Gaussian band matrices and freeness with amalgamation
Dimitri Shlyakhtenko. Random Gaussian band matrices and freeness with amalgamation. International Mathematics Research Notices, 1996(20):1013–1025, 1996
1996
-
[32]
Franz J. Wegner. Disordered system withnorbitals per site:n=∞limit.Physical Review B, 19(2):783–792, 1979. 51
1979
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.