REVIEW 2 major objections 3 minor 1 cited by
On the resolvent convergence of discrete Dirac operators on 3D cubic lattices
T0 review · 2 major / 3 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Lattice Dirac resolvents converge strongly, not in norm
desk verdict Clean 3D proof of strong-but-not-norm resolvent convergence for unmodified lattice Dirac operators with an explicit lower bound; the printed second lower bound has a typo-level error that is easy to fix. 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 proof rests on the step-function embedding $J_h$ that maps lattice functions into $L^2(\mathbb{R}^3)$, its orthogonal projection $P_h$, and the discrete Fourier transform $F_h$ on the torus $[-\pi/h,\pi/h]^3$. Working in the Fourier representation turns the discrete Dirac operator into a multiplication operator whose symbol $\widehat{D}_{m,h}(\xi)$ approaches the continuum symbol $\widehat{D}_m(\xi)$ pointwise, but not uniformly. The failure of uniform convergence is detected by a family of Gaussian test functions $u_h$ concentrated near the boundary of the Brillouin zone; explicit asymptotic estimates show that the projected vector $(D_{m,h}-z)^{-1}P_h u_h$ has norm approaching $\|u_h\|/|m-z|$, proving the lower bound.
What would settle it
Numerically compute the operator norm of $(D_{m,h}\oplus 0_h-z)^{-1}-(D_m-z)^{-1}$ on a finite but large cubic lattice for a fixed $m$ and non-real $z$, taking $h$ smaller at each step; if the norm is ever found to drop below $\max(1/|m-z|,1/|m+z|)$ by a non-vanishing amount, the lower bound of Theorem 2.1(ii) is false.
Extended reading notes
Core claim
The central result is Theorem 2.1: for any non-real $z$ in the resolvent set, the embedded discrete resolvent $(D_{m,h}\oplus 0_h-z)^{-1}$ converges strongly to the continuum resolvent $(D_m-z)^{-1}$ in $L^2(\mathbb{R}^3)^4$ as $h\to 0$, yet the operator norm of their difference satisfies $\liminf_{h\to 0}\|(D_{m,h}\oplus 0_h-z)^{-1}-(D_m-z)^{-1}\| \ge \max(1/|m-z|,1/|m+z|)$. The lower bound is constructed by test functions whose Fourier transform is localised near the corners of the Brillouin zone, where the discrete symbol differs from the continuum symbol by an error that does not vanish in norm.
Load-bearing premise
The proof of strong convergence assumes, without re-proving, that three key lemmas on the discrete Fourier transform and the projection $P_h$ from the two-dimensional paper [3] remain valid verbatim in three dimensions; if any of these lemmas fails in $\mathbb{R}^3$, the strong convergence argument has a gap.
Editorial extensions
If this is right
- Strong resolvent convergence implies that the spectra of the embedded discrete operators converge to the spectrum of the continuum Dirac operator in the Hausdorff distance as $h\to 0$.
- Because norm resolvent convergence fails, functions of the resolvent such as spectral projections and semigroups do not converge uniformly, so uniform approximation of time evolution on the whole space is not guaranteed.
- The explicit lower bound $\max(1/|m-z|,1/|m+z|)$ quantifies the unavoidable error near the mass-shell singularities of the continuum resolvent.
- The paper's remark that central difference operators do not change the fundamental issue indicates that simply altering the finite-difference scheme cannot restore norm convergence.
Reading between the lines
- The construction of the test functions suggests that the non-convergence is a high-frequency phenomenon at the edges of the Brillouin zone, so any discretisation that damps or filters these high frequencies (for example, using a smoother embedding or a non-uniform lattice) might achieve norm resolvent convergence.
- The proof technique, relying only on the symbol comparison and projection estimates, may extend to other lattice geometries and to higher dimensions, yielding the same strong-but-not-norm dichotomy.
- The lower bound depends on the specific step-function embedding; it would be informative to test whether other natural embeddings that also 'preserve the discrete operator' give different limiting behaviour.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the continuum limit of the free Dirac operator discretized on a three-dimensional cubic lattice with forward and backward difference operators. The discrete operator is embedded into L2(R^3)^4 by extending lattice functions to step functions, and the main theorem (Theorem 2.1) states that the resolvent of the embedded discrete Dirac operator, (D_{m,h} ⊕ 0_h - z)^{-1}, converges strongly to (D_m - z)^{-1} as h → 0, but not in the norm resolvent sense, with an explicit lower bound liminf_{h→0} ||(D_{m,h} ⊕ 0_h - z)^{-1} - (D_m - z)^{-1}||_B ≥ max(1/|m-z|, 1/|m+z|). The proof uses the discrete Fourier transform on the torus T^3_{1/h}, compares the discrete and continuum Fourier symbols, and constructs oscillatory Gaussian test functions to establish both the strong convergence and the lower bound on the norm discrepancy.
Significance. If the result is correct, it provides a clean three-dimensional analogue of the authors' two-dimensional result and sharpens the known picture by giving an explicit quantitative obstruction to norm resolvent convergence for the natural step-function embedding. The paper contains detailed, checkable estimates, including explicit constants in the lower bound, and it explains why modifications such as those in [1] and [2] are needed to achieve norm convergence. These are genuine strengths. The main limitations are the reliance on unstated three-dimensional versions of technical lemmas from the authors' prior two-dimensional paper [3], and a local but load-bearing error in the second lower bound of Theorem 2.1(ii). Both are repairable without changing the main claim.
major comments (2)
- [Section 2, proof of Theorem 2.1(ii)] The claimed second lower bound 1/|m+z| is not proved by the stated 'analogously' argument with u_h = (0,0,y_h,0)^T. The Fourier transform of y_h is concentrated at ξ0 = (π/(2h), −π/(2h), 0); at this point, using (2.4) and the matrix (1.2), bD_{m,h}(ξ0)e_3 = (0, 2(1−i)/h, −m, 0)^T, because i∂*_1 − ∂*_2 has symbol 2(1−i)/h there. Thus e_3 is not an eigenvector of the discrete symbol, and the O(1/h) off-diagonal entry prevents (bD_{m,h}−z)^{-1}P_h u_h from being close to (−m−z)^{-1}P_h u_h. The correct vector is u_h = (0,0,0,y_h)^T: at the same ξ0, i∂*_1 + ∂*_2 has symbol 0 and bD_{m,h}(ξ0)e_4 = −m e_4. With this replacement, the Step 3 estimates carry over and give the claimed lower bound. Please make the replacement explicit and repeat the Step 1–3 estimates for the fourth component.
- [Section 2, equations (2.9)–(2.12)] The strong convergence argument uses Lemmas 3.3, 3.5 and 3.6 of [3], which are stated in [3] for two-dimensional square lattices. The current paper is three-dimensional, and these lemmas control the strong convergence of P_h to the identity, of F_h to F, and of the difference F_h − F on Schwartz functions. Since these properties are load-bearing for the convergence proof in (2.9)–(2.12), the authors should either state the three-dimensional versions with proofs or explicitly indicate that the proofs in [3] are dimension-independent and carry over verbatim to R^3.
minor comments (3)
- [Section 2, after equation (2.17)] The displayed formula for the norm of (8/π^2)(y_h)_h has a typo: the leading constant should be 8/(π^{5/4} 2^{3/4}), not 8/π^{5/4} 2^{3/4}; the printed version differs by a factor of 2^{3/2}. This does not affect the subsequent argument, which only needs the norm to tend to a nonzero constant.
- [Section 2, final displayed estimate of Theorem 2.1(ii)] In the line containing '∥(Dmh − z)^{-1}Ph uh∥', the operator 'Dmh' should be 'D_{m,h}'.
- [Section 2, equation (2.9)] The notation (bD_{m,h} − z)^{-1} ⊕ 0_h is used after the Fourier extension of F_h is introduced, but the extension of 0_h to the Fourier side is not defined. Please define the action of this operator on L2(T^3_{1/h}) and its orthogonal complement explicitly.
Circularity Check
No significant circularity: the strong and norm resolvent convergence results are proved by explicit Fourier analysis and test-vector estimates, not by defining the conclusion into the hypotheses.
full rationale
Theorems 2.1(i)-(ii) are derived from direct computations with the discrete Fourier transform, the explicit symbol bD_{m,h}(xi), and carefully normalized test functions; no fitted parameter is renamed as a prediction and no uniqueness theorem is invoked to force a choice. The self-citations to the authors' 2D paper [3] (Lemmas 3.3, 3.5, 3.6 and the projection formula) provide auxiliary Fourier-transform and projection facts whose stated assumptions do not include the 3D resolvent-convergence conclusion, so they are independent technical support rather than circular premises. The paper's claim that the second lower bound in Theorem 2.1(ii) follows 'analogously' from u_h=(0,0,y_h,0)^T may be questionable (the natural eigenvector at the concentration point is e_4 rather than e_3), but that is a correctness or repairability issue in a proof detail, not a circular definitional or self-citation reduction. The central derivation chain is self-contained apart from cited technical lemmas, so the circularity score is 0.
Assumptions & free parameters
assumptions (5)
- standard math Self-adjointness of D_m,h and the resolvent bound ||(D_m,h-z)^{-1}|| ≤ 1/|Im z|
- standard math Density of Schwartz space S(R^3)^4 in L^2(R^3)^4
- domain assumption Lemmas 3.3, 3.5 and 3.6 of [3] on convergence of Ph, F_h and Fh
- standard math The matrix identity (bD_m(ξ)-z)^{-1} = (|ξ|^2+m^2-z^2)^{-1}(bD_m(ξ)+z)
- domain assumption The step-function embedding Jh is the chosen identification of lattice and continuum Hilbert spaces
Cite this review
Pith. "Pith review of On the resolvent convergence of discrete Dirac operators on 3D cubic lattices." pith.science (2026). https://pith.science/paper/PZYXAQUS
@misc{pith2026250701443,
author = {Pith},
title = {Pith review of: On the resolvent convergence of discrete Dirac operators on 3D cubic lattices},
year = {2026},
howpublished = {\url{https://pith.science/paper/PZYXAQUS}},
note = {Machine review of arXiv:2507.01443}
}
read the original abstract
We prove that the discrete Dirac operators in three dimensions converge to the continuum Dirac operators in the strong resolvent sense, but not in the norm resolvent sense.
Forward citations
Cited by 1 Pith paper
-
Continuum limit of discretized matrix-valued Fourier multipliers
Block-matrix Fourier multipliers admit h^{min(2,2γ−β−1)} generalized norm-resolvent convergence after a Wilson-type correction term, while the uncorrected symmetric-difference scheme converges only strongly.
Reference graph
Works this paper leans on
-
[1]
H. Cornean, H. Garde, A. Jensen Discrete approximations to Dirac operators and norm resolvent convergence.J. Spect. Theory 12 (2022), 1589–1622
work page 2022
-
[2]
Nakamura Remarks on discrete Dirac operators and their continuum limits.J
S. Nakamura Remarks on discrete Dirac operators and their continuum limits.J. Spect. Theory 14 (2024) 255–269
work page 2024
-
[3]
K.M. Schmidt, T. Umeda, Continuum limits for discrete Dirac operators on 2D square lat- tices. Analysis and Mathematical Physics 13 (2023), art. no. 46 Karl Michael Schmidt: School of Mathematics, Cardiff University, Senghennydd Road, Cardiff CF24 4AG, W ales, UK Email address: schmidtkm@cardiff.ac.uk Tomio Umeda: Department of Mathematical Sciences, Univ...
work page 2023
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.