{"id":"717f5042-b9c1-4913-a27d-c4f71262e3b7","arxiv_id":"2501.10287","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In the cuprate ladder Sr14Cu24O41, the magnetic response of doped holes is far weaker than the Hubbard model predicts, indicating a large nearest-neighbor attraction that enhances d-wave-like hole pairing.","lead":"Researchers measured magnetic excitations in a copper-oxide ladder using resonant X-rays and compared the spectra to quantum many-body simulations. The data disagree with the plain Hubbard model and point to a strong nearest-neighbor attraction that pairs holes, a potential missing ingredient for high-temperature superconductivity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central V inference rests on the single-ion Rspin normalization of SM Section 2; if the QP spin-flip channel has a different RIXS matrix element than assumed, the observed suppression is a normalization artifact and V need not be attractive.","rationale":"The reader's weakest_assumption listed both the RIXS-to-S(q,omega) normalization/matrix-element uncertainty and the zero-temperature versus 260 K comparison. I focus on the matrix-element issue because it directly targets the interpretation of the key observable: the weak QP spin-flip branch. The paper's main conclusion that an attractive V of order -1.0t to -1.25t is needed follows almost entirely from the intensity ratio between the QP spin-flip branch and the two-triplon continuum. If that ratio is distorted by an incorrect RIXS cross-section for the QP channel, the fitted V loses its empirical basis. This is not an accusation of error; it is a request for a quantitative check that the single-ion normalization is adequate. The paper's qualitative observation that a low-energy branch is present above TCO and absent below is robust, and the DMRG methodology is standard and reproducible. However, the decisive quantitative step, converting RIXS intensity into absolute S(q,omega), relies on a single-ion approximation whose validity for the QP spin-flip branch is asserted rather than demonstrated. The proposed multi-site RIXS calculation or a polarization-dependent measurement would settle this. If the test shows the single-ion normalization is adequate, the paper's conclusion is substantially strengthened; if not, the central claim fails regardless of whether the temperature comparison is handled. Thus the reader's conditional verdict is appropriate, and my concern does not change the overall disposition but sharpens the condition that must be met.","tokens_in":16701,"tokens_out":10616,"duration_ms":107329,"concrete_test":"Compute the RIXS cross-section for the doped two-leg ladder using a multi-site cluster (e.g., 8x2 or 16x2) with the fitted parameters (t = 0.38 eV, t_perp = 0.84t, t' = -0.3t, U = 8t, V = 0) and the Kramers-Heisenberg formula, and compare the resulting QP spin-flip and two-triplon RIXS intensities with the corresponding bare S(q,omega) features. If the RIXS QP-to-two-triplon ratio is suppressed by a factor comparable to the observed suppression (roughly an order of magnitude) even at V = 0, the single-ion Rspin normalization is the culprit and the inferred V is an artifact. Alternatively, repeat the RIXS measurement at the same momentum with sigma and pi incident polarizations; if the normalized S(q,omega) for the QP branch is polarization-dependent, the single-ion spin-flip cross-section is insufficient.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that a large attractive nearest-neighbor interaction V (about -1.0t to -1.25t) is required to explain the faint quasiparticle spin-flip branch at 260 K. The experimental S(q,omega) is obtained by dividing the RIXS intensity by a single-ion spin-flip cross-section Rspin(epsilon, epsilon', Omega_i) computed for an isolated Cu2+ ion (SM Section 2, Ispin proportional to Rspin times S(q,omega)). This procedure assumes that every magnetic feature, including the QP spin-flip branch, has the same RIXS matrix element as a local single-ion spin flip. In a doped ladder, however, the QP spin-flip excitation involves a hole quasiparticle moving in a singlet background; the local electronic configuration at the flippable site is modified by the presence of the hole and by strong correlations. The single-ion Rspin cannot capture this and may differ substantially for the QP branch relative to the two-triplon continuum. If the intrinsic RIXS cross-section for the QP spin-flip is reduced, or has a different energy or angle dependence, dividing by Rspin would artificially suppress the extracted S(q,omega) for that branch, mimicking the effect of an attractive V. The authors cite reference [12] for the validity of the correspondence in doped systems, but the QP spin-flip branch is precisely the feature whose matrix element is least protected by the single-ion mapping. This is the most load-bearing assumption because the entire inference of V and the enhanced pairing energy depends on comparing the intensity of this branch with the two-triplon continuum. A secondary but related concern is the zero-temperature DMRG compared to 260 K data; however, the matrix-element issue alone would invalidate the central claim even if the temperature comparison were exact.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports Cu L-edge RIXS measurements of the spin-ladder compound Sr14Cu24O41 and compares the extracted magnetic dynamical structure factor with DMRG calculations for a two-leg Hubbard ladder and an extended Hubbard ladder with nearest-neighbor interaction V. The central experimental observation is a weak, dispersive quasiparticle spin-flip branch above the charge-order temperature, whose intensity is strongly suppressed relative to the single-band Hubbard model. The authors show that introducing an attractive nearest-neighbor Coulomb interaction V in the range -1.0t to -1.25t, with U-V held fixed to preserve the exchange couplings, reproduces the suppression of this branch, and they further compute that this V enhances the hole-pair binding energy by almost an order of magnitude while retaining a d-wave-like pairing symmetry. The paper also presents orbiton dispersion data suggesting a crossover between one- and two-dimensional physics. The authors conclude that the Hubbard model is insufficient for cuprate ladders and that an attractive intersite interaction may be a universal ingredient for d-wave pairing in cuprates.","tokens_in":17017,"tokens_out":4258,"duration_ms":47057,"significance":"If the central inference is correct, the paper provides a rare experimental signature of beyond-Hubbard physics in a cuprate material and quantifies a large attractive nearest-neighbor interaction that could resolve the ongoing debate about whether the pure Hubbard model supports robust superconductivity. The experimental data are of high quality: the RIXS measurements have 35 meV resolution, the two-triplon dispersion is well characterized, the doping level is verified by XAS, and the charge-order transition is independently determined. The DMRG calculations are state-of-the-art, with clear documentation of the model Hamiltonian, parameter extraction via Bayesian optimization, and the dynamical spin structure factor computed on 64×2 clusters. The paper also explicitly discusses boundary effects and pair-pair correlations, which is commendable. However, as detailed below, the quantitative claims about V and the order-of-magnitude binding-energy enhancement rest on a single-model comparison and a normalization procedure whose validity for the quasiparticle spin-flip branch is not established.","major_comments":[{"comment":"The conversion of raw RIXS intensity into S(q,ω) uses a single-ion spin-flip cross-section Rspin(ε,ε′,Ω_i) to normalize all magnetic features, including the quasiparticle spin-flip branch. The central inference of an attractive V depends on this branch being intrinsically weak after normalization. However, the RIXS matrix element for a quasiparticle spin flip in a doped ladder may differ substantially from that of an isolated Cu2+ ion: the local electronic configuration around the hole is strongly modified by doping, and the cited reference [12] primarily validates the single-ion mapping for the two-triplon continuum, not for the quasiparticle spin-flip channel. If the intrinsic cross-section for this branch is suppressed by the same many-body effects that the authors attribute to V, then dividing by the single-ion Rspin would artificially lower the extracted S(q,ω), mimicking the effect of V<0. This is a load-bearing assumption. The authors should either compute the RIXS cross-section for the doped ladder within their DMRG framework (e.g., along the lines of Ref. [30]) to verify that the matrix element for the quasiparticle spin-flip branch is unchanged from the single-ion value, or provide an independent calibration that rules out an energy- or momentum-dependent matrix-element artifact.","section":"SM Section 2 (Extracting S(q,ω)) and Fig. 2"},{"comment":"The quantitative conclusion that V is -1.0t to -1.25t is obtained by tuning V to reproduce the suppression of the quasiparticle spin-flip intensity in the same data set. The subsequent calculation that this V enhances the hole-pair binding energy by 'almost an order of magnitude' is then a consequence of the fitted V, not an independent prediction or test of enhanced pairing. This is a circularity in the presentation: the binding-energy enhancement is an output of the model whose key parameter was already fit to the magnetic response. The authors should state this explicitly, provide a goodness-of-fit metric or uncertainty estimate for the inferred V range, and avoid presenting the binding-energy enhancement as independent corroboration. If possible, they should identify a distinct, a priori prediction (for example, the momentum dependence of the suppression or the behavior under additional doping) that could be tested against new data.","section":"Fig. 3 and Appendix (DMRG calculations, Binding energy)"},{"comment":"The key experimental spectra were measured at 260 K, above TCO, while the DMRG calculations are zero-temperature. Thermal broadening, thermally excited quasiparticles, and fluctuation effects can independently reduce the intensity of the quasiparticle spin-flip branch relative to the T=0 Hubbard-model calculation, even for V=0. The paper does not assess finite-temperature effects, and the 40 K data are obscured by the charge-order gap, so there is no direct low-temperature check of the Hubbard-model prediction. Without a quantitative estimate of thermal effects on the quasiparticle spin-flip intensity, the comparison is not fully controlled. The authors should either estimate the temperature dependence (e.g., using finite-temperature DMRG or a simplified phenomenological model) or restrict the quantitative claim to a statement of qualitative suppression, deferring the V inference to a lower-temperature measurement where charge order does not interfere.","section":"Comparison of zero-temperature DMRG to 260 K data (Figs. 2 and 3)"}],"minor_comments":[{"comment":"The definition of binding energy B.E. = 2E_GS(N-1) - E_GS(N) - E_GS(N-2) is standard, but for a 64×2 cluster with open boundary conditions the finite-size corrections to this quantity can be as large as the binding energy itself. Please provide an estimate of finite-size errors (e.g., by repeating the calculation for different cluster lengths) or a discussion of why they are negligible.","section":"Appendix Eq. (5)"},{"comment":"The choice to hold U-V fixed at 8t when varying V is an important modeling assumption. The text states that this keeps J constant, but the strong-coupling expression J = -4t^2/(U-V) is only exact in the strong-coupling limit and for a specific charge-transfer process. Please justify this choice more carefully and discuss whether alternative ways of introducing V (e.g., fixing U and varying V) would change the conclusions.","section":"Appendix, model parameters"},{"comment":"In the description of the normalization procedure, the paper states that the orbital excitation integrated intensity is matched to the single-ion ED calculation. The figure shows the θ-dependent scaling factors, but the text does not describe how the normalization is propagated to the magnetic excitations or what uncertainty is introduced by the elastic and phonon subtraction. A short error-propagation analysis, or at least a statement of the typical scale of the normalization correction, would make the quantitative comparison more transparent.","section":"SM Section 2 and Fig. S7"},{"comment":"The orbiton dispersion and dimensional-crossover discussion in Fig. 4 is interesting but somewhat disconnected from the central pairing claim. The connection to the suppression of quasiparticle spin flips is not made explicit, and the statement that the results 'may be directly relevant' to two-dimensional cuprates is speculative. Consider either providing a more direct link (e.g., through the effect of dimensionality on V) or clearly labeling this as a separate observation.","section":"Fig. 4 and main text, final paragraph"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and important question, and the experimental data are of high quality. The qualitative observation of a suppressed quasiparticle spin-flip branch over the Hubbard-model expectation is likely to be of interest to the community. However, the two main quantitative conclusions—the range of V and the order-of-magnitude binding-energy enhancement—are currently supported by a single modeling comparison that suffers from a circularity (V is fit to the data, then used to compute the binding energy) and from a normalization assumption (single-ion Rspin) that is least protected for the very feature used to infer V. These issues are fixable in principle: the authors could compute the RIXS cross-section for the doped ladder, assess thermal effects, and present the V inference with clearer statistics. If they choose to do so, the paper could become a strong contribution. If not, the central claims as stated are not yet established. I recommend major revision rather than rejection, because the core experimental finding is credible and the theoretical framework is appropriate, but the quantitative claims need substantial additional support or significant tempering."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is experimental: above the charge-order transition, the quasiparticle spin-flip branch in Sr14Cu24O41 is much weaker than a doped single-band Hubbard ladder predicts. That qualitative observation looks solid. The RIXS data are high resolution, the background subtraction is careful, and the branch is absent below T_CO and faint above, which is a nice internal control. The DMRG calculations are also well done: parameters are extracted from the undoped two-triplon dispersion via Bayesian optimization, and the comparison is made at the same doping. The orbiton dispersion and the dimensional-crossover argument in Fig. 4 are a useful bonus, not just padding. Credit where it is due: this is a well-executed measurement that identifies a real discrepancy with the Hubbard model.\n\nThe soft spots are in the quantitative chain, not in the observation. The conversion of RIXS intensity into S(q,omega) divides by a single-ion spin-flip cross-section Rspin. The authors acknowledge this ignores doping and cite ref [12] for validity in doped systems, but that reference does not specifically test the quasiparticle spin-flip branch—the one feature whose matrix element is least protected by the single-ion mapping. If that branch has a different RIXS matrix element or a different energy/angle dependence, the apparent suppression could be a normalization artifact. The stress-test note makes this point well, and I think it is the main load-bearing assumption. Second, the DMRG spectra are zero-temperature while the key data are at 260 K; thermal broadening or precursor charge-order fluctuations could weaken the Hubbard-model branch on their own. That is a secondary but real concern. Third, V is tuned to reproduce the suppression and then used to compute the binding-energy enhancement, so the factor-of-ten binding energy is not an independent prediction. There are no error bars on the fitted V or the weak peak amplitudes, and the extrapolation to 2D cuprates is speculative.\n\nNone of these are fatal. The qualitative discrepancy is there, and the attractive-V explanation is plausible and consistent with earlier theory. But the paper overstates what the evidence supports. The abstract's \"nearly an order of magnitude\" binding-energy enhancement should be softened, and the matrix-element and finite-temperature questions should be addressed or explicitly left as caveats. I would send this to serious peer review: the experimental observation is important enough that a good referee can help sharpen the interpretation. I'd also bring it to our reading group—it will generate a useful discussion about what RIXS intensities do and do not tell you.\n\nRecommendation: engage with it. Require revision, but do not desk reject.","headline":"A clear new RIXS observation of a suppressed quasiparticle spin-flip branch in a cuprate ladder, with a plausible but quantitatively conditional attribution to an attractive nearest-neighbor V.","tokens_in":17735,"tokens_out":1517,"would_cite":true,"duration_ms":17118,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In the cuprate ladder Sr14Cu24O41, the magnetic signal from doped holes is far weaker than the single-band Hubbard model predicts, and the paper traces the suppression to a strong nearest-neighbor attraction that enhances d-wave-like hole…","keywords":["cuprate ladder","Hubbard model","nearest-neighbor attraction","resonant inelastic x-ray scattering","dynamical spin structure factor","d-wave pairing","density matrix renormalization group","Sr14Cu24O41"],"falsifier":"Measure the quasiparticle spin-flip branch in a doped two-leg ladder where the nearest-neighbor interaction is known to be small, or compute the 260 K spectra with finite-temperature methods that include charge-order fluctuations; if the branch then appears with full Hubbard-model intensity, the inferred attractive $V$ is not required to explain the suppression.","tokens_in":16430,"feed_emoji":"🧲","tokens_out":6829,"duration_ms":60771,"temperature":0.7,"pith_summary":"The paper uses high-resolution resonant inelastic x-ray scattering on the self-doped cuprate ladder Sr14Cu24O41 to isolate a low-energy branch of magnetic excitations coming from the doped holes. That branch is nearly absent at low temperature and only faintly present above the charge-order transition, in sharp contrast to the intense quasiparticle spin-flip branch predicted by density-matrix renormalization group calculations for a single-band Hubbard ladder. The authors argue the discrepancy is not disorder or charge order but an additional nearest-neighbor attractive interaction $V$ of order $-1.0t$ to $-1.25t$, which binds holes on the same rung and suppresses hole spin-flip scattering. In their model this attraction raises the hole-pair binding energy by almost an order of magnitude and preserves d-wave-like pairing correlations. If correct, the result identifies a missing ingredient beyond the Hubbard model that may be needed for robust d-wave superconductivity in the cuprates.","feed_headline":"Faint spin branch in cuprate ladder reveals strong hole pairing","feed_subtitle":"A nearest-neighbor attraction of about -1t to -1.25t boosts d-wave-like pairing by nearly an order of magnitude.","key_machinery":"The central object is the extended Hubbard model on a two-leg ladder, with leg hopping $t$, rung hopping $t_\\perp=0.84t$, diagonal hopping $t'=-0.3t$, on-site repulsion $U=8t$, and a nearest-neighbor interaction $V$; the paper compares its dynamical spin structure factor $S(q,\\omega)$, computed with the Krylov-space correction-vector DMRG method, against RIXS spectra normalized to a single-ion spin-flip cross section. The key mechanism is the response of $V$: a repulsive $V$ barely changes the spectra, while an attractive $V$ makes it energetically favorable for doped holes to occupy the same rung, restoring intact rung singlets that can be excited into the two-triplon continuum and suppressing the quasiparticle spin-flip branch. The pairing strength is quantified by the hole-pair binding energy $2E_{GS}(N-1)-E_{GS}(N)-E_{GS}(N-2)$, and the pairing symmetry by the relative sign of rung-rung and rung-leg singlet pair correlations.","core_discovery":"The central claim is that the magnetic response of the doped holes in Sr14Cu24O41, observed as a quasiparticle spin-flip branch in the dynamical spin structure factor $S(q,\\omega)$, is strongly suppressed relative to the single-band Hubbard model, and that the suppression is a direct spectroscopic signature of enhanced hole pairing caused by a large nearest-neighbor attractive interaction. Using model parameters fixed by the undoped two-triplon dispersion, the paper shows the Hubbard model produces an intense spin-flip branch that is absent in the measured spectra. Introducing $V \\approx -1.0t$ to $-1.25t$ in the extended Hubbard ladder reproduces the data: holes bind into rung pairs, the pair binding energy grows to about $0.074t$, and the quasiparticle spin-flip intensity decreases monotonically while the two-triplon continuum sharpens. The pair correlations retain a d-wave-like structure, with rung and leg pair correlations of opposite sign, and the measured orbital dynamics place the ladder in a crossover regime between one and two dimensions, supporting the extension of the conclusion to two-dimensional cuprates.","pith_inferences":["If the claim holds, high-resolution RIXS on lightly doped two-dimensional cuprates should reveal a similarly suppressed quasiparticle spin-flip weight relative to Hubbard-model predictions, providing a direct test in the 2D case.","The inferred interaction is consistent with an electron-phonon origin; a sharper test would be an isotope-exchange experiment that shifts phonon energies and checks whether the quasiparticle spin-flip intensity recovers.","The paper's ladder geometry has no local $C_4$ symmetry, yet the pairing retains d-wave-like sign structure; an interesting extension is whether the same robustness persists in models with explicit symmetry-breaking disorder or anisotropic $V$.","Because the suppression sets in monotonically with $|V|$, a quantitative map of the quasiparticle spin-flip intensity across different ladder compounds could serve as a spectroscopic thermometer for the effective intersite attraction."],"forward_implications":["The single-band Hubbard model, with only on-site repulsion, does not fully account for the magnetic excitations of doped cuprate ladders; an attractive nearest-neighbor interaction is required to match the measured spectra.","An attractive $V$ of about $-1.0t$ to $-1.25t$ increases the hole-pair binding energy by nearly an order of magnitude and leaves the pairing symmetry d-wave-like.","Because the suppression is observed above the charge-order transition, the enhanced pairing is intrinsic to the doped ladder and not a byproduct of the charge-order gap.","The combination of a gapped triplon continuum and a confined orbiton places this ladder in a dimensional crossover regime, supporting the relevance of the result to two-dimensional cuprates.","Similar values of $V$ inferred in one-dimensional chain cuprates suggest the attractive interaction is shared across cuprate families and may be a common ingredient for robust superconductivity."],"supporting_citations":[{"why":"Predicts the quasiparticle spin-flip branch and weak pairing of lightly doped two-leg ladders that the experiment tests.","marker":"[28]"},{"why":"Provides the DMRG spin and charge dynamics of two-leg ladders showing the intense quasiparticle spin-flip branch that the measured spectra contradict.","marker":"[30]"},{"why":"Supplies the single-ion model used to normalize RIXS intensity into an absolute dynamical spin structure factor.","marker":"[33]"},{"why":"Gives the disordered-Hubbard-ladder result used to rule out disorder as the cause of the suppression.","marker":"[36]"},{"why":"Fixes the ladder exchange parameters through the one- and two-triplon spectra.","marker":"[37]"},{"why":"Provides the Krylov-space correction-vector DMRG method used for all theoretical spectra.","marker":"[38]"},{"why":"Reports a similarly strong nearest-neighbor attraction in one-dimensional cuprate chains, supporting the generic value of $V$.","marker":"[40]"},{"why":"Shows d-wave-like pair correlations in a two-leg ladder, the symmetry the paper finds persists with $V$.","marker":"[46]"},{"why":"Shows nearest-neighbor attraction enhances superconductivity in the extended Hubbard model, the broader implication of the measured pairing.","marker":"[16]"}],"fun_headline_variants":["Hubbard model fails to capture cuprate ladder pairing","Nearest-neighbor attraction boosts cuprate hole pairing","Suppressed spin branch reveals strong pairing in cuprate ladder","Cuprate ladder spin response points beyond Hubbard model","Strong hole pairing arises from nearest-neighbor attraction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inference relies on the conversion of raw RIXS intensity into an absolute $S(q,\\omega)$ using a single-ion spin-flip cross-section, and on the assumption that zero-temperature DMRG on an isolated 64x2 ladder represents the 260 K thermal state without additional broadening or charge-order precursors.","fun_headline_variants_meta":{"raw":{"variants":["Hubbard model fails to capture cuprate ladder pairing","Nearest-neighbor attraction boosts cuprate hole pairing","Suppressed spin branch reveals strong pairing in cuprate ladder","Cuprate ladder spin response points beyond Hubbard model","Strong hole pairing arises from nearest-neighbor attraction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000517,"raw_usage":{"total_tokens":2495,"prompt_tokens":922,"completion_tokens":1573,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":1493}},"tokens_in":538,"tokens_out":1573,"duration_ms":11273,"temperature":1.0,"reasoning_tokens":1493,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:15:33.903942+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the quasiparticle spin-flip branch in a doped two-leg ladder where the nearest-neighbor interaction is known to be small, or compute the 260 K spectra with finite-temperature methods that include charge-order fluctuations; if the branch then appears with full Hubbard-model intensity, the inferred attractive $V$ is not required to explain the suppression.","supporting_citations":[{"cited_title":"Propagation of an orbiton in the antiferromagnets: theory and experimental verification","cited_arxiv_id":"1912.11363","evidence_quote":"Shows nearest-neighbor attraction enhances superconductivity in the extended Hubbard model, the broader implication of the measured pairing."}],"review_version":1}