{"id":"95d31c43-9ca4-43f7-b542-c0c02159cfaa","arxiv_id":"2412.14773","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The spin-1/2 XXZ pyrochlore model exhibits a wide low-temperature non-magnetic phase between two ferromagnetic domes, including quantum spin ice, with a possible lattice nematic ground state near the antiferromagnetic XY point.","lead":"Using a recent functional renormalization group method, this paper maps the full temperature-versus-anisotropy phase diagram of a spin-1/2 XXZ magnet on the pyrochlore lattice, finding a broad low-temperature non-magnetic region that contains quantum spin ice and antiferromagnetic Heisenberg and XY behaviors. It also reports enhanced lattice-symmetry-breaking responses near the antiferromagnetic XY model, suggesting a possible nematic ground state.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Phase boundaries rely on two-size xi/L crossings in an uncontrolled low-T regime; the upper boundary is unbenchmarked and disagrees with CMFT.","rationale":"After reading the manuscript and the reader's verdict, I find the principal risk to the central claim is the uncontrolled low-temperature extrapolation of the one-loop T-flow PMFRG, combined with a weak finite-size analysis. The authors explicitly state in Sec. III A that quantitative accuracy is expected to degrade below T ~ J, yet most of the phase boundaries are located at T < 0.1J. The QSI0-FM_perp benchmark is reassuring but is a single point; the upper boundary conflicts with the CMFT result (0.554pi vs 0.613pi) and is not externally checked. The use of only two cluster sizes for the xi/L collapse means the 'collapse' is just a crossing, not a validation of scaling. A three-size scaling fit or a convergence check with more Matsubara frequencies would either support or refute the reported boundaries. I do not see an internal mathematical contradiction; the concern is about numerical control. Since the reader already marks CONDITIONAL with MODERATE confidence, my read does not change the verdict. I would add the concrete test above as a prerequisite for acceptance.","tokens_in":21703,"tokens_out":4258,"duration_ms":31303,"concrete_test":"Re-analyze the xi/L data in Fig. 3 with a three-size finite-size scaling fit (N=459, 1029, 1941) to extract Tc and nu for each theta boundary. If the best-fit Tc deviates from the reported two-size crossing by more than 10%, or if the fit quality is poor, the phase boundaries are not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central phase diagram is load-bearing on the claim that one-loop T-flow PMFRG remains quantitatively accurate at T << J, a regime the authors themselves state is beyond perturbative control (Sec. III A: 'the PMFRG is perturbatively controlled only at T >> J'). The only external benchmark is the QSI0-FM_perp boundary (theta=-0.027pi vs QMC -0.033pi); the upper boundary theta=0.554pi has no independent check and differs from CMFT (0.613pi). Moreover, critical temperatures are determined from a collapse of xi/L using only the two largest system sizes, N=1029 and 1941, as stated in Sec. III A and Fig. 3. With two sizes, a 'collapse' reduces to a curve crossing: it cannot validate the scaling form, the universality class, or provide an error estimate. A systematic truncation error or an unreliable crossing would shift the reported boundaries and could reverse the claimed quantum order-by-disorder at the upper boundary.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the temperature-coupling phase diagram of the spin-1/2 XXZ pyrochlore model using the pseudo-Majorana functional renormalization group in its temperature-flow variant. It reports a large low-temperature non-magnetic regime between θ≈-0.027π and θ≈0.554π, containing the antiferromagnetic Ising, Heisenberg, and XY limits, surrounded by FM⊥ and FMz ferromagnetic phases. It also presents spin structure factors showing pinch-point patterns, and linear-response calculations suggesting a lattice nematic candidate for the antiferromagnetic XY model. The QSI0-to-FM⊥ transition agrees with quantum Monte Carlo, while the upper boundary disagrees with cluster mean-field theory.","tokens_in":21859,"tokens_out":4653,"duration_ms":30608,"significance":"If the phase boundaries are accurate, this is the first unified finite-temperature phase diagram for the full XXZ pyrochlore family, and the identification of a candidate lattice nematic state in the XY model would be an important step for a heavily studied frustrated magnet. The paper is commendably explicit about the method's regime of perturbative control and about the limitations of the nematic response functions; the agreement with QMC at the lower boundary and the comparison with classical Monte Carlo give useful external anchors. However, the central diagram currently lacks error estimates and the upper boundary rests on a single approximate method, so the quantitative claims should be treated with caution until additional benchmarks are provided.","major_comments":[{"comment":"Critical temperatures are read off from a collapse of ξµ/L using only the two largest system sizes (N=1029 and 1941). With only two sizes, the collapse criterion is a crossing of two curves; it does not test the scaling form, determine the universality class, or yield an error estimate. Since every phase boundary in Fig. 1(a) is obtained this way, the phase diagram has no reported uncertainty. I ask the authors to add at least one more system size, perform a systematic finite-size scaling analysis, or explicitly quantify the resulting uncertainty in θ and T_c; without this, the precision implied by the phase diagram (e.g., θ=-0.027π vs QMC -0.033π) is not established.","section":"Sec. III A and Fig. 3"},{"comment":"The upper boundary of the non-magnetic phase, θ=0.554π, is not benchmarked by any independent quantum method and deviates from the cluster mean-field result 0.613π by roughly 0.06π. The text interprets this as 'a larger stability of the FMz phase within PMFRG,' but given the paper's own statement in Sec. III A that the method is perturbatively controlled only for T≫J and in Sec. V that it 'is expected to lose quantitative accuracy at T<J,' a method error of this magnitude is equally plausible. The quantum order-by-disorder claim at the upper boundary needs either an independent cross-check (e.g., QMC where sign-problem free, or a different FRG truncation) or a demonstrated convergence with respect to the Matsubara cutoff and system size.","section":"Sec. IV A and Table II"},{"comment":"The proposal of a lattice nematic ground state for the antiferromagnetic XY model is based exclusively on the linear response χ⊥R to seed fields. As the manuscript correctly notes in Sec. III A, large χ_R is not rigorously connected to nematic long-range order, and the method cannot detect the nematic transition. Moreover, the response functions cannot distinguish the proposed C3-nematic from the in-plane spin nematic of Ref. [23]. The conclusion in Sec. IV C that 'the pyrochlore XY model shows strong tendencies for realizing a nematic ground state' therefore overstates what the current data can establish; please reframe as a candidate consistent with the data and explicitly list the alternative orders.","section":"Sec. IV C and Fig. 5"}],"minor_comments":[{"comment":"The word 'antiferromagntic' is a typo for 'antiferromagnetic'.","section":"Sec. II"},{"comment":"The phrase 'steepest descend' should be 'steepest descent'.","section":"Eq. (4) and surrounding text"},{"comment":"The numerical specification sentence leaves unclear whether the 48 positive Matsubara frequencies are used for all three system sizes; please clarify the numerical setup.","section":"Sec. III A"},{"comment":"The phase boundary values are given to three decimals (e.g., -0.027π, 0.554π) but no error bars are shown; please add error bars or a statement that the last digit is not significant.","section":"Fig. 1(a) and Table II"},{"comment":"The red crosses marking the locations of the structure factors in Fig. 1(a) are mentioned in the caption but are not visible in the figure as provided; please ensure all structure-factor points are clearly marked.","section":"Sec. IV B"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely a good fit for the journal, but the phase-diagram claims need strengthening before publication. I would ask the editor to weigh whether the lack of error bars and the single-method upper boundary are acceptable given the paper's emphasis on benchmarking. The authors have been transparent about the method's limitations, which is a positive sign."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline: this paper gives the first T-flow PMFRG phase diagram for the spin-1/2 XXZ pyrochlore model, and its most convincing result is the lower boundary of the non-magnetic region, where the predicted FM_perp-to-QSI0 transition at theta = -0.027 pi sits close to the QMC value -0.033 pi. That benchmark is real. The rest of the phase map is plausible but methodologically less certain.\n\nWhat is genuinely new: the complete theta-T map, the claim of quantum order-by-disorder at the upper disorder-to-FM_z boundary (theta = 0.554 pi vs classical 0.60 pi), and the enhanced C3 and combined I/C3 nematic responses in the AF XY region. The structure factor analysis is a solid extra: broadened pinch points and the QSI0/QSIpi distinction come out nicely, and the comparison with classical Monte Carlo is useful.\n\nThe soft spots are the usual ones for PMFRG pushed to T << J. The authors are honest that the method is perturbatively controlled only for T >> J, and the two largest system sizes (N = 1029, 1941) are used to define a \"scaling collapse\" that is really just a crossing of xi/L curves. Two sizes cannot validate the scaling form or give error bars, so the boundaries in Fig. 1(a) are effectively unquantified. That is a real limitation, but the one QMC benchmark at the lower boundary keeps the method from looking arbitrary. The upper boundary is the bigger worry: it disagrees with the CMFT value 0.613 pi, and the claimed order-by-disorder effect there rests on that unbenchmarked boundary. The XY nematic proposal is also indirect - it is a linear response to seed fields, not a direct four-spin correlator, and the authors themselves note they cannot rule out other interpretations. Minor issues: the Fig. 1(b) caption says the Ising S_z(Q) is a constant 1/4 while the text says it has pinch points; Table I has a dangling \"4.21 [81]\" entry that should be keyed to the discrete-Ising normalization. No code or data are released.\n\nWho is this for: practitioners in frustrated pyrochlore magnetism and fRG method developers. It is not a slam dunk, but it is a serious, publishable map with at least one externally checked boundary. For peer review, I would send it out, not desk reject it. The referee should ask for three system sizes (or a consistency check), error estimates on the critical values, and ideally some additional benchmark near the upper boundary. Artifact release would help but is not a blocker.\n\nNet: worth the field's attention; expect heavy revision.","headline":"First full T-flow PMFRG phase diagram for the XXZ pyrochlore magnet; the lower boundary is QMC-validated, but the upper boundary and nematic proposal rest on weaker evidence.","tokens_in":22439,"tokens_out":4828,"would_cite":true,"duration_ms":33324,"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":"A temperature-flow renormalization group maps the full phase diagram of the spin-1/2 XXZ pyrochlore model and argues that the antiferromagnetic XY limit favors a lattice nematic ground state.","keywords":["XXZ pyrochlore model","pseudo-Majorana functional renormalization group","quantum spin ice","temperature-flow formalism","pyrochlore lattice","spin structure factor","lattice nematic order","frustrated magnetism"],"falsifier":"A definitive check would be an unbiased low-temperature calculation of the same boundaries, for example sign-free quantum Monte Carlo in the unfrustrated transverse sector or a tensor-network or high-order series calculation near $\\theta=0.554\\pi$, looking for whether the upper non-magnetic boundary sits near the PMFRG value or the cluster mean-field value $0.613\\pi$, and whether the QSI$_0$/FM$_\\perp$ boundary remains at $\\theta\\approx-0.033\\pi$. Within PMFRG itself, rerunning the scaling analysis with a fourth system size near $N=4000$ would reveal whether the $\\xi/L$ collapse of the two largest sizes continues or drifts.","tokens_in":21470,"feed_emoji":"🧲","tokens_out":7132,"duration_ms":54625,"temperature":0.7,"pith_summary":"The paper sets out the full temperature–coupling phase diagram of the spin-1/2 nearest-neighbour XXZ model on the pyrochlore lattice, a central frustrated three-dimensional magnet. Using the temperature-flow pseudo-Majorana functional renormalization group, it finds a large low-temperature non-magnetic region from $\\theta=-0.027\\pi$ to $\\theta=0.554\\pi$ that contains the zero-flux quantum spin ice phase near the antiferromagnetic Ising limit, the antiferromagnetic Heisenberg point at $\\theta=\\pi/4$, and the antiferromagnetic XY point at $\\theta=\\pi/2$. This matters because it places several long-studied models into one phase diagram and shows where quantum fluctuations destroy magnetic order. The lower boundary agrees with quantum Monte Carlo, and the study proposes a lattice nematic ground state for the antiferromagnetic XY model based on enhanced symmetry-breaking responses.","feed_headline":"Pyrochlore spin model stays non-magnetic all the way from Ising to XY","feed_subtitle":"One phase diagram covers quantum spin ice, Heisenberg, and XY limits down to T = 0.01J.","key_machinery":"The workhorse is the pseudo-Majorana functional renormalization group in the temperature-flow variant (T-flow PMFRG), in which spins are rewritten as three Majorana fermions and the renormalization group cutoff is the physical temperature, so a single numerical run covers the whole cooling history. The flow equations for the self-energy and two-particle vertex are solved under the one-loop truncation. Magnetic transitions are located by a finite-size scaling collapse of the correlation length ratio $\\xi^\\mu/L$ computed from the peak of the spin structure factor $S^\\mu(Q)$ for system sizes $N=459,1029,1941$. Nematic tendencies are probed by linear response functions $\\chi^\\mu_R$ to perturbations that strengthen and weaken bonds in patterns breaking $C_3$, inversion, or both.","core_discovery":"The central claim is that the spin-1/2 XXZ pyrochlore model, with couplings parametrized as $J_\\perp=J\\sin\\theta$ and $J_z=J\\cos\\theta$, has a broad low-temperature non-magnetic phase for $\\theta\\in[-0.027\\pi,\\,0.554\\pi]$. Inside this phase sit the zero-flux quantum spin ice (QSI$_0$) regime near the antiferromagnetic Ising limit, the antiferromagnetic Heisenberg model at $\\theta=\\pi/4$, and the antiferromagnetic XY model at $\\theta=\\pi/2$. The phase is bounded on one side by ferromagnetic order in the $xy$-plane and on the other by Ising-type ferromagnetic order along $z$; the QSI$_0$-to-FM$_\\perp$ transition at $\\theta=-0.027\\pi$ agrees with the quantum Monte Carlo value $\\theta=-0.033\\pi$. The paper further argues, from growing linear responses to $C_3$ and combined $C_3$/inversion symmetry-breaking perturbations, that the antiferromagnetic XY model tends to a lattice nematic ground state. In the disordered region the longitudinal and transverse structure factors show broadened pinch points that sharpen as magnetic order is approached.","pith_inferences":["Beyond the paper, the same response-function logic could be extended to test in-plane spin-nematic order ($J_x\\ne J_y$) once the method is generalized to XYZ couplings; the authors already identify this as the natural next step.","If the lattice-nematic proposal survives, it would unify the classical thermal order-by-disorder selection of collinear states with a quantum analogue in the XY limit, giving a single mechanism across classical and quantum versions.","The discrepancy at the upper boundary ($0.554\\pi$ vs $0.613\\pi$) is the sharpest place to discriminate between approximations; a future unbiased calculation there would also calibrate how much the one-loop truncation underestimates quantum fluctuations.","The observed leakage of transverse spin correlations into the longitudinal structure factor at the XY point suggests that any experimental probe of $S^z$ in an XY-like material would see quantum-fluctuation-induced signal even with $J_z=0$, which may be a useful diagnostic."],"forward_implications":["If the phase diagram is correct, the zero-flux quantum spin ice phase occupies only a narrow window near the antiferromagnetic Ising limit, and antiferromagnetic Heisenberg and XY pyrochlore magnets stay non-magnetic down to $T\\approx0.01J$.","The agreement at the lower boundary would show that T-flow PMFRG can locate phase transitions well beyond its perturbatively controlled regime, making it a usable tool for other strongly frustrated three-dimensional magnets.","The identified quantum order-by-disorder at both boundaries means quantum fluctuations shrink the non-magnetic region relative to the classical model, reversing earlier variational expectations at the upper boundary.","Broadened pinch points in the structure factors give momentum-resolved fingerprints that neutron scattering or other probes could look for in pyrochlore materials.","A lattice nematic ground state of the antiferromagnetic XY model would resolve the competition between spin-nematic and QSI$_\\pi$ proposals for this limit."],"supporting_citations":[{"why":"Quantum Monte Carlo locating the QSI$_0$/FM$_\\perp$ transition, the benchmark for the lower phase boundary.","marker":"[33]"},{"why":"Quantum Monte Carlo value $\\theta=-0.033\\pi$ used as the reference for the same transition.","marker":"[34]"},{"why":"Cluster mean-field and series study providing the upper boundary $\\theta=0.613\\pi$ and the spin-nematic proposal that the PMFRG result is compared against.","marker":"[23]"},{"why":"Gauge mean-field theory predicting distinct transverse structure-factor maxima for QSI$_0$ and QSI$_\\pi$, used to identify both regimes.","marker":"[61]"},{"why":"PFFRG calculation of lattice-symmetry-breaking response functions in the pyrochlore Heisenberg model, the baseline for the nematic-response analysis.","marker":"[48]"},{"why":"Introduces the pseudo-Majorana functional renormalization group at finite temperature and its asymptotic exactness at high $T/|J|$.","marker":"[37]"},{"why":"Introduces the temperature-flow scheme that lets a single run cover the full cooling range down to small $T$.","marker":"[39]"},{"why":"Quantitative functional renormalization for three-dimensional Heisenberg models, source of the correlation-length scaling-collapse transition detection.","marker":"[30]"},{"why":"Classical XXZ pyrochlore phase diagram and order-by-disorder selection, used for the quantum-versus-classical comparison.","marker":"[24]"},{"why":"Quantum Monte Carlo critical temperature for the ferromagnetic Heisenberg model used in the benchmark table.","marker":"[82]"}],"fun_headline_variants":["Pyrochlore spin model stays non-magnetic from Ising to XY","One phase diagram for quantum spin ice, Heisenberg, and XY","Low-temperature pyrochlore order: spin liquid to possible nematic","FRG maps pyrochlore quantum phases down to T=0.01J"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the one-loop truncated, temperature-flow equations remain quantitatively accurate down to $T\\approx0.01J$, far below the $T\\gg J$ regime where the expansion is controlled, and that a scaling collapse using only the two largest system sizes identifies genuine second-order transitions.","fun_headline_variants_meta":{"raw":{"variants":["Pyrochlore spin model stays non-magnetic from Ising to XY","One phase diagram for quantum spin ice, Heisenberg, and XY","Low-temperature pyrochlore order: spin liquid to possible nematic","FRG maps pyrochlore quantum phases down to T=0.01J"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000277,"raw_usage":{"total_tokens":1701,"prompt_tokens":1048,"completion_tokens":653,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":664,"completion_tokens_details":{"reasoning_tokens":585}},"tokens_in":664,"tokens_out":653,"duration_ms":5946,"temperature":1.0,"reasoning_tokens":585,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:55:55.723720+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A definitive check would be an unbiased low-temperature calculation of the same boundaries, for example sign-free quantum Monte Carlo in the unfrustrated transverse sector or a tensor-network or high-order series calculation near $\\theta=0.554\\pi$, looking for whether the upper non-magnetic boundary sits near the PMFRG value or the cluster mean-field value $0.613\\pi$, and whether the QSI$_0$/FM$_\\perp$ boundary remains at $\\theta\\approx-0.033\\pi$. Within PMFRG itself, rerunning the scaling analysis with a fourth system size near $N=4000$ would reveal whether the $\\xi/L$ collapse of the two largest sizes continues or drifts.","supporting_citations":[{"cited_title":"Kato and S","cited_arxiv_id":null,"evidence_quote":"Quantum Monte Carlo locating the QSI$_0$/FM$_\\perp$ transition, the benchmark for the lower phase boundary."},{"cited_title":"Banerjee, S","cited_arxiv_id":null,"evidence_quote":"Quantum Monte Carlo value $\\theta=-0.033\\pi$ used as the reference for the same transition."},{"cited_title":"Benton, L","cited_arxiv_id":null,"evidence_quote":"Cluster mean-field and series study providing the upper boundary $\\theta=0.613\\pi$ and the spin-nematic proposal that the PMFRG result is compared against."},{"cited_title":"Desrochers and Y","cited_arxiv_id":null,"evidence_quote":"Gauge mean-field theory predicting distinct transverse structure-factor maxima for QSI$_0$ and QSI$_\\pi$, used to identify both regimes."},{"cited_title":"Hering, V","cited_arxiv_id":null,"evidence_quote":"PFFRG calculation of lattice-symmetry-breaking response functions in the pyrochlore Heisenberg model, the baseline for the nematic-response analysis."},{"cited_title":"Niggemann, B","cited_arxiv_id":null,"evidence_quote":"Introduces the pseudo-Majorana functional renormalization group at finite temperature and its asymptotic exactness at high $T/|J|$."},{"cited_title":"Schneider, J","cited_arxiv_id":null,"evidence_quote":"Introduces the temperature-flow scheme that lets a single run cover the full cooling range down to small $T$."},{"cited_title":"Niggemann, J","cited_arxiv_id":null,"evidence_quote":"Quantitative functional renormalization for three-dimensional Heisenberg models, source of the correlation-length scaling-collapse transition detection."},{"cited_title":"M¨ uller, A","cited_arxiv_id":null,"evidence_quote":"Quantum Monte Carlo critical temperature for the ferromagnetic Heisenberg model used in the benchmark table."}],"review_version":1}