{"id":"472ce60a-5849-4c57-ac43-33ad7b39bf23","arxiv_id":"2411.16526","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In square ice with domain-wall boundary conditions, the interior spin liquid is radially inhomogeneous and shows coexistence of pinch points and Bragg peaks, making it an unconventional Coulomb phase.","lead":"This paper simulates the ground states of square ice under special boundary conditions and finds that the disordered interior is not a uniform spin liquid: its vertex populations and magnetic correlations vary with distance from the center, with magnetic order growing toward the edge. It reports that algebraic correlation signatures (pinch points) coexist with magnetic Bragg peaks, which would make the system an unconventional Coulomb phase.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Coulomb-phase claim rests on an unquantified 'no broadening' of a pinch point buried under a Bragg peak; a Bragg tail could mimic the same width, so direct real-space scaling of connected correlations is needed.","rationale":"The reader's weakest assumption and my stress-test concern coincide: the apparent unpinned, unbroadened pinch point is the load-bearing evidence for the unconventional Coulomb phase, and it is not sufficient because the singularity itself is unresolved and a Bragg peak occupies the same reciprocal-space location. The paper is careful and honest, explicitly acknowledging in Sec. IV that Henley's sufficient conditions for a Coulomb phase are not fulfilled, so the burden falls entirely on direct evidence of algebraic correlations. The current evidence is a null result: 'no broadening' is consistent both with a true algebraic singularity and with a featureless fluctuating component hidden under a finite-size-broadened Bragg peak. A direct real-space connected-correlation scaling test would settle which interpretation is correct. The paper's numerical characterization of vertex populations, radial magnetization, and inhomogeneous correlations is credible and useful regardless of the outcome, so the reader's CONDITIONAL verdict remains appropriate; no change in verdict is needed, only the proposed check to convert the condition into a definitive statement.","tokens_in":10326,"tokens_out":4401,"duration_ms":64727,"concrete_test":"Recompute, from the stored 1000 configurations, the connected correlation C(r)=<δs_i·δs_j> with δs_i = s_i - <s_i>, averaged over pairs inside each annular region of Fig. 4(a) with separation r along the local radial and azimuthal directions, for N=30, 60, and 101. If C(r) is consistent with A(θ)/r^2 over at least a decade in r and stable with N, the Coulomb-phase claim is supported; if C(r) decays faster than a power law or exhibits a finite correlation length, the apparent pinch-point width was dominated by the Bragg peak and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and Sec. IV assert that algebraic spin-spin correlations coexist with Bragg peaks, i.e., an unconventional Coulomb phase. The only direct evidence is Sec. III.D and Fig. 5: the pinch-point feature at q=(2,2) is 'not broadened' relative to the square ice, although the authors state that 'the pinch point singularity cannot be resolved' because a magnetic Bragg peak sits at the same wavevector. This is a null measurement. The white-contour width in Fig. 5 can be set by the Bragg peak's finite-size and truncation line shape (the authors themselves attribute surrounding ripples to the annular form factor), not by the underlying correlation exponent, so an absence of extra broadening does not establish an algebraic singularity. The concern is amplified by the paper's own Sec. IV admission that Henley's sufficient conditions fail: coarse-grained regions are strongly correlated and the ordered component cannot be separated by a Helmholtz decomposition. A smooth, short-range-correlated fluctuating component with a strong ordered background could produce the same 'unbroadened' appearance. Thus the central claim is under-determined by the data presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the ground-state manifold of the square ice model under domain-wall boundary conditions (DWBC), which produces an arctic-circle phase separation into an interior disordered region and ordered corner regions. Using a loop-flip algorithm, the authors compute local magnetizations, vertex populations, and magnetic structure factors (MSFs) in annular regions inside the arctic circle. They report that all quantities vary radially, that the liquid resembles the conventional square ice near the center but orders near the perimeter, and that diffuse intensity coexists with Bragg peaks. The central claim, stated in the abstract and Sec. IV, is that this coexistence, together with the apparent absence of broadening of a pinch-point feature in the MSF, implies that the arctic square ice is an unconventional Coulomb phase with algebraic spin correlations coexisting with Bragg peaks. The paper explicitly notes that Henley's sufficient conditions for a Coulomb phase are not fulfilled.","tokens_in":10559,"tokens_out":3552,"duration_ms":38156,"significance":"If the central claim were established, the paper would document a distinctive boundary-condition-driven spin liquid in which algebraic fluctuations coexist with an inhomogeneous ordered texture on a charge-neutral vacuum, connecting arctic-circle phenomena to fragmented spin-liquid physics. The numerical protocol is straightforward and parameter-free: ground-state configurations are generated by loop flips with no fitted parameters, and the open-boundary square ice serves as an external benchmark. The qualitative observations of a radial vertex-population variation, a radial average magnetization profile, and the coexistence of Bragg peaks with diffuse intensity in annular MSFs are well supported by the presented maps and cuts. However, the paper's defining claim—that the fluctuating component is algebraic—is not established by the evidence shown; this is the load-bearing weakness that requires additional analysis.","major_comments":[{"comment":"The central claim that algebraic spin correlations coexist with Bragg peaks is supported only by the observation that the pinch-point feature in the annular MSF is \"not broadened\" relative to the square ice, even though the authors state that the pinch-point singularity cannot be resolved because a magnetic Bragg peak sits at the same wavevector. Since the white-contour width in Fig. 5 can be set by the Bragg peak line shape and by the form factor of the annular region (the surrounding ripples are attributed to the form factor in Sec. III.C), an absence of extra broadening is a null measurement and does not by itself establish an algebraic singularity. Please provide a direct real-space diagnostic, such as the connected spin-spin correlator within annular regions, with an explicit test for power-law decay and a quantitative comparison with the square ice.","section":"Section III.D, Fig. 5"},{"comment":"The annular MSFs and intensity cuts are presented without statistical uncertainties or convergence diagnostics: no error bars are shown in Figs. 4(c,d) or Fig. 5, and there is no discussion of how many statistically independent configurations were used after the loop-flip updates, nor of how the results depend on the number of loop updates n and on lattice size N. Because the core pinch-point claim is an absence of broadening, it is essential to show that the measured width is stable under increased sampling and system size; otherwise finite-size or sampling artifacts cannot be excluded. Please report error bars and convergence tests for the quantities used to support the Coulomb-phase claim.","section":"Section III.C, Eq. (2)"},{"comment":"The paper explicitly acknowledges that Henley's sufficient conditions for a Coulomb phase are not fulfilled: coarse-grained regions are strongly correlated, and the ordered component cannot be separated from the fluctuating component by a Helmholtz decomposition. This admission is appropriate, but it increases the burden on the direct evidence for algebraic correlations. As it stands, the MSF analysis does not separate the fluctuating contribution from the Bragg contribution, so the assertion in Sec. III.D and the abstract that the correlations are algebraic is underdetermined. Please either separate the fluctuating component (e.g., by subtracting a fitted Bragg/background model or by computing the connected correlator after removing the local average spin texture) or soften the claim to the coexistence of Bragg peaks with a diffuse, structured component whose algebraic nature remains an open question.","section":"Section IV"},{"comment":"The main text states that \"the pinch point singularity cannot be resolved,\" while the Fig. 6 caption states that pinch points are \"unambiguously revealed.\" These statements need to be reconciled. Please specify quantitatively what feature in the sector MSFs is identified as a pinch point (for instance, the behavior of the intensity minimum at q=(2,2) along the two orthogonal axes) and how its width is extracted from the data.","section":"Section III.D, Fig. 6 caption"}],"minor_comments":[{"comment":"In the second paragraph, \"Typically,n = N^2 loops\" is missing a space after the comma, and \"constraint\" is misspelled as \"contraint\"; these should be corrected.","section":"Section II"},{"comment":"Equation (2) uses N as the total number of spins, but when the MSF is computed in annular or sector regions the normalization becomes ambiguous; please specify how N is defined for each region and whether the normalization affects the relative intensities shown in Figs. 4 and 5.","section":"Section III.C, Eq. (2)"},{"comment":"The text contains typos \"corse-grained\" for \"coarse-grained\" and \"there is not reason a priori\" for \"there is no reason a priori\"; please correct them.","section":"Section IV"},{"comment":"The abstract and main text call the features at the zone center \"magnetic Bragg peaks,\" while the Fig. 3 caption uses \"Bragg-like features\"; please use consistent terminology throughout.","section":"Abstract and Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is an exploratory numerical characterization with no code or data repository and no machine-checked proofs. The main claim is fixable in principle by adding a real-space scaling analysis of connected correlations and proper statistical uncertainties; for this reason I recommend major revision rather than rejection. The fit to the journal's scope is reasonable, but the paper would be stronger if the Coulomb-phase language were calibrated to what the data actually demonstrate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid numerical characterization of the square ice model under domain-wall boundary conditions, and I think that part is worth engaging with. The radial dependence of vertex populations and the gradual onset of order as the arctic curve is approached are documented carefully, with a healthy respect for what the exact solutions already tell us. The MSF analysis of annular regions is a reasonable way to probe spatial variation, and the figures are informative. Credit where it is due: the authors are explicit that Henley's sufficient conditions fail, and they do not pretend the ordered component can be separated by a Helmholtz decomposition. That honesty is a real strength.\n\nThe soft spot is the central claim. The paper says pinch points signaling algebraic correlations coexist with Bragg peaks, but the only direct evidence is that the pinch-point feature is \"not broadened\" relative to the square ice. The issue is that the singularity itself cannot be resolved because a Bragg peak sits on top of it, and the width of the feature could easily be controlled by the Bragg peak's line shape, the finite system size, or the annular form factor. The authors themselves attribute the surrounding ripples to the form factor, so they cannot then lean on an absence of extra broadening to establish an algebraic singularity. A direct real-space scaling of the connected correlation function, fitted in a region away from the ordered component, would settle this. Without that, the claim is under-determined. The paper should either soften the conclusion to \"compatible with\" a Coulomb phase or supply the quantitative test.\n\nAlso, the analysis would be much stronger with error bars or convergence diagnostics, and the lack of code or data is a minor but real weakness. The radial density variation itself is not new, but the MSF angle-resolved and annular results are a useful addition to the literature.\n\nWho is this for? Groups working on artificial spin ice and Coulomb phases in vertex models. I think it deserves a serious referee, not a desk reject, because the observations are credible and the question is interesting. My recommendation to the editor: send it to review, and make sure the referees push on the pinch-point evidence. If the authors add a quantitative correlation scaling and temper the abstract, this could become a solid contribution.","headline":"A careful numerical study of the six-vertex model under DWBC, but the headline Coulomb-phase claim is not actually established by the data presented.","tokens_in":11051,"tokens_out":1380,"would_cite":false,"duration_ms":15480,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B20","82B26"],"pacs":[],"model":"deepseek-v4-flash","headline":"Arctic square ice keeps algebraic spin correlations even as magnetic order sets in.","keywords":["square ice","Coulomb phase","spin liquid","arctic circle","domain wall boundary conditions","pinch points","magnetic structure factor","vertex model"],"falsifier":"Compute the radial correlation function $C(r) = \\langle \\mathbf{S}_i \\cdot \\mathbf{S}_j \\rangle$ inside the arctic disk for a sequence of lattice sizes $N = 30, 60, 101, 200$: if, at fixed separation normalized by the disk radius, $C(r)$ decays exponentially rather than as a power law, the algebraic Coulomb-phase claim fails. Equivalently, measure the full width at half maximum of the diffuse intensity around $\\mathbf{q} = (2,2)$ in the annular structure factor for increasing $N$; a width that grows as $1/N$ would show that the apparent pinch point is only a finite-size effect.","tokens_in":10143,"feed_emoji":"🧊","tokens_out":3893,"duration_ms":37949,"temperature":0.7,"pith_summary":"The paper claims that the square ice model under domain wall boundary conditions hosts an unusual Coulomb phase in its disordered interior: the spin liquid is confined inside an arctic circle, surrounded by magnetically ordered regions, and the liquid's properties vary radially. Near the center the spins behave like the conventional square ice, while approaching the arctic curve magnetic order gradually develops. Remarkably, the magnetic structure factor computed inside the disk still shows pinch points whose width is no broader than in the pure square ice, even though Bragg peaks sit on top of them. The authors conclude that algebraic spin-spin correlations and magnetic order coexist, making the arctic square ice an unconventional Coulomb phase related to fragmented spin liquids but with no magnetic charge injection.","feed_headline":"Arctic square ice is a Coulomb phase with Bragg peaks","feed_subtitle":"Ground-state numerics show sharp pinch points and radial magnetic order coexist inside the arctic disk.","key_machinery":"The central object is the arctic square ice, the spin-liquid state confined inside the arctic circle when the square ice model obeys domain wall boundary conditions. The key diagnostic is the magnetic structure factor computed in annular rings within the disk, whose pinch points are the reciprocal-space signature of algebraic spin correlations. A loop flip algorithm generates a statistical ensemble of $10^{3}$ ice-rule-obeying configurations, and radial profiles of vertex populations and average spin magnitude track how order sets in toward the perimeter.","core_discovery":"In the ground state of square ice with domain wall boundary conditions, the spin liquid occupying the arctic disk is not a homogeneous Coulomb phase. All local quantities—vertex populations, average spin texture, and pairwise correlations—depend on the distance r from the disk center: they match the conventional square ice for r/d < 1/5, change continuously in the intermediate regime, and saturate into ordered type-II tiling beyond the arctic curve. When the magnetic structure factor is computed in annular regions of fixed radius, the diffuse pinch point at reciprocal lattice vectors remains sharp and un-broadened relative to the open-boundary square ice, despite the growing Bragg intensity at the zone center. Since these ordered correlations cannot be separated from the fluctuating component by a Helmholtz decomposition, the authors argue the arctic square ice is a Coulomb phase of a new kind: ordered and disordered at once, within a charge-neutral vacuum.","pith_inferences":["A natural extension is to check whether the same coexistence of sharp pinch points and boundary-induced Bragg peaks appears for other boundary conditions in the six-vertex model, which would make the radial order a generic feature of arctic-type phase separation rather than a square-ice accident.","The result suggests a broader class of 'boundary Coulomb phases' in which long-range order is induced purely by constraints at the sample edge while the bulk retains algebraic correlations; testing this on kagome or other ice models would probe its generality.","One could test the unconventional nature directly by measuring real-space spin correlations in the arctic disk as a function of system size: a correlation length that grows with N would support a true algebraic phase, whereas saturation would indicate only a dressed paramagnet dressed by the boundary texture."],"forward_implications":["If correct, the arctic square ice provides a concrete spin model in which algebraic spin correlations coexist with magnetic Bragg peaks in the same ground state manifold.","The radial dependence of vertex populations means experimental imaging of arctic square ice should expect local ice-rule statistics to drift from the 38/62 square-ice fractions near the center toward pure type-II tiling near the arctic curve.","Because the ordered component cannot be extracted by Helmholtz decomposition, the result distinguishes boundary-driven order from the fragmentation of magnetic moments into charge-crystal and Coulomb channels.","The un-broadened pinch point offers a testable reciprocal-space signature that could be searched for in scattering experiments on artificial square ice arrays with domain wall boundary conditions.","The fine structure of the Bragg peaks, with sub-peaks at the zone center, encodes the quadrant-wise type-II tiling and may serve as a fingerprint for identifying the arctic phase in experiments."],"supporting_citations":[{"why":"Supplies the experimental realization of the arctic square ice in a programmable qubit lattice and the first observation of the arctic circle, motivating the numerical ground-state study.","marker":"[18]"},{"why":"Establishes the square ice as a Coulomb phase with algebraic correlations and magnetic monopole excitations, serving as the reference state for all comparisons.","marker":"[5]"},{"why":"Provides the phenomenological description of the square ice magnetic structure factor used as the open-boundary benchmark for the diffuse and pinch-point intensity.","marker":"[10]"},{"why":"States Henley's sufficient conditions for a Coulomb phase, which the paper invokes to argue that the arctic square ice is unconventional because those conditions are not met.","marker":"[16]"},{"why":"Introduces the concept of magnetic fragmentation, the framework the paper compares with the arctic square ice's coexistence of ordered and fluctuating components.","marker":"[28]"},{"why":"Describes the loop flip algorithm used to shuffle the ground state manifold, which the paper adapts to the domain wall boundary conditions.","marker":"[29]"}],"fun_headline_variants":["Arctic square ice: a Coulomb phase that orders radially","Pinch points and Bragg peaks coexist in arctic square ice","Arctic square ice: ordered and disordered at once","Arctic square ice's spin liquid shows radial order","Unconventional Coulomb phase: arctic square ice revealed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that the arctic square ice is a Coulomb phase rests on the assumption that the sharp, un-broadened feature seen at the pinch-point location inside the annular structure factor is a genuine algebraic singularity and not a finite-size artifact or a leftover of the Bragg peak that sits at the same wavevector.","fun_headline_variants_meta":{"raw":{"variants":["Arctic square ice: a Coulomb phase that orders radially","Pinch points and Bragg peaks coexist in arctic square ice","Arctic square ice: ordered and disordered at once","Arctic square ice's spin liquid shows radial order","Unconventional Coulomb phase: arctic square ice revealed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000674,"raw_usage":{"total_tokens":3044,"prompt_tokens":894,"completion_tokens":2150,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":2080}},"tokens_in":510,"tokens_out":2150,"duration_ms":16568,"temperature":1.0,"reasoning_tokens":2080,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:00:14.439090+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the radial correlation function $C(r) = \\langle \\mathbf{S}_i \\cdot \\mathbf{S}_j \\rangle$ inside the arctic disk for a sequence of lattice sizes $N = 30, 60, 101, 200$: if, at fixed separation normalized by the disk radius, $C(r)$ decays exponentially rather than as a power law, the algebraic Coulomb-phase claim fails. Equivalently, measure the full width at half maximum of the diffuse intensity around $\\mathbf{q} = (2,2)$ in the annular structure factor for increasing $N$; a width that grows as $1/N$ would show that the apparent pinch point is only a finite-size effect.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental realization of the arctic square ice in a programmable qubit lattice and the first observation of the arctic circle, motivating the numerical ground-state study."},{"cited_title":"Perrin, B","cited_arxiv_id":null,"evidence_quote":"Establishes the square ice as a Coulomb phase with algebraic correlations and magnetic monopole excitations, serving as the reference state for all comparisons."},{"cited_title":"Rougemaille and B","cited_arxiv_id":null,"evidence_quote":"Provides the phenomenological description of the square ice magnetic structure factor used as the open-boundary benchmark for the diffuse and pinch-point intensity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"States Henley's sufficient conditions for a Coulomb phase, which the paper invokes to argue that the arctic square ice is unconventional because those conditions are not met."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the concept of magnetic fragmentation, the framework the paper compares with the arctic square ice's coexistence of ordered and fluctuating components."},{"cited_title":"Coraux, N","cited_arxiv_id":null,"evidence_quote":"Describes the loop flip algorithm used to shuffle the ground state manifold, which the paper adapts to the domain wall boundary conditions."}],"review_version":1}