{"id":"f0a76b17-3b0c-4134-b40b-89ff050abe81","arxiv_id":"2504.12964","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Including a confining background field in Dyson-Schwinger equations makes Lee-Yang edge singularities terminate at the Roberge-Weiss line, mu_i = pi T / 3, above T = 0.235 GeV, while leaving the critical exponent unchanged.","lead":"This paper computes QCD phase transitions in a model with a confining gluon background field, and finds that Lee-Yang edge singularities stop at the Roberge-Weiss line at imaginary chemical potential. It matters because it shows where such singularities can still be used to locate the QCD critical endpoint.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The RW-induced LYES cutoff at T_e = 0.235 GeV is computed entirely from a one-loop background-field potential with vacuum-fitted gluon/ghost/vertex inputs; absent a stability check, the termination temperature is a truncation-dependent prediction rather than a robust QCD statement.","rationale":"Read in good faith, the paper does what it claims: it couples a constant A4 background to the quark DSE, solves the quark gap equation with A4 determined from the minimum of the Polyakov-loop potential, and studies Lee-Yang edge singularities in the complex chemical potential plane. The RW symmetry analysis is standard, and the qualitative mechanism that a first-order RW transition line acts as a branch cut stopping the chiral LYES is physically plausible and consistent with the broader literature. The computation is laid out in enough detail to follow, and the claims in Section II.B are honest about the one-loop truncation, although the phrase 'fully coupled, self-consistent' in Section II.A is stronger than what is solved since the gluon, ghost, and vertex inputs are fixed. I do not see an internal inconsistency in the derivation. The load-bearing numerical output is T_e = 0.235 GeV, and it is controlled by the one-loop background-field potential with vacuum-fitted propagators. The reader identified exactly this as the weakest assumption, and I agree. Since there is no quantitative estimate of how the one-loop/fixed-input approximation shifts T_e, the result cannot be treated as a robust quantitative prediction, but there is also no evidence that the qualitative termination mechanism fails. The reader's CONDITIONAL verdict is therefore appropriate, and no verdict change is needed.","tokens_in":16312,"tokens_out":7155,"duration_ms":79227,"concrete_test":"Recompute dV/dA4 including the next-order (two-loop) contribution to the background-field effective potential, keeping all other inputs identical, and re-extract the temperature T_e at which Im mu_LYE = pi T/3. If T_e shifts by more than ~20 MeV, or if the LYES trajectory no longer meets the RW line, the central termination claim is not robust to truncation. A cheaper cross-check is to repeat the same one-loop calculation with the gluon parameters {a,b,c,d,e,f,g} varied within their vacuum-fit uncertainties and track T_e; the quoted 0.235 GeV should be reported with the resulting spread.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II.B evaluates V'(phi) from the one-loop term only (Eq. (12) and Fig. 3), with the dressed gluon propagator (Eq. (10)), ghost propagator (Eq. (16)), and quark-gluon vertex (Eq. (11)) frozen at vacuum fits. The paper states this explicitly: 'we analyze exclusively the one-loop term in Fig. 3'. All subsequent outputs — the RW transition temperature T_c = 0.155 GeV, the A4 minima, and especially the LYES endpoint T_e = 0.235 GeV — are determined by this fixed-input, one-loop effective potential. The central claim is not protected by RW symmetry alone: the C-plus-Z3 reflection symmetry around mu_i = pi T/3 only says the free energy is symmetric there; it does not specify whether, or at what temperature, the chiral LYES branch point meets the RW branch cut. That meeting is a dynamical condition, and it is computed from the one-loop potential. If higher-loop terms or finite-T/mu refitting of the gluon and ghost propagators shift the A4 minima, T_e and the termination behavior can move. No systematic variation, error bar, or independent lattice-constrained check of the RW potential is provided, so the numerical cutoff is currently a truncation-dependent model prediction rather than a robust quantitative QCD statement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends a previous Dyson-Schwinger equation (DSE) study of QCD phase transitions to include a constant gluonic background field A4. The authors solve the quark gap equation and the background-field equation self-consistently, obtaining the Polyakov loop potential and the chiral condensate. They then continue the calculation to complex chemical potential, identify the Roberge-Weiss (RW) transition, and compute the trajectory of Lee-Yang edge singularities (LYES) of the chiral transition. The central new result is that the LYES trajectory terminates at the RW line mu_i = pi*T/3 at a temperature T_e = 0.235 GeV, which the authors attribute to RW symmetry. They also report that the LYES scaling exponent beta*delta is unchanged from their previous work (beta*delta = 1.38) and extract a slope parameter c2 = 5.3 GeV^(1-beta*delta).","tokens_in":16528,"tokens_out":11020,"duration_ms":107231,"significance":"If the computed cutoff is robust, the paper makes a nontrivial prediction: above T_e approximately 0.235 GeV, the chiral LYES are absent below the RW line, which would restrict the use of LYES extrapolations for locating the QCD critical endpoint. The paper also demonstrates a self-consistent coupling between the A4 background field and the chiral condensate within the DSE framework and explicitly exhibits RW symmetry in the Polyakov loop potential. These are useful contributions. However, the main quantitative result rests on a one-loop truncation of the background-field potential with gluon, ghost, and vertex inputs frozen at vacuum fits, so the significance is conditional on the stability of that truncation.","major_comments":[{"comment":"The central quantitative prediction T_e = 0.235 GeV is obtained from the one-loop background-field potential V' in Eq. (12), with the gluon propagator (Eq. (10)), the quark-gluon vertex (Eq. (11)), and the ghost dressing (Eq. (16)) fixed at vacuum fits, as stated in Section II.B: 'we analyze exclusively the one-loop term in Fig. 3.' Because the location of the A4 minima controls both the RW transition temperature and the LYES termination temperature, this truncation is load-bearing: higher-loop contributions, or a finite-T/mu refitting of the propagators and vertex, could shift the minima and hence T_e. The manuscript provides no stability check, such as a parameter variation, a two-loop estimate, or a lattice benchmark of the RW potential. I ask the authors to assess the sensitivity of T_e and of the termination behavior to the truncation and to state the resulting uncertainty on the central claim.","section":"II.B, Eq. (12)"},{"comment":"Section IV.B reports a directly computed CEP at (T, mu_B) = (0.095, 0.735) GeV but later states that the scaling analysis yields an extracted CEP location at (118, 608) MeV. These values differ by roughly 23 MeV in T and 127 MeV in mu_B, and the text does not reconcile them. Since the final paragraph recommends the use of LYES scaling for T <= 160 MeV as a plausible way to determine the CEP, this discrepancy is directly relevant to the paper's conclusions. Please clarify which quantity (118, 608) MeV refers to and, if it is an extrapolated value distinct from the direct CEP, discuss the origin of the difference and its implications for the extrapolation method.","section":"IV.B"},{"comment":"The claim that the LYES trajectory terminates because of RW symmetry is asserted rather than demonstrated. The text states that RW symmetry entails a reflection of the singularity to mu_i < 1/3 pi T, but no explicit Riemann-surface argument or continuity proof is given to show that the LYES branch point must meet the RW branch cut at T_e = 0.235 GeV rather than continuing onto another sheet. Please provide a more formal description of the analytic structure, or at least a numerical demonstration that the singularity cannot be continued past the RW line, to substantiate the causal interpretation that the cutoff is caused by RW symmetry.","section":"IV.B (termination mechanism)"}],"minor_comments":[{"comment":"There are several typographical and grammatical issues, e.g., 'singulairties' in the Introduction, 'the imaginary of phi8' in Section IV.A, and 'has represented a cut-off temperature' in Section IV.B; the text should be carefully proofread.","section":"Throughout"},{"comment":"In Eq. (12), the symbol P is used for the momentum sum/integral without definition; please define it explicitly (e.g., T sum_n integral d^3q/(2 pi)^3).","section":"II.B, Eq. (12)"},{"comment":"The name 'functional-lattice' for the gluon propagator in Eq. (10) is not explained; a brief description or a reference to the origin of this parameterization would help the reader.","section":"II.A, Eq. (10)"},{"comment":"The phrase 'demonstrate the emergence of RW symmetry' is imprecise: RW symmetry is an exact property of the theory, and the paper demonstrates its realization or spontaneous breaking in the computed potential. Please rephrase.","section":"Abstract and III.B"},{"comment":"In the scaling analysis, beta*delta = 1.38 is fixed from previous work rather than extracted from the data. The authors do state 'taking beta*delta = 1.38', but they should explicitly note that the critical exponent is an input and comment on how the conclusions would change if beta*delta were left free in the fit.","section":"IV.B, Eq. (34)"},{"comment":"The RW transition below T_c is called a 'crossover'; for imaginary chemical potential this term is nonstandard, since there is no order parameter and the transition is analytic. Please rephrase to avoid confusion.","section":"IV.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a follow-up of the authors' earlier work [1], and the novelty is the cutoff of the LYES trajectory at the RW line. The main concern is truncation sensitivity of the one-loop background-field potential. I would like the authors to provide quantitative evidence that T_e is not an artifact of this truncation; if they can do so, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing you should know: this paper reports a genuinely new numerical observation—the Lee-Yang edge singularity trajectory of the chiral transition cuts off at the Roberge-Weiss line mu_i = pi T/3, at T_e = 0.235 GeV, and the authors argue this is due to the RW symmetry, not the chiral dynamics alone. That is absent from their earlier no-background-field work and from the cited lattice/DSE literature. The paper also demonstrates RW symmetry emerging in a DSE framework with a gluonic background field solved self-consistently with the quark gap equation.\n\nWhat it does well: the calculation is laid out clearly. The equations, parameter tables, and symmetry analysis are sufficiently detailed to follow. The comparison between the with-A4 and without-A4 cases is explicit, and the scaling fit to Im mu_LYE = i c2 (T - T_CEP)^{beta delta} is transparent. The qualitative mechanism—that the RW transition line provides a branch cut which stops the LYES trajectory—is physically plausible and consistent with what we know about first-order transition lines.\n\nThe soft spots are real but not fatal. The biggest one is the truncation: the background-field potential V' is evaluated at one loop only, with the gluon propagator, ghost propagator, and quark-gluon vertex frozen at vacuum fits. The paper says this explicitly. That means T_e = 0.235 GeV and the RW transition temperature T_c = 0.155 GeV are all predictions of that fixed-input, one-loop potential. There is no stability check, no parameter variation, no estimate of how higher-loop terms or finite-T/mu refitting would move the numbers. So the exact cutoff temperature is a truncation-dependent model prediction, not a robust QCD statement. The qualitative cutoff, though, is likely robust because RW symmetry itself is exact, and the branch cut reasoning does not depend on the details of the potential—only the temperature at which the LYES hits it does.\n\nAlso, the phrase 'fully coupled, self-consistent' overstates what is actually solved: the gluon and vertex are inputs, not solved for. And the explanation of why the LYES cannot cross the RW line is argued in words rather than shown analytically—the statement about singularities on one side mapping to the other side deserves a more explicit treatment. Minor: beta*delta = 1.38 is imported from the authors' own previous paper and the slope c2 is fit in the same framework, so there is some circularity, but it is mild and the RW symmetry check is independent.\n\nWho is this for: specialists in functional QCD (DSE/fRG) and anyone working on locating the QCD critical endpoint via LYES extrapolations. It deserves a serious referee; I would send it to review. A good referee will ask for a truncation-dependence check, but the paper is coherent, honest, and the new observation is worth reporting.","headline":"A serious DSE calculation with a new and plausible observation—LYES terminating at the Roberge-Weiss line—but the reported cutoff temperature is only as trustworthy as the one-loop, vacuum-fitted truncation behind it.","tokens_in":17182,"tokens_out":4178,"would_cite":true,"duration_ms":36688,"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 Lee-Yang edge singularity trajectory of QCD terminates at the Roberge-Weiss line, $\\mu_i = \\frac{1}{3}\\pi T$, at $T_e = 0.235$ GeV, because the Roberge-Weiss symmetry cuts it off.","keywords":["Lee-Yang edge singularities","Roberge-Weiss symmetry","chiral phase transition","deconfinement","gluonic background field","Polyakov loop","Dyson-Schwinger equations","QCD phase diagram"],"falsifier":"Compute the two-loop contribution to $V'(\\bar A_4)$ with the same vacuum-fitted inputs, or refit the gluon propagator and vertex at finite temperature and complex chemical potential, and track the minimum of $\\bar A_4$ near $T_e$: if a chiral Lee-Yang edge singularity with $\\operatorname{Im}\\mu > \\pi T/3$ appears for $T > 0.235$ GeV, the cut-off claim fails. On the lattice, an analytic continuation of the chiral susceptibility that locates a singularity above the Roberge-Weiss line in that temperature range would likewise falsify the termination.","tokens_in":2077,"feed_emoji":"⚛️","tokens_out":2018,"duration_ms":79435,"temperature":0.7,"pith_summary":"This paper claims that QCD's Lee-Yang edge singularity, the complex-chemical-potential singularity used to locate the chiral critical endpoint, is cut off by Roberge-Weiss symmetry before it can serve as an extrapolation tool at high temperature. Working with Dyson-Schwinger equations coupled to a constant gluonic background field $\\bar A_4$ that represents confinement, the authors solve the quark gap equation and the background-field equation self-consistently and find that the Lee-Yang edge singularity trajectory terminates exactly when its imaginary part reaches $\\mu_i = \\frac{1}{3}\\pi T$, at $T_e = 0.235$ GeV. This cut-off, they argue, is caused by the Roberge-Weiss symmetry, not by the chiral dynamics alone. If correct, this means that above $T_e$ there is no Lee-Yang edge singularity left to extrapolate toward the QCD critical endpoint, while below it the critical scaling exponent $\\beta\\delta = 1.38$ is unchanged by the presence of the confining background field.","feed_headline":"Lee-Yang singularity stops at the Roberge-Weiss line at 0.235 GeV","feed_subtitle":"The chiral singularity used to locate QCD's critical endpoint disappears; the critical exponent survives.","key_machinery":"The central object is the constant gluonic background field $\\bar A_4$, diagonalized as $\\bar A_4 = \\frac{1}{g} 2\\pi T(\\phi_3 t_3 + \\phi_8 t_8)$, whose physical value is fixed by the minimum of the Polyakov-loop potential. The machinery is the coupled set of Dyson-Schwinger equations: the quark gap equation for the three color components, with $\\bar A_4$ appearing as a color-dependent shift of the fermionic Matsubara frequencies $\\tilde\\omega_n = \\omega_n + i\\mu + g\\bar A_4$, together with the one-loop background-field equation $V'(\\bar A_4) = 0$ whose propagators and vertex are the same 'functional-lattice' vacuum-fit ans\\\"atze used in the authors' earlier work. The Roberge-Weiss symmetry, realized as invariance under $\\mu_i = \\frac{2\\pi T}{3} n$ combined with center transformations, selects the physical minimum and is what confines the Lee-Yang edge singularity trajectory to $\\mu_i \\le \\frac{1}{3}\\pi T$.","core_discovery":"The central discovery is that a self-consistently determined gluonic background field changes the fate of the Lee-Yang edge singularity trajectory at high temperature while leaving its critical behavior intact. In the coupled scheme, the background field shifts the Matsubara frequencies of the three color components differently, and the one-loop background-field potential, minimized together with the chiral order parameter, produces a first-order Roberge-Weiss transition at imaginary chemical potential $\\mu_i = (1+2n)\\pi T/3$ above $T \\simeq 0.155$ GeV. As temperature rises, the imaginary part of the chiral Lee-Yang edge singularity grows toward this line; at $T_e = 0.235$ GeV it reaches $\\mu_i = \\frac{1}{3}\\pi T$, and no singularity exists beyond it, because the Roberge-Weiss symmetry forces the physical $\\bar A_4$ configuration to be mirrored across that axis and the singularity is absorbed into the branch cut of the first-order Roberge-Weiss transition. Below the cut, the trajectory follows the scaling law $\\operatorname{Im}\\mu_{\\mathrm{LYE}} = i c_2 (T - T_{\\mathrm{CEP}})^{\\beta\\delta}$ with $\\beta\\delta = 1.38$, and the implied endpoint location is consistent with the previous calculation without the background field.","pith_inferences":["Inference: If this termination is generic, any QCD approach that respects Roberge-Weiss periodicity, including lattice simulations and functional renormalization group studies, should also find no chiral Lee-Yang edge singularity beyond $\\mu_i = \\pi T/3$; the Roberge-Weiss branch cut would act as a universal wall for singularity-based extrapolations.","Inference: The same wall may apply to singularities in other conserved-charge directions, meaning estimates of critical points based on singularity locations in the baryon plane are confined to temperatures below $T_e$.","Inference: Because the quantitative value $T_e = 0.235$ GeV inherits the one-loop truncation, the symmetry protects only the existence and location of the cut-off line, not the precise temperature; recomputing with two-loop terms or a finite-temperature refit of the gluon is the natural check."],"forward_implications":["For $T > T_e = 0.235$ GeV, no chiral Lee-Yang edge singularity exists, so Lee-Yang-based extrapolation to the QCD critical endpoint is invalid in that temperature range.","For $T \\lesssim 0.160$ GeV the trajectory still follows the chiral scaling law, so low-temperature extrapolation to the critical endpoint remains plausible; the extracted endpoint is $(T, \\mu_B) = (118, 608)$ MeV.","The critical exponent $\\beta\\delta = 1.38$ is insensitive to the confining background field, indicating that chiral dynamics dominates the Lee-Yang edge singularity critical behavior in the region $T - T_{\\mathrm{CEP}} \\lesssim 0.04$ GeV.","The Roberge-Weiss transition becomes first order above $T \\simeq 0.155$ GeV at $\\mu_i = (1+2n)\\pi T/3$, and its branch cut is what absorbs the Lee-Yang edge singularity at higher temperature.","Including the background field raises the chiral pseudo-critical temperature from $0.155$ GeV to $0.162$ GeV, showing that confinement back-coupling strengthens chiral symmetry breaking."],"supporting_citations":[{"why":"Supplies the previous DSE treatment of Lee-Yang edge singularities without the background field, the baseline trajectory and scaling analysis that this paper extends.","marker":"[1]"},{"why":"Establishes the one-loop background-field equation for the gluon one-point function and demonstrates its physical viability in the present truncation.","marker":"[19]"},{"why":"Defines Roberge-Weiss symmetry and the Roberge-Weiss transition line used to explain the termination of the Lee-Yang edge singularity trajectory.","marker":"[28]"},{"why":"Provides the earlier DSE treatment of the background field and Roberge-Weiss symmetry on which the coupled framework is built.","marker":"[31]"},{"why":"Supplies lattice QCD Lee-Yang edge singularity points from multipoint Pad\\'e that the present trajectory is compared against.","marker":"[33]"},{"why":"Provides the 'functional-lattice' gluon propagator parameterization used in the quark gap equation and Polyakov-loop potential.","marker":"[59]"},{"why":"Supplies the perturbative one-loop Weiss potential term entering the background-field potential.","marker":"[60]"},{"why":"Supplies the [4,4] Pad\\'e lattice Lee-Yang edge singularity extraction used for comparison and scaling fits.","marker":"[66]"}],"fun_headline_variants":["Lee-Yang singularity hits wall at Roberge-Weisse line","Chiral singularity vanishes at 0.235 GeV, scaling survives","RW symmetry truncates Lee-Yang edge singularity trajectory","Background field kills Lee-Yang singularity above 0.235 GeV"],"cache_read_input_tokens":19072,"weakest_assumption_plain":"The calculation's load-bearing premise is that a single-loop approximation for the background-field potential, with the gluon propagator and quark-gluon vertex fixed by vacuum fits, correctly decides where the $\\bar A_4$ condensate sits at finite temperature and complex chemical potential.","fun_headline_variants_meta":{"raw":{"variants":["Lee-Yang singularity hits wall at Roberge-Weisse line","Chiral singularity vanishes at 0.235 GeV, scaling survives","RW symmetry truncates Lee-Yang edge singularity trajectory","Background field kills Lee-Yang singularity above 0.235 GeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1643,"prompt_tokens":1016,"completion_tokens":627,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":556}},"tokens_in":632,"tokens_out":627,"duration_ms":6629,"temperature":1.0,"reasoning_tokens":556,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:18:33.671848+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two-loop contribution to $V'(\\bar A_4)$ with the same vacuum-fitted inputs, or refit the gluon propagator and vertex at finite temperature and complex chemical potential, and track the minimum of $\\bar A_4$ near $T_e$: if a chiral Lee-Yang edge singularity with $\\operatorname{Im}\\mu > \\pi T/3$ appears for $T > 0.235$ GeV, the cut-off claim fails. On the lattice, an analytic continuation of the chiral susceptibility that locates a singularity above the Roberge-Weiss line in that temperature range would likewise falsify the termination.","supporting_citations":[{"cited_title":"Roberge and N","cited_arxiv_id":null,"evidence_quote":"Defines Roberge-Weiss symmetry and the Roberge-Weiss transition line used to explain the termination of the Lee-Yang edge singularity trajectory."}],"review_version":1}