{"id":"390d1645-acf1-4d43-a9be-0d9e151da2bf","arxiv_id":"1908.08335","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Pyroxene chains with comparable direct and oxygen-mediated hopping can host a spin-orbital liquid with low-energy spectral weight, and NaRuSi2O6 is proposed as a concrete candidate.","lead":"This paper studies magnetic chain materials called pyroxenes and shows when their spin and orbital excitations can stay mobile instead of freezing into ordered patterns. It identifies the compound NaRuSi2O6 as a promising candidate for a quantum liquid with strong orbital fluctuations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ru/Os pyroxenes have strong spin-orbit coupling, absent from the model; with λ_SOC≈0.1 eV (Ru) and larger (Os), the low-energy t2g^5 space collapses to a J_eff=1/2 Kramers doublet, so the proposed NaRuSi2O6 orbital-liquid claim is not supported.","rationale":"The single most load-bearing condition for the central claim is that the candidate material's low-energy Hilbert space is the six-dimensional spin-times-t2g manifold of Eq. (3). For NaRuSi2O6 this fails before any question of structural decoupling or the value of Δ, because Ru3+ has strong spin-orbit coupling that is absent from the model. The reader's weakest assumption about ideal orbital relabeling and free Δ is valid and contributes to uncertainty, but it is a quantitative idealization within the same Hilbert space; SOC changes the Hilbert space itself. I therefore retain the CONDITIONAL verdict rather than rejecting the paper, since the abstract mechanism for Ti pyroxenes and the SU(6) model remains internally coherent, and the SOC issue could in principle be addressed by adding the term and re-evaluating the candidate.","tokens_in":11930,"tokens_out":16279,"duration_ms":173454,"concrete_test":"Perform DFT+U including spin-orbit coupling for NaRuSi2O6, or add an on-site λ L·S term to the Wannier-interpolated t2g Hamiltonian, and compute the single-ion t2g^5 spectrum. Then compare the J_eff=1/2-to-J_eff=3/2 splitting with J≈2t1^2/U≈10 meV and with the orbital gaps of Eq. (3). If the SOC-induced gap is several times J, the six-dimensional spin-orbital Hilbert space used in the DMRG is not the low-energy sector of the proposed material, so the central material claim would need to be replaced by a spin-orbit-entangled description rather than an SU(6) orbital liquid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (3) and the DMRG phase diagram assume a six-dimensional local Hilbert space with independent SU(2)-spin and SU(3)-type orbital operators and no on-site L·S term. The central material claim extends this to Ru3+/Os3+ pyroxenes, but those ions have large spin-orbit coupling. For octahedral t2g^5, L·S splits the six states into a J_eff=1/2 Kramers doublet and a J_eff=3/2 quartet separated by about 3λ/2; with λ_Ru≈0.1 eV and λ_Os≈0.3–0.5 eV, this splitting is at least 0.15 eV. The reported NaRu hopping parameters are t1=0.132 eV and t2=0.085 eV, and the exchange scale is J≈2t1^2/U≈10 meV. Thus SOC is comparable to or larger than the bandwidth and an order of magnitude larger than J; orbital degrees of freedom are quenched at low energies instead of being liberated. The first-principles estimates in Table I were obtained without SOC, so even t2/t1=0.64 and the orbital mixing are not reliable low-energy parameters for 4d/5d pyroxenes. The general t2/t1 mechanism may survive for 3d Ti, but the specific NaRuSi2O6 candidate is not supported by the presented model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies quasi-one-dimensional pyroxene Mott insulators with one electron or hole in the t2g orbitals. The authors construct a three-orbital Hubbard model, relabel orbitals along the zigzag chain to decouple the kinetic terms into three identical channels with direct hopping t1 and oxygen-mediated hopping t2, and derive an effective spin-orbital exchange Hamiltonian (Eq. (3)). Using DMRG, they compute the dynamical spin structure factor and find that for moderate crystal-field splitting Δ and t2/t1~1 the spectral weight shifts to low energies, indicating liberated orbital fluctuations, whereas large Δ yields dimerized spin singlets. They report DFT/Wannier parameters: NaTiSi2O6 has t2/t1=0.21 and NaRuSi2O6 (hypothetical, using the NaTi structure) has t2/t1=0.64, and they suggest NaRuSi2O6 as a candidate for a liquid with tightly bound spin and orbital excitations.","tokens_in":12295,"tokens_out":7371,"duration_ms":76127,"significance":"The paper contains a solid analysis of the three-orbital Hubbard chain: the SU(6)-symmetric Sutherland point is exactly solvable, the conformal-field-theory scaling dimension of the 4kF perturbation is computed (5/3), and the DMRG spectral functions are obtained with stated truncation errors. The identification of the t2/t1 ratio as a control parameter for orbital liberation is physically transparent, and the DFT/Wannier hopping parameters are a constructive step toward material-specific modeling. However, the material-specific conclusion rests on the absence of spin-orbit coupling, which is not justifiable for the proposed 4d/5d t2g^5 ions; this undermines the NaRuSi2O6 candidate and, with it, the paper's central claim as stated.","major_comments":[{"comment":"The six-dimensional local Hilbert space of Eq. (3) contains independent spin and orbital operators and no L·S term. For a t2g^5 ion in an octahedral environment, spin-orbit coupling λ splits the six states into a J_eff=1/2 Kramers doublet and a J_eff=3/2 quartet separated by approximately 3λ/2. With λ_Ru≈0.1 eV and λ_Os≈0.3–0.5 eV, this splitting is at least 0.15 eV, comparable to the hoppings t1=0.132 eV and t2=0.085 eV in Table I and an order of magnitude larger than the exchange scale J≈2t1^2/U≈10 meV. The orbital degrees of freedom are therefore quenched at low energies instead of being liberated. The GGA/Wannier parameters, obtained without SOC, are not valid low-energy parameters for Ru/Os pyroxenes. The suggestion that NaRuSi2O6 is a liquid with tightly bound spin and orbital excitations is not supported by the presented model.","section":"Main text, final paragraph; Supplemental Table I and Fig. 7"},{"comment":"The predicted low-energy spectral weight and the crossover to the orbital-liquid regime depend on Δ/J, yet no independent estimate of Δ is provided for NaRuSi2O6. The DFT band structures in Supplemental Fig. 7 show the yz orbital moving higher in energy, but the actual Δ entering Eq. (3) is not quantified. As a result, the material candidate is placed in the desired regime by a free parameter, and the central material prediction is not falsifiable as presented. The authors should either compute or constrain Δ, or specify a measurable prediction (for example, the field or temperature dependence of the spin gap) that distinguishes the liquid regime from the dimerized regime.","section":"Main text, 'the value of Δ is uncertain and is taken as a free parameter'; Figs. 2 and 4"},{"comment":"The exact decoupling into three degenerate hopping paths assumes an ideal undistorted edge-sharing geometry. Real pyroxene chains have octahedral tilts and, in NaTiSi2O6, dimerization; Table I itself contains nonzero inter-orbital hoppings, e.g., the yz-zx element of 0.0204 eV, that are discarded by the relabeling. The statement that the band is exactly three times degenerate is therefore an idealization whose accuracy for real compounds is not quantified. This weakens the quantitative mapping of t2/t1 and Δ from DFT to the model, although it does not invalidate the model as a theoretical construct.","section":"Supplemental Eq. (8) and Fig. 3(c)"}],"minor_comments":[{"comment":"The PACS numbers '74.72.-h, 74.72. Gh' contain a typographical error (a space before 'Gh') and appear to be incorrect; the authors should verify them.","section":"PACS numbers"},{"comment":"The caption does not specify the values of t2/t1 and Δ/J for every panel (a)–(g); please add this information so the reader can connect each panel to the parameter regimes discussed in the text.","section":"Figure 2"},{"comment":"In the limit of large Δ, the sentence 'For J = 0 each site has two degenerate orbitals' is confusing because J elsewhere denotes the exchange constant; the intended meaning (vanishing exchange coupling) should be stated explicitly.","section":"Large-Δ discussion"},{"comment":"Equation (3) and the subsequent δVk expression are difficult to follow because the relation between the relabeled orbitals a, b, c and the original xy, yz, zx orbitals is only given through Fig. 3; a short table or explicit relabeling definition in the text would improve readability.","section":"Equation (3)"},{"comment":"The abstract and introduction refer to 'the pyroxene family' broadly, but only NaTiSi2O6 and a hypothetical NaRuSi2O6 are discussed; a sentence clarifying the intended scope would avoid overgeneralization.","section":"Abstract and introduction"},{"comment":"The comparison of the DMRG spectral weight at about J/2 to the 10 meV heat-capacity gap uses J=2t1^2/U fixed by the 53 meV spin gap; this is a consistency check rather than a fitted prediction, and stating this explicitly would increase clarity.","section":"Heat-capacity comparison"}],"recommendation":"reject","confidential_remarks":"The theoretical analysis of the three-orbital chain is competent and the DMRG results are useful, but the paper's central material claim is contradicted by the basic spin-orbit physics of 4d/5d t2g^5 ions. If the authors were to resubmit a version that either includes SOC in the model or restricts the claim to the model itself rather than to NaRuSi2O6/NaOsSi2O6, the contribution could be publishable; as it stands, the proposed material candidate is not supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The tripled-degeneracy relabeling is genuinely elegant: for the ideal zigzag chain, the three t2g orbitals can be redefined so the kinetic term decouples into three identical hopping paths, giving an SU(6)-symmetric Sutherland point when t1=t2 and Delta=0. That is a fresh observation and it makes the paper worth reading. The DMRG spectral functions and the flow from gapless to dimerized behavior as Delta/J grows are carefully presented, and the DFT Wannier parameters for NaTiSi2O6 (t2/t1=0.21) and NaRuSi2O6 (0.64) are useful numbers.\n\nThe soft spots are concentrated in the material claim. The model has no spin-orbit coupling. For a t2g^5 ion in an octahedral field, SOC splits the six states into a J=1/2 Kramers doublet and a J=3/2 quartet with a separation of order 3 lambda/2. With lambda_Ru ~0.1 eV, that is about 0.15 eV, an order of magnitude bigger than the estimated exchange scale J ~10 meV. So the low-energy Hilbert space is not the six-dimensional one the paper uses; the orbital fluctuations are quenched. The DFT hoppings were computed without SOC, so t2/t1=0.64 is not a reliable low-energy parameter for Ru. The same objection holds for Os, where lambda is larger. The general t2-mechanism might survive in 3d Ti systems with weak SOC, but the NaTiSi2O6 numbers give t2/t1=0.21, which the paper itself says is too small to liberate orbitals. So the specific NaRuSi2O6 candidate is not supported by the current model.\n\nOther issues are minor. Delta is a free parameter, and the heat capacity comparison is a post-hoc consistency check rather than a prediction. The Hund's coupling is set to zero; less important for one hole, but it does affect the exchange anisotropy. The ideal-geometry decoupling is broken by real tilts and dimerization, and the supplement has a handful of proofreading errors.\n\nWho is this for? People working on SU(N) magnetism, DMRG spectral functions, and pyroxene materials. The tripled-degeneracy mapping deserves serious refereeing and could be cited. The material prediction should be revised to account for SOC or restricted to systems where SOC is negligible. I'd recommend sending it out, with a request to address the SOC issue head-on.","headline":"A clean SU(6) mapping for pyroxene chains, but the Ru-based material claim needs SOC input before it can fly.","tokens_in":12790,"tokens_out":3262,"would_cite":true,"duration_ms":31989,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.72.-h","74.72.Gh"],"model":"deepseek-v4-flash","headline":"Pyroxene chains with comparable direct and oxygen-mediated hopping amplitudes can host a gapless spin-orbital quantum liquid, and the ruthenium compound NaRuSi2O6 is a candidate.","keywords":["spin-orbital liquid","pyroxene","t2g orbitals","SU(6) Sutherland model","Mott insulator","dimerization","NaRuSi2O6","DMRG"],"falsifier":"Synthesize a single crystal of NaRuSi$_2$O$_6$ and measure its low-energy magnetic excitation spectrum by inelastic neutron scattering or its specific heat: observation of a gapless, dispersive continuum with weight concentrated near $q = \\pm \\pi/3$ would support the spin-orbital liquid, while a flat gapped mode near $2J$ would indicate dimerization. Alternatively, a more accurate ab initio calculation that puts $t_2/t_1$ well below 0.3 for NaRuSi$_2$O$_6$ would place it firmly in the dimerized regime and falsify the paper's main proposal.","tokens_in":11755,"feed_emoji":"⚛️","tokens_out":10753,"duration_ms":90203,"temperature":0.7,"pith_summary":"The paper argues that in pyroxene quasi-one-dimensional Mott insulators, the competition between two hopping amplitudes—direct metal–metal hopping $t_1$ and oxygen-mediated hopping $t_2$—controls whether orbital degrees of freedom stay frozen or become active. When $t_2$ is close to $t_1$, the effective low-energy model has SU(6) symmetry and is integrable, with gapless excitations that carry both spin and orbital quantum numbers. When $t_2 \\ll t_1$, orbital fluctuations are quenched and the system dimerizes into spin singlets, as in NaTiSi$_2$O$_6$. First-principles calculations give $t_2/t_1 = 0.64$ for NaRuSi$_2$O$_6$, leading the authors to propose this material as a candidate for a spin-orbital quantum liquid.","feed_headline":"Ruthenium pyroxene is a candidate spin-orbital quantum liquid","feed_subtitle":"When the two hopping amplitudes nearly match, the third t2g orbital switches on and spin-orbital excitations stay gapless.","key_machinery":"The load-bearing object is the three-orbital Hubbard model of a zigzag chain of edge-sharing MO$_6$ octahedra, with the $t_{2g}$ orbitals ($xy, yz, zx$) occupied by one electron or hole per site. The key simplification is a site-dependent relabeling of orbitals that makes the kinetic terms diagonal and splits the chain into three identical hopping channels, leaving only two hopping amplitudes $t_1$ and $t_2$, plus a single crystal-field parameter $\\Delta$. In the strong-coupling limit $U \\gg t$, this maps onto the SU(6)-symmetric Sutherland Hamiltonian, a permutation-exchange model with spin and orbital (isospin) operators, solvable by Bethe ansatz; the paper then uses DMRG time evolution to compute the spin dynamical structure factor, which tracks how spectral weight moves from the gapless continuum to localized dimers as $t_2/t_1$ and $\\Delta$ vary.","core_discovery":"The central discovery is that the third $t_{2g}$ orbital—usually neglected in minimal two-orbital descriptions—becomes dynamically relevant when the oxygen-mediated hopping integral $t_2$ is comparable to the direct hopping $t_1$. In that regime the three-orbital Hubbard model on the zigzag chain can be relabeled so that the kinetic terms decouple into three degenerate hopping paths, and in the strong-coupling limit the model maps onto the SU(6)-symmetric Sutherland Hamiltonian $H = J \\sum_k P^o_{k,k+1} P^s_{k,k+1}$ with $J = 2t^2/U$. This integrable point has gapless, fractionalized spin-orbital excitations. Any anisotropy—$t_2 \\neq t_1$ or crystal field $\\Delta \\neq 0$—generates a relevant perturbation that opens spectral gaps and drives dimerization, but the perturbation's high scaling dimension ($5/3$) means the SU(6)-like degeneracy survives for small anisotropy. The paper uses DMRG spectral functions to show that as anisotropy grows, spectral weight shifts from the gapless continuum to the local dimer-breaking scale $2J$, while a residue of soft weight near $J/2$ explains the ~10 meV heat-capacity gap of NaTiSi$_2$O$_6$. For NaRuSi$_2$O$_6$, first-principles calculations yield $t_2/t_1 = 0.64$ and no energetic preference for the dimerized structure, so the authors suggest this material as a candidate for a liquid with tightly bound spin and orbital excitations.","pith_inferences":["The same orbital-relabeling decoupling may apply to other edge-sharing octahedral chains (e.g., osmates or iridates with $d^5$ configurations), suggesting a broader family of candidate spin-orbital liquids where $t_2/t_1$ is naturally large.","The gapless SU(6) point inherited from the Sutherland model is a quantum critical point; even for moderate anisotropy the scaling dimension $5/3$ of the relevant perturbation implies that specific heat and spin correlations should show pronounced power-law precursors, which could be searched for in the proposed Ru compound before full synthesis.","A direct experimental dial: hydrostatic pressure on NaTiSi$_2$O$_6$ should increase $t_2/t_1$ by shortening Ti–O–Ti bonds, continuously reducing the dimer gap toward zero if the paper's phase diagram is correct."],"forward_implications":["If NaRuSi$_2$O$_6$ realizes the predicted liquid, its low-energy spin and orbital excitations will be gapless, dispersive, and locked together, a rare quasi-1D example of an SU(N)-type quantum disordered state.","The paper resolves an existing puzzle in NaTiSi$_2$O$_6$: the heat-capacity gap (~10 meV) is much smaller than the singlet–triplet gap ($2J \\approx 53$ meV) because residual spectral weight at ~$J/2$ survives even in the dimerized regime.","The theory provides a control parameter, $t_2/t_1$, for engineering pyroxene materials: substituting Ti by Ru or Os, or applying pressure to modify M–O–M angles, should tune a family of compounds across the dimerization-to-liquid transition.","In the dimerized phase, orbital order is expected to accompany spin singlet formation, yielding a specific pattern of occupied orbitals that can be checked by resonant x-ray scattering."],"supporting_citations":[{"why":"Supplies the integrable SU(6)-symmetric Sutherland Hamiltonian, the gapless model that anchors the symmetric limit of the three-orbital chain.","marker":"[15]"},{"why":"Establishes the $t_{2g}$ hopping paths in pyroxenes, including the oxygen-mediated shoulder-to-shoulder hoppings that define $t_2$.","marker":"[1]"},{"why":"Provides the measured spin gap of about 53 meV in NaTiSi$_2$O$_6$, the dimerized reference point for the theoretical phase diagram.","marker":"[6]"},{"why":"Reports the ~10 meV heat-capacity gap that the theory's soft spectral weight near $J/2$ is invoked to explain.","marker":"[5]"},{"why":"Shows that oxygen-mediated hopping can dominate in 4d/5d transition-metal compounds, motivating the analogous $t_2$ enhancement predicted for ruthenium pyroxene.","marker":"[10]"},{"why":"Derives the indirect hopping paths in edge-sharing chains, identifying the geometry that yields the $t_2$ channel.","marker":"[9]"},{"why":"Provides the DMRG time-evolution method used to compute the spin spectral functions that map the dimerization transition.","marker":"[16]"}],"fun_headline_variants":["Pyroxene model reveals SU(6) spin-orbital liquid state","Ruthenium pyroxene shows integrable spin-orbital liquid","Orbital liberation in pyroxenes: route to quantum liquid","Three-orbital pyroxene chain supports gapless spin-orbital excitations","SU(6) symmetry unlocks quantum liquid in pyroxene magnets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole analysis assumes the real pyroxene chain behaves as an ideal three-orbital Hubbard model in which, after relabeling orbitals, the kinetic terms separate into three identical hopping paths with only two amplitudes $t_1$, $t_2$ and one crystal field $\\Delta$; real octahedral tilts, distortions, and dimerization break this exact decoupling, and $\\Delta$ is treated as a free parameter.","fun_headline_variants_meta":{"raw":{"variants":["Pyroxene model reveals SU(6) spin-orbital liquid state","Ruthenium pyroxene shows integrable spin-orbital liquid","Orbital liberation in pyroxenes: route to quantum liquid","Three-orbital pyroxene chain supports gapless spin-orbital excitations","SU(6) symmetry unlocks quantum liquid in pyroxene magnets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00021,"raw_usage":{"total_tokens":1407,"prompt_tokens":935,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":373}},"tokens_in":551,"tokens_out":472,"duration_ms":4147,"temperature":1.0,"reasoning_tokens":373,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:41:05.637644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Synthesize a single crystal of NaRuSi$_2$O$_6$ and measure its low-energy magnetic excitation spectrum by inelastic neutron scattering or its specific heat: observation of a gapless, dispersive continuum with weight concentrated near $q = \\pm \\pi/3$ would support the spin-orbital liquid, while a flat gapped mode near $2J$ would indicate dimerization. Alternatively, a more accurate ab initio calculation that puts $t_2/t_1$ well below 0.3 for NaRuSi$_2$O$_6$ would place it firmly in the dimerized regime and falsify the paper's main proposal.","supporting_citations":[{"cited_title":"Manmana, Kaden Hazzard, Gang Chen, Adrian E","cited_arxiv_id":null,"evidence_quote":"Supplies the integrable SU(6)-symmetric Sutherland Hamiltonian, the gapless model that anchors the symmetric limit of the three-orbital chain."},{"cited_title":"(2) Att1 =t2, ∆ = 0 the solution is ϵ = 2t cosk","cited_arxiv_id":null,"evidence_quote":"Establishes the $t_{2g}$ hopping paths in pyroxenes, including the oxygen-mediated shoulder-to-shoulder hoppings that define $t_2$."},{"cited_title":"Isobe, E","cited_arxiv_id":null,"evidence_quote":"Provides the measured spin gap of about 53 meV in NaTiSi$_2$O$_6$, the dimerized reference point for the theoretical phase diagram."},{"cited_title":"Hikihara and Y","cited_arxiv_id":null,"evidence_quote":"Reports the ~10 meV heat-capacity gap that the theory's soft spectral weight near $J/2$ is invoked to explain."},{"cited_title":"Mila and F.-C","cited_arxiv_id":null,"evidence_quote":"Shows that oxygen-mediated hopping can dominate in 4d/5d transition-metal compounds, motivating the analogous $t_2$ enhancement predicted for ruthenium pyroxene."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the indirect hopping paths in edge-sharing chains, identifying the geometry that yields the $t_2$ channel."},{"cited_title":"Sutherland, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the DMRG time-evolution method used to compute the spin spectral functions that map the dimerization transition."}],"review_version":1}