{"id":"6d0855e0-5cb1-41ba-8723-741e0e409ed0","arxiv_id":"2501.13322","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"An N-layer Ising model in the infinite interlayer coupling limit reduces to a single layer with N times the coupling, which the authors argue stabilizes the charge-density wave and explains d-wave-like experimental signals in cuprates.","lead":"A simple Ising-model calculation claims that charge-density waves in multilayer cuprates are strongly stabilized by layer coupling, and that the resulting pattern may explain why some experiments mistook this charge order for d-wave superconductivity. The mathematical result is shown only in an extreme limit and is not tied to measured material parameters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stabilization is proven only at J'/J → ∞; real cuprate parameters give J'/J of order unity, so the effective coupling NJ and enhanced CCDW temperature are not established for Bi2212/Bi2223.","rationale":"The reader's weakest assumption correctly identifies the infinite-coupling limit as the load-bearing premise. The paper's stabilization conclusion follows from the exact dimensional reduction at J'/J → ∞, but real Bi2212 and Bi2223 have finite interlayer/intralayer oxygen distances, so J'/J is of order unity. The paper explicitly acknowledges in the Supplementary Material that it can only take the K' → ∞ limit, equivalent to J'/J → ∞, and provides no finite-ratio analysis. The ferromagnetic sign issue flagged by the reader is also real and would misidentify the ground state as a uniform alignment rather than the alternating ±δe CCDW pattern of Fig. 13; however, on a bipartite lattice a sublattice gauge transformation maps the ferromagnetic Ising model to the antiferromagnetic one, so the partition-function and dimensional-reduction results would survive such a reinterpretation. The finite-J'/J problem is not repairable by a transformation and therefore is the more decisive weakness. Since this is the same concern the reader raised, the verdict remains REJECT; no adjustment is needed.","tokens_in":19684,"tokens_out":4568,"duration_ms":48100,"concrete_test":"Perform a transfer-matrix or Monte Carlo study of the N = 2 and N = 3 layer Ising model of Eqs. (4)-(6) for the staggered order parameter appropriate to the CCDW, at finite ratios J'/J = 1.2, 1.5, 2, and 5 on lattices up to at least 64×64 sites per layer. For each ratio, extract the ordering temperature T_c(N, J'/J) and compare it with N times the single-layer Onsager critical temperature. If T_c(N, J'/J) is substantially below N·T_c(1) for J'/J ≈ d/d' (a ratio plausibly between 1 and 2 in Bi2212/Bi2223), then the stabilization mechanism does not apply to real cuprates and the paper's central physical conclusion fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that multilayer structure multiplies the intralayer coupling by N and thereby strongly stabilizes the CCDW—is derived exclusively in the K' → ±∞ limit, i.e. J'/J → ∞. The paper's own parameterization in Eq. (7) gives J'/J = d/d', where d and d' are near-neighbor intralayer and interlayer oxygen distances. Even if d' < d as the authors assert for Bi2212/Bi2223, this ratio is a finite number of order unity, not infinite. No finite-J' calculation, estimate, or perturbative argument is provided anywhere in the main text or Supplementary Material. At finite J'/J, spins in adjacent layers are not rigidly locked; interlayer fluctuations reduce the effective in-plane coupling below NJ, and the partition function does not reduce to a single layer. Consequently, the quantitative basis for the claimed N-fold enhancement of the CCDW ordering temperature—and hence for the explanations of enhanced THz emission and the apparent d-wave symmetry of many experiments—is absent for the actual materials discussed. The mathematical limit result itself may be correct, but it is a singular limit whose physical relevance has not been established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes an N-layer Ising model on the oxygen sublattice of the CuO2 planes as a description of the commensurate charge-density wave (CCDW) in underdoped Bi2201 (N=1), Bi2212 (N=2), and Bi2223 (N=3). The Hamiltonian in Eqs. (4)-(6) has near-neighbor intralayer and interlayer couplings J and J' with J,J'>0. The central result, established by finite-size partition-function evaluations in the Supplementary Material, is that in the K'→±∞ limit the ratio Z_LMN(K,K')/[Z_112(K')]^{LM} reduces to Z_LM1(NK)/2^{LM}, so the N-layer system behaves as a single-layer Ising model with coupling NJ. From this the authors conclude that multilayer stacking strongly stabilizes the CCDW and offer this as an explanation of enhanced THz emission from Bi2212 mesas and of experiments that have been interpreted as evidence for d_{x^2-y^2}-wave superconductivity.","tokens_in":19792,"tokens_out":19487,"duration_ms":184456,"significance":"The partition-function reduction is a clean and internally consistent result: the finite-size checks for strips and small two-dimensional sections are reproducible with the stated methods, no parameters are fitted to data, and the reduction is a parameter-free statement about the Hamiltonian as written. If it applied at realistic couplings, the argument would give a simple mechanism for multilayer enhancement of charge order. However, the model's physical identification is currently incorrect (ferromagnetic ground state versus the claimed alternating charge pattern), and the enhancement is proven only in the singular J'/J→∞ limit, whereas the paper's own Coulomb parameterization gives J'/J of order unity for Bi2212 and Bi2223. The interpretive claims about THz emission and the superconducting order-parameter symmetry go beyond what is derived.","major_comments":[{"comment":"The Hamiltonian in Eqs. (4)-(6) with J, J' > 0 has negative signs, so it is the ferromagnetic Ising model: a bond with σ_i = σ_j contributes −J, while a bond with σ_i = −σ_j contributes +J. Its ground state is therefore uniform (all σ equal), whereas the CCDW described in the text and in Figs. 13-16 is an alternating checkerboard of ±δe charges with zero net charge. The passage after Eq. (7) stating that the negative signs incorporate the Coulomb preference for opposite charge deviations is backwards for this Hamiltonian; the intended physics requires the opposite sign convention (H ∝ +J Σ σσ). On the bipartite square lattice the partition functions of the two conventions are equal under a staggered gauge transformation, so the reduction result may survive, but the manuscript as written misidentifies the model's ground state and its relationship to the physical CCDW, and this must be corrected explicitly.","section":"§V, Eqs. (4)-(7) and Figs. 13-16"},{"comment":"The reduction to an effective single layer with coupling NJ is derived and verified only in the K'→±∞ limit at fixed K, equivalently J'/J→∞, a limitation the Supplementary Material itself states after Eq. (22) ('we can only take the limit K'→∞, which is equivalent to the J'/J→∞ limit'). The paper's own parameterization in Eq. (7) gives J'/J = d/d', which is of order unity for Bi2212 and Bi2223 even if d' < d as asserted. No finite-J' computation, estimate, or perturbation argument is supplied, so the claims in the abstract that the CCDW is 'strongly enhanced and stabilized' in the multilayer compounds and that the coupling becomes NJ are unsupported extrapolations from a singular limit. To make the central claim load-bearing for real materials, the authors need to provide finite-J' results (for example transfer-matrix, Monte Carlo, or mean-field calculations of the bilayer and trilayer Ising model at physical J'/J, including the resulting ordering temperature), or restrict the conclusions explicitly to the singular limit.","section":"§V and Supplementary §I"},{"comment":"The abstract and Section IV attribute the enhanced THz emission from Bi2212 mesas and the d_{x^2-y^2}-wave interpretation of many phase-sensitive experiments to the CCDW stabilization described by the model. These statements are not derived in the manuscript: there is no calculation connecting the Ising partition function to Josephson emission linewidths or powers, and no calculation showing that a static checkerboard charge pattern would produce the interference signatures (tricrystal rings, SQUIDs, twist junctions) attributed to d-wave pairing. These are interpretive assertions that should either be substantiated by explicit calculations or removed from the abstract and conclusions, or clearly labeled as speculation.","section":"Abstract and §IV"}],"minor_comments":[{"comment":"The manuscript contains several typos and misspellings, including 'It has has long been known' (Section I), 'In ths model' (Section V), and 'easly found' (Supplementary Eq. (1)); these should be corrected in a revision.","section":"Throughout"},{"comment":"The footnote to Ref. [40] recounts a private conversation at a conference and is not appropriate material for a reference footnote; it should be removed or replaced with a neutral statement about the published work.","section":"Ref. [40] footnote"},{"comment":"The closing paragraph of the Supplementary Material contains two unresolved citation placeholders ('[?]') for the statements about the YBCO gap anisotropy, and these need to be supplied.","section":"Supplementary closing paragraph"},{"comment":"The abstract calls J and J' 'repulsive interactions,' but with the negative signs in Eqs. (4)-(6) they are ferromagnetic couplings in the spin language; the terminology should be reconciled with the corrected sign convention.","section":"Abstract"},{"comment":"The claim that the checkerboard charge pattern has d_{x^2-y^2} symmetry is asserted without a derivation; the precise relation between the real-space pattern in Fig. 13 and the polar form in Fig. 14 should be defined explicitly, since it supports the discussion in Section IV.","section":"Figs. 13 and 14"}],"recommendation":"major_revision","confidential_remarks":"The core mathematical observation is modest but appears correct and checkable; the publication decision should hinge on whether the finite-J' gap can be filled and whether the interpretive claims are scaled back. The manuscript is heavily self-cited (roughly fifteen of fifty-four references are by the corresponding author), which is worth editorial awareness but not, by itself, disqualifying. The paper's position on the s-wave versus d-wave controversy is outside the current consensus; I have based my comments only on internal consistency and support, not on that disagreement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does one clean thing: it shows by explicit finite-lattice partition function evaluations that for L×M×N Ising sections with interlayer coupling J'→∞, the partition function reduces to that of a single layer with J replaced by NJ. The Mathematica results look internally consistent, and the reduction is what one expects from the standard dimensional-reduction property of layered Ising models. That part is solid, as far as it goes. The supplement gives reproducible formulas, which is real evidence.\n\nThe problems start when the paper moves from that limit to the cuprates. The stabilization argument requires J'/J→∞, but the paper's own Eq. (7) gives J'/J = d/d', a ratio of order unity for Bi2212 and Bi2223. No finite-J' calculation, estimate, or perturbative argument is provided. The stress-test note has this right: the N-fold enhancement is not established for the materials discussed. The abstract's language ('therefore strongly enhanced and stabilized') overstates what the math supports.\n\nSecond, the Hamiltonian as written is ferromagnetic (H = -JΣσσ with J,J'>0). Its ground state is uniform spin alignment, not the alternating ±δe pattern of Fig. 13. The authors say the minus sign accounts for the opposite charge deviations, but a repulsive Coulomb interaction between like deviations would require +JΣσσ. On a bipartite lattice a sublattice gauge transformation maps one to the other, but the paper never makes that mapping explicit, so the claimed ground-state interpretation is muddled. This is fixable—state that the physical charges are staggered spins—but as written it is a genuine inconsistency.\n\nThird, the connections to THz emission and to the d-wave experiments are asserted, not derived. The logic is conditional ('likely accounting for'), and those sections read more like a research program than a result.\n\nTo the paper's credit, it is transparent about what it proves: only the J'→∞ limit for finite sections, with no general proof for arbitrary L,M,N. The finite-lattice formulas are reproducible, and the paper does not hide the limitation.\n\nWho is this for? Someone asking whether multilayer coupling can raise CDW ordering temperatures. The paper gives a clear Yes in the singular limit, and an Unproven for real materials. The question is legitimate and the math, though narrow, appears correct. I would send it to a serious referee, but the referee should demand (i) a fix or explicit gauge treatment of the Hamiltonian sign, (ii) finite-J' results or a scaling argument for why order-unity J'/J is enough, and (iii) a substantial toning down of the claims about THz and d-wave. Without those, reject.","headline":"Careful finite-lattice calculation of the J'→∞ limit, but the physical claims for cuprates rest on a singular limit with no finite-J' support, plus an unaddressed sign issue in the Hamiltonian.","tokens_in":20408,"tokens_out":2518,"would_cite":false,"duration_ms":41078,"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":"In the strong interlayer-coupling limit, an N-layer stack of CuO2 planes behaves as a single layer with coupling NJ, stabilizing commensurate charge-density waves in bilayer and trilayer cuprates.","keywords":["commensurate charge-density wave","pseudogap","cuprate superconductors","multilayer Ising model","Bi2212","Bi2223","terahertz emission","order parameter symmetry"],"falsifier":"Compare the onset temperature of the commensurate charge order in Bi2201, Bi2212, and Bi2223 at matched hole doping: the model predicts a systematic rise with layer number, governed by the effective coupling $NJ$ in the strong-coupling limit; if the onset does not rise with $N$, or rises far less than the $N$-fold enhancement, the central claim is falsified. A second check is to measure or compute $J'/J$ from the oxygen-oxygen distances and screening: if $J'/J$ is not large compared with 1, the strong-coupling identity is not the governing regime for these materials.","tokens_in":19354,"feed_emoji":"","tokens_out":17368,"duration_ms":158719,"temperature":0.7,"pith_summary":"The paper proposes that the pseudogap observed above the superconducting transition in underdoped bismuth-based cuprates is a commensurate charge-density wave (CCDW): within each CuO2 plane, alternating oxygen sites carry excess charges $\\pm\\delta e$, forming a pattern below a temperature $T_p$ but above $T_c$. It models this order with an $N$-layer Ising model on a square lattice, with repulsive intralayer coupling $J$ and interlayer coupling $J'$ between neighboring oxygen sites. The central result is that, for finite $L\\times M\\times N$ sections, the partition function in the $J'\\to\\pm\\infty$ limit reduces exactly to that of a single $L\\times M$ layer with $J$ replaced by $NJ$. Two- and three-layer compounds therefore behave like monolayers with doubled or tripled coupling, which stabilizes the CCDW and raises its ordering temperature. If this picture is correct, the nodal 'd-wave-looking' features in pseudogap measurements are charge-order features rather than superconducting pairing signatures, and the strong terahertz emission from underdoped Bi2212 mesas follows from the rigidity of the commensurate pattern.","feed_headline":"Stacking CuO2 layers multiplies charge-order coupling by N","feed_subtitle":"Stacking CuO2 planes turns a monolayer's coupling J into NJ, stabilizing the charge order behind the pseudogap.","key_machinery":"The central object is the $N$-layer Ising partition function $Z_{LMN}(K,K')$ on an $L\\times M\\times N$ lattice of oxygen-site charge deviations $\\sigma_{i,j,k}=\\pm1$, with $K=J/k_BT$ and $K'=J'/k_BT$. The mechanism is a strong-coupling reduction: for $K'\\to\\pm\\infty$, the interlayer bonds force $\\sigma_{i,j,1}=\\cdots=\\sigma_{i,j,N}$ in every column, so the intralayer couplings of the $N$ layers add, and the partition function ratio $Z_{LMN}(K,K')/[Z_{112}(K')]^{LM}$ tends to $Z_{LM1}(NK)/2^{LM}$. This 'effective single layer with coupling $NJ$' identity carries the stabilization argument: it converts a multilayer charge-order problem into the single-layer 2D Ising problem at $N$ times the monolayer coupling.","core_discovery":"The paper's central claim is that multilayer stacking multiplies the effective intralayer coupling by the number of layers. For an Ising model on an $L\\times M\\times N$ lattice with couplings $K=J/k_BT$ and $K'=J'/k_BT$, the $J'\\to\\pm\\infty$ limit locks each column of $N$ spins into a single configuration, so the intralayer bonds of all $N$ layers contribute coherently. Explicit finite-lattice evaluations for $N=2$ and $N=3$ give $Z_{LMN}(K,K')/[Z_{112}(K')]^{LM}\\to Z_{LM1}(NK)/2^{LM}$, i.e. an effective monolayer with coupling $NJ$. The paper states this reduction as the mechanism by which the CCDW is enhanced and stabilized in Bi2212 and Bi2223 relative to monolayer Bi2201, and uses it to argue that the pseudogap's $d_{x^2-y^2}$-shaped nodal pattern is a charge-order property. It further claims this accounts for the enhanced THz emission from underdoped Bi2212 mesas and for experiments that were interpreted as evidence for a $d_{x^2-y^2}$ superconducting order parameter.","pith_inferences":["If the dimensional reduction persists at finite $J'$, the model predicts a monotonic rise of the CCDW onset temperature with layer number even at realistic coupling ratios; comparing Bi2201, Bi2212, and Bi2223 at matched hole doping would test this directly.","The model's rigid-lattice, classical-charge description omits lattice fluctuations and quenched disorder; a testable consequence is that the sharpest CCDW onset and the largest THz emission enhancement should occur in the most stoichiometric underdoped samples.","Because $J'$ is set by the interlayer oxygen spacing $d'$, c-axis pressure or strain should tune the CCDW stabilization: reducing $d'$ raises $J'$ and should move the charge-order onset upward."],"forward_implications":["Bilayer Bi2212 and trilayer Bi2223 should exhibit commensurate charge order that is substantially more stable than in monolayer Bi2201, with effective intralayer couplings $2J$ and $3J$, respectively.","The nodal, V-shaped density of states that defines the pseudogap should be understood as the CCDW's $d_{x^2-y^2}$-symmetric charge pattern, not as evidence for a $d_{x^2-y^2}$ superconducting gap.","Narrow-linewidth terahertz emission from underdoped Bi2212 mesas is consistent with a pinned commensurate charge order that cannot drift and disrupt synchronized Josephson emission, unlike an incommensurate CDW.","Phase-sensitive experiments that report a $d_{x^2-y^2}$-looking pattern in Bi2212 need re-examination: if they probe the CCDW's nodal pattern, they constrain the charge order rather than the superconducting order parameter.","In the $J'\\to\\infty$ limit, the charge-order transition of an $N$-layer stack is the 2D Ising transition at coupling $NJ$, so the ordering temperature should rise with layer number when interlayer coupling dominates."],"supporting_citations":[{"why":"Supplies the layered-superconductor background and the structural fact that the interlayer oxygen distance $d'$ is smaller than the intralayer distance $d$, so $J'>J$.","marker":"[1]"},{"why":"Provides the transition-metal-dichalcogenide precedent of a commensurate CDW on chalcogen sites stabilized by interlayer coupling, motivating the oxygen-site CCDW.","marker":"[4]"},{"why":"Provides the STM data showing an isotropic nodeless superconducting gap and a separate nodal pseudogap persisting above $T_c$ in Bi2212.","marker":"[28]"},{"why":"Provides the 45-degree twisted-flake Josephson experiment across the superconducting dome, used to argue the superconducting order parameter is not $d_{x^2-y^2}$.","marker":"[37]"},{"why":"Supplies the twisted single-layer Bi2201 Josephson result with a substantial s-wave component, extending the phase-sensitive evidence to the monolayer compound.","marker":"[41]"},{"why":"Gives the high-$T_c$ intrinsic Josephson junction emitter data whose narrow linewidths and underdoped enhancement the CCDW argument is invoked to explain.","marker":"[43]"},{"why":"Shows the pseudogap and superconducting gap in Bi2201 scale differently with temperature, supporting the pseudogap as a distinct ordered state above $T_c$.","marker":"[50]"},{"why":"Provides the exact solution of the single-layer 2D Ising model, the target that the effective single-layer reduction reproduces with coupling $NJ$.","marker":"[53]"}],"fun_headline_variants":["Layer count N multiplies charge-order coupling","Stacking planes boosts charge order by factor N","CuO2 stacking turns coupling J into NJ","Multilayer cuprates: charge coupling scales as N","Charge-density waves get N-fold boost from layer stacking"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The stabilization argument is derived in the $J'\\to\\infty$ limit with $J$ held fixed, but real Bi2212 and Bi2223 have interlayer and intralayer oxygen-oxygen spacings that differ only modestly, so $J'/J$ is of order unity; the paper does not provide a finite-$J'$ calculation showing that the effective coupling $NJ$ and the enhanced stabilization survive at realistic coupling ratios.","fun_headline_variants_meta":{"raw":{"variants":["Layer count N multiplies charge-order coupling","Stacking planes boosts charge order by factor N","CuO2 stacking turns coupling J into NJ","Multilayer cuprates: charge coupling scales as N","Charge-density waves get N-fold boost from layer stacking"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000733,"raw_usage":{"total_tokens":3379,"prompt_tokens":1148,"completion_tokens":2231,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":764,"completion_tokens_details":{"reasoning_tokens":2157}},"tokens_in":764,"tokens_out":2231,"duration_ms":17510,"temperature":1.0,"reasoning_tokens":2157,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:16:42.408140+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the onset temperature of the commensurate charge order in Bi2201, Bi2212, and Bi2223 at matched hole doping: the model predicts a systematic rise with layer number, governed by the effective coupling $NJ$ in the strong-coupling limit; if the onset does not rise with $N$, or rises far less than the $N$-fold enhancement, the central claim is falsified. A second check is to measure or compute $J'/J$ from the oxygen-oxygen distances and screening: if $J'/J$ is not large compared with 1, the strong-coupling identity is not the governing regime for these materials.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the layered-superconductor background and the structural fact that the interlayer oxygen distance $d'$ is smaller than the intralayer distance $d$, so $J'>J$."},{"cited_title":"Zhong, Y","cited_arxiv_id":null,"evidence_quote":"Provides the STM data showing an isotropic nodeless superconducting gap and a separate nodal pseudogap persisting above $T_c$ in Bi2212."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 45-degree twisted-flake Josephson experiment across the superconducting dome, used to argue the superconducting order parameter is not $d_{x^2-y^2}$."},{"cited_title":"We note that all overdoped samples studied to date have excess oxygen in their grain boundaries","cited_arxiv_id":null,"evidence_quote":"Supplies the twisted single-layer Bi2201 Josephson result with a substantial s-wave component, extending the phase-sensitive evidence to the monolayer compound."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the high-$T_c$ intrinsic Josephson junction emitter data whose narrow linewidths and underdoped enhancement the CCDW argument is invoked to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows the pseudogap and superconducting gap in Bi2201 scale differently with temperature, supporting the pseudogap as a distinct ordered state above $T_c$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the exact solution of the single-layer 2D Ising model, the target that the effective single-layer reduction reproduces with coupling $NJ$."}],"review_version":1}