A note on the spectral distribution of non-Hermitian block matrices with Toeplitz blocks
Pith reviewed 2026-05-13 16:51 UTC · model grok-4.3
The pith
The limiting spectral distribution of non-Hermitian block matrices with Toeplitz blocks is given by the geometric mean of the associated GLT symbols.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For non-Hermitian block matrix sequences whose blocks are Toeplitz, the empirical spectral distribution converges to the distribution induced by the geometric mean of the matrix-valued GLT symbol of the sequence.
What carries the argument
The geometric mean of matrices applied to the GLT symbol of the block matrix sequence.
Load-bearing premise
The block matrix sequences admit a GLT symbol whose geometric mean correctly describes the limiting spectral distribution for the non-Hermitian case.
What would settle it
A numerical experiment on a sequence of large non-Hermitian block Toeplitz matrices in which the empirical spectral measure deviates from the measure predicted by the geometric mean of the GLT symbol.
Figures
read the original abstract
In the present paper, we are concerned with the study of the spectral distribution of matrix-sequences showing a non-Hermitian block structure with Toeplitz blocks. We use the notion of geometric mean of matrices and the theory of Generalized Locally Toeplitz (GLT) sequences to perform our analysis and produce some numerical tests and visualizations to confirm our theoretical derivations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the limiting spectral distribution of sequences of non-Hermitian block matrices with Toeplitz blocks. It combines the geometric mean of matrices with Generalized Locally Toeplitz (GLT) theory to derive an effective symbol for the block structure and supports the claim with numerical tests and visualizations.
Significance. If the geometric-mean construction can be rigorously justified for non-Hermitian symbols, the result would extend GLT techniques to a new class of block-structured non-Hermitian matrices, offering a compact description of their eigenvalue distributions with possible utility in numerical linear algebra.
major comments (1)
- The central construction relies on the geometric mean of the GLT symbols of the individual Toeplitz blocks to obtain the limiting eigenvalue distribution of the non-Hermitian block matrix sequence. The classical matrix geometric mean A#B is defined via the principal square root and requires both arguments to be positive definite Hermitian; no alternative definition (e.g., via holomorphic functional calculus on the joint spectrum or via the GLT symbol directly) is supplied for general complex-valued symbols. Because the blocks are arbitrary Toeplitz (hence their symbols need not be Hermitian or positive), the step that equates the geometric mean of the symbols to the effective GLT symbol of the block matrix is formally unsupported. This is load-bearing: without it, the claimed spectral distribution formula does not follow from standard GLT theory.
minor comments (1)
- The abstract states that numerical tests 'confirm our theoretical derivations' but provides no details on matrix dimensions, sampling of the symbol, or quantitative error measures used in the visualizations.
Simulated Author's Rebuttal
We thank the referee for the thorough review and for identifying the need to clarify the definition of the geometric mean in our non-Hermitian setting. We provide a point-by-point response below and will update the manuscript to include the missing definition.
read point-by-point responses
-
Referee: The central construction relies on the geometric mean of the GLT symbols of the individual Toeplitz blocks to obtain the limiting eigenvalue distribution of the non-Hermitian block matrix sequence. The classical matrix geometric mean A#B is defined via the principal square root and requires both arguments to be positive definite Hermitian; no alternative definition (e.g., via holomorphic functional calculus on the joint spectrum or via the GLT symbol directly) is supplied for general complex-valued symbols. Because the blocks are arbitrary Toeplitz (hence their symbols need not be Hermitian or positive), the step that equates the geometric mean of the symbols to the effective GLT symbol of the block matrix is formally unsupported. This is load-bearing: without it, the claimed spectral distribution formula does not follow from standard GLT theory.
Authors: We agree with the referee that the classical definition of the matrix geometric mean requires positive definite Hermitian matrices and that our manuscript does not explicitly provide an alternative definition for general complex symbols. This is indeed a gap in the current presentation. In the revised version of the paper, we will introduce a section on matrix functions where we define the geometric mean of two matrices A and B (with appropriate spectra) as exp( (log A + log B)/2 ), using the principal branch of the matrix logarithm via holomorphic functional calculus. This definition extends naturally to the pointwise application on the GLT symbols, which are continuous matrix-valued functions. We will also discuss the conditions under which this is well-defined for the symbols of Toeplitz blocks. With this addition, the application of GLT theory to derive the spectral distribution follows as stated. We thank the referee for this suggestion, which will improve the rigor of the work. revision: yes
Circularity Check
No circularity: derivation applies external GLT theory and matrix geometric mean without self-referential reduction
full rationale
The paper's abstract states that the analysis of the spectral distribution for non-Hermitian block matrices with Toeplitz blocks relies on the geometric mean of matrices combined with Generalized Locally Toeplitz (GLT) sequence theory, followed by numerical confirmation. No equations or steps are exhibited that reduce a claimed prediction or limiting distribution to a fitted parameter or self-defined quantity by construction. The central construction invokes prior GLT results as an external framework rather than deriving them internally, and the geometric mean is presented as a standard tool applied to the symbols without evidence of ansatz smuggling or uniqueness imported via self-citation chains. The derivation chain remains self-contained against external benchmarks, with numerical tests serving as independent verification rather than tautological confirmation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption GLT theory extends to non-Hermitian block matrices with Toeplitz blocks
Reference graph
Works this paper leans on
-
[1]
A. Adriani, I. Furci, C. Garoni, S. Serra-Capizzano,Spectral and singular value distribution of se- quences of block matrices with rectangular Toeplitz blocks. Part I: Asymptotically rational block size ratios, Journal of Numerical Mathematics (2026)
work page 2026
-
[2]
A. Adriani, C. Garoni, S. Serra-Capizzano,Spectral and singular value distribution of sequences of block matrices with rectangular Toeplitz blocks. Part II: Asymptotically irrational block size ratios, Journal of Numerical Mathematics (2026)
work page 2026
-
[3]
T. Ando, C.-K. Li, R. Mathias,Geometric means, Linear Algebra Appl., 385 (2004), pp. 305–334
work page 2004
-
[4]
G. Barbarino, C. Garoni, S. Serra-Capizzano,Block generalized locally Toeplitz sequences: theory and applications in the unidimensional case, Electron. Trans. Numer. Anal. 53 (2020), pp. 28–112
work page 2020
-
[5]
G. Barbarino, C. Garoni, S. Serra-Capizzano,Block generalized locally Toeplitz sequences: theory and applications in the multidimensional case, Electron. Trans. Numer. Anal. 53 (2020), pp. 113–216
work page 2020
-
[6]
G. Barbarino, C. Garoni, M. Mazza, S. Serra-Capizzano,Rectangular GLT sequences, Electron. Trans. Numer. Anal., 55 (2022)
work page 2022
-
[7]
G. Barbarino, S. Serra-Capizzano,Non-Hermitian perturbation of Hermitian matrix-sequences and applications to the spectral analysis of numerical approximation of partial differential equations, Nu- merical Linear Algebra with Applications, 27 (3) (2020)
work page 2020
-
[8]
P. Benedusi, A.J. Ellingsrud, H. Herlyng, M.E. Rognes,Scalable approximation and solvers for ionic electrodiffusion in cellular geometries, SIAM J. on Scientific Computing, 46(5) (2024)
work page 2024
-
[9]
A. B¨ ottcher , B. Silbermann,Introduction to Large Truncated Toeplitz Matrices, Springer-Verlag, New York (1999)
work page 1999
-
[10]
D. Bertaccini, M. Donatelli, F. Durastante, S. Serra-Capizzano,Optimizing a multigrid Runge–Kutta smoother for variable-coefficient convection-diffusion equations, Linear Algebra Appl., 533 (2017), pp. 507-535
work page 2017
-
[11]
C. Garoni,Introduction to the theory of generalized locally Toeplitz sequences and its applications, Mathematical Notes (2026)
work page 2026
- [12]
- [13]
- [14]
-
[15]
A. Ilyas, M.F. Khan, V. Loi, S. Serra-Capizzano,A maximal distribution result and numerical tests for geometric means of HPD GLT matrix-sequences with degenerate commuting/non-commuting GLT symbols, Linear Algebra Appl. (2026). https://doi.org/10.1016/j.laa.2025.08.004
-
[16]
A. Ilyas, M.F. Khan, R.L. Sormani, G. Tento, S. Serra-Capizzano,Determining the space dependent coefficients in space-time fractional diffusion equations via Krylov preconditioning, Linear Algebra Appl. (2026). https://doi.org/10.1016/j.laa.2026.02.024
-
[17]
F. Kubo, T. Ando,Means of positive linear operators, Math. Ann., 246(3) (1980), pp. 205–224. 18
work page 1980
-
[18]
M. Moakher,A differential geometric approach to the geometric mean of symmetric positive-definite matrices, SIAM J. Matrix Anal. Appl., 26(3) (2005), pp. 735–747
work page 2005
-
[19]
Nakamura,Geometric means of positive operators, Kyungpook Math
N. Nakamura,Geometric means of positive operators, Kyungpook Math. J., 49(1) (2009), pp. 167– 181
work page 2009
-
[20]
S. Serra-Capizzano,Spectral theory of matrix-sequences: perspectives of the GLT analysis and beyond, arXiv:2509.19544 (2025)
-
[21]
S. Serra-Capizzano,Generalized locally Toeplitz sequences: Spectral analysis and applications to discretized partial differential equations, Linear Algebra Appl. (2003)
work page 2003
-
[22]
P. Tilli,Locally Toeplitz sequences: Spectral properties and applications, Linear Algebra Appl., 278 (1998), pp. 91–120
work page 1998
- [23]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.