REVIEW 4 major objections 7 minor 10 references
Universal estimates for the density of states for aperiodic block subwavelength resonator systems
T0 review · 4 major / 7 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper proves that for one-dimensional random block subwavelength resonator arrays, the integrated density of states converges almost surely to a non-random continuous function, and that the fractal-like hybridisation-region spectrum…
desk verdict A useful, clearly written paper that proves an ergodic limit for the discrete capacitance model and introduces an effective meta-atom algorithm; the gap to the physical resonances is unproven but probably fixable. 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 load-bearing object is the infinite Jacobi operator $J=V^{1/2}CV^{1/2}$ obtained as the $N\to\infty$ limit of the symmetrised generalised capacitance matrix of the resonator chain, with off-diagonal bands $s(i)=v_{i-1}v_i s_{i-1}^{-1}(\ell_{i-1}\ell_i)^{-1/2}$ and diagonal entries $q(i)=v_i^2\ell_i^{-1}(s_{i-1}^{-1}+s_i^{-1})$. Because i.i.d. block sampling makes the resonator sequence a bi-infinite Markov chain, the shift group acts metrically transitively on $J$, yielding ergodicity of the spectrum and convergence of finite-size integrated densities of states. The second mechanism is the propagation-matrix formalism: each block has a $2\times2$ transfer matrix, and frequencies with $|\operatorname{tr}P_{B_d}(\lambda)|>2$ lie in a bandgap for that block, so in a hybridisation region at least one block type is gapped and the corresponding eigenmodes decay exponentially. That decay justifies replacing the array by a catalogue of finite meta-atoms—local sequences beginning and ending with the active block—whose defect eigenfrequencies, precomputed from small capacitance matrices, reproduce the hybridisation-region density of states.
What would settle it
Solve the full Helmholtz resonance problem (2.2) at a fixed small contrast $\delta>0$ for random block chains of increasing length $M$, and compare the empirical integrated density of states with the capacitance-matrix limit $N(J,\lambda)$; if the Wasserstein distance does not tend to zero as $M\to\infty$ and $\delta\to0$, the deterministic density is an artefact of the discrete approximation. A sharper test: for $\delta$ small enough that the $O(\delta)$ eigenvalue error is below the peak spacing, the exact subwavelength resonances should reproduce the same meta-atom peak positions in the hybridisation region.
Extended reading notes
Core claim
On its own terms, the central discovery is that the integrated density of states of an i.i.d. block-disordered one-dimensional resonator system converges almost surely as the number of blocks $M\to\infty$ to a non-random measure $N(J,\mathrm{d}\lambda)$, and the distribution function $\lambda\mapsto N(J,\lambda)$ is continuous: Theorem 3.9, imported from metric-transitivity theory. The spectrum is partitioned by the pass bands of the constituent blocks: frequencies gapped for every block carry no states; frequencies passed by every block form a smooth band; and frequencies passed by some but not all blocks form hybridisation regions with non-zero density concentrated on self-similar peaks that are only weakly smoothed as the system grows. The paper further claims that the peaks are the defect modes of finite local block arrangements, called meta-atoms, so the density of states can be predicted by enumerating meta-atoms, computing their defect eigenvalues once, and scanning the block sequence in linear time. The same meta-atom procedure is demonstrated for bound-length, hyperuniform chunk and softmax, and Fibonacci quasiperiodic sampling, where it typically performs as well as or better than for i.i.d. sampling.
Load-bearing premise
All theorems are proved for the discrete capacitance matrix, whose eigenvalues match the true subwavelength resonances only up to an $O(\delta)$ error, and the paper does not show this error is uniform as the number of blocks tends to infinity.
Editorial extensions
If this is right
- For any single realization of a large random block chain, the empirical cumulative eigenvalue count approaches one fixed continuous non-random curve, so one large finite sample is a statistically representative proxy for the infinite system.
- The density of states is zero on the intersection of all constituent blocks' bandgaps and positive in both shared pass bands and hybridisation regions, which classifies the observable spectrum using only the block propagation matrices.
- The apparent fractal jumps in the empirical cumulative density are finite-size signatures: the infinite-limit integrated density is continuous, so the jagged structure is progressively smoothed as $M$ grows.
- The meta-atom algorithm reconstructs the hybridisation-region cumulative density in $O(M)$ time (or $O(ML)$ when the catalogue is scaled with length), making the deterministic limit computable for very long chains.
- The same linear-time meta-atom estimation works, and often converges faster, for bound-length, hyperuniform chunk and softmax, and Fibonacci quasiperiodic sampling, provided the meta-atom catalogue is tailored to the sampling rule.
Reading between the lines
- A direct consequence of Theorem 2.1 is an open uniformity question the paper leaves implicit: the capacitance eigenvalues match the physical resonances only up to $O(\delta)$, and if that error is not uniform in $M$, the proven deterministic density of states describes the discrete tight-binding model rather than the Helmholtz resonators.
- The decay mechanism is generic: in any one-dimensional aperiodic wave system where one constituent is opaque in a frequency band, the global density of states should be assemblable from a weighted catalogue of finite defect patterns, so the meta-atom idea could be tried on layered dielectrics, phononic chains, or non-Hermitian arrays.
- The Section 4.2 caveat that the meta-atom estimator ignores edge effects means its advertised linear-time accuracy holds strictly away from boundaries; the authors argue this is acceptable because small arrays can be diagonalised directly, but it sets a precise domain of validity for the estimator.
- A testable refinement suggested by the sampling comparison is to adapt the meta-atom set to the sampling rule, discarding impossible patterns and weighting by occurrence probability; this should remove the accuracy reversal observed for Fibonacci sequences at large meta-atom length.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the density of states (DoS) of one-dimensional aperiodic block subwavelength resonator systems. It introduces a discrete capacitance model, forms the infinite-volume Jacobi operator, and uses Pastur's theory of metrically transitive operators to prove that, for i.i.d. random block sequences, the integrated density of states converges to a non-random, continuous function as the number of blocks goes to infinity. The authors also propose a tripartite spectral decomposition into bandgaps, shared pass bands, and hybridisation regions, and they develop a fast 'meta-atom' algorithm to estimate the DoS in the hybridisation region. Numerical experiments are presented for i.i.d., bound-length, hyperuniform, and quasiperiodic samplings.
Significance. If the claims hold, the paper gives a rigorous ergodic-theoretic foundation for the DoS of disordered subwavelength resonator arrays and offers a linear-time numerical method. The metric transitivity proof (Proposition 3.6) is clean, the application of Pastur's theorem is transparent, and the open code and detailed numerical comparisons are strengths. However, the physical relevance of the main theorem is limited by an unproven identification between continuous and discrete resonances, and the tripartite decomposition and the meta-atom algorithm are largely heuristic. The paper is therefore best seen as a rigorous study of the discrete capacitance model with suggestive numerical evidence for the physical system.
major comments (4)
- [Section 2.1, Theorem 2.1] The approximation ω_i(δ) = √δ λ_i + O(δ) is used to identify the physical resonant frequencies with the discrete eigenvalues, and the paper states 'we will often use λ_i and ω_i interchangeably'. However, all subsequent theorems and figures concern the eigenvalues λ_i of the generalised capacitance matrix (or of the Jacobi operator J). The O(δ) error is not shown to be uniform in the number of resonators N, so the empirical measure of ω_i/√δ could differ from that of λ_i in the thermodynamic limit M→∞ at fixed δ, and then δ→0. This gap directly affects the title's claim of a density of states for subwavelength resonator systems. Either provide a uniform (in N and i) error estimate or clearly restrict the claims to the discrete model.
- [Section 5, §§5.1–5.3] The paper extends the convergence and determinism results of Section 3 to dependent samplings. For bound-length sampling (§5.1) it asserts 'all the convergence results from Section 3 continue to hold also under this sampling' without proof; for softmax sampling (§5.2, Eq. (5.2)) the process is defined by a feedback rule and metric transitivity is not verified; for the Fibonacci tiling (§5.3) no probability space is even specified. Since Theorem 3.9 is quoted from Pastur for metrically transitive operators, each of these samplings requires explicit verification of the hypotheses. Without this, the universal determinism of the DoS for these cases is an assumption rather than a theorem.
- [Section 2.2, Figure 2; Section 4] The tripartite decomposition of the spectrum into shared pass band, bandgap, and hybridisation region is presented as a 'complete description'. The bandgap part is rigorously supported by Theorem 2.6, but the smoothness of the DoS in shared pass bands and the fractal-like behaviour in hybridisation regions are only demonstrated numerically. The analytical argument in §4.1 (Proposition 4.1) is a perturbation bound for a single eigenpair and does not imply the universal structure of the density of states. Please either provide formal statements (e.g., asymptotics or scaling laws for the IDS in the hybridisation region) or explicitly label this trichotomy as a numerical observation and adjust the abstract.
- [Section 4.2, Algorithm 1] The meta-atom algorithm is a central contribution, but it is presented without any theoretical error estimate. The Wasserstein-distance convergence shown in Figure 5 is empirical and depends on the chosen blocks and sampling parameters; no proved bound relates the output of Algorithm 1 to the infinite-volume IDS of Theorem 3.9. If the algorithm is intended as a heuristic numerical tool, the paper should state this; if it is intended as a 'universal estimate', a convergence theorem or error bound is needed.
minor comments (7)
- [Definition 2.5] The phrase 'from the left edge of the resonator x_i^L to the left edge of the following resonator x_i^L' is internally inconsistent; the second symbol should be x_{i+1}^L.
- [Figures 2–7] Labels like '10□18' appear to be rendering artifacts of negative exponents; please regenerate the figures so that all annotations are legible.
- [Remark 3.4] The truncation size is denoted by N (with (2N+1) resonators), while N was earlier used for the total number of resonators in a finite system. This overloading is confusing and should be disambiguated.
- [References [4], [6]] The core results rely heavily on the companion preprint [4] and the forthcoming book [6]. Please state the dependence explicitly and, where possible, include the relevant statements in the text or make the references available.
- [Appendix B, Eq. (B.1)] The notation 'C_α := V C_α' reuses the same symbol for the matrix and its generalised version; a distinct notation would improve clarity.
- [Figure 3 caption] The caption says 'M = 105 realisations', which presumably means 10^5 realisations; please correct.
- [Abstract and §4] The term 'fractal-like' is used without a precise definition; since the paper itself notes that hybridisation smooths the density, please clarify the intended meaning (e.g., self-similarity of peaks at a fixed resolution).
Circularity Check
No direct circularity: the DoS convergence rests on Pastur's external ergodic theorem, and the meta-atom estimates are non-fitted approximations; self-citations are present but not load-bearing.
full rationale
The paper's central convergence claim, Theorem 3.9, is quoted from Pastur's external monograph [27]; it is not derived from the authors' prior work. The only inputs are the i.i.d./Markov block-sampling construction and the verification that the induced resonator sequences give a shift-invariant Jacobi operator. Proposition 3.6 checks the metric-transitivity hypothesis directly, so the application of [27] is a genuine import of an external theorem rather than a restatement of the desired conclusion. The meta-atom algorithm likewise does not fit parameters to the target density of states: it enumerates local block arrangements, computes their spectra from the same generalized capacitance model, and compares the resulting CDF to the full empirical CDF via the Wasserstein distance, with Proposition 4.1 providing an a priori perturbation estimate. The tripartite spectral classification uses Theorem 2.6 quoted from [4], a self-citation, but the DoS convergence does not depend on that theorem, and the theorem is a checkable mathematical statement rather than a fitted output. The only substantive caveat is the identification of lambda_i with omega_i via Theorem 2.1's O(delta) error, whose uniformity in N is not established; this is an approximation gap and a physical-modeling risk, not a circular definition. No equation in the paper reduces to its own output by construction, so there is no significant circularity.
Assumptions & free parameters
free parameters (4)
- R: number of attached single-resonator buffers in meta-atom evaluation =
4
- L: maximum meta-atom length =
varied up to 14
- P: maximum number of single resonators in a meta-atom =
scaled as L/2
- beta: softmax sampling temperature =
not stated
assumptions (6)
- domain assumption The generalized capacitance matrix C=VC accurately describes the N subwavelength resonant frequencies in leading order as delta tends to 0 (Theorem 2.1 of [21]), and this approximation is stable as N tends to infinity.
- standard math Pastur's ergodic theorems for metrically transitive operators apply to the Jacobi operator J(mu) built from the resonator sequence mu (Theorem 3.9).
- ad hoc to paper The integrated density of states lambda maps to N(J, lambda) is continuous for the block-disordered Jacobi operators considered (Theorem 3.9(iii)).
- domain assumption With i.i.d. sampling, any frequency in the pass band of at least one constituent block is approached by eigenvalues of the full random system, so the spectral trichotomy (shared pass band, bandgap, hybridisation region) is exhaustive.
- standard math Eigenmodes at frequencies in the single-resonator bandgap decay exponentially through single-resonator blocks, so local dimer arrangements are only weakly coupled.
- domain assumption The finite capacitance matrix eigenvalue count converges to the IDS of the infinite Jacobi operator as N tends to infinity, with edge effects negligible.
invented entities (1)
-
Meta-atoms: finite block sequences constrained by M^P_L used to compute defect modes
Cite this review
Pith. "Pith review of Universal estimates for the density of states for aperiodic block subwavelength resonator systems." pith.science (2026). https://pith.science/paper/5ANOISFP
@misc{pith2026250516677,
author = {Pith},
title = {Pith review of: Universal estimates for the density of states for aperiodic block subwavelength resonator systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/5ANOISFP}},
note = {Machine review of arXiv:2505.16677}
}
abstract
We consider the spectral properties of aperiodic block subwavelength resonator systems in one dimension, with a primary focus on the density of states. We prove that for random block configurations, as the number of blocks $M\to \infty$, the integrated density of states converges to a non-random, continuous function. We show both analytically and numerically that the density of states exhibits a tripartite decomposition: it vanishes identically within bandgaps; it forms smooth, band-like distributions in shared pass bands (a consequence of constructive eigenmode interactions); and, most notably, it exhibits a distinct fractal-like character in hybridisation regions. We demonstrate that this fractal-like behaviour stems from the limited interaction between eigenmodes within these hybridisation regions. Capitalising on this insight, we introduce an efficient meta-atom approach that enables rapid and accurate prediction of the density of states in these hybridisation regions. This approach is shown to extend to systems with quasiperiodic and hyperuniform arrangements of blocks.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[28]
Computation and visualization of photonic quasicrystal spectra via Bloch’s theorem
Rodriguez Alejandro W, McCauley Alexander P, A vnielYehuda, andJohnson Steven G, “Computation and visualization of photonic quasicrystal spectra via Bloch’s theorem”, in:Physical Review B77.10 (2008), p. 104201
work page 2008
-
[29]
Some electronic properties of a one-dimensional crystal model
Saxon D. S. andHutner R. A., “Some electronic properties of a one-dimensional crystal model”, in:Philips J. Res.4 (1949), pp. 81–122
work page 1949
-
[30]
Almost periodic functions and partial differential operators
Shubin Mikhail Aleksandrovich, “Almost periodic functions and partial differential operators”, in:Russian Mathematical Surveys33.2 (1978), p. 1
work page 1978
-
[31]
Almost periodic Schrödinger operators: a review
Simon Barry, “Almost periodic Schrödinger operators: a review”, in:Advances in Applied Mathematics3.4 (1982), pp. 463–490
work page 1982
-
[32]
Electrons in Disordered Systems and the Theory of Localization
Thouless D.J., “Electrons in Disordered Systems and the Theory of Localization”, in: Phys. Rep.13.3 (1974), pp. 93–142
work page 1974
-
[33]
Torquato Salvatore, “Hyperuniform states of matter”, in:Phys. Rep.745 (2018), pp. 1–95
work page 2018
-
[34]
Light in correlated disordered media
Vynck Kevin, Pierrat Romain, Carminati Rémi, Froufe-Pérez Luis S., Schef- fold Frank, Sapienza Riccardo, Vignolini Silvia, and Sáenz Juan José, “Light in correlated disordered media”, in:Rev. Mod. Phys.95.4 (2023), p. 045003
work page 2023
-
[35]
Engineered disorder in photonics
Yu Sunkyu, Qiu Cheng-Wei, Chong Yidong, Torquato Salvatore, and Park Namkyoo, “Engineered disorder in photonics”, in:Nature Reviews Materials6.3 (2021), pp. 226–243
work page 2021
Show all 10 references
-
[36]
Crystalline Metamaterials for Topological Properties at Subwavelength Scales
Yves Simon, FleuryRomain, BerthelotThomas, Fink Mathias, LemoultFabrice, and Lerosey Geoffroy, “Crystalline Metamaterials for Topological Properties at Subwavelength Scales”, in:Nat. Commun.8.1 (2017), p. 16023
2017
-
[37]
A Review on Non-Hermitian Skin Effect
Zhang Xiujuan, Zhang Tian, Lu Ming-Hui, and Chen Yan-Feng, “A Review on Non-Hermitian Skin Effect”, in:Advances in Physics: X7.1 (2022), p. 2109431,doi: 10.1080/23746149.2022.2109431. Habib Ammari ETH Zürich, Department of Mathematics, Rämistrasse 101, 8092 Zürich, Switzerland...
2022
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.