{"id":"012bb587-a3c0-4fdb-9eb3-cb0be643905c","arxiv_id":"2608.02720","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Using the 1/N expansion, the authors show that the O(N)×Z2 biconical model in 3≤d<4 has a critical fixed point whose Z2 symmetry remains spontaneously broken at arbitrarily high temperature for large finite N.","lead":"A large-N calculation shows that the biconical scalar model in 2+1 dimensions keeps its Z2 symmetry broken no matter how high the temperature rises. The paper unifies two previously separate constructions of this 'thermal order' effect into one analytic framework.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-temperature NLO corrections are not computed; the d=3 PSSB minimum is fixed at LO with NLO-determined couplings, and the 3<d<4 flat-direction lift lacks error-controlled numerics.","rationale":"The reader's CONDITIONAL verdict already identifies the main structural gaps, and my read does not move that verdict. I agree with the reader that the d=3 NLO thermal corrections are not computed and that the 3<d<4 flat-direction lift relies on a numerical positivity check without error bars. I differ from the reader's emphasis on the field-space-localization assumption; in my reading, restricting to cluster-decomposing states is a standard and legitimate way to define spontaneous symmetry breaking, and the convex-envelope objection does not remove the localized minima. The more load-bearing gap is the order mixing: the thermal minimum is determined at LO in 1/N, while the coupling that creates it is determined by the NLO zero-temperature analysis. Without a finite-T NLO computation, the claim that the minimum survives for arbitrarily high temperatures is an assertion about a regime the paper does not explicitly calculate. A concrete finite-T NLO evaluation would settle whether the O(1/N) corrections are uniformly small and whether the Hessian remains positive. If that computation revealed a T-dependent 1/N correction that can overcome the LO Hessian, the central claim would fail; if it confirms uniform smallness, the conditional concerns would be resolved. I also noticed an apparent inconsistency between equation (2.7) for γφ and Table 1, but since that affects subleading corrections rather than the leading PSSB mechanism, I did not make it the headline concern.","tokens_in":31123,"tokens_out":40515,"duration_ms":405085,"concrete_test":"Compute the finite-temperature NLO effective potential in d=3 at the LO minimum, including all loops generated by the χ6 vertices listed in (3.26), with the T=0 counterterms already fixed by the NLO analysis. Check that the O(1/N) shift in χ̄2 and in the Hessian eigenvalues is bounded by C/N with a T-independent C over T/Λ ranging from 10^-4 to 10^-1. If the Hessian can become negative for some N≥N0 or the shift grows with ln(Λ/T), the Section 5.2 stability argument fails and the PSSB claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central d=3 claim in Section 5.2 is that the LO thermal effective potential has a Z2-breaking minimum for every T>0 because the χ6 term with g=15360 stabilizes the negative-ρ region. However, g=15360 and α=-1 are fixed by the zero-temperature NLO scale-invariance condition, and only the logarithmically divergent part of the T=0 potential is used in that step. The finite-temperature NLO effective potential is never computed: the expectation values in (5.11) and (5.18) and the Hessian in (5.19) are leading order in 1/N, with NLO-fixed couplings substituted in. The paper asserts that O(1/N) corrections cannot make the vev vanish for sufficiently large N, but it provides no uniform-in-T bound on those corrections. In particular, the finite-T NLO contribution can in principle introduce T-dependent terms of the same power T^3, possibly with ln(Λ/T) factors, that shift or destabilize the LO minimum. In 3<d<4, the analogous issue appears in the flat-direction lift, which depends on the numerical positivity check in Appendix F performed at a single cutoff Λ/T=2000 with no error bars. The concern is not that a contradiction is shown, but that the central claim is established at leading order in 1/N using couplings determined at next order, leaving the order that fixes the couplings unchecked at finite temperature.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the O(N)×Z2 biconical scalar model in 3≤d<4 spacetime dimensions using the 1/N expansion. It computes the zero-temperature effective potential at leading and next-to-leading order, imposes scale invariance to fix the critical couplings, derives beta functions and relevant scaling exponents, and checks the d→4 limit against the epsilon expansion. At finite temperature, using the leading-order thermal effective potential with the next-to-leading-order-fixed couplings, the authors find a Z2-breaking minimum for every T>0 in 3≤d<4, and conclude that for large finite N the Z2 symmetry remains spontaneously broken at arbitrarily high temperatures. The paper also claims to establish the equivalence between the critical biconical model and the construction of [13].","tokens_in":31422,"tokens_out":9371,"duration_ms":89612,"significance":"If the central claim holds, the paper provides a unified analytic large-N treatment of persistent spontaneous symmetry breaking in local conformal field theories, reproducing and extending results from the epsilon expansion and the functional RG. The zero-temperature analysis is careful and internally consistent: the NLO effective potential is computed explicitly, the logarithmic divergences are isolated, and the d→4 limit is verified on two loci. The beta functions and scaling exponents agree with FRG numerics and the epsilon expansion. The main weakness is that the finite-temperature PSSB claim is established only at leading order in 1/N for the thermal potential, with no explicit control of NLO thermal corrections; the numerical positivity check in Appendix F also lacks error estimates. These gaps leave the 'large finite N' claim under-supported rather than refuted.","major_comments":[{"comment":"The thermal expectation values (5.11) and (5.18) are derived from the LO thermal effective potential, while the couplings α and g are fixed by the NLO zero-temperature analysis. The NLO finite-temperature effective potential is never computed, so the statement in Section 5.2 that 'O(N^{-1}) corrections ... cannot make it vanish for sufficiently large N' is an assertion without a uniform-in-T estimate. If the NLO correction to the vev contains a factor growing with T, no finite N would suffice at arbitrarily high T. Please either compute the NLO thermal corrections or provide a bound that is uniform in T.","section":"Section 5.2 (and 5.1)"},{"comment":"The lifting of the flat direction at finite T for 3<d<4 rests on the positivity of F(0) and F'(0), established numerically in Figure 1 at a single cutoff Λ/T=2000 with no error bars or convergence study. Since this positivity is the mechanism that selects φ=0 and yields (5.11), please provide a quantitative cutoff-dependence analysis or an analytic argument; the current evidence is numerical only.","section":"Appendix F / Section 5.1"},{"comment":"At T=0 the NLO effective potential is defined only for ρ=φ²+αχ²≥0 (see Section 4 and the footnote on p.9), whereas the finite-temperature d=3 analysis accesses ρ<0. The paper does not discuss how the 1/N expansion behaves in this region, nor whether NLO terms can destabilize the LO minimum found at χ²≈0.0596T. The claim of PSSB for large finite N therefore requires at least an argument that the NLO corrections remain bounded uniformly in T in the ρ<0 region.","section":"Section 5.2 / Section 4"}],"minor_comments":[{"comment":"The restriction to field-space-localized states is a structural assumption; please add a short discussion in Section 7 of whether mixed (non-clustering) states could change the PSSB conclusion.","section":"Footnote p.9"},{"comment":"The claimed agreement with [13] 'up to a factor of 1/2' should be made precise; specify the corresponding formula in [13] and whether the factor is d-dependent.","section":"Eq. (5.11)"},{"comment":"The axes are not labeled in the text; add axis labels and describe the numerical subtraction procedure in more detail.","section":"Figure 1"},{"comment":"The sentence 'setting m=0' in the NLO correction (5.9) is an approximation; state the order of the neglected terms explicitly.","section":"Section 5.1"},{"comment":"The phrase 'large finite N' is not quantified; consider stating the condition, e.g., N sufficiently large so that the leading-order vev dominates, with corrections bounded by a T-independent constant.","section":"Abstract / Section 5"},{"comment":"Small typos in diagrammatic notation, such as 'spines' for spine insertions, should be corrected.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The paper's main new result—the finite-T PSSB claim—rests on the same leading-order thermal analysis as [13] plus an NLO-fixed coupling choice. The authors should be asked to clarify the novelty relative to [13], especially for d=3, and to provide the missing NLO thermal computation. The numerical check in Appendix F, while suggestive, is not rigorous enough for the central claim. Overall, the paper is a solid T=0 large-N analysis with an incomplete finite-T extension; major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. This is the first explicit large-N treatment of the biconical model in 3≤d<4. The NLO effective potential, the critical couplings α=-(d-1)(d-2)/2 and g=15360 in d=3, and the beta functions all check out, and the d→4 limit agrees on two loci. The match with FRG numerics is reasonable, and the identification of the biconical fixed point with the Komargodski–Popov construction is a genuine consolidation. The sharp prediction χ̄²≈0.0596 T in d=3 is worth taking seriously.\n\nThe soft spots are real but not disqualifying. The finite-temperature effective potential is computed only at LO; the NLO corrections are never evaluated. The d=3 minimum uses NLO-fixed couplings in the LO thermal potential, and the claim that O(1/N) corrections cannot make the vev vanish rests on an assertion without a uniform bound. In 3<d<4 the flat-direction lift depends on a numerical positivity check at a single cutoff with no error bars. These are gaps, not contradictions. A referee should ask whether the NLO thermal contributions can shift or destabilize the minimum; the authors should either answer it or explicitly label it as an open problem.\n\nThe paper is honest about its structural assumption: the effective potential is restricted to field-space-localized states, and the footnote on p.9 says so plainly. That premise is standard for this kind of analysis. The self-citations are appropriate given the prior literature.\n\nThis is a paper for people working on thermal order, large-N vector models, and RG flows in 2+1 dimensions. It deserves a serious referee. I would send it out, and I would expect the review to ask for the NLO finite-T calculation or a clear statement of its absence. As is, the central claim is plausible and well-supported at LO, but the key order that fixes the couplings is unchecked at finite temperature.","headline":"The paper earns its place: first explicit large-N treatment of the biconical model, but the d=3 PSSB claim is established only at LO with NLO-fixed couplings and the NLO finite-T corrections are not computed.","tokens_in":31937,"tokens_out":2378,"would_cite":true,"duration_ms":21874,"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 paper claims that in the O(N)×Z2 biconical model with 3 ≤ d < 4, the Z2 symmetry remains spontaneously broken at arbitrarily high temperatures for large finite N, and that this fixed point unifies earlier epsilon-expansion, FRG, and…","keywords":["persistent spontaneous symmetry breaking","biconical model","large-N expansion","effective potential","finite-temperature field theory","Z2 symmetry","renormalization group","conformal field theory"],"falsifier":"A lattice Monte Carlo simulation of the $O(N)\\times\\mathbb{Z}_2$ biconical model in $2+1$ dimensions at the critical couplings, with large $N$, measuring $\\langle\\chi\\rangle$ as a function of $T$: the paper predicts $\\langle\\chi\\rangle^2\\simeq0.0596\\,T$ for all $T>0$ (up to $1/N$ corrections), so observing $\\langle\\chi\\rangle\\to0$ at any finite temperature, or a restoration transition at some $T_c$, would falsify the claim.","tokens_in":30873,"feed_emoji":"🌡️","tokens_out":9418,"duration_ms":78028,"temperature":0.7,"pith_summary":"This paper argues that in the $O(N)\\times\\mathbb{Z}_2$ biconical scalar model in $3\\le d<4$ spacetime dimensions, and in particular at $d=3$, the $\\mathbb{Z}_2$ symmetry stays spontaneously broken at arbitrarily high temperature for large but finite $N$, contrary to the usual expectation that heating restores symmetry. Using the $1/N$ expansion, the authors compute the effective potential at leading and next-to-leading order, identify the infrared fixed point (the critical biconical model), and find stable thermal minima in which the $\\chi$ field has a nonzero expectation value for every $T>0$. They also show that this fixed point is the same as a previously studied large-$N$ conformal construction, providing a common analytic framework for earlier $\\epsilon$-expansion and functional-RG evidence. If correct, the result gives local, unitary, UV-complete relativistic field theories in $2+1$ dimensions whose discrete symmetry is never restored by heating.","feed_headline":"Z2 stays broken at all temperatures in the biconical model","feed_subtitle":"Large-N calculation finds χ's thermal expectation value never vanishes for 3 ≤ d < 4.","key_machinery":"The central object is the effective potential restricted to field-space-localized states, computed by shifting the fields $\\phi\\to N^{1/2}\\phi+\\eta$, $\\chi\\to N^{1/2}\\chi+\\xi$ and resumming the cactus (self-energy insertion) diagrams. Its LO input is the self-energy $m^2=\\Sigma_{\\eta_2}$ satisfying $m^2=\\frac{\\lambda_\\phi}{2}\\sum_K\\frac{1}{K^2+m^2}+\\frac{\\lambda_\\phi}{2}(\\phi^2+\\alpha\\chi^2)$, whose IR solution $m^2=\\frac{2d}{d-2}\\nu_d(\\phi^2+\\alpha\\chi^2)^{2/(d-2)}$ carries the scale-invariance condition $\\lambda_{\\phi\\chi}^2=\\lambda_\\phi\\lambda_\\chi$. At NLO, scale invariance of the potential fixes $\\alpha$ to an isolated negative value and, in $d=3$, requires a dynamically generated $\\chi^6$ term; the NLO thermal corrections then lift the LO flat direction and select $\\bar\\phi=0$, $\\bar\\chi\\neq0$. The same machinery yields the $\\beta$ functions and the three relevant scaling exponents of the fixed point.","core_discovery":"At the critical biconical fixed point of the $O(N)\\times\\mathbb{Z}_2$ model, for $3\\le d<4$ and large finite $N$, the thermal effective potential has two degenerate minima related by $\\mathbb{Z}_2$, with $\\bar\\phi=0$ and $\\bar\\chi\\neq0$ at every temperature. In $3<d<4$ the expectation value is $\\bar\\chi^2=\\frac{\\zeta(d-2)\\Gamma((d-2)/2)}{\\pi^{d/2}(d-2)(d-1)}T^{d-2}+O(N^{-1})$; at $d=3$ it is $\\bar\\chi^2\\simeq 0.0596\\,T+O(N^{-1})$, and the $O(N^{-1})$ corrections cannot cancel it for sufficiently large $N$. The same analysis fixes the critical couplings: scale invariance at NLO selects $\\alpha=\\lambda_{\\phi\\chi}/\\lambda_\\phi=-(d-1)(d-2)/2+O(N^{-1})$, and in $d=3$ the RG flow generates a $\\chi^6$ interaction with critical coupling $g=15360+O(N^{-1})$. The paper further establishes that this critical biconical model coincides with the earlier large-$N$ construction of a temperature-resistant $O(N)\\times\\mathbb{Z}_2$ CFT, thereby unifying the $\\epsilon$-expansion, FRG, and large-$N$ approaches to persistent spontaneous symmetry breaking.","pith_inferences":["If the equivalence with the earlier large-$N$ construction is exact, then new observables computed in one formulation—for instance the stress-tensor two-point function or OPE coefficients—should agree in the other; checking one at NLO would test whether the identification is genuine universality or an accidental match at the computed orders.","The mechanism suggests that other multicritical models with two competing order parameters, such as the $O(N)\\times O(M)$ biconical model, may exhibit persistent symmetry breaking through the same entropic stabilization; a large-$N$ analysis of that model would be a direct test.","Because the $d=3$ fixed point requires an emergent $\\chi^6$ interaction, a lattice model in the same universality class should display effective sextic couplings in its infrared action; measuring critical exponents on the lattice and comparing with the large-$N$ predictions would probe the mechanism beyond perturbation theory."],"forward_implications":["For large finite $N$, the $\\mathbb{Z}_2$ symmetry of the critical biconical model is spontaneously broken at every temperature, with $\\langle\\chi\\rangle_T = a_\\chi T^{\\Delta_\\chi}$ and $a_\\chi\\neq0$.","The critical biconical model and the earlier large-$N$ construction are the same fixed point, so results such as anomalous dimensions and scaling exponents transfer between the two formulations.","At $d=3$, the $\\chi^6$ interaction is an unavoidable part of the IR theory: it is generated by the RG flow, and scale invariance fixes its coupling to $g=15360+O(N^{-1})$.","Reaching the critical theory requires tuning a codimension-three surface: the relevant deformations are the mass-like couplings $v$, $r$ and the quartic combination $\\kappa$, while $\\alpha$ and $g$ are weakly irrelevant.","The resulting anomalous dimensions and scaling exponents agree with numerical FRG results and reduce to the $\\epsilon$-expansion near $d=4$, providing a cross-check on the fixed-point identification."],"supporting_citations":[{"why":"Introduced the biconical model as a candidate for persistent symmetry breaking in 2+1 dimensions using the epsilon expansion; the paper's d→4 limit reproduces and extends its beta functions and fixed point.","marker":"[4]"},{"why":"Established symmetry breaking at all temperatures in the same class of conformal theories; the paper recovers its leading d=4−ϵ result for the thermal expectation value.","marker":"[5]"},{"why":"Provided numerical FRG evidence for the critical biconical fixed point in 2+1 dimensions; the paper compares its scaling exponents and anomalous dimensions against these numerics.","marker":"[12]"},{"why":"Constructed a closely related O(N)×Z2 large-N conformal theory with temperature-resistant order; the paper's central result is that this construction and the critical biconical model are the same fixed point.","marker":"[13]"},{"why":"Introduced the biconical model in the context of multicritical phenomena and supplies the model's original definition.","marker":"[27]"},{"why":"Gave the known large-N anomalous dimension of the critical O(N) model used to check the φ-field anomalous dimension.","marker":"[29]"}],"fun_headline_variants":["Z2 breaking survives infinite heat in biconical model","Biconical model: Z2 order defies all temperatures","Thermal order never melts in biconical model","Persistent Z2 breaking at any temperature, 2+1D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the true thermal vacuum is a field-space-localized state; if the physical vacuum is instead a mixture of macroscopically separated configurations, the computed minima—and the persistent symmetry breaking—could disappear.","fun_headline_variants_meta":{"raw":{"variants":["Z2 breaking survives infinite heat in biconical model","Biconical model: Z2 order defies all temperatures","Thermal order never melts in biconical model","Persistent Z2 breaking at any temperature, 2+1D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000264,"raw_usage":{"total_tokens":1658,"prompt_tokens":1052,"completion_tokens":606,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":536}},"tokens_in":668,"tokens_out":606,"duration_ms":5958,"temperature":1.0,"reasoning_tokens":536,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:00:33.079702+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice Monte Carlo simulation of the $O(N)\\times\\mathbb{Z}_2$ biconical model in $2+1$ dimensions at the critical couplings, with large $N$, measuring $\\langle\\chi\\rangle$ as a function of $T$: the paper predicts $\\langle\\chi\\rangle^2\\simeq0.0596\\,T$ for all $T>0$ (up to $1/N$ corrections), so observing $\\langle\\chi\\rangle\\to0$ at any finite temperature, or a restoration transition at some $T_c$, would falsify the claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Established symmetry breaking at all temperatures in the same class of conformal theories; the paper recovers its leading d=4−ϵ result for the thermal expectation value."},{"cited_title":"Hawashin, J","cited_arxiv_id":null,"evidence_quote":"Provided numerical FRG evidence for the critical biconical fixed point in 2+1 dimensions; the paper compares its scaling exponents and anomalous dimensions against these numerics."},{"cited_title":"Nelson, J.M","cited_arxiv_id":null,"evidence_quote":"Introduced the biconical model in the context of multicritical phenomena and supplies the model's original definition."}],"review_version":2}