{"id":"02d89118-473a-4bc4-94d1-313859884e4f","arxiv_id":"2511.22706","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Strains in two orthogonal directions produce dissipationless momentum flow in the third direction, with magnitude set by multiband quantum geometry and integer topological invariants.","lead":"A theoretical paper predicts that certain three-dimensional magnetic or quantum-geometric crystals, when squeezed in two directions at once, develop a dissipationless flow of momentum in the third direction. The effect is a nonlinear, odd 'viscoelastic' response whose strength can be tuned by integer topological invariants, offering a mechanical probe of quantum geometry.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"O(w^2) truncation of the strained Hamiltonian is unresolved and could invalidate Eqs. (4)-(5) as the full second-order odd response; the Delta/ln Delta issue is a chain-rule artifact.","rationale":"The reader's leading concern is the strongest one: the truncation of H(w) at linear order in strain is precisely the point on which the physical interpretation of η as the response to real biaxial strain depends, and the manuscript contains an explicit sentence saying that higher-order derivative terms survive. This is unresolved and must be fixed before Eqs. (4)-(5) can be taken as the complete second-order odd viscoelastic response. However, the reader's secondary point about Delta_ba in Eq. (4) versus ln Delta_ba in SM Eq. (32) is not a defect, since they are equivalent by the chain rule. I found no more severe internal inconsistency: the Kubo derivation is standard in structure, the two- and three-band decompositions are gauge-invariant, and the model calculations agree with the symmetry analysis. The abstract's 'nonmetricity' promise is undefined, but it is not central to the tensor formulas. Thus the appropriate outcome remains the same conditional verdict: the truncation issue must be settled, the nonmetricity wording clarified, and the numerics made reproducible, while the core proposal remains plausible if the check above comes out favourably.","tokens_in":19001,"tokens_out":14776,"duration_ms":143247,"concrete_test":"For the perturbed MRW Hopf model (End Matter), expand H(w) exactly to second order in w, extract U_ij,kl = ∂^2H/∂w_ij∂w_kl|_0, and compute the first-order-in-U contribution to the odd combination defined in Eq. (12) (the T-U retarded susceptibility). If this contribution is nonzero for the η_xx;yy,zz component, Eqs. (4)-(5) are incomplete; if it vanishes identically under the odd projector, or is canceled by a Ward identity, the truncation is vindicated. A practical version: evaluate the full finite-strain response of the model at small w by exact diagonalization and compare its antisymmetrized part with the truncated Kubo result on a 32^3 k-grid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the replacement of the strained Hamiltonian by H' = -T_ij w_ij, with T_ij = ∂H/∂w_ij|_0. For a second-order response, the quadratic part U_ij,kl = ∂^2H/∂w_ij∂w_kl|_0 can contribute at first order in perturbation theory and hence at O(w^2) to <T_mn>, unless it is annihilated by the odd projector. The main-text Discussion explicitly states that 'the higher-order derivative terms survive in the w_ab -> 0 limit', while SM Sec. II only justifies the form of the first-derivative stress operator at w=0; it does not show that the quadratic derivatives cannot enter the antisymmetrized second-order response. If U survives the odd combination in Eq. (12), the numerical magnitudes in Figs. 2-3 and the identification of η with the response to biaxial strain are incomplete. The reader's secondary 'Delta vs ln Delta' discrepancy is not real: ∂ ln Δ = (∂Δ)/Δ, so main-text Eq. (4) and SM Eq. (32) coincide. The abstract's 'nonmetricity' is undefined elsewhere, but it is not load-bearing for Eqs. (4)-(5).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a class of nonlinear odd viscoelastic effects in three-dimensional electronic systems. It defines a second-order stress-stress response tensor η_{ij;kl,mn} via a nested Kubo-like commutator of stress operators, claims the response is dissipationless, and derives a momentum-resolved decomposition into a two-band term built from symplectic Berry-type forms and a three-band term built from three-state quantum geometric tensors. The effect is demonstrated numerically in a perturbed Moore-Ren-Wen Hopf insulator, where the response grows with the Hopf invariant under k_z rescaling, and in a flattened chiral three-band model, where it tracks the chiral invariant and the three-state geometry. A real-space Wannier-center calculation is presented as a complementary check. The central physical claim is that biaxial strain produces a transverse momentum current whose magnitude is controlled by quantum geometry and integer topological invariants.","tokens_in":19331,"tokens_out":14793,"duration_ms":129232,"significance":"If the computed tensor is the complete second-order odd viscoelastic response, the paper introduces a genuinely new class of dissipationless mechanical responses and connects them to multiband quantum geometry and topological invariants in a falsifiable way. The strengths of the manuscript are its detailed Kubo derivation in the Supplemental Material, the gauge-invariant two-band/three-band decomposition, the independent numerical computation of Hopf and chiral invariants, and the complementary Wannier-center calculation. The proposed effect is experimentally addressable in magnetic topological insulators. However, the completeness of Eqs. (4)-(5) is conditional on resolving the strain-Hamiltonian truncation issue, so the significance is currently prospective rather than established.","major_comments":[{"comment":"The response is computed after truncating the strained Hamiltonian to H(w) ≈ H(0) - w_ab T_ab, with T_ab = ∂H/∂w_ab|0. For a second-order response this truncation is not automatically justified. The quadratic terms (1/2)w_ab w_cd U_ab,cd, where U_ab,cd = ∂²H/∂w_ab∂w_cd|0, contribute to ⟨T_ij⟩ at first order in perturbation theory, i.e., at O(w²), through kernels of the form ⟨[T_ij, U_ab,cd]⟩; additionally, the strain dependence T_ij(w) = T_ij + w_cd U_ij,cd produces mixed first-order contributions. SM Sec. II only proves that the first strain derivative at w=0 equals the stress tensor; it does not show that U is annihilated by the odd projector in SM Eq. (12). The Discussion's statement that 'the higher-order derivative terms survive in the w_ab → 0 limit' directly contradicts the preceding assertion that they can be neglected. Unless a proof of cancellation is supplied, or the U terms a","section":"Discussion (strain-implementation paragraph); SM Sec. II"},{"comment":"Eq. (1) defines η_{ij;kl,mn} as a single nested commutator, (i/ω)⟨[T_mn,[T_kl,T_ij]]⟩. The derivation in the SM instead defines the odd response as the antisymmetrized combination σ^odd_{A;B1,B2} = (2i/ω)χ_{A;B1,B2} - (i/ω)χ_{B1;A,B2} - (i/ω)χ_{B2;B1,A} in SM Eq. (13). These two objects are not equal; the Jacobi identity does not reduce the antisymmetrized combination to the single commutator. Consequently, Eq. (1) is not the quantity actually computed in Eqs. (4)-(5), and the relation J_ij = η w_kl w_mn in Eq. (2) does not follow from Eq. (1) as written. Please replace Eq. (1) by the SM definition, or explicitly state that Eq. (1) is a schematic representative of the antisymmetrized correlator.","section":"Main text Eq. (1); SM Eqs. (11)-(13)"}],"minor_comments":[{"comment":"The two-band formula is written as (1/Δ_ba) Ω (∂Δ_ba) in the main text and as Ω ∂ ln Δ_ba in the SM. These are equivalent by the chain rule, but the notation should be harmonized to avoid apparent inconsistency.","section":"Main text Eq. (4); SM Eq. (32)"},{"comment":"The abstract mentions 'nonmetricity,' but this term is never defined or used in the main text or SM. Please either define it and connect it to the formulas, or remove it.","section":"Abstract"},{"comment":"The definition of Q^{abc} differs in band-index order: main text has A^{ac}A^{cb}A^{ba}, while the SM has A^{ab}A^{bc}A^{ca}. Since a real part is taken, the expressions may coincide up to conjugation and relabeling, but the notation should be made consistent.","section":"Main text Eq. (5); SM Eq. (26)"},{"comment":"The caption lists χ values (0, p, −2p, ...) but does not explain how the plotted curves correspond to those values. Please add explicit labels per curve or panel.","section":"Fig. 2 caption"},{"comment":"The quantity d x_c/dw_kl is a static strain derivative of a Wannier center, not a time-dependent current. The text calls it a 'net current' without deriving the relation to the momentum current J_ij. Please clarify this connection.","section":"Nonlinear viscoelasticity in Wannier basis"}],"recommendation":"major_revision","confidential_remarks":"The strain-truncation issue is the main technical risk: the Discussion itself acknowledges that higher-order derivative terms survive, which undermines the claim that Eqs. (4)-(5) are complete. This is potentially fixable by proving that the quadratic strain terms are annihilated by the odd projector, or by including them, but it must be addressed explicitly. The Eq. (1) definition mismatch is a separate but easily fixable issue. I do not recommend rejection if the authors can resolve these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core claim is new and worth taking seriously: a nonlinear (second-order) odd viscoelastic tensor in 3D, decomposed into two-band and three-band geometric pieces, with model calculations showing scaling with Hopf and chiral invariants. The mapping to nonlinear Hall physics is apt, and the symmetry analysis is careful — they identify which symmetries kill the response and check the same numerically. The sum rule and the equivalence between normal-stress and shear-stress components are also nice, concrete results.\n\nThe soft spots are real but addressable. The main one is the treatment of O(w^2) terms in the strained Hamiltonian. The Kubo formula they use assumes the perturbation is linear in strain, H' = -T_ij w_ij. But the full H(w) also has a quadratic part, and that quadratic part contributes to the second-order response at first order in perturbation theory. The main-text Discussion says these higher-derivative terms 'survive in the w_ab -> 0 limit' while the SM claims they vanish because w=0 at the point where T_ij is defined. That is not a resolution — it only fixes the stress operator, not the perturbation Hamiltonian. The authors need to compute or bound the quadratic-derivative contribution, or show it is annihilated by the odd projector. As written, Eqs. (4)-(5) may be incomplete.\n\nThe other issue flagged in the reading note is the Delta vs ln Delta discrepancy between main-text Eq. (4) and SM Eq. (32). I checked: these are equivalent by the chain rule, since d(ln Δ) = dΔ/Δ. Not a real problem. The abstract's 'nonmetricity' is never defined in the main text, which is careless but minor. The numerical work is standard and re-implementable; no code is shipped, but the models are simple enough.\n\nThis paper deserves a serious referee. The effect is new, the derivation is serious, and the models are concrete. But the truncation issue is load-bearing, so I would not accept it as is. Send it to review, and let the referees press on the O(w^2) contribution.","headline":"A genuinely new nonlinear odd viscosity tensor with topological fingerprints, but the second-order strain truncation needs to be fixed before the central formulas can be trusted.","tokens_in":19824,"tokens_out":3687,"would_cite":false,"duration_ms":35218,"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 paper predicts that biaxial strain on a three-dimensional topological crystal produces a dissipationless momentum current in the third direction, controlled by quantum geometry and integer topological invariants.","keywords":["nonlinear odd viscosity","viscoelasticity","quantum geometry","Hopf invariant","three-state quantum geometric tensor","strain response","topological phase transition","momentum current"],"falsifier":"Compute the second-order stress-stress correlator using the exact strain-dependent Hamiltonian including the O(w^2) terms: if the result differs from Eqs. (4)–(5), the predicted topological scaling is incomplete. Alternatively, apply biaxial strain to a candidate Hopf insulator and look for a transverse momentum current scaling linearly with the Hopf invariant.","tokens_in":18875,"feed_emoji":"🌀","tokens_out":8051,"duration_ms":61559,"temperature":0.7,"pith_summary":"The paper predicts a class of nonlinear odd viscoelastic effects in three-dimensional crystalline electron fluids: when a crystal is strained in two orthogonal directions, a dissipationless momentum current flows in the third direction. The central claim is that this second-order stress-stress response is entirely governed by the Hilbert-space geometry of the occupied Bloch states, decomposing into a two-band term built from symplectic Berry-type connection forms and a three-band term built from three-state quantum geometric tensors. The effect scales with integer topological invariants — the Hopf invariant in two-band models, and a chiral invariant in three-band models — and changes sharply across topological phase transitions. A sympathetic reader would care because the effect offers a bulk experimental fingerprint of multiband quantum geometry and of topological invariants that are otherwise hard to access, and it extends odd-viscosity physics from linear to nonlinear order.","feed_headline":"Biaxial strain drives a momentum flow in the third direction","feed_subtitle":"Quantum geometry and topological invariants set the strength of a new nonlinear odd viscosity in 3D crystals.","key_machinery":"The objects doing the work are strain-dressed geometric tensors: the non-Abelian Berry connection A^{ab}_{ij} = i⟨u_a|∂_{w_ij}u_b⟩, the symplectic form Ω^{ab}_{ij,kl} = -2 Im(A^{ab}_{ij}A^{ba}_{kl}), and the three-state quantum geometric tensor Q^{abc}_{ij,kl,mn} = A^{ab}_{ij}A^{bc}_{kl}A^{ca}_{mn}. The two-band response is essentially Σ (1/Δ_{ba}) Ω^{ba}_{ij,(kl}∂_{mn)}Δ_{ba}, and the three-band response is Σ (Δ_{cb}/Δ_{ba}Δ_{ac}) Re[Δ_{cb}Q^{abc} + Δ_{ba}Q^{cab} + Δ_{ac}Q^{bca}], with normalized symmetrization of index pairs. These tensors convert geometric phase structure of the Bloch states into a prediction for the nonlinear momentum current, and their symmetry properties determine when","core_discovery":"The paper's central claim is that the dc nonlinear viscoelastic tensor η_{ij;kl,mn} — the second-order response of a momentum current to two static strains — is nonzero in three-dimensional magnetic and geometrically nontrivial phases, and decomposes over the Brillouin zone into two gauge-invariant pieces. The two-band piece is built from the strain-dressed symplectic Berry form Ω^{ba}_{ij,kl} weighted by the inverse gap; the three-band piece is built from products of three non-Abelian Berry connections forming the three-state quantum geometric tensor Q^{abc}. Numerically, the two-band response grows in integer steps as the Hopf invariant is scaled in a two-band Hopf insulator, and the three","pith_inferences":["Editorial: the same two-band/three-band decomposition likely extends to higher-order viscoelastic responses, where the paper's linear-in-strain truncation of the Hamiltonian would break down; those higher orders are where the dropped O(w^2) terms might enter.","Editorial: a biaxial-strain experiment on a candidate Hopf insulator could serve as a bulk probe of the Hopf invariant itself, complementing surface-state probes, provided the O(w^2) terms do not spoil the scaling.","Editorial: because the three-state quantum geometric tensor also controls other second-order geometric responses (shift-current-like and photogalvanic), the NOVE may share a common geometric origin with those transport channels, suggesting correlated signatures.","Editorial: a clean internal check is to recompute η in the same toy models using the full strain-dependent Hamiltonian, including the O(w^2) terms the paper drops; if the result changes, the topological scaling law would need refinement."],"forward_implications":["Magnetic topological insulators in three dimensions, including candidate axion insulators, should exhibit a nonlinear momentum current under biaxial strain, with the response strength tied to integer topological invariants.","The effect gives a bulk experimental fingerprint of multiband Hilbert-space geometry — specifically the skewness and chirality of Wannier functions encoded in three-state quantum geometric tensors.","Because the three-band term survives in flattened bands and after the protecting symmetry is broken, it separates geometric contributions from dissipative ones, allowing a clean readout of quantum geometry.","The tensor obeys the sum rule η_xx;yy,zz + η_yy;zz,xx + η_zz;xx,yy = 0, so only two independent components need to be measured to characterize the effect.","The formal mapping between nonlinear viscoelastic and nonlinear Hall responses means biaxial strain protocols can be adapted from existing Hall viscosity experiments, with the current response deviating from linear scaling in the presence of nonzero correlators."],"fun_headline_variants":["Quantum geometry bends momentum into a third axis","Topological invariants amplify nonlinear odd viscosity","Two strains, one flow: quantum geometry directs momentum","Dissipationless odd viscoelasticity from quantum geometry"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole derivation rests on taking the stress operator as T_ij = ∂H/∂w_ij evaluated at zero strain, thereby dropping O(w^2) terms in the strain-dependent Hamiltonian; the paper's own Discussion states these higher-order derivative terms survive in the w→0 limit, which, if true, would make Eqs. (4)–(5) incomplete at second order.","fun_headline_variants_meta":{"raw":{"variants":["Quantum geometry bends momentum into a third axis","Topological invariants amplify nonlinear odd viscosity","Two strains, one flow: quantum geometry directs momentum","Dissipationless odd viscoelasticity from quantum geometry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001949,"raw_usage":{"total_tokens":7407,"prompt_tokens":640,"completion_tokens":6767,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":384,"completion_tokens_details":{"reasoning_tokens":6707}},"tokens_in":384,"tokens_out":6767,"duration_ms":44810,"temperature":1.0,"reasoning_tokens":6707,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T19:43:19.110952+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the second-order stress-stress correlator using the exact strain-dependent Hamiltonian including the O(w^2) terms: if the result differs from Eqs. (4)–(5), the predicted topological scaling is incomplete. Alternatively, apply biaxial strain to a candidate Hopf insulator and look for a transverse momentum current scaling linearly with the Hopf invariant.","supporting_citations":[],"review_version":1}