{"id":"a4203f2e-163b-494c-bf79-22f3008d6b6d","arxiv_id":"2607.26008","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Strain-controlled loops around tensor-space degeneracies produce a measurable sign-reversal parity of the EFG principal axis in TiO2, SnO2, and MgO.","lead":"This paper shows that the electric-field-gradient tensor at a single atomic site can be steered through a closed strain loop, making a principal axis return flipped — a binary topological invariant called return parity. It gives first-principles evidence that rutile TiO2, SnO2, and MgO realize these singular geometries and that existing hyperfine spectroscopy could measure them.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"TiO2 degeneracy is supported only by quadratic fits; direct DFT at Q* and convergence checks are needed to rule out a fitting artifact.","rationale":"The reader's verdict of CONDITIONAL rests on the same weakest assumption: the PBE DFT EFG tensors and their quadratic strain fits are quantitatively faithful, and the common zero of dz and dx at Q* is a real isolated degeneracy. My stress-test reading independently identifies this as the single most load-bearing concern. The paper provides no direct DFT calculation at Q*, no analysis of fit residuals, no convergence tests with respect to cutoff, k-mesh, or pseudopotential choice, and no deposited data. The loop parity demonstration is mathematically robust but is executed on interpolated/fitted tensors, so it does not by itself certify the degeneracy. This concern is addressable: direct DFT at Q*, raw-data loop transport, and pseudopotential/convergence checks would settle it. Since the reader already assigned CONDITIONAL and my analysis does not change that assessment, I recommend UNCHANGED rather than ACCEPT, REJECT, or UNVERDICTED. There is no evidence of internal mathematical inconsistency or fraud; the issue is missing empirical verification of a key input to an otherwise coherent argument.","tokens_in":21132,"tokens_out":6198,"duration_ms":95964,"concrete_test":"Run self-consistent PBE DFT at the claimed degeneracy point Q*=(0,-0.0120387) and at 16-32 points on a small loop around it, computing dz and dx directly from the raw EFG tensors. Accept the claim only if |dz| and |dx| are below ~1e-4 of the tensor norm and the loop winding computed from the raw DFT tensors gives w1=1. Then repeat with a different pseudopotential (e.g., PAW without semicore, or an all-electron method) and with doubled cutoff/k-mesh; if the degeneracy disappears or shifts by more than ~0.1% strain, the claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central material claim is the isolated degeneracy in rutile TiO2 at Q*=(0, -1.204%) in the QB1g-QB2g strain plane (Methods E4, Eq. M47). This is established solely from a six-parameter quadratic fit of each Cartesian EFG component (Eq. M46). The Jacobian determinant in Eqs. M48-M49 is computed from the fitted polynomial derivatives, not from direct numerical differentiation of DFT data. If the true strain response contains cubic or higher-order terms—plausible at ~1.2% strain—the common zero of dz and dx could be an artifact: the fitted model can vanish where the actual DFT tensor remains gapped. The loop parity w1=1 (Fig. 3c) is then obtained by transporting eigenvectors along loops built from the same fitted tensors, so it does not independently confirm the degeneracy. No convergence tests (cutoff, k-mesh, PAW vs all-electron), no error bars, and no deposited data are provided; the EFG is known to be highly sensitive to the near-nucleus density (Ref. 48). The mathematical framework is sound, but the load-bearing empirical fact—that a true degeneracy exists and is not a fitting artifact—is not yet verified beyond the model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a spectral-geometric and topological framework for real symmetric traceless rank-2 tensor observables, applied to electric-field gradients at a single probe site as functions of symmetry-adapted strain. It introduces invariant coordinates I2, I3 and the normalized shape parameter p=√6 I3/I2^{3/2}, shows that the standard magnitude-ordered (|V33|,η) chart folds at the degeneracy stratum Σ0 (η=0) and at the determinant-zero ordering seam Δ0 (η=1), and defines a Z2 return parity w1 for a transported principal axis around a loop, identified with the first Stiefel–Whitney invariant. Using PBE/PAW DFT, the authors report an isolated in-plane EFG degeneracy in rutile TiO2 at approximately (0,−1.204%) in the B1g–B2g strain plane, with w1=1 for encircling loops; in SnO2 an isolated degeneracy plus a symmetry-protected extended degeneracy branch; and in cubic MgO a full-rank 5×5 linear response matrix allowing local control of the traceless EFG. They propose orientation-resolved single-crystal TDPAC/NMR/Mössbauer measurements to observe the return parity.","tokens_in":21509,"tokens_out":17013,"duration_ms":234638,"significance":"I found the mathematical core sound and the exposition mostly clear. The invariant parametrization and the two-level class-AI reduction with winding/transport are correct as far as I checked; the connection to the Czjzek–Evenson unfolding is a nice addition. The novelty lies in applying tensor-field eigenframe topology to a local observable under external strain control, rather than to a field over physical space. The proposed spectroscopic route is falsifiable. There is no circularity: the classification is applied to independently computed DFT tensors and the invariant/winding identities used are standard. The main weakness is empirical verification: the central TiO2 claim (and the analogous SnO2 point degeneracy) rests on a quadratic fit without residuals or direct DFT at the degeneracy, no convergence tests are reported, and the data are not yet available. These are fixable within a revision. If the degeneracies are confirmed, the paper would be a valuable demonstration of a long-known topological structure in a new, experimentally accessible setting.","major_comments":[{"comment":"The isolated TiO2 degeneracy at Q*≈(0,−1.204%) is the load-bearing fact for the w1=1 loop in Fig. 3. It is established only by six-parameter quadratic fits (M46) to the Cartesian EFG components, with no reported residuals, sample-point distribution, strain range, or fit quality. The Jacobian (M48) and its determinant (M49) are computed from the fitted polynomial derivatives, not from direct numerical differentiation of DFT data. At ~1.2% strain, cubic or higher-order terms are plausible, so the common zero of dz and dx could be an artifact of second-order truncation. Please add direct DFT calculations at Q* and at nearby points bracketing the zero, provide fit residuals/cross-validation and an uncertainty estimate for Q*, and report convergence tests with respect to cutoff, k-mesh, PAW potential (Ti_sv vs a harder/all-electron treatment), and LASPH/ADDGRID settings for TiO2 and SnO2. The","section":"Section 4 / Methods E4, Eqs. M46–M49"},{"comment":"The text states that the transported arrows are shown 'directly' in Fig. 3c, but it is not specified whether the loop data are DFT calculations at discrete points or evaluations of the fitted quadratic model. This matters because only the former independently confirms the degeneracy. Please clarify, report the number of loop points and the discretization error, and give the numerical values of ⃗v2(θ)·⃗v2(0) at the final point for C0 and C1, including any uncertainty. The same clarification is needed for the SnO2 loop in Fig. 5.","section":"Section 4, Fig. 3c,d"},{"comment":"The extended λ1=λ2 degeneracy branch in SnO2 is convincingly explained by the xy⊕z block structure preserved by the B1g/B2g controls, so that part is robust. However, the isolated point degeneracy in SnO2 is presumably located by the same quadratic-fit procedure as TiO2 and needs the same verification: residuals, direct DFT at/near the point, and convergence checks. The caption of Fig. 5 states that the isolated point acts as a winding centre; please confirm that the winding/parity computation uses DFT tensors rather than the fitted surface.","section":"Section 5, Figs. 4 and 5"},{"comment":"The claim of complete local control of Sym0(3,R) in MgO rests on the statement that the 5×5 linear response matrix is full rank to numerical accuracy. The matrix, its singular values/condition number, and the rank tolerance are not reported. Please provide the matrix or its singular values (or the two Schur-lemma scalars for the Eg and T2g blocks) and state the numerical precision. This is needed to support the headline result that all five EFG components are accessible.","section":"Section 6, Fig. 6d"}],"minor_comments":[{"comment":"'insufficient' should be 'insufficient'.","section":"Section 3, page 6"},{"comment":"'choosen' should be 'chosen'; also define θ in the caption.","section":"Fig. 3c caption"},{"comment":"Please specify the units of the Jacobian entries in Eq. (M48) and the strain sampling range used for the fits.","section":"Methods E4"},{"comment":"LASPH and ADDGRID are specified for MgO but not for TiO2/SnO2; clarify whether they were used, since EFG output is sensitive to these settings (Ref. 48).","section":"Methods E2/E3"},{"comment":"The statement that data 'will be deposited' is prospective; for a computational paper with load-bearing numerical claims, providing the key data (or a DOI) during review would allow verification.","section":"Data availability"},{"comment":"The first-order Taylor expansion is not used in the subsequent analysis; consider removing it or connecting it to the fits.","section":"Section 4, Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The main risk is verification, not circularity or scholarship. I do not see citation problems; the self-citation [16] is clearly a related in-press paper and not used as evidence for the topological claims. I recommend major revision because the requested checks are standard and within scope. The paper should not be rejected on the basis of the mathematical framework, which is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper is worth reading for its clean assembly of eigenframe topology for local rank-2 tensors and for switching the problem from physical-space tensor fields to control-space loops, with EFG as the observable. That framing and the DFT-driven examples are genuinely new. But the central empirical claim — an isolated Σ0 degeneracy in rutile TiO2 with w1=1 — is supported only by quadratic fits to DFT data, not by a direct calculation at the degeneracy point, and the loop parity is computed from the same fitted field. Right now that's a model-based prediction, not a demonstrated material fact.\n\nWhat it does well: the math is standard (degeneracy stratum Σ0, ordering seam Δ0, Stiefel-Whitney parity) but the paper's presentation is careful and the connection to reporting conventions (|V33|, η) is well explained. The Czjzek/Evenson prior work is credited properly, and the MgO 5×5 full-rank response is a nice, clean result. The proposal to measure return parity via single-crystal TDPAC or Mössbauer is concrete, and the SnO2 point-vs-line degeneracy distinction is a good illustration.\n\nThe soft spots are in proportion: the TiO2 degeneracy is the load-bearing piece, and Methods E4 confirms the stress-test concern. The quadratic fit (Eq. M46) is used to locate Q*; the Jacobian in Eq. M48 is the derivative of the fitted polynomial, not of the DFT data. No residuals, no convergence tests, no direct DFT at Q*, and the data are not yet deposited. EFG is known to be sensitive to near-nucleus density, so a pseudopotential or cutoff effect could plausibly move the degeneracy. None of this is fatal to the framework, and the authors may well be right — the zero is at a symmetry-allowed location and the block structure makes a crossing natural. But the paper currently asserts a first-principles discovery that is actually a fit-extrapolated prediction. That gap is addressable: compute EFG directly at Q* and a small neighborhood, check cutoff and k-mesh convergence, deposit the data. The loop parity should then be recomputed from the raw data rather than the polynomial.\n\nFor peer review: yes, send it out. A competent referee can force the needed checks; the framework and proposed experiment deserve attention. But I would not cite the TiO2 holonomy as established fact until the direct calculation appears.","headline":"Good framework, standard math, but the key TiO2 degeneracy is only a fitted zero — needs direct DFT and deposited data before the central claim lands.","tokens_in":21981,"tokens_out":3559,"would_cite":false,"duration_ms":49321,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The electric-field-gradient tensor at a probe site can be steered by strain around an isolated degeneracy in tensor space; a transported principal axis returns reversed, defining a Z2 return parity that hyperfine spectroscopy can in princip","keywords":["electric field gradient","eigenframe topology","Stiefel-Whitney invariant","tensor singularities","strain control","hyperfine spectroscopy","return parity","rutile TiO2"],"falsifier":"Calculate the EFG on a denser strain grid around Q*=(0,−1.204%) in rutile, or with higher-order strain expansions, and check whether the common zero of dz and dx persists with a nonzero Jacobian. Alternatively, run a strain-resolved single-crystal hyperfine experiment (for example, 119mSn-based TDPAC in SnO2 or a comparable probe in rutile) and reconstruct the principal-axis projector around the predicted encircling loop: observing w1=0 where the paper predicts w1=1 would falsify the central claim.","tokens_in":21092,"feed_emoji":"🔄","tokens_out":4387,"duration_ms":66876,"temperature":0.7,"pith_summary":"This paper argues that the standard way of reporting symmetric rank-2 tensor observables—through magnitude-ordered principal values and axes—folds tensor space and hides a genuinely global effect. For the electric-field-gradient (EFG) tensor at a single site, the authors show by first-principles calculations that symmetry-adapted strain can steer the tensor through its five-dimensional space, and that a closed strain loop encircling an isolated degeneracy makes a continuously tracked principal axis return reversed. That reversal is a binary return parity, the first Stiefel–Whitney invariant, and it cannot be changed by smooth deformations of the loop as long as no degeneracy is crossed. The paper establishes rutile TiO2 as a prototype with an isolated degeneracy and nontrivial loop parity, shows SnO2 has point- or line-like degeneracies depending on the control slice, and demonstrates that deviatoric strain modes in MgO give full local control of all five EFG components. If the argument is right, topological structure familiar from tensor fields over space becomes a measurable property of a single local observable under external control.","feed_headline":"Strain loop flips a tensor axis in rutile","feed_subtitle":"A site's EFG axis returns reversed after a closed strain loop, a topological invariant hyperfine spectroscopy can reconstruct.","key_machinery":"The central object is the EFG tensor as a traceless real symmetric 3×3 matrix, viewed as a map from symmetry-adapted strain amplitudes into the five-dimensional tensor space Sym0(3,R). When two principal values approach each other, the problem reduces to a real 2×2 block with coordinates (dz,dx) in a Pauli basis; the winding of this two-vector around the degeneracy point defines the Z2 Stiefel–Whitney return parity w1. This two-level real-symmetric reduction, together with the invariant shape parameter p=√6 I3/I2^{3/2}, carries the argument: p=±1 marks the degeneracy stratum Σ0 and p=0 marks the ordering seam Δ0, separating local chart-wall artifacts from global eigenframe holonomy.","core_discovery":"The central claim is that the eigenframe of a local rank-2 tensor observable carries a Z2 topological charge in control-parameter space, and this charge is physically accessible. For the traceless EFG tensor at a probe nucleus, the degeneracy stratum Σ0 (where two principal values coincide, η=0) is the organizing singularity: a loop in strain space that links Σ0 makes a chosen arrow on a principal axis return reversed, w1=1, even though the tensor is smooth and nondegenerate on the loop. In rutile TiO2, first-principles strain trajectories find an isolated transverse crossing of Σ0 at the strain point Q*≈(0,−1.204%), with local conical-intersection structure of the principal-value branches a","pith_inferences":["The same strain-controlled eigenframe holonomy should appear in other quadrupolar probe nuclei and crystals wherever an isolated EFG degeneracy is encircled; the predicted w1=1 is a concrete signature to search for in existing hyperfine datasets.","Because w1 is invariant under smooth loop deformation, it provides a noise-robust topological classification of strain paths: any experimentally realized loop that keeps an eigenvalue gap and links Σ0 should exhibit the same sign reversal, even if the loop geometry drifts.","The MgO full-rank control result suggests a broader principle: when the symmetry-adapted strain basis and the traceless-tensor basis match, strain becomes a complete local handle on the observable's tensor space, which could be used to engineer loops with prescribed winding rather than relying on accidental degeneracies.","The connection to real symmetric two-band topology hints that local tensor observables could serve as a classical, tabletop realization of Stiefel–Whitney-type physics, distinct from electronic band theory and not requiring spatially extended states."],"forward_implications":["For rutile TiO2, a closed strain loop in the QB1g–QB2g plane that encircles the isolated degeneracy must return the tracked principal axis reversed, with w1=1, independent of the loop's precise shape or radius.","For SnO2, an isolated degeneracy point supports a gapped encircling loop with nontrivial parity, whereas the extended degeneracy branch is symmetry-protected within the chosen slice and requires an additional Eg-type strain mode to lift the degeneracy and permit a transverse gapped loop.","For MgO, the full-rank linear response means arbitrary local EFG tensors can be generated by deviatoric strain, so singular features and control loops can be designed rather than only discovered.","The invariant parameters I2, I3, and p give a chart-independent description of the EFG spectrum, so the cusp-like features seen in (|V33|,η) can be recognized as crossing or approaching the two distinguished loci Σ0 and Δ0.","Orientation-resolved hyperfine methods—single-crystal TDPAC, NMR rotation patterns, or Mössbauer line-intensity analysis—can in principle reconstruct the principal-axis projectors needed to observe the return parity."],"fun_headline_variants":["Strain loop flips principal axis in rutile","Tensor eigenframe has topological parity","Closed strain loop reveals tensor flip","EFG axis returns reversed after strain loop"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The DFT-computed EFG tensors and their quadratic strain fits are faithful: the common zero of dz and dx at Q* must be a true isolated degeneracy of the tensor response, not an artifact of second-order strain truncation or of the pseudopotential choice.","fun_headline_variants_meta":{"raw":{"variants":["Strain loop flips principal axis in rutile","Tensor eigenframe has topological parity","Closed strain loop reveals tensor flip","EFG axis returns reversed after strain loop"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000165,"raw_usage":{"total_tokens":1102,"prompt_tokens":776,"completion_tokens":326,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":273}},"tokens_in":520,"tokens_out":326,"duration_ms":5635,"temperature":1.0,"reasoning_tokens":273,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T03:15:01.583269+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Calculate the EFG on a denser strain grid around Q*=(0,−1.204%) in rutile, or with higher-order strain expansions, and check whether the common zero of dz and dx persists with a nonzero Jacobian. Alternatively, run a strain-resolved single-crystal hyperfine experiment (for example, 119mSn-based TDPAC in SnO2 or a comparable probe in rutile) and reconstruct the principal-axis projector around the predicted encircling loop: observing w1=0 where the paper predicts w1=1 would falsify the central claim.","supporting_citations":[],"review_version":2}