{"id":"ac15a869-7a2e-4309-96fc-2449efbbaba8","arxiv_id":"2501.15352","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"TRG computations show that Lüscher's admissibility condition moves the first-order transition of the 2D U(1) gauge-Higgs model at theta = pi closer to theta = pi under the field-theoretical definition, and locate the critical endpoint at M_c = 2.9974765(14).","lead":"This paper uses tensor renormalization group (TRG) to simulate a 2D lattice gauge-Higgs model with a topological theta term, avoiding both the Monte Carlo sign problem and topological freezing. It shows that Lüscher's admissibility condition places the phase transition at theta = pi even with the field-theoretical definition of the theta term.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The crucial epsilon=0.1 result lacks convergence checks; with K_g=K_h=20 and D_BTRG=160 fixed, the admissible-region width ~epsilon may be undersampled, so the claimed transition near theta=pi could be a numerical artifact.","rationale":"The paper's strongest and most interesting claim is that smaller epsilon moves the field-theoretical transition to theta=pi, showing an advantage of the Lüscher action. The numerical evidence for this claim is Fig. 6 and Fig. 7, but the entire calculation is carried out at fixed algorithmic parameters with no convergence demonstration in this manuscript. The reader's weakest_assumption identified exactly this issue: K_g=K_h=20 and D_BTRG=160 are taken to be converged, with all checks deferred to Ref. [14]. I agree with that identification. My stress-test refines it: the danger is not generic truncation error but specifically the quadrature resolution of the epsilon=0.1 admissible region, whose width is an order of magnitude smaller than the typical spacing of 20 Gauss-Legendre nodes in the link-angle integration. Because the field-theoretical definition (2.4) and the admissibility constraint (1.2) together make the integrand sharply peaked around admissible configurations, the finite K_g approximation is especially delicate. The paper also contains an unexplained sign of the topological susceptibility relative to Eq. (3.3): since the charge density jumps upward at the transition, the susceptibility defined as its derivative should have a positive peak, while the plotted curves appear negative. I do not make this the primary attack because the central claim rests mainly on the charge density and ground-state degeneracy, but it remains an unresolved internal inconsistency that should be corrected. The independent checks—agreement with Ref. [22] for the Wilson action and with Ref. [21] for M_c—are real positive evidence, yet they do not cover the epsilon=0.1 regime that carries the headline conclusion. Therefore a conditional verdict remains appropriate: the claim is plausible and well-motivated but not numerically secured in the present proceedings paper.","tokens_in":571,"tokens_out":5826,"duration_ms":154286,"concrete_test":"Repeat the field-theoretical computation at beta=10, M=4, epsilon=0.1 (the setup of Figs. 6 and 7) with K_g=K_h=32 and 48, and with D_BTRG=200 and 240. Extract the transition point from the discontinuity in the topological charge density or from the peak of d<Q>/dθ for each parameter set. If the transition point shifts by more than about 0.02π, or if the topological charge density values change by more than a few percent, the epsilon=0.1 conclusion is not numerically established. As a secondary check, recompute the critical endpoint M_c=2.9974765(14) with K=32 and D_BTRG=200 to verify that the Section 3 result is stable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the advantage of Lüscher's admissibility condition for the field-theoretical theta term, demonstrated by Figs. 6 and 7: for epsilon=0.1, the transition appears almost exactly at theta=pi. This conclusion rests on the numerical accuracy of the TRG calculation at beta=10, M=4, with all algorithmic parameters fixed to (K_g, K_h, D_BTRG)=(20,20,160) in Section 3. No convergence checks or error estimates are shown in this paper; the footnote only defers to the companion paper Ref. [14]. This is load-bearing because the epsilon=0.1 integrand is extremely narrow: Eq. (1.2) restricts each plaquette to ||1-P||<epsilon, i.e. plaquette phase |φ| ≲ 0.1 rad. The link-angle integration in Eq. (2.5) is replaced by a 20-point Gauss-Legendre quadrature on [-pi,pi], whose node spacing near zero is substantially larger than this admissible-region width. A product quadrature with only 20 points per link may not resolve the sharp admissible shell, and the bond truncation at D_BTRG=160 may introduce further systematic bias. If the quadrature or truncation is not converged at these parameters, the extracted transition point could shift away from pi, directly undermining the claimed advantage. The consistency of the Wilson-action computation with Ref. [22] and the M_c comparison with Ref. [21] provide independent support, but neither validates the untested epsilon=0.1 regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper studies the two-dimensional U(1) gauge-Higgs model with a topological theta term under Lüscher's admissibility condition, using the tensor renormalization group (TRG). The authors construct a tensor network representation via Gauss-Legendre and Gauss-Laguerre quadratures, then contract the network with bond-weighted TRG. They compute the topological charge density, topological susceptibility, ground-state degeneracy, and transfer-matrix scaling dimensions. For the logarithmic theta term, they locate the critical endpoint at M_c=2.9974765(14) and report evidence for 2D Ising universality. For the field-theoretical theta term, they find that Lüscher's action moves the first-order transition point closer to theta=pi, and for epsilon=0.1 the transition lies almost exactly at theta=pi. The paper concludes that TRG simultaneously resolves the complex-action problem and the topological-freezing problem for this model.","tokens_in":8781,"tokens_out":9821,"duration_ms":89696,"significance":"If the results are correct, the paper provides a useful demonstration that TRG can handle admissible lattice gauge theories with theta terms without suffering from the sign problem or topological freezing. The central quantitative claim, that the Lüscher action with small epsilon moves the field-theoretical theta-term transition to theta=pi, is a falsifiable prediction that can be checked by independent methods. Strengths of the paper include the absence of fitting parameters in the path-integral evaluation, the benchmarks against independent dual Monte Carlo results, the high-precision level-spectroscopy determination of the critical endpoint and central charge, and the explicit treatment of the topological-freezing problem. The main weakness is that the numerical convergence for the new epsilon=0.1 regime is not established within the paper, which makes the central claim dependent on an unverified choice of algorithmic parameters.","major_comments":[{"comment":"The central claim that the transition point for the field-theoretical theta term lies almost at theta=pi for epsilon=0.1 rests entirely on BTRG calculations at the fixed algorithmic parameters (K_g, K_h, D_BTRG) = (20, 20, 160). No convergence check or systematic error estimate is presented for this regime. The footnote in Section 3 defers the algorithmic parameter dependence of the critical endpoint to Ref. [14], but that reference concerns the logarithmic definition and does not cover the field-theoretical epsilon=0.1 observables. This matters because the admissible shell has width ||1-P|| < 0.1 while the 20-point Gauss-Legendre grid for the link angles has spacing of order 0.5 rad near zero; the product quadrature may not resolve the narrow shell, and the bond truncation at D_BTRG=160 may introduce additional bias. Please provide at least one convergence sequence in K_g, K_h, and D_BTRG for the epsilon=0.1 topological charge density and susceptibility, or give a quantitative estimate of the resulting systematic uncertainty in the extracted transition point.","section":"Section 3, Figs. 6 and 7"},{"comment":"The topological susceptibility defined in Eq. (3.3) is the variance of the topological charge and is non-negative, but the plotted \"Topological Susceptibility\" curves in Figs. 4-6 are negative in the displayed range (for example, the inset of Fig. 4 reaches about -0.06 and the L=2^7 curve in Fig. 4 is negative around theta/pi=1). The equation and the figures are therefore inconsistent. Please correct the sign convention in one of them and confirm that the comparison with Ref. [22] is made for the same quantity.","section":"Section 3, Eq. (3.3) and Figs. 4-6"},{"comment":"The paper does not discuss how the field-theoretical theta term (2.4) behaves under the admissibility condition. Since Im P(n) = sin(arg P(n)) approximates Im ln P(n) = arg P(n) up to corrections of order (arg P)^3, taking epsilon small makes the field-theoretical action approach the logarithmic definition (2.3). This provides a simple explanation for the near-pi transition at epsilon=0.1. To make the claimed \"advantage of the admissibility condition\" precise, the authors should state this relation and ideally compare the epsilon=0.1 field-theoretical results with a logarithmic-definition run at the same beta and M; otherwise the reader cannot tell whether the observed effect is specific to Lüscher's action or just a convergence of Eq. (2.4) to Eq. (2.3).","section":"Section 2, Eq. (2.4) and Section 4"}],"minor_comments":[{"comment":"The word \"grantees\" should be \"guarantees\".","section":"Section 2, after Eq. (2.3)"},{"comment":"The inset graphs do not have visible axis labels or legends; please add them or describe the inset content fully in the caption.","section":"Section 3, Figs. 4 and 6"},{"comment":"The caption lists two groups of lattice sizes without a single unified legend, making it hard to distinguish the L=2^3 through L=2^6 curves from the L=2^7 through L=2^10 curves; a single legend would improve readability.","section":"Section 3, Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The main unresolved issue is the lack of algorithmic-parameter convergence checks for the epsilon=0.1 results, which are the basis of the paper's central claim. The companion paper Ref. [14] may contain the needed checks, but the present manuscript should either include them or clearly state which of its results are covered there. The sign inconsistency in the susceptibility figures should also be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The piece worth your attention is the epsilon=0.1 comparison in Figs. 6 and 7: it claims that with a sufficiently small admissibility parameter, the field-theoretical theta term places the transition almost exactly at theta=pi. That is new and, if correct, a nice demonstration that Lüscher's action cures a lattice artifact for this definition of the topological term. The rest of the paper—the M_c determination, level spectroscopy, central charge—is a summary of the authors' companion paper (Ref. [14]), and it is consistent with independent Monte Carlo results.\n\nThe paper does several things well. The TRG setup is straightforward, and the Wilson-action results are directly compared with existing dual simulations. The qualitative trend (smaller epsilon moves the transition toward pi) is clearly visible and physically plausible. The M_c value 2.9974765(14) matches the dual Monte Carlo within errors, which gives confidence that the method is under control.\n\nThe soft spots are real, though none are fatal on their own. First, the epsilon=0.1 regime is exactly where the numerics need the most care, and no convergence checks are shown. The admissible shell is roughly |plaquette phase| < 0.1 rad, while a 20-point Gauss-Legendre quadrature on [-pi,pi] has node spacing around 0.3 rad near zero. A product quadrature may not resolve that narrow shell, and the bond truncation D=160 could also bias the result. The paper fixes these parameters and defers to the companion paper, but that companion paper does not cover this field-theoretical definition. So the central claim 'almost at theta=pi' is backed by a single parameter choice. The stress-test note lands on this, and I agree it is the load-bearing issue.\n\nSecond, the topological susceptibility plots show negative values, while Eq. (3.3) defines chi_Q = - (1/V) d^2 ln Z / dtheta^2, which should be positive (it is the variance of the topological charge). Either the plots follow a different sign convention or the definition in the text is missing a minus; either way, it needs fixing. Third, no error bars appear on any of the field-theoretical results, and no code or data are provided. For a proceedings contribution that is not unusual, but it makes independent verification harder.\n\nThe sign issue and missing error bars are minor. The quadrature concern is not; it directly affects whether the epsilon=0.1 result is real. Still, the consistency with Ref. [22] for the Wilson action and the M_c consistency with Ref. [21] suggest the TRG machinery is working. I would send this to peer review: a referee can reasonably ask for a convergence study in K_g and D at epsilon=0.1, and that is a well-posed request rather than a request for a new paper.\n\nBottom line: if you work on TRG or theta terms, read this together with the companion paper. The new result is promising but provisional; do not quote the epsilon=0.1 transition location as settled until the convergence check appears.","headline":"The epsilon=0.1 comparison is the only genuinely new result, and it currently rests on fixed TRG parameters with no convergence check; the rest is a competent summary of the authors' earlier work.","tokens_in":9371,"tokens_out":3294,"would_cite":false,"duration_ms":30343,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The tensor renormalization group resolves both the complex action and topological freezing problems and shows that Lüscher's admissibility condition places the field-theoretical θ-term transition almost exactly at θ = π.","keywords":["tensor renormalization group","lattice gauge theory","admissibility condition","theta term","U(1) gauge-Higgs model","topological freezing","complex action problem","phase transition"],"falsifier":"Recompute the topological charge density and ground-state degeneracy with larger algorithmic parameters (for example, $D_{\\rm BTRG}=240$ or $320$, and $K_g=K_h=30$) and with $\\epsilon$ between $1$ and $0.1$; if the critical endpoint $M_c=2.9974765(14)$ or the location of the first-order transition at $\\epsilon=0.1$ shifts beyond the quoted errors, the central claim fails.","tokens_in":8213,"feed_emoji":"⚛️","tokens_out":9762,"duration_ms":76186,"temperature":0.7,"pith_summary":"Monte Carlo simulations of lattice gauge theories with a topological $\\theta$ term are blocked twice: the term makes the action complex, and Lüscher's admissibility condition makes topology change very slowly. This paper evaluates the path integral with the tensor renormalization group (TRG), which sidesteps both problems because it sums over configurations directly rather than sampling them. For the two-dimensional U(1) gauge-Higgs model, the paper finds that the field-theoretical definition of the $\\theta$ term, whose transition point with the standard Wilson action deviates from $\\theta=\\pi$ at finite coupling, moves essentially to $\\theta=\\pi$ when the admissibility parameter $\\epsilon$ is small. For the logarithmic definition, the transition sits at $\\theta=\\pi$ and the critical endpoint is $M_c = 2.9974765(14)$ with two-dimensional Ising critical behavior. If these results hold, the TRG opens a practical route to topological terms in lattice gauge theories that standard Monte Carlo cannot reach.","feed_headline":"Topological freezing defeated: transition lands at θ = π","feed_subtitle":"TRG solves the complex-action and freezing problems, pinning the U(1) gauge-Higgs transition at θ = π.","key_machinery":"Lüscher's admissibility condition restricts each plaquette $P_{\\mu\\nu}(n)$ to satisfy $\\|1 - P_{\\mu\\nu}(n)\\| < \\epsilon$, and the associated gauge action is infinite outside that set, splitting configuration space into topological sectors. The tensor network machinery is the bond-weighted tensor renormalization group (BTRG), which coarse-grains the network built from Gauss–Legendre quadrature over link angles and Gauss–Laguerre quadrature over Higgs radial variables; the path integral is then a tensor contraction that includes every topological sector in one computation. The field-theoretical $\\theta$ term enters as $S_\\theta = i\\theta/(2\\pi)\\sum_n \\mathrm{Im}\\,P_{12}(n)$, and the machinery's key role is showing how the truncation imposed by $\\epsilon$ controls the lattice artifacts of that term.","core_discovery":"The central discovery is that the TRG, applied to a tensor network built from Gauss quadrature, computes the full path integral of the U(1) gauge-Higgs model with a $\\theta$ term under the admissibility condition, including all topological sectors, so neither the complex action problem nor the topological freezing problem appears. With the logarithmic $\\theta$ term, the model has a first-order transition at $\\theta=\\pi$ for large enough mass, and tensor-network level spectroscopy locates the endpoint at $M_c=2.9974765(14)$, with combined scaling dimension $x_{\\rm cmb}=3/16$ and central charge $c=0.50(7)$, indicating the two-dimensional Ising universality class. With the field-theoretical $\\theta$ term, the paper claims a clear advantage of the admissible Lüscher action over the Wilson action: at $\\beta=10$, $M=4$, the transition point moves from noticeably below $\\theta=\\pi$ for the Wilson action to almost exactly $\\theta=\\pi$ for $\\epsilon=0.1$, and the two-fold ground-state degeneracy signals spontaneous $\\mathbb{Z}_2$ breaking. The paper thereby argues that the admissibility condition restores the continuum topological interpretation for the field-theoretical definition at finite lattice spacing.","pith_inferences":["The paper leaves implicit that the same TRG strategy transfers to four-dimensional gauge theories, where the field-theoretical $\\theta$ term is the standard lattice definition; the admissibility condition may serve as a general device to control its lattice artifacts.","A testable extension would be to map the critical endpoint for the field-theoretical $\\theta$ term as a function of $\\beta$ and $\\epsilon$, which the paper announces will appear elsewhere; if the endpoint extrapolates to a universal value as $\\epsilon\\to 0$, that would strengthen the case that the field-theoretical definition becomes continuum-like.","The observed enhancement of finite-volume effects at smaller $\\epsilon$ suggests a practical trade-off: smaller $\\epsilon$ improves the $\\theta$ periodicity but demands larger volumes, which TRG can afford because its cost grows only logarithmically with volume."],"forward_implications":["All topological sectors are summed in a single TRG contraction, so no fixed-sector simulation or topological-charge update is needed under the admissibility condition.","For the logarithmic $\\theta$ definition, the first-order transition occurs exactly at $\\theta=\\pi$ for $M \\ge M_c$, and the critical endpoint $M_c = 2.9974765(14)$ is compatible with the dual-simulation result $M_c=2.989(2)$.","For the field-theoretical definition, shrinking the admissibility parameter $\\epsilon$ pushes the transition point toward $\\theta=\\pi$ at finite $\\beta$, demonstrating that the admissibility condition suppresses the lattice artifact of this discretization.","The tensor-network level spectroscopy and central-charge measurement ($c=0.50(7)$) establish that the endpoint is in the two-dimensional Ising universality class."],"supporting_citations":[{"why":"Introduces the admissibility condition and the gauge action that defines the model and separates topological sectors.","marker":"[1]"},{"why":"Reports the topological freezing problem in Monte Carlo simulations with the admissibility action, the obstacle the TRG removes.","marker":"[2]"},{"why":"Companion paper giving the tensor network representation and the algorithmic parameter dependence, including convergence of the critical endpoint.","marker":"[14]"},{"why":"Supplies the bond-weighted TRG algorithm used for the tensor contractions.","marker":"[15]"},{"why":"Provides the expected phase structure: first-order transition at $\\theta=\\pi$ with spontaneous $\\mathbb{Z}_2$ breaking for large mass.","marker":"[17]"},{"why":"Ground-state degeneracy criterion used to bound the critical mass $M_c$.","marker":"[18]"},{"why":"Tensor-network-based level spectroscopy method used to locate the critical endpoint.","marker":"[19]"},{"why":"Refines level spectroscopy with the combined scaling dimension used to identify the two-dimensional Ising universality class.","marker":"[20]"},{"why":"Dual simulation result $M_c=2.989(2)$ against which the TRG critical endpoint is compared.","marker":"[21]"},{"why":"Dual simulation with the field-theoretical $\\theta$ term and Wilson action against which the TRG results are quantitatively compared.","marker":"[22]"}],"fun_headline_variants":["Admissible Lüscher action shifts transition exactly to θ=π","TRG defeats complex action and freezing, finds θ=π transition","Lattice U(1) gauge-Higgs: TRG clears all obstacles, pinpoints θ=π","No sign problem, no freezing: TRG lands on θ=π","Admissibility wins: U(1) transition at θ=π"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the finite algorithmic truncations ($K_g=K_h=20$ quadrature points and bond dimension $D_{\\rm BTRG}=160$) are converged, since no systematic convergence check appears in this paper and the companion paper is cited instead.","fun_headline_variants_meta":{"raw":{"variants":["Admissible Lüscher action shifts transition exactly to θ=π","TRG defeats complex action and freezing, finds θ=π transition","Lattice U(1) gauge-Higgs: TRG clears all obstacles, pinpoints θ=π","No sign problem, no freezing: TRG lands on θ=π","Admissibility wins: U(1) transition at θ=π"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001673,"raw_usage":{"total_tokens":6619,"prompt_tokens":909,"completion_tokens":5710,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":5610}},"tokens_in":525,"tokens_out":5710,"duration_ms":35639,"temperature":1.0,"reasoning_tokens":5610,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:22:30.911819+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the topological charge density and ground-state degeneracy with larger algorithmic parameters (for example, $D_{\\rm BTRG}=240$ or $320$, and $K_g=K_h=30$) and with $\\epsilon$ between $1$ and $0.1$; if the critical endpoint $M_c=2.9974765(14)$ or the location of the first-order transition at $\\epsilon=0.1$ shifts beyond the quoted errors, the central claim fails.","supporting_citations":[{"cited_title":"Gu and X.-G","cited_arxiv_id":null,"evidence_quote":"Ground-state degeneracy criterion used to bound the critical mass $M_c$."},{"cited_title":"Ueda and M","cited_arxiv_id":null,"evidence_quote":"Refines level spectroscopy with the combined scaling dimension used to identify the two-dimensional Ising universality class."},{"cited_title":"Dual simulation of the 2-dimensional lattice U(1) gauge-Higgs model with a topological term","cited_arxiv_id":"1508.00681","evidence_quote":"Dual simulation with the field-theoretical $\\theta$ term and Wilson action against which the TRG results are quantitatively compared."}],"review_version":1}