{"id":"20b658c3-30c3-4e8c-8126-0948d94a6220","arxiv_id":"2505.16677","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For random block subwavelength resonator chains, the integrated density of states converges to a non-random continuous limit, and a meta-atom expansion predicts its fractal-like structure in hybridisation regions.","lead":"This paper analyzes long chains of subwavelength resonators built from random arrangements of resonator blocks, and proves that the distribution of resonant frequencies becomes a predictable, non-random curve as the chains grow. It also introduces a fast meta-atom sampling method that accurately reproduces the density of states in the most complex spectral regions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper proves DoS convergence for the discrete capacitance model, but Theorem 2.1's O(δ) approximation to the true subwavelength resonances is not shown uniform in N; the physical DoS may differ.","rationale":"The reader's weakest assumption is the same one I would flag: the paper's universal DoS theorem concerns the discrete Jacobi operator/capacitance matrix, not the actual subwavelength resonance frequencies of the Helmholtz problem (2.2). Theorem 2.1 provides an O(δ) approximation of each ω_i by √δ λ_i, but no estimate of the remainder uniform in the number of blocks M is given. Because the paper's central claim is about physical resonator systems and all figures use the discrete eigenvalues, this gap is load-bearing. If the O(δ) constant grows with M, the thermodynamic limit of the physical IDS may differ from the discrete IDS, and the deterministic continuity theorem would not apply to the physical problem. I agree with the reader's assessment; the concern does not overturn the paper's contribution but should be addressed. The proposed numerical test would settle it: compare the exact (or higher-order) subwavelength frequencies with the discrete eigenvalues for increasing M and fixed small δ. The paper's metric transitivity proof and the meta-atom algorithm with released code are independent strengths, but they do not close the discrete-continuum gap.","tokens_in":22691,"tokens_out":9017,"duration_ms":73288,"concrete_test":"Compute the exact subwavelength eigenfrequencies ω_i(δ) for the two-block system of Example 2.2 at δ=10^-6 and M=2^p, p=5,...,12, by solving the resonance problem (2.2) numerically (e.g., via a root-finder on the transfer-matrix transmission condition or via higher-order asymptotic beyond Theorem 2.1). Compare the empirical CDF of ω_i/√δ with the empirical CDF of the eigenvalues of C via the Wasserstein distance. If the distance fails to be O(√δ) uniformly in M (e.g., does not stay bounded by a constant independent of M), then the physical IDS could differ from the discrete IDS, and Theorem 3.9 does not apply to the resonator system.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theorem (Theorem 3.9) is proved for the infinite Jacobi operator J = V^{1/2} C V^{1/2}, and all numerical DoS plots are empirical CDFs of the eigenvalues of the finite generalised capacitance matrix C. However, the physical quantity advertised in the title—the density of states of a subwavelength resonator system—is the density of the actual resonant frequencies ω_i(δ) solving (2.2). Theorem 2.1 gives ω_i(δ) = √δ λ_i + O(δ), with an error term whose dependence on the number of blocks M is never estimated. The paper's blanket statement in Section 2.1 that 'we will often use λ_i and ω_i interchangeably' is therefore an unproven identification. For the weak limit of the empirical measures to be the same for ω_i/√δ and λ_i, one needs the O(δ) error to be uniform in i and M (at least in a counting sense). If the constant in the O(δ) grows with M, e.g. through the condition number of C or through small spectral gaps, then taking the thermodynamic limit M→∞ at fixed δ and then δ→0 can produce a different IDS. Since Theorem 3.9, the tripartite decomposition, and the meta-atom algorithm all operate on the discrete model, this gap means the universal continuity and determinism of the physical DoS is not established by the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":22978,"tokens_out":10118,"duration_ms":86628,"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":[{"comment":"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":"Section 2.1, Theorem 2.1"},{"comment":"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":"Section 5, §§5.1–5.3"},{"comment":"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":"Section 2.2, Figure 2; Section 4"},{"comment":"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.","section":"Section 4.2, Algorithm 1"}],"minor_comments":[{"comment":"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.","section":"Definition 2.5"},{"comment":"Labels like '10□18' appear to be rendering artifacts of negative exponents; please regenerate the figures so that all annotations are legible.","section":"Figures 2–7"},{"comment":"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.","section":"Remark 3.4"},{"comment":"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.","section":"References [4], [6]"},{"comment":"The notation 'C_α := V C_α' reuses the same symbol for the matrix and its generalised version; a distinct notation would improve clarity.","section":"Appendix B, Eq. (B.1)"},{"comment":"The caption says 'M = 105 realisations', which presumably means 10^5 realisations; please correct.","section":"Figure 3 caption"},{"comment":"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).","section":"Abstract and §4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the companion preprint [4] and the forthcoming work [6] by the same group. The editor may wish to ensure that the overlap with [4] is properly delineated. In addition, the central theorem's physical interpretation hinges on an unproven uniform-in-N estimate connecting the discrete and continuous resonance problems; this could be addressed either by adding such an estimate or by reframing the paper's claims as specific to the discrete capacitance model. The tripartite spectral classification and the meta-atom algorithm, while numerically convincing, are not proven with the same rigor as the ergodic convergence theorem and should be labelled accordingly in any revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on arXiv:2505.16677. The paper is worth reading, but read it as a paper about the discrete capacitance model, not literally about the Helmholtz resonances.\n\nThe genuinely new thing is the meta-atom construction for the density of states in hybridisation regions, plus the demonstration that it works for quasiperiodic and hyperuniform sampling. The numerics are convincing and the algorithm is linear-time in the number of blocks. The paper is also honest: it distinguishes the rigorous ergodic part (Pastur applied to metrically transitive Jacobi operators) from the heuristic parts (the trichotomy, the meta-atom accuracy). The proof of metric transitivity for the resonator sequence is clean, and the convergence result follows faithfully from Pastur.\n\nThe biggest soft spot is the gap between the discrete model and the physical resonances. Theorem 2.1 gives ω_i = sqrt(δ) λ_i + O(δ), but the O(δ) is not shown to be uniform in N. Since the IDS theorems and all numerical plots are for the λ_i's, the paper does not actually prove that the density of states of the physical system converges to the stated non-random continuous limit. This is not likely to be wrong in practice—the capacitance approximation is standard—but it is unproven, and the paper's blanket \"we will use λ_i and ω_i interchangeably\" hides the issue. A referee should ask for a uniformity estimate or a careful statement limiting the claims to the discrete model.\n\nTwo smaller things. The \"if and only if\" in the trichotomy (bandgap iff in all block gaps) has a rigorous \"if\" from the Saxon-Hutner theorem and a plausible but heuristic \"only if\". I don't think this is load-bearing, but it should be flagged as a conjecture rather than a theorem. And the meta-atom convergence is empirical; there is no bound on the Wasserstein error as a function of L and P. That's fine for an applied paper, but it means \"universal estimates\" is a bit of an overpromise.\n\nOverall, the paper is solid, clearly written, and a useful contribution to the subwavelength literature. The central ergodic theorem is correct as stated for the discrete operator, and the algorithm is a practical tool. I would send it to peer review. The referee should focus on the uniform O(δ) issue and secondarily on the status of the trichotomy.","headline":"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.","tokens_in":23567,"tokens_out":3247,"would_cite":true,"duration_ms":26150,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35J05","35C20","35P20"],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["density of states","block disordered systems","subwavelength resonators","hybridisation regions","metric transitivity","Jacobi operators","quasiperiodic sampling","hyperuniform sampling"],"falsifier":"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.","tokens_in":22417,"feed_emoji":"📊","tokens_out":12340,"duration_ms":96343,"temperature":0.7,"pith_summary":"The paper studies long chains of subwavelength acoustic resonators built by sampling independently from a small set of building blocks. It sets out to show that, once the chain is long, the distribution of resonant frequencies becomes a single non-random, continuous curve, independent of the particular random sample; numerically it finds that this limiting density of states splits into three regimes: zero in bandgaps, smooth bands where every block type passes, and a fractal-like set of peaks in hybridisation regions where some block types are gapped. It goes on to explain that the fractal-like peaks come from eigenmodes decaying exponentially through the gapped blocks, and uses that insight to build a meta-atom estimator that reconstructs the hybridisation-region spectrum in time linear in the number of blocks. A reader should care because it says a deterministic, computable spectral law governs nominally disordered acoustic metamaterials, and the exotic fractal-looking spectrum is a finite-size signature rather than the infinite limit.","feed_headline":"Aperiodic resonator chains get one universal density of states","feed_subtitle":"For random block arrays the eigenvalue count becomes a fixed continuous curve, computable in linear time via meta-atoms.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the metric-transitivity theorems from which the almost-sure convergence to a non-random continuous integrated density of states is imported.","marker":"[27]"},{"why":"Introduces the block-disordered resonator model, propagation matrices, the Saxon-Hutner-type bandgap theorem, and the Thouless criterion used throughout.","marker":"[4]"},{"why":"Establishes the existence of $N$ subwavelength resonances and the generalised capacitance matrix whose eigenvalues approximate them with an $O(\\delta)$ error.","marker":"[21]"},{"why":"Provides the random word model whose one-dimensional localisation results the paper extends to classical resonator systems.","marker":"[15]"},{"why":"Saxon-Hutner theorem underlying the claim that a frequency gapped in every block is gapped for the whole system.","marker":"[29]"},{"why":"Gives the autocovariance/Fourier characterisation of hyperuniformity used to certify chunk and softmax sampling.","marker":"[33]"}],"fun_headline_variants":["Aperiodic resonator arrays yield a universal density of states","Random block resonators converge to one continuous DOS","Meta-atom trick speeds up fractal DOS prediction","One density of states for all aperiodic resonator chains","Subwavelength resonators: universal DOS via meta-atoms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Aperiodic resonator arrays yield a universal density of states","Random block resonators converge to one continuous DOS","Meta-atom trick speeds up fractal DOS prediction","One density of states for all aperiodic resonator chains","Subwavelength resonators: universal DOS via meta-atoms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000368,"raw_usage":{"total_tokens":1982,"prompt_tokens":958,"completion_tokens":1024,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":946}},"tokens_in":574,"tokens_out":1024,"duration_ms":8477,"temperature":1.0,"reasoning_tokens":946,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:57:37.066273+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Some electronic properties of a one-dimensional crystal model","cited_arxiv_id":null,"evidence_quote":"Saxon-Hutner theorem underlying the claim that a frequency gapped in every block is gapped for the whole system."},{"cited_title":"Hyperuniform states of matter","cited_arxiv_id":null,"evidence_quote":"Gives the autocovariance/Fourier characterisation of hyperuniformity used to certify chunk and softmax sampling."}],"review_version":1}