{"id":"3d764814-d92d-4b8a-ad37-fb1a03d82a26","arxiv_id":"2510.15369","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a 1D incommensurate Schrödinger toy model, the second-order semiclassical DoS matches momentum-space numerics at small ε, and the residual oscillations at larger ε are traced to harmonic-oscillator levels at band-structure critical points.","lead":"This paper compares two ways of computing the density of states of incommensurate (moiré-like) Schrödinger operators on a 1D toy model: a momentum-space truncation and a semiclassical expansion in the incommensurability ratio. It shows the two agree for small ratio, and explains the mismatch at larger ratio as oscillations near moiré-scale 'Van Hove' points, captured by a harmonic oscillator model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Momentum-space reference is not demonstrably converged in the critical σ=0.08, ε=0.01 regime; negative DoS values indicate artifacts that could explain part of the discrepancy.","rationale":"The reader's weakest assumption—that the reference momentum-space computation is converged—is indeed the most load-bearing assumption in the paper. Every comparison that establishes the central claim (agreement at ε=0.001, discrepancy at ε=0.01, and the quantitative fit of the harmonic oscillator model) uses the momentum-space result as ground truth. The negative values in Fig. 5 are a concrete, checkable red flag that this assumption may fail in exactly the regime where the discrepancy is explained. Because the paper does not report the KPM expansion order, the reader cannot determine whether the negative values are due to under-resolution of the Chebyshev expansion, insufficient truncation, or some other numerical issue. The proposed test—increasing the KPM order and checking both nonnegativity and convergence—would settle the point in a straightforward way. If the test shows the reference is converged, the paper's conclusions stand; if not, the attribution of the discrepancy to the semiclassical truncation is weakened. Since this concern is exactly the reader's weakest_assumption and the reader already issued a CONDITIONAL verdict, my review does not move the verdict; it reinforces the conditionality with a concrete, falsifiable check.","tokens_in":26290,"tokens_out":2755,"duration_ms":24079,"concrete_test":"Recompute the momentum-space ν_ε,σ for the two critical parameter sets—(σ=0.08, ε=0.01) and (σ=0.04, ε=0.01)—with increasing KPM/Chebyshev expansion order N (e.g., N = 2^13, 2^15, 2^17) while keeping W, L, h fixed. For each N, (i) check that min_E ν_ε,σ(E) ≥ -tol for a small tolerance (e.g., 10^-6 or 10^-8); (ii) compute the sup-norm difference between successive orders, δ_N = ||ν_ε,σ^(N) - ν_ε,σ^(N/2)||_∞; and (iii) compare the momentum-space result at the highest N to the semiclassical approximation (2.29). If the negative values vanish, δ_N is below a small fraction of the observed discrepancy, and the difference from the semiclassical result is unchanged, then the reference is converged and the paper's attribution survives. If instead the negative values persist or δ_N is comparable to the discrepancy, the comparison must be repeated with a better-converged reference before drawing c","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that the second-order semiclassical expansion (2.29) matches the true DoS at small ε and fails at ε=0.01 due to missing higher-order corrections—is validated entirely against the momentum-space computation (2.10). But in the very regime where the failure is analyzed, σ=0.08 and ε=0.01 (Fig. 5, right), the plotted momentum-space ν_ε,σ(E) takes negative values. Since ν_ε,σ is the convolution of a positive measure with a Gaussian, it must be nonnegative. A negative value is a rigorous sign of numerical artifact, not a property of the true DoS. The paper reports the truncation parameters W=80, L=5000, h=0.01 for σ=0.08 but never reports the KPM/Chebyshev expansion order, so the reference accuracy is unquantified. In the σ=0.04 experiments (Figs. 6–11) the smearing is even narrower, requiring a higher expansion order; the same issue may affect those results. If the momentum-space artifacts are comparable in magnitude to the method-to-method differences in Fig. 5, then attributing the discrepancy entirely to the semiclassical truncation is not cleanly established, and the harmonic-oscillator validation in Figs. 7, 10, and 11 is also called into question because it uses the momentum-space curves as ground truth.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares two numerical approaches for the density of states of one-dimensional incommensurate Schrödinger operators of the form H_ε = -1/2 d²/dx² + V(x,(1+ε)x): a momentum-space method with KPM evaluation, and a semiclassical expansion in ε. The main theoretical output is Theorem 2.5, which gives explicit sum-over-states formulae (2.25)–(2.27) for the first three coefficients L0, L1, L2 of the expansion Tr[f(H_ε)] = L0 + εL1 + ε²L2 + ..., with a detailed derivation in Appendix A. Numerically, the authors compare L0, L1, L2 with finite differences of the momentum-space DoS, compare the full regularized DoS for ε=0.001 and ε=0.01, and attribute discrepancies at ε=0.01 to oscillations at Van Hove critical points captured by a harmonic-oscillator effective model (3.5)–(3.6).","tokens_in":26505,"tokens_out":10324,"duration_ms":82639,"significance":"If the numerical convergence questions are resolved, this is a valuable contribution. The paper provides a practical, parameter-free second-order semiclassical formula for the DoS of two-scale Schrödinger operators, with a self-contained algebraic derivation in Appendix A, and a predictive harmonic-oscillator description of DoS oscillations at band critical points. The comparisons at small ε (Figures 1–4) show impressive agreement, and the authors report their discretization parameters, which aids reproducibility. However, the momentum-space reference computation is not fully characterized in the very regime used to validate the failure analysis, and some quantitative claims lack supporting derivations.","major_comments":[{"comment":"For σ=0.08 and ε=0.01, the momentum-space approximation of ν_{ε,σ}(E) (solid blue) takes negative values in the energy range [9,18]. Since ν_{ε,σ} is the convolution of the positive DoS measure with a positive Gaussian, it must be nonnegative. A negative value is a rigorous indicator of numerical artifact, not a property of the true DoS. The manuscript does not discuss this, and the KPM/Chebyshev expansion order is not reported, so the reference accuracy is unquantified. This is exactly the regime in which the paper claims the second-order semiclassical method fails due to missing higher-order corrections; the reference must be demonstrated converged (e.g., by increasing Chebyshev order and showing stability) before that attribution is clean. The same issue affects the σ=0.04 results in §3.3, which use the momentum-space curves as ground truth.","section":"§3.2, Figure 5 (right)"},{"comment":"The momentum-space approximations of L1 and L2 are obtained by finite differences in ε, but the step size and stencil are not reported. The reported agreement (errors up to 10^-5 for L1 and 10^-3 for L2 at σ=0.08) could be limited by finite-difference truncation rather than by the methods themselves. Without the step size and a convergence check in the step, the consistency of the expansion coefficients is not fully established. Please provide this information.","section":"§3.2, Figures 1–3"},{"comment":"The harmonic approximation formula (3.6) is stated without derivation. In particular, the prefactor ε/|Ω| is not derived from the semiclassical quantization, and the calculation is restricted to a single band with Bp assumed zero. A derivation, or at least a precise statement of the approximations leading to (3.6), is necessary to assess the quantitative claims in Figures 7, 10, and 11. Currently the prefactor and the neglect of interband coupling are not justified, weakening the claim that the harmonic model is quantitative.","section":"§3.3, Eq. (3.6)"}],"minor_comments":[{"comment":"The abstract and body contain several typographical artifacts ('Schr¨odinger', 'DoS', etc.). Please proofread.","section":"General"},{"comment":"References [16] and [17] are identical; one should be removed or replaced.","section":"References"},{"comment":"The theory is stated for irrational ε, but all numerical experiments use rational ε=0.001 and 0.01. Please explain whether the momentum-space formula (2.8) remains valid for rational ε, or why these rational approximations are representative of the incommensurate limit.","section":"§2.2.2 / §3.2"},{"comment":"The constant C_{d,m} in (2.28) is not stated to depend on V and f; clarify the dependence explicitly.","section":"Theorem 2.5"},{"comment":"The Wigner transform formula introduces ℏ=1 without explaining the ε-scaling that is used in the subsequent localization argument; a short clarification would improve readability.","section":"§3.3, Wigner transform"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid formal core (Appendix A) and strong numerical agreement at small ε. The main risk is that the momentum-space reference is not demonstrably converged in the very regime used to validate the failure analysis; this is fixable and should be required. In particular, the negative values in Fig. 5 and the unreported Chebyshev order are substantial. The rational-ε issue and the missing derivation of (3.6) also need attention. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this as a careful numerical comparison of two DoS methods for a 1D incommensurate toy model, plus a physical explanation of the observed discrepancies. The main new content is real: the explicit second-order sum-over-states formulas L1 and L2 in Theorem 2.5 (Eqs. (2.26)–(2.27)) are not in the cited prior work, and the harmonic-oscillator model for the oscillations near band critical points is a genuinely cheap and partly predictive tool. The parameters in Table 1 come from the band structure, not from fitting the oscillations, and the agreement for p=2 is quite good. That is worth a serious referee's time.\n\nThe validation at large smearing (σ=0.4) is clean and impressive: errors down to 1e-14 for L0, 1e-11 for L1, 1e-6 for L2. These are not fitted numbers; the formulas are algebraically derived and the comparison uses an independent discretization. I agree with you that this part is solid.\n\nThe soft spots are real but not fatal. The most important is the one you flagged: in Fig. 5, at σ=0.08 and ε=0.01, the momentum-space ν_{ε,σ} goes negative. Since the exact quantity is a positive measure convolved with a Gaussian, negative values are a numerical artifact. The paper does not report the KPM expansion order, so we cannot tell how converged that reference is. That weakens the clean attribution of the ε=0.01 discrepancies to the truncated semiclassical expansion, and it also puts some pressure on the harmonic-model validation that uses the momentum-space curves as ground truth. But the issue is addressable: the authors could report the KPM order, show convergence in that order at σ=0.08, and explain the negative values. I would not call this a refutation; the broad outline (good agreement at small ε, oscillations at less small ε, harmonic explanation for small enough ε) is probably right.\n\nA second, more minor issue is that the error bound in Theorem 2.5 is technically vacuous for the small-σ Gaussians used here, because the derivative norms blow up as σ^{-(2m+6d+8)}. So the numerical agreement at σ=0.04–0.08 is empirical, not certified. The paper does not hide this, but it could say it more explicitly.\n\nFinally, the harmonic model is quantitative for only one of the six critical points (p=2); for p=1 it is only qualitative, and p=3–6 only for very small ε. The authors acknowledge this, and the Wigner-transform explanation for why is reasonable, but it does limit the scope of the headline claim.\n\nWho is this for? People working on semiclassical asymptotics for incommensurate systems, and experimental/theoretical groups that want a quick predictor for moiré flat-band-like features. It is a solid numerical paper with a useful new formula and an interesting physical hypothesis; it deserves peer review. My recommendation: send it out, but ask the authors to (1) report KPM order and convergence checks, especially at σ=0.08, (2) address the negative DoS directly, and (3) consider releasing the code and data. If those are fixed in revision, I would be happy to cite it.","headline":"Solid, genuinely useful numerical study with a real but addressable caveat: the momentum-space reference looks under-converged exactly where the paper's main discrepancy story plays out.","tokens_in":27179,"tokens_out":2228,"would_cite":true,"duration_ms":21066,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q20","35P20"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a 1D incommensurate Schrödinger operator, the density of states matches a second-order semiclassical expansion at small mismatch, and residual oscillations are harmonic-oscillator levels at band critical points.","keywords":["density of states","incommensurate systems","semiclassical expansion","moiré materials","Van Hove singularities","momentum-space method","harmonic oscillator approximation","Schrödinger operators"],"falsifier":"Recompute the momentum-space/KPM approximation for ε = 0.01, σ = 0.08 with a systematically increased Chebyshev expansion order (reporting the order) until the DoS approximation is nonnegative everywhere to numerical precision. If the converged reference still disagrees with the second-order semiclassical curve in the energy range [9,18], the paper's attribution of the mismatch to missing higher-order semiclassical corrections is supported; if the reference shifts, the mismatch is partly a reference artifact.","tokens_in":26010,"feed_emoji":"⚛️","tokens_out":6292,"duration_ms":54522,"temperature":0.7,"pith_summary":"This paper compares two computational approaches to the density of states (DoS) of a Schrödinger operator with two incommensurate periodicities—the setting behind moiré materials with a small twist or lattice mismatch. The first approach is a momentum-space planewave calculation; the second expands the DoS in powers of the small mismatch parameter ε. The paper supplies explicit formulas for the zeroth, first, and second-order terms of that expansion as integrals over the band structure of a commensurate two-parameter symbol, and verifies them against the momentum-space method. It then argues that when the full DoS curves disagree at larger ε, the disagreement is caused by oscillations at the semiclassical analogues of Van Hove singularities, and that these oscillations are described—quantitatively for small ε—by harmonic oscillator energy levels attached to each band critical point. If correct, this provides an inexpensive, parameter-free explanation of sharp oscillatory features in moiré spectra.","feed_headline":"Moiré density-of-states ripples are quantized oscillator levels","feed_subtitle":"A two-term semiclassical expansion matches the exact density of states at small mismatch; the leftover ripples fit harmonic energy levels.","key_machinery":"The central object is the operator-valued Weyl symbol h(k,X), the Bloch Hamiltonian of the commensurate approximation at local disregistry X. Its band structure (Ej(k,X)) and spectral matrix elements carry the argument: the semiclassical coefficients L0, L1, and L2 are integrals of divided differences of the test function over these bands, and the critical points of Ej drive an effective harmonic oscillator model Heff_p = E(p) − (Ap/2)d²/dx² + (ε²Cp/2)x². The Wigner transform of the oscillator eigenstates is used to diagnose whether the harmonic model is valid in a given region of phase space.","core_discovery":"For Hε = −½Δ + V(x,(1+ε)x) with smooth, periodic V, Theorem 2.5 establishes that Tr[f(Hε)] = L0(f) + εL1(f) + ε²L2(f) + O(ε³), where each Lj is an explicit sum-over-states integral over the eigenvalues λn(k,X) of the operator-valued symbol h(k,X) = ½(−i∇+k)² + V(x,X) on the torus, together with its spectral matrix elements. The paper checks these coefficients by finite differences against the momentum-space method and finds full-DoS agreement for ε = 0.001. At ε = 0.01, the mismatch is localized near critical points (k0,X0) of the band functions, and the paper shows that near such a point the effective Hamiltonian is a quantum harmonic oscillator with eigenvalues E(p) + εωp(n+½). This harmon","pith_inferences":["Editorial inference: In the ε = 0.01, σ = 0.08 comparison, the plotted momentum-space DoS dips below zero in the energy range [9,18], even though the exact Gaussian-smeared DoS must be nonnegative; this suggests the reference calculation is not fully converged there, so part of the reported mismatch may reflect reference error rather than missing semiclassical terms.","Editorial inference: The harmonic-oscillator mechanism offers a practical diagnostic: from the second derivatives of any band critical point one can predict where sharp DoS peaks will appear in a moiré spectrum, without solving the full aperiodic problem.","Editorial inference: In nearly degenerate band regions, such as the second/third band crossing near E ≈ 10 in this model, the single-band harmonic model breaks down even for small ε; a multi-band effective model would be a natural next test.","Editorial inference: Because L0, L1, and L2 depend only on the potential and not on ε, they could be reused as inexpensive surrogates in parameter sweeps over twist angles or lattice mismatches."],"forward_implications":["For sufficiently small ε, the DoS of an incommensurate system can be computed from the band structure of a single commensurate symbol, avoiding direct simulation of the aperiodic system.","The explicit formulas for L0, L1, and L2 make the semiclassical expansion directly usable once the spectral decomposition of h(k,X) is known.","The harmonic-oscillator analysis predicts that sharp DoS oscillations appear at energies determined by band critical points and their second derivatives, explaining what would otherwise look like numerical artifacts.","At less-small ε, missing higher-order semiclassical corrections are the reason the truncated expansion fails, pointing to where improved approximations are needed."],"fun_headline_variants":["Quantum oscillator explains moiré DoS ripples","DoS ripples in moiré are harmonic oscillator levels","Small-mismatch density ripples match oscillator energies","Moiré density ripples trace harmonic levels","Semiclassical ripples fit harmonic oscillator eigenvalues"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the momentum-space reference calculation is converged wherever it is used as ground truth—especially in the ε = 0.01, σ = 0.08 case, where the plotted exact DoS dips below zero, which the true Gaussian-smeared density of states cannot.","fun_headline_variants_meta":{"raw":{"variants":["Quantum oscillator explains moiré DoS ripples","DoS ripples in moiré are harmonic oscillator levels","Small-mismatch density ripples match oscillator energies","Moiré density ripples trace harmonic levels","Semiclassical ripples fit harmonic oscillator eigenvalues"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000276,"raw_usage":{"total_tokens":1516,"prompt_tokens":808,"completion_tokens":708,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":632}},"tokens_in":552,"tokens_out":708,"duration_ms":6938,"temperature":1.0,"reasoning_tokens":632,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T09:25:46.437092+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the momentum-space/KPM approximation for ε = 0.01, σ = 0.08 with a systematically increased Chebyshev expansion order (reporting the order) until the DoS approximation is nonnegative everywhere to numerical precision. If the converged reference still disagrees with the second-order semiclassical curve in the energy range [9,18], the paper's attribution of the mismatch to missing higher-order semiclassical corrections is supported; if the reference shifts, the mismatch is partly a reference artifact.","supporting_citations":[],"review_version":1}