{"id":"7c445f81-67dc-40aa-86f4-0aeef6d36274","arxiv_id":"2607.06462","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":6,"one_line_summary":"CDMFT calculations on the Emery model show that charge-transfer gap size and oxygen hole content are two independent mechanisms controlling the maximum superconducting Tc, with oxygen hole content being the dominant driver.","lead":"This paper uses a computational method (CDMFT) to study the Emery model of copper-oxygen planes and finds that two quantities — the charge-transfer gap size and the oxygen hole content — independently control the maximum superconducting temperature. A smart generalist would read this because it identifies which microscopic parameters to tune to raise Tc in cuprate superconductors and in future cold-atom simulators.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The 'dominant variable' claim rests on mean-field Tc; the superfluid stiffness that sets the true 2D Tc may be more sensitive to Δ_I than to 2p_p, potentially rotating the gradient in Fig. 12b.","rationale":"The reader correctly identified the KT physics limitation as the weakest assumption. My analysis confirms this is the most load-bearing concern, and I have made it more specific: the issue is not just that KT fluctuations reduce Tc generically, but that the correction is likely Δ_I-dependent because the superfluid stiffness is directly tied to the charge-transfer gap through quasiparticle coherence. This specifically threatens the 'dominant variable' part of the central claim while leaving the 'two independent mechanisms' part relatively secure. The reader's CONDITIONAL verdict with MODERATE confidence already accounts for this concern appropriately. Other issues (8 data points, narrow parameter range, sign problem blocking the most relevant regime, quantitative disagreements with experiment) are real but secondary: they weaken the precision of all conclusions equally, whereas the KT/stiffness concern selectively targets the dominance ranking. The paper's strengths — honest acknowledgment of limitations, standard methodology, internally consistent trends across 8 points, and the non-trivial observation that Δ_I and 2p_p contour lines are not parallel in the Zaanen-Sawatzky-Allen diagram — support a conditional acceptance. The concrete test I propose is computationally feasible within the existing CDMFT framework and would directly settle whether the dominance conclusion survives KT corrections.","tokens_in":35190,"tokens_out":6417,"duration_ms":335074,"concrete_test":"Compute the superfluid stiffness ρ_s at the CDMFT level (via the static current-current correlator or the second derivative of the free energy with respect to a uniform vector potential) for the same 8 parameter points. Construct the gradient of ρ_s(Δ_I, 2p_p) and compare its direction to the gradient of Tc^CDMFT(Δ_I, 2p_p) shown in Fig. 12(b). If the two gradient directions differ by more than ~30°, the KT-corrected Tc would likely reorder the dominance conclusion, weakening the claim that 2p_p is the dominant variable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that oxygen hole content 2p_p is the dominant variable (Sec. V.B, Fig. 12b) is derived from the gradient of Tc^CDMFT, which measures the onset of nonzero pairing order parameter Φ. In two dimensions, the physical transition temperature is set by the Kosterlitz-Thouless transition, T_KT ∝ ρ_s(T_KT), where ρ_s is the superfluid stiffness. The superfluid stiffness depends on quasiparticle coherence and effective bandwidth, both of which are directly controlled by the charge-transfer gap Δ_I: increasing Δ_I deepens the charge-transfer regime, reducing the quasiparticle weight Z and hence ρ_s. This means the KT correction to Tc^CDMFT is likely Δ_I-dependent — larger for larger Δ_I — which would steepen the effective Tc vs Δ_I gradient relative to the Tc vs 2p_p gradient. If this correction is large enough, the gradient direction in Fig. 12(b) could rotate toward the Δ_I axis, changing the conclusion about which variable dominates. The authors acknowledge that KT fluctuations 'can alter the Tc^max trends' (Sec. V.C) but do not assess whether the correction is preferentially Δ_I-dependent in a way that would specifically threaten the dominance claim. This is the softest point: the independence of the two mechanisms (part 1 of the claim) is established from parent-state properties and is robust to KT corrections, but the dominance ranking (part 2) is not.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This paper uses cellular dynamical mean-field theory (CDMFT) to study the maximum superconducting critical temperature $T_c^{max}$ in the Emery model for hole-doped cuprates. Using the Zaanen-Sawatzky-Allen diagram as a framework, the authors fix the hopping parameters ($t_{pd}=1.5$, $t'_{pp}=1$, $t_{pp}=1$) and systematically vary the bare charge-transfer energy $U_d$ and $d-p$ energy distance. They compute $T_c^{CDMFT}$ for eight parent insulating states corresponding to two charge gap sizes ($Delta_I = 0.6, 0.9$). The central claim is that $Delta_I$ and the oxygen hole content $2p_p$ are two independent mechanisms controlling $T_c^{max}$, with $T_c^{max}$ monotonically increasing as $Delta_I$ decreases and $2p_p$ increases. The authors further claim that $2p_p$ is the dominant variable, as inferred from the gradient of the interpolated $T_c^{max}$ surface in the $(Delta_I, 2p_p)$ plane (Fig. 12).","tokens_in":35409,"tokens_out":1428,"duration_ms":430668,"significance":"The paper tackles an important problem: disentangling the roles of charge-transfer gap size and oxygen hole content in determining $T_c$ in cuprates. The CDMFT methodology is standard and well-implemented, and the authors are transparent about limitations (Sec. V.C). A clear strength is the direct computation of $T_c^{CDMFT}$ rather than relying solely on the superconducting order parameter $Phi$ as a proxy, allowing the authors to verify (Fig. 9b, Fig. 10b) that $Phi^{max}$ and $T_c^{max}$ can exhibit different rates of change. The identification of two distinct, non-proportional mechanisms for controlling $T_c$ is a valuable conceptual contribution. The falsifiable predictions for cold-atom simulators are timely.","major_comments":[{"comment":"Sec. V.B, Fig. 12: The central claim that $2p_p$ is the 'dominant variable' rests on the gradient of the interpolated $T_c^{max}$ surface. However, the gradient direction is not invariant under rescaling of the axes. The text states that $Delta_I$ varies in [0.6, 0.9] (absolute change 0.3) while $2p_p$ varies in [0.26, 0.33] (absolute change 0.07). Because the axes have 'equal aspect ratio' but the physical quantities have different units and ranges, the gradient direction is an artifact of the chosen axis scaling. If $Delta_I$ were plotted in eV (as in the secondary axis of Fig. 10, where $t_{pp}=0.4$ eV gives a range of ~0.24 to 0.36 eV), the aspect ratio and hence the gradient direction would change. The dominance claim requires a parameter-free metric, such as comparing $partial T_c^{max}/partial Delta_I$ and $partial T_c^{max}/partial (2p_p)$ in physical units, or a regression-based","section":null},{"comment":"sensitivity analysis. As presented, the gradient in Fig. 12(b) is not a physically meaningful measure of relative importance.","section":null},{"comment":"Sec. V.C and Fig. 12: The $T_c^{max}$ surface in Fig. 12 is an interpolation over only 8 data points (5 at $Delta_I=0.6$, 3 at $Delta_I=0.9$). The gradient field and contour lines are thus derived from a very sparse grid, and the monotonicity claim is supported by limited sampling. The authors should explicitly state the interpolation method used and add uncertainty estimates or error bars to the data points in Fig. 12, so that the robustness of the gradient direction and the monotonicity claim can be assessed.","section":null}],"minor_comments":[{"comment":"Sec. V.C: The acknowledgment that KT fluctuations 'can alter the $T_c^{max}$ trends' is appreciated, but the specific risk to the dominance claim is not addressed. The superfluid stiffness $rho_s$ depends on quasiparticle coherence and bandwidth, which are directly controlled by $Delta_I$. If the KT correction is preferentially $Delta_I$-dependent, it could rotate the effective gradient in Fig. 12(b). A brief discussion of this specific risk would strengthen the paper.","section":null},{"comment":"Fig. 2(a): The labels for constant $2p_p$ lines list values from 0.24 to 0.28, but the text in Sec. III.A.2 states the step is 0.1, which should presumably be 0.01.","section":null},{"comment":"Fig. 9(a): The $T_c^{max}$ values are converted to Kelvin using $t_{pp}=0.4$ eV, yielding values around 96-108 K. However, Sec. V.C notes that this is approximately twice the experimental $T_c$ for LSCO. This discrepancy should be noted on the figure or in its caption to avoid confusion.","section":null},{"comment":"Sec. III.A.2: The oxygen hole content is defined as $2p_p = 2(2-n_p)$. It would help to clarify early on whether $p_p$ refers to holes per oxygen orbital or per CuO$_2$ unit cell, as the factor of 2 for orbital degeneracy can be a source of confusion.","section":null},{"comment":"Appendix B, Fig. 14: The data points from Ref. [30] use different hopping amplitudes, as noted in the text. It would be helpful to indicate which specific parameters differ so the reader can assess comparability.","section":null}],"recommendation":"major_revision","confidential_remarks":"The core issue with the gradient direction in Fig. 12 is load-bearing for the paper's central claim of dominance. The authors need to either provide a scaling-invariant metric for relative importance or significantly soften the dominance claim. The reliance on a companion paper [39] for the $Delta_I=0.6$ results is acceptable given that this paper extends to $Delta_I=0.9$ and performs the 2D optimization, but the interpolation sparsity remains a concern. The KT physics concern raised in the stress-test is valid but secondary to the axis-scaling problem, which is a more immediate and fixable issue within the manuscript's scope."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful reading and for identifying two important issues with our analysis in Sec. V.B. Both points are well-taken and will be addressed in the revised manuscript.","responses":[{"response":"The referee is correct that the gradient direction of an interpolated surface is not invariant under rescaling of the axes, and that our claim of 2p_p being the 'dominant variable' as inferred from the gradient in Fig. 12(b) is therefore not well-justified as presented. We acknowledge this as a genuine weakness in our argumentation. We will revise the manuscript to address this in two ways. First, we will remove or substantially soften the claim that 2p_p is the 'dominant variable' as inferred from the gradient direction, since the gradient is indeed scale-dependent and not a parameter-free metric. Second, we will provide a quantitative comparison using the actual data points rather than the interpolated surface. Specifically, from the raw data: (i) at fixed Δ_I = 0.6, T_c^max changes by approximately ΔT_c ≈ 0.0015 (in units of t_pp) as 2p_p varies over a range of ~0.05; (ii) at approximately fixed 2p_p, T_c^max changes by approximately ΔT_c ≈ 0.001 as Δ_I varies over a range of 0.3. We will compute the partial derivatives ∂T_c^max/∂Δ_I and ∂T_c^max/∂(2p_p) directly from the data and compare them in physical units (converting to eV using t_pp = 0.4 eV). We will also perform a simple linear regression of T_c^max on both variables to assess relative importance via standardized coefficients. We note that with only 8 data points, any such quantitative assessment will carry significant uncertainty, which we will state explicitly. The qualitative finding that both Δ_I and 2p_p independently control T_c^max (demonstrated by the fact that the two curves in Fig. 11(a) do not overlap) does not depend on the gradient argument and remains valid. However, the quantitative ranking of which variable is 'dominant' requires the more careful analysis we will add.","revision_made":"yes","referee_comment":"Sec. V.B, Fig. 12: The central claim that 2p_p is the 'dominant variable' rests on the gradient of the interpolated T_c^max surface. However, the gradient direction is not invariant under rescaling of the axes... The dominance claim requires a parameter-free metric, such as comparing partial T_c^max/partial Delta_I and partial T_c^max/partial (2p_p) in physical units, or a regression-based sensitivity analysis. As presented, the gradient in Fig. 12(b) is not a physically meaningful measure of relative importance."},{"response":"The referee is correct on all counts. We will make the following revisions. First, we will explicitly state the interpolation method: we used a standard bilinear interpolation on a uniform grid, as implemented in matplotlib's contourf/contour routines. We will state this in the figure caption and in the text. Second, we will add error bars to the data points in Fig. 12. The uncertainty in T_c^max arises from two sources: (a) the finite temperature grid used to locate the superconducting transition (the temperature spacing is Δβ = 0.5–1.0, corresponding to ΔT ≈ 0.0004–0.0008 in units of t_pp), and (b) statistical Monte Carlo errors on the order parameter Φ, which affect the identification of the transition point. We will estimate and display these error bars. Third, we will add an explicit caveat that the monotonicity claim is based on only 8 data points and that additional data points, particularly at intermediate values of Δ_I (e.g., Δ_I = 0.75) and at additional values of 2p_p within each Δ_I slice, would be needed to robustly establish monotonicity. We note that the sign problem, as discussed in Sec. V.C and Appendix C, limits our ability to reach lower temperatures and thus constrains the density of data we can obtain. We will also add a note that the interpolated surface and gradient field in Fig. 12(b) should be interpreted as a guide to the eye rather than a quantitative reconstruction, especially given the sparsity of the grid.","revision_made":"yes","referee_comment":"Sec. V.C and Fig. 12: The T_c^max surface in Fig. 12 is an interpolation over only 8 data points (5 at Delta_I=0.6, 3 at Delta_I=0.9). The gradient field and contour lines are thus derived from a very sparse grid, and the monotonicity claim is supported by limited sampling. The authors should explicitly state the interpolation method used and add uncertainty estimates or error bars to the data points in Fig. 12, so that the robustness of the gradient direction and the monotonicity claim can be assessed."}],"tokens_in":34997,"tokens_out":1452,"duration_ms":71877,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper extends the authors' companion work [39] from one charge-gap value to two, enabling a 2D optimization of Tc^max over both charge-transfer gap size Δ_I and oxygen hole content 2p_p. The genuinely new result is the demonstration that these are two independent mechanisms — Fig. 11(a) shows that the same 2p_p can give different Tc^max at different Δ_I, and vice versa. The gradient analysis in Fig. 12(b) claiming 2p_p is the dominant variable is also new, but this is where I'd put an asterisk (more below). The paper does several things well. The CDMFT methodology is standard and competently executed. The Zaanen-Sawatzky-Allen framework is a clean organizing principle for navigating the parameter space. The authors are unusually transparent about limitations — they openly state the KT caveat, the sign problem blocking access to the most relevant parameter region, the factor-of-2 discrepancy in Tc, the factor-of-3 discrepancy in gap size, and even the wrong sign on the pd-vs-2pp slope compared to experiment. That honesty earns trust. The stress-test concern about KT physics is valid but needs calibration. The concern is that superfluid stiffness ρ_s depends on quasiparticle weight Z, which is Δ_I-dependent, so the KT correction could be preferentially larger for bigger Δ_I, potentially rotating the gradient in Fig. 12(b) toward the Δ_I axis. This is a real concern, but it attacks only the dominance ranking (part 2 of the claim), not the independence of the two mechanisms (part 1). The independence result is derived from parent-state properties and is robust to KT corrections. So the paper's strongest claim survives; its flashiest claim is on shakier ground. The sparsity of data points (8 points interpolated onto a surface) is a legitimate concern but proportionate to what the authors claim — they frame this as an entry point, not a final word. This paper is for researchers working on cuprate theory and cold-atom simulators who want an organizing framework for Tc trends. It deserves a serious referee who can assess whether the KT caveat undermines the dominance claim or merely qualifies it. My read: the independence result is solid, the dominance ranking is suggestive but not yet robust. Recommend peer review.","headline":"Two-mechanism Tc optimization in the Emery model: real but preliminary, with a load-bearing caveat on the dominance claim","tokens_in":36263,"tokens_out":557,"would_cite":true,"duration_ms":72991,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.72.-h","74.20.Mn","71.27.+a","71.30.+h"],"model":"glm-5.2","headline":"Oxygen holes drive superconducting temperature in cuprate model","keywords":["Emery model","cuprate superconductivity","charge-transfer gap","oxygen hole content","Zaanen-Sawatzky-Allen diagram","cellular dynamical mean-field theory","d-wave superconductivity"],"falsifier":"If including Kosterlitz-Thouless phase fluctuations were to shift the optimal doping differently for different charge gap sizes, the ranking of oxygen hole content as the dominant variable could change.","tokens_in":35235,"feed_emoji":"🔑","tokens_out":1273,"duration_ms":150840,"temperature":0.7,"pith_summary":"This paper studies the Emery model — the standard description of electrons in the copper-oxygen planes of cuprate superconductors — using a computational method called cellular dynamical mean-field theory. The authors systematically vary two microscopic parameters — the energy difference between copper and oxygen orbitals and the strength of electron repulsion on copper — and track how the maximum superconducting transition temperature changes. They then translate these model parameters into two physically measurable quantities: the charge-transfer gap size and the oxygen hole content, which measures how much the hole is shared between copper and oxygen orbitals. The central claim is that these two quantities are independent mechanisms controlling the maximum superconducting temperature. The superconducting temperature rises monotonically as the charge gap shrinks and as the oxygen hole content grows, with oxygen hole content being the dominant driver. The authors also find that the highest temperatures are reached not only near the metal-insulator boundary, as previously known, but also deep in the charge-transfer regime where the copper repulsion is much larger than the copper-oxygen energy difference — a regime that had not been previously identified as favorable for superconductivity.","feed_headline":"Oxygen holes, not just charge gap, set superconducting temperature ceiling","feed_subtitle":"A two-variable map of the Emery model shows oxygen hole content is the dominant driver of maximum Tc, offering a design rule for cuprates.","key_machinery":"The Zaanen-Sawatzky-Allen diagram serves as the navigational map. It plots the copper onsite repulsion U_d against the bare charge-transfer energy Delta (the energy distance between the upper Hubbard band and the oxygen orbital energy). Contour lines of constant charge gap size and constant oxygen hole content are overlaid on this diagram. The key structural observation is that these two sets of contour lines are not parallel, proving that charge gap size and oxygen hole content are not proportional to each other — they are genuinely independent parameters. This independence is what allows the authors to disentangle their separate effects on the superconducting temperature.","core_discovery":"The charge-transfer gap size and the oxygen hole content are two independent mechanisms controlling the maximum superconducting transition temperature in the Emery model. The temperature increases monotonically as the charge gap decreases and as the oxygen hole content increases, and the oxygen hole content is the dominant variable. This is established by computing superconducting domes at eight points across the Zaanen-Sawatzky-Allen diagram — which classifies insulating states by the relative sizes of the copper repulsion and the copper-oxygen energy difference — at two fixed charge gap sizes, and then mapping the results onto the two-dimensional space of physical observables. The gradient","pith_inferences":["If oxygen hole content is indeed the dominant driver, then material design strategies that increase copper-oxygen covalency — for instance through chemical pressure, strain, or apical oxygen manipulation — should be more effective at raising the superconducting temperature than strategies that only narrow the charge gap.","The paper's quantitative limitations (oxygen hole content range matching only La2-xSrxCuO4, predicted temperatures about twice experimental values, optimal doping underestimated by a factor of three) suggest that extending the calculation to include Kosterlitz-Thouless physics and variable hopping parameters could shift the quantitative predictions while preserving the qualitative trend that oxyge","The observation that the parent insulating state's properties predict superconducting trends upon doping hints at a deeper organizing principle: the mixed copper-oxygen character of the ground-state hole may be a more fundamental indicator of pairing strength than any single energy scale, which would connect to the broader question of why charge-transfer insulators host higher superconducting temp"],"forward_implications":["For proposed ultracold-atom realizations of the Emery model, the paper predicts that maximizing superconducting temperature requires tuning the system deep into the charge-transfer regime (large copper repulsion relative to the copper-oxygen energy difference) while keeping the charge gap small, which corresponds to staying near the metal-insulator boundary.","The finding that oxygen hole content is the dominant variable suggests that experimental correlations between superconducting temperature and oxygen hole content in real cuprates should be robust against small variations in the charge-transfer gap, providing a theoretical basis for interpreting those measurements.","The paper provides a two-dimensional optimization framework that can be extended: varying the hopping parameters (currently fixed) would add further dimensions and could reveal whether the dominance of oxygen hole content persists across a broader parameter space.","The identification of a second favorable region deep in the charge-transfer regime, beyond the previously known optimum near the metal-insulator boundary, suggests a design principle for engineering new superconducting materials: large copper-oxygen covalency combined with a narrow charge gap."],"fun_headline_variants":["Oxygen holes dominate maximum Tc in the Emery model","Charge gap and oxygen holes set maximum Tc in the Emery model","Two mechanisms control maximum superconducting Tc in the Emery model","Oxygen hole content is the primary driver of maximum Tc","Mapping maximum Tc by charge gap and oxygen holes"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The superconducting temperature is computed at a mean-field level that neglects Kosterlitz-Thouless phase fluctuations, which are essential in two dimensions. These fluctuations could reduce the maximum temperature and shift the optimal doping, potentially reordering which variable — oxygen hole content or charge gap size — dominates.","fun_headline_variants_meta":{"raw":{"variants":["Oxygen holes dominate maximum Tc in the Emery model","Charge gap and oxygen holes set maximum Tc in the Emery model","Two mechanisms control maximum superconducting Tc in the Emery model","Oxygen hole content is the primary driver of maximum Tc","Mapping maximum Tc by charge gap and oxygen holes"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1676,"prompt_tokens":577,"completion_tokens":1099,"prompt_tokens_details":null},"tokens_in":577,"tokens_out":1099,"duration_ms":41155,"temperature":1.0,"reasoning_tokens":1053,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T04:53:35.037208+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If including Kosterlitz-Thouless phase fluctuations were to shift the optimal doping differently for different charge gap sizes, the ranking of oxygen hole content as the dominant variable could change.","supporting_citations":[],"review_version":1}