{"id":"d430eaa6-a46c-4c7a-af6b-5ad06e16c49e","arxiv_id":"2506.19637","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In a strictly short-range XY model made of a plane intersected by parallel planes, true long-range order appears along the intersection lines when the parallel planes enter a Berezinskii-Kosterlitz-Thouless critical phase.","lead":"The paper reports simulations of a new lattice geometry: a flat sheet of XY spins crossed by many parallel sheets, all coupled only to nearest neighbors. Above a critical coupling the crossing lines develop true long-range magnetic order, which ordinary 2D models forbid. The result suggests a new way to stabilize low-dimensional order using critical fluctuations, possibly relevant for future superfluid experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The y-line LR order claim rests on fits to G_y = a + b L^{-q}; a critical phase with logarithmic decay G_y ~ (ln L)^{-p} would fit equally well, so the positive intercept is not yet conclusive evidence of true long-range order.","rationale":"The reader's weakest assumption focuses on whether each P plane remains a standard 2D XY model despite the subextensive coupling to the V plane. That assumption is supported by the paper's own checks: G_P and G_z scale as BKT with eta = 1/4 at K_2, so it is not the most fragile point. The most load-bearing concern is different: even granting the P planes are critical, the evidence that the y-lines enter true LR order is a positive intercept in fits to a + b L^{-q}, and that positive intercept can be mimicked by a slowly varying logarithmic decay with a small exponent, which the paper itself already invokes for the x-lines. Since the central claim is LR order, this model-selection ambiguity directly threatens the headline result. The concrete test is feasible with the publicly deposited Monte Carlo data and would distinguish the two scenarios. The verdict should remain CONDITIONAL: the paper should not be accepted as establishing LR order until this alternative is excluded. The reader's CONDITIONAL verdict already requires conditions, but not specifically this one, so my read does not change the verdict category.","tokens_in":47470,"tokens_out":12662,"duration_ms":136487,"concrete_test":"Using the deposited data (Zenodo 10.5281/zenodo.19901475), refit G_y and <M_y^2> at P1-P4 to (i) the power-law-plus-constant form a + b L^{-q}, (ii) the pure logarithmic-decay form A (ln L + C)^{-p}, and (iii) a mixed form a + b (ln L)^{-p}, over identical L ranges, and compare chi-squared per degree of freedom. If the logarithmic forms yield comparable chi-square and an extrapolated thermodynamic-limit value consistent with zero (or with no significant positive intercept after including corrections), then the claim of true LR order along y is not supported. If the power-law-plus-constant form is decisively preferred with a > 0 at the 3-sigma level across all four parameter sets, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that for K >= K_2 the y-lines enter true LR order with G_y -> a > 0, based on fits of G_y and <M_y^2> to Eqs. (5)-(6), i.e., a + b L^{-q}. The load-bearing step is the assertion that a > 0 in the thermodynamic limit. Over the simulated range L = 24-384, ln L changes only from about 3.2 to 6.0, so a critical phase with logarithmic decay G_y = A (ln L)^{-p} with small p (the x-lines have p about 0.04 in Eq. (8)) would also be well fitted by a + b L^{-q} with a positive intercept: the slow logarithmic variation is absorbed into the constant a plus a power-law correction. The paper never fits G_y or <M_y^2> to a logarithmic form, so the positive intercept does not discriminate between LR order and a critical phase with log decay. The SM argument (Section III.A, Fig. 6) that the slope of g_y(r) versus ln r 'does not increase in magnitude' excludes power-law decay but not a constant slope, which is exactly what a pure logarithmic decay would produce. Thus the numerical evidence as presented does not settle whether the thermodynamic limit has a > 0 or a = 0 with log decay. The authors themselves call for a field-theoretical derivation, conceding there is no analytic support.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter studies a quasi-2D XY model consisting of a vertical (V) plane intersected by L parallel (P) planes, with strictly nearest-neighbor couplings. The authors perform Monte Carlo simulations up to L=384 and report a phase diagram in which, for 0<W<J_BKT, the V plane undergoes two BKT-type transitions at K1 and K2=J_BKT, while for K>=K2 the y-lines of the V plane develop true long-range order described by G_y=a+bL^{-q} with a>0 and q≈0.51(2); the x-lines are instead claimed to remain critical, with G_x~[ln(L/l0)]^{-qhat}. The proposed mechanism is that critical fluctuations of the P planes mediate effective ordering interactions along the intersection lines.","tokens_in":47851,"tokens_out":6749,"duration_ms":75303,"significance":"If the central claim holds, this is a striking result: a strictly short-range classical model in a quasi-2D geometry would host true long-range order along a one-dimensional subspace at finite temperature, mediated by BKT criticality of the surrounding planes, in contrast to conventional Mermin-Wagner-type expectations. The manuscript is technically careful in several respects: it uses multiple independent observables (G_y, <M_y^2>, <M_yk^2>, xi_y, R), reports chi-squared values and systematic L_min studies, and makes the data openly available. The phase-transition analysis at K1 and K2 is credible. The main weakness is that the evidence for a positive intercept in the y-line correlations is not yet decisive against a critical phase with logarithmic decay; given the counterintuitive nature of the claim, this discrimination is essential.","major_comments":[{"comment":"The central claim that G_y tends to a>0 in the thermodynamic limit is not fully supported because the fitted form G_y=a+bL^{-q} is not discriminated from a critical logarithmic decay G_y=A[ln(L/l0)]^{-p} over the simulated range L=24-384, where ln L grows only from about 3.2 to 6.0. The x-line fits in Eq. (8) yield exponents qhat as small as 0.026-0.049, so a similarly slow logarithmic decay for the y-lines would be absorbed into an apparent constant plus a power-law correction. The SM Fig. 6 slope analysis excludes a power-law decay of g_y(r), but a pure logarithmic decay corresponds to a constant slope in that plot and is not excluded. The authors should directly fit G_y and <M_y^2> to logarithmic forms, report the resulting chi-squared and stability with L_min, and, if possible, use a discriminating scaling collapse or an effective-exponent extrapolation that separates a>0 from a=0 with logarithmic decay.","section":"Long-range Ordered Phase, Eq. (5); SM III.A, Fig. 6"},{"comment":"The abstract states that the perpendicular direction exhibits quasi-long-range order, but the body (Eq. (8) and End Matter) describes the x-lines as a critical phase with logarithmic decay G_x~[ln(L/l0)]^{-qhat}. In standard usage, quasi-long-range order denotes power-law decay; logarithmic decay is a different critical behavior. This terminology should be reconciled: if the x-direction is logarithmically critical for K>=K2, the abstract and the phase-diagram labels should say so, and 'QLRO' should be reserved for the K1<K<K2 regime.","section":"Abstract and Main Results"}],"minor_comments":[{"comment":"For G_P at K2 with W=0.8, the free-eta fit gives eta=0.2676(8) with L_min=16, which is not close to 1/4 at the quoted precision; the text says the estimates are again close to 1/4. Please comment on this deviation and on whether the fixed-eta=1/4 fit is preferred despite its larger chi-squared.","section":"SM Table XXVI"},{"comment":"In the K2(L) fits, the free fit returns K2=1.07(2) while the fixed fit uses K2=1.11996; the text then treats 1.11996 as the thermodynamic limit. The choice is reasonable, but the manuscript should state more explicitly that the free fit is statistically consistent with the fixed value and that the fixed-value fits are used only to reduce uncertainty.","section":"Successive Phase Transitions, Eq. (3)"},{"comment":"The final estimate q=0.51(2) is quoted without a transparent combination of the results in Table I; a weighted average or an explicit statement of the spread across P1-P4 would make the procedure reproducible.","section":"Final estimate of q"},{"comment":"The phrase 'complete phase diagram' is stronger than the presented data: the K1 line is determined at only five W values, and the K2 line is verified for a limited set of couplings. A phrase such as 'phase diagram for the studied parameter range' would be more accurate.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is potentially significant and the numerical study is substantial, but the central claim of true y-direction long-range order requires a model-selection test against logarithmic decay. The authors should be encouraged to provide the requested log-fit comparison and, if the two forms are statistically indistinguishable, to soften the claim accordingly. No concerns about citation practice or authorship attribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The model is new and worth knowing about: a V plane intersected by L P planes with subextensive coupling, where critical BKT fluctuations in the P planes are claimed to stabilize true LR order along the intersection lines. If correct, that is a real mechanism, distinct from the 3D E-Log surface scenario and from dissipation-induced order. The MC work is extensive: four parameter sets, sizes up to L=384, multiple observables (G_y, <M_y^2>, <M_yk^2>, xi_y, R) all fitting consistently to Goldstone-mode scaling with q around 0.5, and the data is deposited on Zenodo. The K_2 transition is checked against the known J_BKT = 1.11996 with free fits returning close values before fixing, which is honest practice.\n\nThe soft spots are real but mostly moderate. The abstract overstates the x-direction as 'quasi-long-range order' when the body leaves open a logarithmic decay; that should be fixed. The stress-test concern about the y-direction has genuine bite: over L=24-384, ln L changes only from about 3.2 to 6.0, so a critical phase with G_y ~ (ln L)^{-p}, p small, would also be well fit by a + b L^{-q} with a positive intercept. The SM slope argument (g_y(r) vs ln r does not steepen) does not exclude a pure log decay; a constant slope is exactly what log decay produces. So the positive intercept is not a clean discriminator on its own. However, the divergence of xi_y/L as L^0.25 is harder to reconcile with a critical phase with log decay, and the residual magnetization fits give well-determined nonzero limits. The central claim is plausible but not airtight.\n\nAlso, no source code is released, only data. The authors concede the lack of an analytic derivation and call for field-theoretic work, which is appropriate.\n\nWho it is for: people working on low-dimensional order, boundary criticality, and the E-Log/TQF comparison. It deserves a serious referee. I would send it out, asking the authors to explicitly fit log-decay forms for G_y, <M_y^2>, and xi_y/L, and to tone down the abstract's x-direction language.","headline":"A genuinely novel geometry for evading Mermin-Wagner, with thorough Monte Carlo evidence that is strongly suggestive but not quite conclusive against a log-decay critical phase.","tokens_in":48387,"tokens_out":5130,"would_cite":true,"duration_ms":57958,"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":"The authors report that a strictly short-range quasi-2D XY model develops true long-range order along its intersection lines, with the order switching on at the BKT transition of the intersecting planes.","keywords":["quasi-2D XY model","long-range order","Berezinskii-Kosterlitz-Thouless transition","critical fluctuations","Mermin-Wagner theorem","Goldstone mode","Monte Carlo simulation","anisotropic superfluidity"],"falsifier":"At $K = J_{\\rm BKT}$, measure the y-line correlation $G_y$ on lattice sizes well beyond $L = 384$ and test whether the intercept $a$ in $G_y = a + b L^{-q}$ extrapolates to a positive, stable value; in the same runs, check whether the BKT transition of an isolated P plane shifts when $W$ is varied. If the intercept tends to zero or the P-plane transition moves with $W$, the proposed mechanism fails.","tokens_in":47262,"feed_emoji":"🧲","tokens_out":19252,"duration_ms":167392,"temperature":0.7,"pith_summary":"This paper reports a strictly short-range quasi-two-dimensional XY model in which a single vertical plane of spins is crossed perpendicularly by a stack of parallel XY planes, all coupled by nearest-neighbor ferromagnetic interactions. The authors argue that when the parallel planes are in their Berezinskii-Kosterlitz-Thouless (BKT) critical phase, their critical fluctuations mediate an effective long-range coupling along the intersection lines, giving the vertical plane true long-range order along those lines while the perpendicular direction stays quasi-long-range ordered. Large-scale Monte Carlo simulations and finite-size scaling locate the onset of this order at exactly the BKT coupling of the intersecting planes and find a universal Goldstone-mode exponent $q \\approx 0.51(2)$. If correct, this is a classical, finite-temperature mechanism that stabilizes directional superfluid order in a quasi-2D system with only short-range couplings, and it gives a concrete route to test such order in optical-lattice emulators.","feed_headline":"Force a short-range XY model into true long-range order","feed_subtitle":"The mechanism turns critical 2D fluctuations into directional 1D order, with a universal Goldstone exponent.","key_machinery":"The central object is the orthogonal intersection geometry: a vertical V plane crossed by $L$ parallel P planes, with coupling strength $W$ inside the V plane and $K$ on all other nearest-neighbor bonds, so each P plane touches the V plane only along a line of $L$ sites and remains macroscopically equivalent to a standard 2D XY model. The load-bearing mechanism is the BKT critical phase of the P planes—the scale-invariant quasi-long-range ordered phase of the 2D XY model—which at the transition point $K = J_{\\rm BKT}$ has correlation exponent $\\eta = 1/4$ and for larger $K$ has a continuously varying exponent. These critical fluctuations mediate an effective interaction along the y-lines of the V plane. The long-range order is diagnosed through the finite-size scaling forms $G_y = a + b L^{-q}$ and $\\langle M_y^2 \\rangle = a + b L^{-q}$, whose nonzero intercept $a$ in the thermodynamic limit is the signature of true order, together with the companion forms $\\langle M_{yk}^2 \\rangle = L^{-q}(a + b L^{-\\omega})$ and $(\\xi_y/L)^2 = L^q(a + b L^{-\\omega}) + c$; all four quantities give the same universal $q \\approx 0.51(2)$, the Goldstone-mode exponent of the emergent ordered direction.","core_discovery":"For couplings $W$ inside the vertical (V) plane below the 2D XY BKT value $J_{\\rm BKT} \\approx 1.11996$, the model undergoes successive transitions as the coupling $K$ within the intersecting (P) planes is increased. A first BKT transition at $K_1 \\approx 0.75$ (for $W = 0.8$) takes the V plane from disorder into a quasi-long-range ordered phase. A second transition occurs at $K_2 = J_{\\rm BKT}$, inherited from the simultaneous BKT transition of every P plane: the y-lines of the V plane, i.e., the intersection lines, enter a true long-range ordered phase in the thermodynamic limit, with spin correlation $G_y = a + b L^{-q}$ and $a > 0$, while the x-lines remain critical, with $G_x$ decaying as $(\\ln L)^{-\\hat q}$ (the data do not fully exclude a very weak power law). The long-range order is anisotropic and displays Goldstone-mode physics, with $q \\approx 0.51(2)$ independent of $W$ and $K$ over the studied range. The paper interprets this as the critical fluctuations of two-dimensional P planes—which on their own cannot order at finite temperature—mediating the effective interaction that stabilizes one-dimensional long-range order along their intersection with the V plane.","pith_inferences":["If $q$ is exactly $1/2$, the effective interaction mediated along the intersection lines is likely an inverse-square ($1/r^2$) interaction in one dimension, which is marginal for $O(2)$ order; a field-theoretic derivation would presumably show how the BKT critical plane generates exactly this interaction with a universal amplitude.","The mechanism should extend to other continuous symmetries (for example $O(3)$ Heisenberg spins) provided the environment has a critical phase rather than an isolated critical point, since the P-plane low-temperature BKT phase provides the tunable slow decay that drives the order.","A sharper numerical test than the paper gives would be to measure the spin stiffness of the y-lines in the thermodynamic limit: it should be nonzero in the long-range ordered phase and zero along x, directly confirming the anisotropic Goldstone physics.","The same geometry in a quantum setting—a one-dimensional bosonic chain coupled transversely to a critical two-dimensional bath—should show analogous dissipation-free long-range order, connecting this classical mechanism to impurity and comb-lattice problems."],"forward_implications":["The onset of the long-range ordered phase is pinned to the P-plane BKT coupling $K_2 = J_{\\rm BKT}$ for all $0 < W < J_{\\rm BKT}$, so the ordering is a sharply tunable transition rather than a crossover.","In the ordered phase the system is a directional superfluid: phase coherence along the intersection lines is true long-range order, while the perpendicular direction is quasi-long-range; in a cold-atom realization this should appear as size-independent interference contrast along y and decaying contrast along x.","The exponent $q \\approx 0.51(2)$ is independent of $W$ and $K$ across the four parameter sets studied, marking the phase as a universality class of 1D order mediated by a 2D critical environment.","Because the P planes remain ordinary 2D XY models, the critical environment is continuously tunable by $K$, which may allow experimental control of the emergent order by simply changing the lattice depth of the intersecting planes."],"supporting_citations":[{"why":"Supplies the precise BKT coupling $J_{\\rm BKT} = 1.11996$ of the 2D XY model, the value to which $K_2$ is pinned.","marker":"[41]"},{"why":"Provides the Wolff single-cluster Monte Carlo algorithm used for all simulations in the paper.","marker":"[42]"},{"why":"Supplies the BKT finite-size scaling formula $K_n(L) = K_n + a[\\ln(L/l_0)]^{-2}$ used to extrapolate $K_1$ and $K_2$.","marker":"[45]"},{"why":"Provide the multiplicative logarithmic correction $(\\ln L)^{1/8}$ at the BKT transition used to collapse the correlation data.","marker":"[43, 44]"},{"why":"Supplies the finite-size scaling forms $G = a + b L^{-q}$ and the associated Goldstone-mode analysis used to establish the long-range ordered phase.","marker":"[6]"},{"why":"Establish the extraordinary-log phase in the 3D XY surface and interface problem, the contrasting q2D environment against which the present long-range order is measured.","marker":"[15-19]"},{"why":"Provides the standard 2D XY low-temperature exponents used to verify that P-plane correlations match a genuine 2D XY critical phase.","marker":"[48]"}],"fun_headline_variants":["Critical fluctuations unlock true order in short-range XY","XY model's critical phase sparks 1D long-range order","Short-range spins, long-range order: a fluctuation trick","How 2D criticality births 1D order in XY model","New order from critical noise: XY model defies Mermin-Wagner"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument hinges on each intersecting plane remaining a standard two-dimensional XY model—so that the large slow fluctuations that appear at the transition coupling $K = J_{\\rm BKT}$ are unchanged—even though it is joined to the vertical plane along a line of sites; the paper checks this numerically but does not prove it.","fun_headline_variants_meta":{"raw":{"variants":["Critical fluctuations unlock true order in short-range XY","XY model's critical phase sparks 1D long-range order","Short-range spins, long-range order: a fluctuation trick","How 2D criticality births 1D order in XY model","New order from critical noise: XY model defies Mermin-Wagner"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000275,"raw_usage":{"total_tokens":1704,"prompt_tokens":1066,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":682,"completion_tokens_details":{"reasoning_tokens":551}},"tokens_in":682,"tokens_out":638,"duration_ms":6520,"temperature":1.0,"reasoning_tokens":551,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:31:08.891370+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At $K = J_{\\rm BKT}$, measure the y-line correlation $G_y$ on lattice sizes well beyond $L = 384$ and test whether the intercept $a$ in $G_y = a + b L^{-q}$ extrapolates to a positive, stable value; in the same runs, check whether the BKT transition of an isolated P plane shifts when $W$ is varied. If the intercept tends to zero or the P-plane transition moves with $W$, the proposed mechanism fails.","supporting_citations":[{"cited_title":"Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs","cited_arxiv_id":"1210.6116","evidence_quote":"Supplies the precise BKT coupling $J_{\\rm BKT} = 1.11996$ of the 2D XY model, the value to which $K_2$ is pinned."},{"cited_title":"Wolff, Collective monte carlo updating for spin systems, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the Wolff single-cluster Monte Carlo algorithm used for all simulations in the paper."},{"cited_title":"Probability-Changing Cluster Algorithm for Two-Dimensional XY and Clock Models","cited_arxiv_id":"cond-mat/0202161","evidence_quote":"Supplies the BKT finite-size scaling formula $K_n(L) = K_n + a[\\ln(L/l_0)]^{-2}$ used to extrapolate $K_1$ and $K_2$."},{"cited_title":"Bulk and surface properties in the critical phase of the two-dimensional XY model","cited_arxiv_id":"cond-mat/0211584","evidence_quote":"Provides the standard 2D XY low-temperature exponents used to verify that P-plane correlations match a genuine 2D XY critical phase."}],"review_version":1}