{"id":"9d8970fc-7bbb-4a3d-9a8e-885b8cb7bebf","arxiv_id":"2411.15906","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper proves explicit supercell convergence rates for quasiperiodic operator spectra and a spectrum-equality theorem for the superspace lift, though the numerical discretisation of the superspace method still shows spurious eigenvalues without proven error bounds.","lead":"Two standard shortcuts for computing spectra of quasiperiodic operators, the supercell and superspace methods, are shown to converge to the true spectrum, with explicit error rates for the supercell approach. The results put the simulation practice behind quasicrystal metamaterial design on a rigorous footing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Superspace Fourier-mode mismatch in Theorem 4.1 invalidates the proof as printed; the lift in (4.6) uses y-frequency θn but D_θ requires y-frequency n to match the slice.","rationale":"The reader's weakest assumption—the Fourier-mode mismatch in Theorem 4.1—is the most load-bearing concern because it invalidates the proof of a central claim as printed, not merely a rate estimate. I verified algebraically that replacing θn by n in (4.6) makes the action of B_θ on plane waves match that of A_θ on generalized Fourier coefficients, so the theorem is likely true and the flaw is a typo. Nevertheless, the preprint as written does not establish σ(A_θ)=σ(B_θ). The supercell theorem also contains index and exponent typos (Lemma 3.5 states |θ−θ_l|^{-1/2} but the proof yields q_l^{-1}; the recursion in (3.27)-(3.28) has q_n^{-1} instead of q_{n-1}^{-1}), but these are cosmetic and the underlying argument is sound. Therefore the reader's CONDITIONAL verdict is appropriate; no verdict change is needed.","tokens_in":17170,"tokens_out":28725,"duration_ms":234949,"concrete_test":"Repair the lift in (4.6) to ~F(x,y)=e^{ikx}\\sum f_{mn}e^{2πi(mx+ny)} and expand B_θ~F in the Fourier basis of T^2. Verify that the coefficient of e^{2πi(m'x+n'y)} equals the coefficient with which A_θ maps the generalized Fourier coefficient f_{m'n'} to frequency m'+θn'. If the identity holds for all integers m',n', the proof is repairable; if it fails, Theorem 4.1 is unsupported. Independently, compute σ(B_θ) by a fine rectangular-mesh discretisation for the Schrödinger example (3.32) and compare against the supercell limit from Theorem 3.4, checking for extra spectral values not attributable to discretisation pollution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Theorem 4.1 rests on the claim that the lifted operator B_θ acts on the Fourier coefficients of ~F in (4.6) exactly as A_θ acts on the generalized Fourier coefficients of ~f in (4.5). This is false as displayed. In (4.6), ~F(x,y)=e^{ikx}\\sum f_{mn}e^{2πi(mx+θny)}. The directional derivative D_θ = ∂_x + θ∂_y sends e^{2πi(mx+θny)} to 2πi(m+θ^2n)e^{2πi(mx+θny)}, so the slice y=θx of ~F has frequencies m+θ^2n, not the m+θn appearing in ~f and in the Fourier expansion of A_θ. The subspaces E_k in the proof have the same defect. Consequently the claimed equality σ(A_θ)=σ(B_θ) is not established by the argument in Section 4.1. The error is likely a typo—replacing the y-frequency by n makes the mode-matching identity hold algebraically—but as printed the proof of a central theorem is invalid. A secondary gap is that B_θ is non-elliptic, so the 'converse inclusion' via Fourier series of approximate eigenvectors of B_θ needs a spectral theory that is not supplied.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops convergence theory for two numerical strategies for spectra of one-dimensional quasiperiodic elliptic operators. For the supercell method, Theorem 3.4 asserts an explicit Hausdorff-type rate of order q_l^{-1} when the quasiperiodic slope theta is replaced by its continued-fraction approximants p_l/q_l, and the authors give a largely self-contained argument for the key localization lemma. Theorem 4.1 claims equality between the spectrum of the quasiperiodic operator A_theta and that of its superspace lift B_theta, and Section 5 applies the supercell estimates to prove existence and decay of localised interface modes. Numerical examples for Schrödinger and Helmholtz-type operators illustrate both methods.","tokens_in":17415,"tokens_out":15142,"duration_ms":129800,"significance":"If fully established, Theorem 3.4 would provide an explicit convergence rate for supercell approximations of quasiperiodic spectra and would rigorously justify the common practice of reading persistent gaps from supercell band diagrams; this is a useful and timely result for the metamaterials literature. The superspace equality of Theorem 4.1 would likewise put a widely used heuristic on a firm footing. The paper’s supercell arguments are mostly elementary and self-contained, and the continued-fraction construction is a clear strength. However, the superspace proof is invalid as printed, and the interface-mode theorem in Section 5 is not supported by the arguments given; these points must be repaired before the advertised claims can be accepted.","major_comments":[{"comment":"Equation (4.6) does not define a lift that matches the slice used in (4.5). The mode e^{2\\pi i(mx+\\theta ny)} is mapped by D_\\theta = \\partial_x + \\theta \\partial_y to 2\\pi i(m+\\theta^2 n)e^{2\\pi i(mx+\\theta ny)}, whereas the generalized Fourier coefficient in (4.5) has frequency m+\\theta n. Thus the claim that B_\\theta acts on the Fourier coefficients of ~F exactly as A_\\theta acts on the generalized Fourier coefficients of ~f is false as printed. Replacing the y-frequency \\theta n by n repairs the algebraic identity, but the converse inclusion also needs a spectral theory for the non-elliptic operator B_\\theta, which is not supplied by Theorem 2.7. As it stands, Theorem 4.1 is not proven.","section":"4.1, Eq. (4.6)"},{"comment":"The interface-mode theorem is not established. The condition (5.2) uses an undefined index n and has a quantifier over lambda in sigma(A_theta) that conflicts with the fixed lambda whose separation is being assumed. More substantively, “interface eigenvalue for A_{theta_{l0}}” is never defined, and the periodic operator A_{theta_{l0}} has no eigenvalues on the whole line; some truncated or defect problem must be specified. The proof then uses Theorem 3.4 to infer the existence of a limiting eigenvalue, but Theorem 3.4 only bounds distances to the spectra sigma(A_{theta_l}) and sigma(A_theta); it does not by itself produce an isolated eigenvalue of A_theta. Finally, the appeal to [1,17,47] for periodic band-gap interface modes does not automatically transfer to the quasiperiodic limit. The authors should either give a precise statement with a complete proof or clearly label this as a conjecture.","section":"5, Theorem 5.1 and Eq. (5.2)"},{"comment":"The proof of the upward bound is not rigorous as written. After estimating a point lambda in sigma(A_{theta_{l+2}}), the text introduces a recursion d_{n-1} \\le C(1+|\\lambda|+d_n) q_n^{-1} + d_n with d_n called the distance from lambda to sigma(A_{theta_n}) and d_0 = 0; if d_n is the distance for a fixed lambda, then the base case d_{l+2}=0 holds only for lambda in sigma(A_{theta_{l+2}}), and the passage to arbitrary lambda in sigma(A_theta) requires the spectral-inclusiveness argument from Proposition 3.1 (cited as “Theorem 3.1”), which is not spelled out. Likewise, the first inequality “follows from Lemma 3.5 by letting m \\to \\infty” needs a no-spectral-pollution statement for the sequence A_{theta_m}. These gaps are likely repairable, but they are load-bearing for the central theorem.","section":"3.1, proof of Theorem 3.4"}],"minor_comments":[{"comment":"The title contains spacing typos (“supersp ace”, “opera tors”); the whole manuscript should be proofread for similar artifacts.","section":"Title"},{"comment":"The sentence “The claim now follows from Theorem 3.1” should presumably refer to Proposition 3.1 rather than Theorem 3.1.","section":"3.1, end of proof of Theorem 3.4"},{"comment":"The displayed inequality compares a_k(x,\\theta x) with a_k(x,\\theta_l), but the surrounding argument compares A_{\\theta_m} and A_{\\theta_l}; the first argument should be \\theta_m x.","section":"3.1, Eq. (3.22)"},{"comment":"The caption says “as the length of the periodic unit cell decreases”, but q_l increases along the sequence; the intended statement is that the unit cell becomes longer.","section":"3.3, Figure 3.1 caption"},{"comment":"No convergence theorem is stated for the finite-difference or plane-wave discretisations used in the numerical examples; the figures are therefore heuristic illustrations rather than rigorous consequences of Theorems 3.4 and 4.1.","section":"4.2"}],"recommendation":"major_revision","confidential_remarks":"The core supercell result is promising, and the continued-fraction construction is a genuine asset. The superspace theorem is likely fixable by correcting one line, but the interface-mode section needs a substantially more careful treatment before the paper can be accepted. I would not recommend rejection, because the main ideas appear sound and the gaps are repairable within the manuscript’s scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the supercell convergence theorem is the real result and it's good; the superspace equality is not proven as printed due to a mode-matching error that looks like a typo.\n\nThe explicit continued-fraction rates in Theorem 3.4 are a genuine step up from the qualitative convergence in Damanik–Goldstein–Lukic. The proof is mostly self-contained, the constants are explicit enough, and the tiling analogue in Section 3.2 is a nice extra. The numerical sections are honest: they show the superspace plane-wave method pollution and point out that the finite-difference choice mitigates it, which is a useful practical warning. The interface-mode theorem is a reasonable application of the supercell estimates.\n\nThe soft spot is in Section 4.1. Equation (4.6) defines ~F with frequencies m in x and θn in y. But D_θ = ∂_x + θ∂_y sends e^{2πi(mx+θny)} to m + θ²n, not the m+θn that appears in the reduced operator's Fourier expansion. So the asserted preservation of subspaces E_k and the claimed action of B_θ on the Fourier coefficients of ~F do not hold as written. The fix is almost certainly to set the y-frequency to n, i.e. use e^{2πi(mx+ny)} in (4.6); then the algebra works. As printed, the proof of σ(A_θ)=σ(B_θ) does not go through, and since B_θ is non-elliptic, the converse inclusion needs more care than the one-line 'similar argument' provides. This is a repairable flaw, not a fatal one, but it is load-bearing for the superspace half of the paper.\n\nMinor quibbles: in Lemma 3.5 the sentence 'by letting δ→∞' is not what you want (you should optimise δ), and Theorem 5.1 leans on external results a bit heavily. Neither undermines the supercell argument.\n\nBottom line: the supercell contribution deserves attention and the paper should go out after a revision that fixes the superspace proof and adds a few lines on the converse inclusion. A serious referee would be worth it.","headline":"The supercell convergence theorem is a solid, citable result with explicit rates, but the superspace spectral equality is not proven as printed due to a Fourier-mode mismatch that looks like a typo.","tokens_in":17934,"tokens_out":2173,"would_cite":true,"duration_ms":18684,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35P05","81Q10"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that continued-fraction supercell approximants converge to the true quasiperiodic spectrum with explicit error bounds, and that the superspace lift has the same spectrum.","keywords":["quasiperiodic operators","supercell method","superspace method","Floquet-Bloch theory","continued fractions","Fibonacci tiling","spectral gaps","interface modes"],"falsifier":"Carry out the displayed mode matching directly: applying $D_\\theta = \\partial_x + \\theta\\,\\partial_y$ to $e^{2\\pi i(mx+\\theta n y)}$ gives frequency $m + \\theta^2 n$, whereas the one-dimensional operator acts on $e^{2\\pi i(m+\\theta n)x}$ with frequency $m + \\theta n$; the two match only if $\\theta\\in\\{0,1\\}$, so this computation would refute the superspace theorem as printed and force a corrected lift.","tokens_in":16931,"feed_emoji":"📐","tokens_out":9806,"duration_ms":73434,"temperature":0.7,"pith_summary":"Quasicrystalline materials have exotic wave properties, but the differential operators that describe them have spectra that are notoriously hard to compute reliably. This paper targets the two most common computational shortcuts: the supercell method, which repeats a finite patch of the material periodically, and the superspace method, which treats the aperiodic material as a slice of a higher-dimensional periodic structure. For the supercell method it proves that, when the approximating periods come from continued fractions of the irrational parameter, the Floquet-Bloch spectra converge to the true quasiperiodic spectrum in Hausdorff distance at an explicit rate, with errors of order $(1+|\\lambda|)q_l^{-1}$. A practical consequence is that band gaps that persist across supercell approximations, known as super band gaps, are guaranteed to be genuine gaps of the limiting operator. For the superspace method the paper proves that the lifted periodic operator has exactly the same spectrum, and shows that the spurious eigenvalues seen in earlier implementations depend on the choice of numerical discretisation.","feed_headline":"Supercell spectra provably converge to quasiperiodic spectrum","feed_subtitle":"Continued-fraction error bounds show persistent band gaps are real gaps of the limiting operator.","key_machinery":"The engine of the supercell result is the continued-fraction approximant $\\theta_l = p_l/q_l$ combined with a localised approximate-eigenfunction construction. A Bloch eigenfunction of the periodic approximant is multiplied by a smooth cut-off supported on a window of length comparable to $q_l$; the commutator with the operator is then $O(q_l^{-1})$, and the coefficient mismatch between approximants is controlled by the continued-fraction bound $|\\theta_m-\\theta_l| \\le (q_l q_{l+1})^{-1}$. For the superspace claim the load-bearing object is the lifted operator $B_\\theta = \\sum_k a_k(x,y) D_\\theta^k$ with $D_\\theta = \\partial_x + \\theta\\,\\partial_y$, whose Fourier modes are matched to the generalised Fourier modes of the one-dimensional operator. A secondary mechanism is the trace-map recursion for Fibonacci tilings, which provides the explicit super-band-gap criterion used to certify persistent gaps from finitely many iterates.","core_discovery":"The paper's central claim is that two widely used computational heuristics for quasiperiodic spectra are convergent in a strong, quantitative sense. Theorem 3.4 states that if $A_\\theta$ is a self-adjoint elliptic operator with quasiperiodic coefficients and $\\theta_l = p_l/q_l$ is the $l$-th continued fraction approximant of $\\theta$, then every point of the supercell spectrum $\\sigma(A_{\\theta_l})$ lies within distance $C(1+|\\lambda|) q_l^{-1}$ of $\\sigma(A_\\theta)$, and conversely every point of $\\sigma(A_\\theta)$ lies within an explicitly bounded distance of $\\sigma(A_{\\theta_l})$ expressed as a convergent series. Hence the supercell spectra converge to the quasiperiodic spectrum in Hausdorff distance, and persistent spectral gaps in a supercell band diagram must be gaps of the limiting operator. Theorem 4.1 asserts the same spectrum for the superspace lift $B_\\theta = \\sum_k a_k(x,y) D_\\theta^k$ on the two-dimensional torus, with $D_\\theta = \\partial_x + \\theta\\partial_y$, thereby validating the superspace method and attributing its observed spectral pollution to the numerical discretisation rather than to the lifting idea.","pith_inferences":["If the cut-off construction carries over to several variables, the same strategy would give supercell convergence for higher-dimensional quasiperiodic operators, with the continued-fraction rate replaced by a multi-dimensional Diophantine approximation error.","The mode-matching inconsistency in the superspace proof suggests a repaired formulation: taking the lifted plane-wave frequency to be $m + n$ rather than $m + \\theta n$ would make the identity consistent, and the equality of spectra would then plausibly hold as stated.","The error bound can be read as a design rule for quasicrystal metamaterials: one can work backwards from a desired gap resolution to the required continued-fraction index, a practical step the paper leaves implicit.","The same localised-mode argument likely applies to other symmetry-induced interfaces, such as phase-shift or anti-reflection defects, provided the defect eigenvalue is sufficiently well separated from the approximant spectrum."],"forward_implications":["Supercell band diagrams built from continued-fraction approximants converge to the true quasiperiodic spectrum at a quantified rate, so persistent spectral gaps seen in finite supercell plots are provably gaps of the infinite system.","The explicit error bound makes supercell computations a controlled numerical method: a target accuracy fixes the continued-fraction depth needed, with no guesswork about supercell size.","For golden-mean Fibonacci tilings the analogous convergence holds with errors controlled by Fibonacci numbers, and the trace-map criterion identifies super band gaps from finitely many iterates.","The superspace spectrum equals the original quasiperiodic spectrum when the lift is discretised consistently, and the spectral pollution previously reported is tied to the discretisation (notably plane-wave truncation) rather than the method itself.","Localised interface modes in reflection-symmetric quasicrystals persist from supercell approximants, with estimates for their eigenfrequencies and exponential decay rates."],"supporting_citations":[{"why":"Supplies the isospectral-torus periodic-approximation framework that the supercell theorem refines with explicit supercell constructions and rates.","marker":"[19]"},{"why":"Provides the almost-periodic function and operator theory, including the spectrum identity used in the superspace proof.","marker":"[44]"},{"why":"Introduced the superspace spectral method whose spectrum and spectral pollution the paper analyses.","marker":"[42]"},{"why":"Established super band gaps for generalised Fibonacci tilings, which the paper connects to the limiting spectrum.","marker":"[22]"},{"why":"Supplies continued-fraction approximation theory giving the $q_l^{-1}$ error rate used throughout the supercell argument.","marker":"[30]"},{"why":"Gives the rigorous existence of quasiperiodic Floquet solutions for analytic Schrödinger potentials, cited in the superspace discussion.","marker":"[24]"},{"why":"Provides the trace-map recursion for Fibonacci tiling transfer matrices used in the super-band-gap criterion.","marker":"[31]"},{"why":"Used with the Liouville transform to show band gaps support only finitely many interface eigenvalues, supporting the interface-mode theorem.","marker":"[17]"}],"fun_headline_variants":["Supercell spectra provably converge to quasiperiodic limit","Explicit error bounds for supercell quasiperiodic spectra","Persistent band gaps in supercell diagrams are real gaps","Superspace method validated: spectral pollution is numerical only","Supercell and superspace spectra both converge to the quasiperiodic spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The superspace part of the paper hinges on a frequency-matching identity between the original operator's modes and the plane-wave modes of its two-dimensional lift; if that identity is wrong, the claimed equality of spectra does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Supercell spectra provably converge to quasiperiodic limit","Explicit error bounds for supercell quasiperiodic spectra","Persistent band gaps in supercell diagrams are real gaps","Superspace method validated: spectral pollution is numerical only","Supercell and superspace spectra both converge to the quasiperiodic spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000442,"raw_usage":{"total_tokens":2192,"prompt_tokens":851,"completion_tokens":1341,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":467,"completion_tokens_details":{"reasoning_tokens":1257}},"tokens_in":467,"tokens_out":1341,"duration_ms":10142,"temperature":1.0,"reasoning_tokens":1257,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:46:32.362031+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Carry out the displayed mode matching directly: applying $D_\\theta = \\partial_x + \\theta\\,\\partial_y$ to $e^{2\\pi i(mx+\\theta n y)}$ gives frequency $m + \\theta^2 n$, whereas the one-dimensional operator acts on $e^{2\\pi i(m+\\theta n)x}$ with frequency $m + \\theta n$; the two match only if $\\theta\\in\\{0,1\\}$, so this computation would refute the superspace theorem as printed and force a corrected lift.","supporting_citations":[],"review_version":1}