{"id":"73e1ad90-5a10-4bcb-9446-72c18abe1e2f","arxiv_id":"2608.06260","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A phase-space electronic Hamiltonian for solids is proposed, together with a q-dependent Nafie equality, and an approximate Gamma operator is tested only on a 1D model.","lead":"The authors extend phase-space electronic structure theory from molecules to periodic solids, where electronic bands depend on both nuclear position and momentum. They prove a q-dependent version of Nafie's equality and propose a local Gamma operator intended to capture electron-phonon interactions without Berry curvature calculations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The practical claim rests on the ad hoc Gamma operator (Eq. 17), validated only by a q=0 sum rule and a 1D model with about 16% average error; its fidelity for finite q is unproven, so the claimed extraction of electron-phonon information without Berry curvature is not established.","rationale":"The reader's weakest_assumption and my concern coincide: the Gamma operator is the only bridge between the rigorous Nafie identity (which concerns exact d) and the practical band-structure calculation. The paper's Fig. 1 is the sole numerical test; an average 16% deviation already signals non-negligible approximation error, and the absence of q-resolved or sigma-dependent data leaves open the possibility that finite-q electron-phonon matrix elements are qualitatively wrong. The check I propose directly measures the approximation that must hold for the method to work, so it would settle the concern. I do not see a flaw in the formal proof; the issue is the validation of the practical operator, so a conditional verdict remains appropriate.","tokens_in":22961,"tokens_out":6136,"duration_ms":73872,"concrete_test":"In the 1D model of SM Section V, compute the exact derivative-coupling matrix elements d_{Aq}^{c,n}(k) from Eq. 7 (or by finite differences of the Bloch states) and the cell-periodic Gamma matrix elements of SM Eq. 150 on the same k,q grid; plot the q-resolved relative error ||Gamma-d||/||d|| averaged over k for sigma = 0.3, 0.5, 1, and 2 bohr. If the error exceeds about 20% for any q away from the zone center, or if it does not shrink as sigma is varied, then Gamma is not a faithful local proxy for the derivative coupling and the central extraction claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The formal q-dependent Nafie equality (Eq. 11) is derived for the exact derivative coupling d_{Aq}, and the SM proof appears internally consistent. The practical claim, however, requires that the Gaussian-partition-of-unity operator Gamma_{Aq} (Eq. 17) reproduce d_{Aq} well enough that Eq. 21 is a valid proxy for Eq. 12. The paper's only evidence is (i) the q=0 sum rule i hbar sum_A Gamma_{A,q=0}=p, which is satisfied by infinitely many operators and says nothing about q != 0, and (ii) a single 1D one-electron model (Fig. 1), whose BZ-averaged LHS/RHS ratio is 1.16 rather than 1, with no q-resolved worst-case or sigma convergence. The exact d_{Aq} is nonlocal in character (matrix elements of the potential response delta W_{Aq} divided by an energy denominator, Eq. 7), while Gamma_{Aq} is a local momentum weighted by a hand-chosen broadening sigma. There is no argument that the matrix elements of Gamma_{Aq} approximate the corresponding d matrix elements except in the q -> 0 limit. Since H_PS (Eq. 19) is built from Gamma and the extracted 'missing momentum' in Fig. 2 is directly proportional to Gamma matrix elements, an unvalidated Gamma makes the central claim of extracting electron-phonon couplings without Berry curvature unsupported for the finite-q phonons most relevant to solids.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a phase-space electronic structure approach for periodic solids, in which electronic bands depend on both nuclear position Q and nuclear momentum Π. The central formal result is a q-dependent analogue of Nafie's equality (Eq. 11), derived in the Supplementary Material, relating the change in a q-dependent electronic position expectation value with respect to a phonon coordinate Q_Aq to the change in electronic momentum with respect to the conjugate nuclear momentum Π_{A,-q}. The manuscript then introduces an approximate operator Γ_{Aq}(r) (Eq. 17), built from a Gaussian partition of unity, to replace the exact derivative coupling d_{Aq}, and uses it to define a phase-space Hamiltonian H_PS (Eq. 19). The approximation is tested in a one-dimensional one-electron model by comparing the exact and approximate sides of Nafie's equality (Fig. 1), and the H_PS Hamiltonian is used to compute band structures and nuclear-induced electronic momentum for three nuclear-motion scenarios (Fig. 2). The paper claims that this enables extracting electron-phonon information without Berry-curvature calculations.","tokens_in":23366,"tokens_out":5308,"duration_ms":68274,"significance":"If the formal q-dependent Nafie equality (Eq. 11) is correct, it is a valuable extension of molecular sum rules to periodic solids, connecting derivative couplings, Berry-curvature-like tensors, and momentum conservation. The proof in the Supplementary Material is detailed and appears internally consistent, and the symmetry relation in Eq. 15 is a clean formal statement. The practical claim, however, rests on the unproven faithfulness of the approximate operator Γ_{Aq}. The manuscript provides a machine-checkable-style derivation for the exact equality, a reproducible one-dimensional model, and explicit expressions for the Γ matrix elements; these are strengths. At the same time, the validation is limited to a single 1D one-electron model with a BZ-averaged ratio of 1.16 rather than 1, no q-resolved worst-case or σ-convergence study, and, as detailed below, the operator actually used in H_PS is not the operator validated in Fig. 1. The potential payoff—estimating non-adiabatic electronic momentum from a band-structure calculation without Berry curvature—is significant if the approximation can be established, but at present the evidence does not support the central applied claim.","major_comments":[{"comment":"There is a disconnect between the validation and the proposed practical Hamiltonian. The validation in Fig. 1 tests Nafie's equality using the finite-q operator Γ_{Aq} in Eq. 21, whereas the phase-space Hamiltonian H_PS in Eq. 19 contains only the q=0 operator Γ_{A,q=0}. Consequently, diagonalizing H_PS cannot produce finite-q electron-phonon matrix elements or q-dependent nuclear-induced electronic momentum. Fig. 2 varies only zone-center nuclear momenta P_A, so it does not demonstrate extraction of finite-q electron-phonon information. This undercuts the abstract's claim that electron-phonon interactions can be extracted from simple band structure calculations without Berry curvature, unless the scope is explicitly restricted to q=0.","section":"§4, Eqs. 19 and 21; Fig. 1"},{"comment":"The approximation Γ_{Aq} is an ad hoc local operator built from a Gaussian partition of unity with a hand-chosen broadening σ. The only exact property proven is the q=0 sum rule −iℏ Σ_A Γ_{A,q=0} = p (Eq. 159), which is satisfied by infinitely many operators and constrains only the total momentum, not the individual matrix elements ⟨u_{m,k+q}|Γ_{Aq}|u_{n,k}⟩ that enter the physical predictions. The numeric validation in Fig. 1 reports an average LHS/RHS ratio of 1.16 rather than 1, with no q-resolved worst-case analysis and no σ-convergence study. Thus the claim that Γ_{Aq} faithfully approximates the exact derivative coupling d_{Aq} at finite q is not established, and the central practical result—using Γ in place of d—is unsupported.","section":"Eq. 17 and SM §IV"},{"comment":"The exact LHS in Eq. 12 has a squared energy denominator, whereas the approximate RHS in Eq. 21 has only a linear denominator. These expressions coincide only if Γ_{Aq} equals the full derivative coupling d_{Aq}, including its energy-denominator structure. The test in Fig. 1 compares these two different expressions and reports an average deviation of about 16%; for near-degenerate bands the linear-versus-quadratic denominator difference will amplify any error in Γ. A direct comparison of the matrix elements of Γ_{Aq} and d_{Aq} (or a q-resolved and σ-converged test) is needed to support the approximation.","section":"Eq. 12 vs. Eq. 21; Fig. 1"},{"comment":"The formal derivation of Nafie's equality assumes non-degenerate bands and uses intermediate normalization to discard diagonal terms (SM Eqs. 63–69). Fig. 2, however, deliberately uses a metallic ground state with degenerate bands, where the energy denominators in Eq. 7 and Eq. 12 can vanish and the perturbation-theoretic derivation requires qualification. The paper does not explain how the exact equality or the approximate Γ-based test extends to the metallic case, so the physical interpretation of the band crossings and 'missing momentum' in Fig. 2 is not grounded in the derived formalism.","section":"SM §III.B and Fig. 2"}],"minor_comments":[{"comment":"The notation ∂_r in Eq. 20 is not defined; it should be written as ∂/∂r for clarity.","section":"Eq. 20"},{"comment":"Eq. 22 has an unmatched parenthesis in the denominator: '⟨ψ_{1k}(P=0)|P_A · Γ_A |ψ_{0k}(P=0⟩' should have the closing bracket after |ψ_{0k}(P=0)⟩.","section":"Eq. 22"},{"comment":"The caption says 'we plot the ratio of the LHS in Eq. 12 and RHS in Eq. 21' but the figure actually plots the ratio of the LHS to the RHS; please rephrase to avoid ambiguity.","section":"Fig. 1 caption"},{"comment":"There is a typo in the section heading: 'eletronic' should be 'electronic'.","section":"§4 heading"},{"comment":"The sentence 'as q increases, electronic orbitals increases' is unclear because q does not appear in the Hamiltonian H_PS of Eq. 19; please clarify what quantity is plotted as a function of q in Fig. 2.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":"The exact formal result in Eq. 11 and the derivation in the SM are likely the strongest contribution. The practical claim, however, is not yet supported: the approximate operator Γ is introduced ad hoc, the validation is limited to a single 1D model with ~16% average error, and the operator in the proposed H_PS (q=0) is not the one validated for finite q. These are fixable with additional numerical evidence (e.g., a q-resolved, σ-converged test in 1D, or better a 2D/3D model), but as written the Letter overstates what is established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: the formal part is in good shape, the practical part is not. Bradbury and Subotnik prove a q-dependent version of Nafie's equality for periodic solids, and the proof in the SM is detailed and internally consistent. That is a real new result. They also show that the equality can be recast as a symmetry of a mixed nuclear-electronic curvature tensor, which is elegant and should be useful.\n\nWhat the paper does well: it takes the molecular phase-space program and extends it to Bloch states with a clean Wigner-Weyl treatment. The definitions of q-dependent expectation values and the derivation of Eq. 11 are careful. The reduction of both sides to the same Berry-curvature-like expression is demonstrated. The SM is self-contained and the notation is manageable. The 1D numerical test at least shows the idea is not vacuous.\n\nWhere it softens: the approximate operator Gamma, Eq. 17, is a heuristic built from a Gaussian partition of unity with a hand-chosen width sigma. The paper validates it with the q=0 sum rule and a single 1D one-electron model. The BZ-averaged ratio is 1.16, not 1, and there is no q-resolved worst-case analysis or sigma convergence study. Gamma is local and momentum-weighted; the exact derivative coupling is nonlocal, with an energy denominator. Nothing in the paper shows that Gamma reproduces d_{Aq} matrix elements at finite q, which is exactly the regime where the claimed extraction of electron-phonon information is supposed to work. So the central practical claim, that you can get electron-phonon momentum from band structures without Berry curvature, is not yet supported.\n\nI disagree with the stress-test note's implied dismissal: the formal proof is not circular, and the authors are transparent that Gamma is an ansatz. But the note is right that the validation is too thin for the weight of the claim.\n\nWho this is for: researchers working on electron-phonon couplings, non-adiabatic effects in solids, and phase-space electronic structure. They will get a new formal relation and a promising but unproven route to avoid Berry curvature computations. It deserves a serious referee, but the referee should require benchmarks against exact derivative couplings in a small or exactly solvable system, sigma and system-size convergence, and at least a 2D demonstration before the practical claim is accepted.\n\nRecommendation: send to peer review, major revision.","headline":"Solid formal result on a q-dependent Nafie equality for periodic solids, but the practical claim of extracting electron-phonon information without Berry curvature rests on an unvalidated local approximation.","tokens_in":23827,"tokens_out":2363,"would_cite":false,"duration_ms":25192,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.38.-k","71.15.-m"],"model":"deepseek-v4-flash","headline":"This paper proves a q-dependent Nafie equality for periodic solids and argues that a phase-space Hamiltonian with a local approximate coupling operator can extract nuclear-induced electronic momentum from plain band-structure…","keywords":["phase space electronic structure","non-Born-Oppenheimer","Nafie equality","electron-phonon coupling","nuclear momentum","band structure","Berry curvature-free","Wigner-Weyl transform"],"falsifier":"In a three-dimensional solid with a known DFPT implementation, evaluate the exact left side of Eq. (11) from the derivative couplings and compare it with the right side computed through $\\Gamma$ at $\\sigma\\approx 1$ bohr across the full Brillouin zone; a ratio far from unity for optical or zone-boundary phonons, or a strong dependence on $\\sigma$, would falsify the claim that $\\Gamma$ replaces the derivative coupling.","tokens_in":22749,"feed_emoji":"⚛️","tokens_out":8663,"duration_ms":85336,"temperature":0.7,"pith_summary":"The paper argues that the standard frozen-nucleus picture of solids silently drops the momentum nuclei impart to electrons, and that a phase-space electronic Hamiltonian, with bands depending on nuclear momentum as well as position, can restore it. Its central result is a proof that a wavevector-resolved version of Nafie's molecular equality holds in periodic systems, linking the q-dependent change of electronic position under a phonon displacement to the change of electronic momentum under the conjugate nuclear momentum. The paper then proposes to approximate the derivative couplings by a local operator built from a Gaussian partition of unity, so that this momentum transfer can be read off by diagonalizing a momentum-dependent Hamiltonian at ordinary band-structure cost, without computing Berry curvature. If correct, this opens a practical route to electron-phonon and inertial effects, such as band renormalization and momentum-induced band crossings, in ab initio solids. A one-dimensional semiconductor model supports the equality with best agreement near q=0, and a metallic two-site model shows nuclear momentum shifting, and sometimes splitting, the bands.","feed_headline":"Band structure alone can reveal phonon-driven electron momentum","feed_subtitle":"A q-resolved Nafie equality ties electron momentum to phonon motion without Berry curvature.","key_machinery":"The argument runs on two linked objects. First, the q-dependent Nafie equality, written as a symmetry of a mixed nuclear-electronic curvature tensor, $$\\$\\Omega$^q_{nk}(Q_{Aq},\\tilde{k}) = -\\frac{M_A\\hbar}{m_e}\\$\\Omega$^q_{nk}(\\Pi_{A,-q},r),$$ which is equivalent to Eq. (11) and is exact when the electronic states are exact. Second, the phase-space Hamiltonian $$$H^{{PS}}$_k(R,P) = \\sum_A \\frac{(P_A - i\\hbar \\Gamma_{A,q=0,k})^2}{2M_A} + \\frac{(\\hat p+\\hbar k)^2}{2m_e}+V(r,R),$$ where $$\\Gamma_{Aq}(r) = \\frac{1}{2i\\hbar}\\sum_\\ell \\left[\\theta_{A\\ell}(r)$e^{{iR_\\ell\\cdot q}}$\\hat p + \\hat p $e^{{iR_\\ell\\cdot q}}$\\theta_{A\\ell}(r) - \\hbar q $e^{{iR_\\ell\\cdot q}}$\\theta_{A\\ell}(r)\\right]$$ is built from a normalized Gaussian partition of unity $\\theta_{A\\ell}$ of width $\\sigma$. The Gamma operator does the load-bearing work: it replaces the exact derivative coupling $d_{Aq}$ from the Born-Huang expansion, obeys the same sum rule $i\\hbar\\sum_A\\Gamma_{A,q=0}=\\hat p$ and the same crystal-momentum selection rule $\\delta_{k',k+q}$, and lets the electronic momentum be obtained by diagonalizing a momentum-dependent Hamiltonian rather than by computing Berry curvature.","core_discovery":"Within a Wigner-Weyl treatment of nuclear degrees of freedom, the authors prove that periodic solids satisfy a wavevector-resolved analogue of Nafie's equality: $$\\frac{\\partial\\langle \\psi_{nk}|\\hat r|\\psi_{nk}\\rangle_q}{\\partial Q_{Aq}} = \\frac{M_A}{m_e}\\frac{\\partial\\langle\\psi_{nk}|\\hat p|\\psi_{nk}\\rangle_q}{\\partial \\Pi_{A,-q}},$$ evaluated at $\\Pi=0$, with the left side expressible in closed form through the standard DFPT perturbation $\\delta W_{Aq}$ and band-structure ingredients. They then claim this relation supports a practical shortcut: replace the exact derivative coupling with the local operator $\\Gamma_{Aq}$ of Eq. (17), so that the nuclear-induced electronic momentum, the piece usually obtained from Berry curvature, comes out of diagonalizing the phase-space Hamiltonian Eq. (19). In a one-dimensional semiconductor model the approximate equality holds with a Brillouin-zone-averaged ratio of 1.16, best near $q=0$, and in a two-site metallic model nuclear momentum shifts the bands, creates crossings, and produces nonzero electronic momentum only when the nuclear motion breaks symmetry.","pith_inferences":["The authors leave implicit that the $\\sigma$-independence of $\\Gamma$-based results could serve as a practical convergence test: if nuclear-induced momentum changes materially when the Gaussian width is varied, the approximation is not converged for that phonon.","If the method transfers to two and three dimensions, electron-phonon inertial effects could become a post-processing step on ordinary band-structure data, skipping response-theory calculations for quantities like the atomic polar tensor.","Because $\\Gamma$ has no energy denominator, it is formally better behaved than exact derivative couplings at band crossings, though the same feature makes its justification via perturbation theory weakest precisely in degenerate systems."],"forward_implications":["Diagonalizing $H^{PS}$ gives both band energies and nuclear-induced electronic momentum for finite $P$, so momentum transfer between nuclei and electrons no longer requires a separate Berry curvature calculation.","Classical dynamics on a single phase-space eigensurface conserve total momentum, in contrast to conventional Born-Oppenheimer dynamics.","For exact derivative couplings, the q-dependent Nafie equality holds exactly; with $\\Gamma$ in place of $d$, the one-dimensional semiconductor model recovers the equality with an average ratio of 1.16 and near-unity quality for q close to zero.","In the metallic two-site model, an optical phonon produces no electronic momentum, an acoustic phonon shifts the Brillouin zone uniformly, and a symmetry-breaking single-atom motion yields an avoided crossing, showing that nuclear momentum can re-shape band structure."],"supporting_citations":[{"why":"Supplies the original molecular Nafie equality that the paper generalizes to periodic solids.","marker":"[31]"},{"why":"Defines the phonon coordinates $Q_{Aq}$, the perturbation $\\delta W_{Aq}$, and the DFPT machinery used throughout the derivation.","marker":"[25]"},{"why":"Provides the modern theory of polarization mapping $\\hat r \\to i\\partial_k$ that lets the left side of the q-dependent equality be evaluated.","marker":"[27]"},{"why":"Supplies the standard electron-phonon matrix element framework and the relation of derivative couplings to electron-phonon couplings.","marker":"[18]"},{"why":"Identifies the q=0 form of the central expression as the atomic polar tensor and Berry curvature expression the paper avoids.","marker":"[28]"},{"why":"Inspires the phase-space, Wigner-transformed nuclear electronic structure ansatz that the paper adapts to periodic boundary conditions.","marker":"[20]"}],"fun_headline_variants":["Phonon-driven electron momentum without Berry curvature","Band structure alone yields phonon-induced electron momentum","New equality links electron momentum to phonons, no Berry curvature","Solid-state shortcut: nuclear-induced electron momentum from bands","Wavevector-resolved Nafie equality unlocks phonon momentum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the local operator $\\Gamma$, built by splitting space with Gaussian blobs of width $\\sigma$, faithfully reproduces the true derivative coupling for real phonons; the paper establishes only the q=0 sum rule and one one-dimensional semiconductor check.","fun_headline_variants_meta":{"raw":{"variants":["Phonon-driven electron momentum without Berry curvature","Band structure alone yields phonon-induced electron momentum","New equality links electron momentum to phonons, no Berry curvature","Solid-state shortcut: nuclear-induced electron momentum from bands","Wavevector-resolved Nafie equality unlocks phonon momentum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000199,"raw_usage":{"total_tokens":1363,"prompt_tokens":929,"completion_tokens":434,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":357}},"tokens_in":545,"tokens_out":434,"duration_ms":4596,"temperature":1.0,"reasoning_tokens":357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:11:19.433446+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a three-dimensional solid with a known DFPT implementation, evaluate the exact left side of Eq. (11) from the derivative couplings and compare it with the right side computed through $\\Gamma$ at $\\sigma\\approx 1$ bohr across the full Brillouin zone; a ratio far from unity for optical or zone-boundary phonons, or a strong dependence on $\\sigma$, would falsify the claim that $\\Gamma$ replaces the derivative coupling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the original molecular Nafie equality that the paper generalizes to periodic solids."},{"cited_title":"Baroni, S","cited_arxiv_id":null,"evidence_quote":"Defines the phonon coordinates $Q_{Aq}$, the perturbation $\\delta W_{Aq}$, and the DFPT machinery used throughout the derivation."},{"cited_title":"Vanderbilt,Berry Phases in Electronic Structure Theory: Electric Polarization, Orbital Magnetiza- tion and Topological Insulators(Cambridge Univer- sity Press, 2018)","cited_arxiv_id":null,"evidence_quote":"Provides the modern theory of polarization mapping $\\hat r \\to i\\partial_k$ that lets the left side of the q-dependent equality be evaluated."},{"cited_title":"Shenvi, Phase-space surface hopping: Nonadia- batic dynamics in a superadiabatic basis, J","cited_arxiv_id":null,"evidence_quote":"Inspires the phase-space, Wigner-transformed nuclear electronic structure ansatz that the paper adapts to periodic boundary conditions."}],"review_version":1}