{"id":"ac771fd9-7ea0-4502-b688-e01f59d7a05c","arxiv_id":"2505.04569","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the Quark-Meson model, analytic continuation from imaginary to real chemical potential reproduces the chiral phase boundary only up to mu_conv about 146 MeV and has about 150 percent error near the critical endpoint.","lead":"This paper tests a common lattice-QCD trick, reconstructing the phase boundary at real chemical potential from simulations at imaginary chemical potential, inside a simplified Quark-Meson model. The trick works up to about 146 MeV but badly misplaces the critical endpoint, so endpoint estimates from this method deserve caution.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model lacks Roberge-Weiss periodicity and the fit uses imaginary-μ data beyond the RW endpoint, so μ_conv≈146 MeV and the 150% CEP discrepancy are not a faithful benchmark for lattice QCD.","rationale":"The paper is a legitimate, internally consistent model study; the FRG-LPA implementation is standard and the qualitative conclusion that continuation from imaginary μ deteriorates near the CEP is robust. The reader's conditional verdict is appropriate. The most load-bearing weakness is not the arbitrary 10% threshold defining μ_conv, which only affects the interpretation of an 'effective' radius, but the model's missing Roberge-Weiss structure. The authors explicitly acknowledge this gap in the Conclusions. Because Eq. (13) has period 2πT in μ rather than 2πT/3, the phase boundary is smooth across the point where QCD has a genuine phase transition, and the fit of Eq. (24) includes imaginary-μ data beyond the RW endpoint—data that lattice QCD cannot supply. This can skew the extracted κ2, κ4 and the headline μ_conv value. The proposed refit within the RW domain would directly test whether the 146 MeV value survives. This concern is specific, falsifiable, and tied to the paper's stated goal of benchmarking the lattice technique. The verdict remains CONDITIONAL: the authors should either restrict the fit to the RW-limited region or explicitly state that the benchmark applies only to a model without RW structure, thereby limiting the transferability of the quantitative results.","tokens_in":12701,"tokens_out":11497,"duration_ms":112348,"concrete_test":"Refit Eq. (24) to the imaginary-μ data restricted to the Roberge-Weiss domain, |μ_I| ≤ π T / 3 (i.e., μ² ≥ −0.035 GeV² at T ≈ 0.18 GeV), and recompute μ_conv from ε_rel using the same 10% threshold. Report the new κ2, κ4, and μ_conv; if μ_conv shifts by more than ~20% or the fit quality degrades substantially, the original value depends on data outside the QCD-accessible region. In parallel, verify the model's periodicity by computing Tc(μ_I) over a full period in imaginary μ: QCD requires period 2πT/3, whereas the present QM model will show period 2πT if no RW transition is present.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative benchmark for lattice-QCD continuation is compromised by the model's incorrect analytic structure in imaginary chemical potential. The authors admit in the Conclusions that a Polyakov-loop extended version is needed for 'a more precise reconstruction of the imaginary chemical potential dynamic and the RW transition.' Without the Polyakov loop, the quark contribution in Eq. (13) is periodic in μ with period 2πT, whereas in QCD the physical period is 2πT/3 and the Roberge-Weiss transition lies at μ_I/T = π/3. At the pseudocritical temperatures shown (Tc ≈ 176–188 MeV), the RW endpoint is at μ_I ≈ 185–197 MeV, i.e., μ_I² ≈ 0.034–0.039 GeV². The imaginary-μ data in Fig. 3 extend to μ² = −0.1 GeV² (μ_I ≈ 316 MeV), well beyond this point. In the lattice procedure, analytic continuation is restricted to below the RW endpoint; data beyond it are not physically meaningful. Fitting Eq. (24) to model data in this region can bias κ2, κ4 (Table I) and hence the computed μ_conv in Eqs. (26)–(27). The claimed μ_conv ≈ 146 MeV and the ≈150% discrepancy near the CEP therefore do not faithfully represent the lattice reconstruction technique, whose accessible imaginary-μ range is smaller.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript tests the imaginary-μ analytic continuation technique used in lattice QCD by applying it inside the two-flavor Quark-Meson model, computed both in the mean-field approximation and in the FRG-LPA. The authors fit the pseudocritical line with Tc(μ)/Tc = 1 − κ2(μ/Tc)^2 − κ4(μ/Tc)^4 using model data at imaginary μ, continue the fit to real μ, and compare it with the model's direct real-μ pseudocritical line. They define a relative error ε_rel, choose a 10% threshold, and obtain an effective convergence radius μ_conv ≈ 146 MeV in both approximations; near the critical endpoint the discrepancy reaches about 150%. They conclude that the reconstruction works at moderate μ but that CEP locations obtained by continuation from imaginary μ should be regarded with caution.","tokens_in":13028,"tokens_out":8407,"duration_ms":79636,"significance":"The paper provides an internally consistent, non-circular test: the real-μ target line is not used to fit κ2 and κ4, so the comparison measures the quality of the continuation procedure itself. The qualitative conclusion that continuation degrades near the CEP is plausible and well illustrated in both approximations. Its main strength is the direct side-by-side comparison made possible by working in a model where both real and imaginary μ are computable. The quantitative values, however, are not directly transferable to QCD: the model's imaginary-μ analytic structure differs from QCD's (no Roberge-Weiss transition at μ_I/T = π/3), and the extracted μ_conv depends on an arbitrary error threshold. With these caveats addressed, the work would be a useful cautionary benchmark.","major_comments":[{"comment":"The imaginary-μ data used in the fit extend to μ² = −0.1 GeV², i.e. μ_I/T ≈ 1.7–1.8 at the pseudocritical temperatures, well beyond the physical Roberge-Weiss endpoint μ_I/T = π/3. In QCD, analytic continuation from imaginary μ is reliable only below that endpoint; the QM model without a Polyakov loop is 2πT-periodic in μ (see the fermionic source term in Eq. (13)) and therefore has no RW non-analyticity at π/3. Fitting Eq. (24) through this region can bias κ2 and κ4 in Table I and hence μ_conv in Eqs. (26)–(27). Please restrict the fit to the QCD-compatible imaginary-μ range or extend the model with a Polyakov loop, and quantify the resulting change in μ_conv.","section":"Section IV, Fig. 3 and Eq. (13)"},{"comment":"μ_conv is defined by an arbitrarily chosen ε_rel threshold of 0.10, and no uncertainty or sensitivity to that choice is reported. The values in Eqs. (26) and (27) are quoted to four significant digits without error bars even though the fit parameters in Table I have uncertainties. Please propagate the parameter uncertainties and show how μ_conv varies when the threshold is changed (for example 0.05 and 0.15).","section":"Section IV, Eqs. (25)–(27)"},{"comment":"The statement that the relative error is 'of order ε_rel ≈ 1.5 in the proximity of the CEP' is not fully supported in the FRG case, where the CEP is only bracketed because of backbending. Please explain how ε_rel near the CEP is computed in the FRG case and how it varies over the bracketed CEP region.","section":"Section IV, CEP discussion"}],"minor_comments":[{"comment":"Fix the typos 'pbtained', 'extrapolatin', and 'then' (which should be 'than').","section":"Abstract"},{"comment":"The FRG κ4 entry appears as '0255 ± 0.0005' without a leading decimal; clarify whether the intended value is 0.0255 or 0.255.","section":"Table I"},{"comment":"The sentence says 'the presence of the factor of 3 in Eq. (5)', but the factor of 3 appears in the flow equation Eq. (13); the cross-reference should be corrected.","section":"After Eq. (22)"},{"comment":"The dashed 'RW limit' line is not defined in the text; specify whether it is placed at μ_I/T = π/3 (the QCD value) or at the 2πT-periodicity endpoint of the model, and explain its relevance for a model without a Roberge-Weiss transition.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about the need for a Polyakov-loop extended version, but the RW-related caveat is stronger than the conclusions suggest: it affects the very data range used for the fit, not just the quantitative transfer to QCD. I would recommend major revision rather than rejection, since the issue can be addressed by restricting the fit range or adding the Polyakov loop, and by adding threshold sensitivity."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is an honest, internally consistent model benchmark: in the QM model, both mean-field and FRG-LPA, they compare the real-μ phase boundary with the one reconstructed from a standard κ2,κ4 fit to imaginary-μ data. They find a 10% relative-error radius μ_conv ≈ 146 MeV and a relative error of ~150% near the model's CEP. The FRG implementation is nontrivial, and the comparison is set up without circularity. Second, the quantitative benchmark is compromised by the model's analytic structure in imaginary μ. The quark term in Eq. (13) is 2πT periodic in μ, whereas QCD has period 2πT/3 with an RW endpoint at μ_I/T = π/3. The fit in Fig. 3 uses data to μ² = -0.1 GeV², which at Tc ≈ 0.18 GeV is μ_I/T ≈ 1.8—far past the physical RW limit. A lattice continuation would stop at μ_I/T ≈ 1.05. So the quoted μ_conv and 150% discrepancy are properties of the QM model's own imaginary axis, not a faithful benchmark for the lattice procedure.\n\nOn the credit side, the authors are open about the two main soft spots: they fix the 10% threshold by hand and acknowledge in the Conclusions that a Polyakov-loop version is needed for a proper treatment of the RW transition. They also flag the FRG backbending that prevents a precise CEP location in that scheme. The qualitative conclusion—analytic continuation from imaginary μ becomes unreliable as the CEP is approached—is almost certainly correct and is consistent with the DSE work [28]; it is not new, but this paper gives a different model-based illustration. The citation pattern looks fine.\n\nThe soft spots beyond the RW issue are the arbitrary threshold defining μ_conv (no sensitivity to it), the absence of code or data to check the FRG numerics, and a likely typo in Table I for the FRG κ4 value. None of these by themselves sink the paper. The RW/periodicity issue is the one that needs to be fixed or at least carefully caveated before the numbers can be used as a lattice benchmark.\n\nBottom line: this deserves a serious referee. I would send it to review but with the expectation that the authors either restrict the imaginary-μ fit to the physical RW domain or switch to a Polyakov-loop extended model, and that they report how μ_conv moves with the threshold. For a reader, treat the qualitative warning as solid and the quantitative figures as model-dependent.","headline":"A clean model calculation whose headline numbers are built on imaginary-μ data in a region that lattice QCD cannot access, so the quantitative benchmark is not faithful to the lattice technique.","tokens_in":13508,"tokens_out":6827,"would_cite":false,"duration_ms":63687,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Aw"],"model":"deepseek-v4-flash","headline":"Continuing the QCD phase boundary from imaginary chemical potential works only up to about 146 MeV and misses the critical endpoint by roughly 150 percent.","keywords":["QCD phase diagram","critical endpoint","imaginary chemical potential","analytic continuation","Quark-Meson model","Functional Renormalization Group","chiral phase transition","convergence radius"],"falsifier":"Run the same reconstruction comparison in a second low-energy effective model whose phase boundary includes deconfinement effects, or replace the polynomial fit with a ratio-of-polynomials resummation: if the effective convergence radius moves far from 146 MeV or the near-endpoint discrepancy moves far from 150%, the numbers are model- or ansatz-dependent; if they persist, the cautionary conclusion is generic.","tokens_in":12509,"feed_emoji":"⚛️","tokens_out":10489,"duration_ms":92992,"temperature":0.7,"pith_summary":"This paper tests a technique used in lattice QCD to reconstruct the chiral phase boundary at real chemical potential by analytic continuation from imaginary chemical potential—a trick that avoids the sign problem. The authors run the Quark-Meson model, a low-energy effective model of chiral symmetry breaking, in both mean-field and functional renormalization group treatments, so they can compute the true phase boundary at real chemical potential and compare it directly with the extrapolated one. They find the reconstruction is reliable only up to an effective convergence radius $\\mu_{\\rm conv} \\approx 146$ MeV; beyond that the extrapolated line diverges from the true line, and near the model's critical endpoint the discrepancy reaches about 150%. The conclusion is that critical-endpoint locations obtained from imaginary-$\\mu$ continuation should be regarded with caution, even though the technique works well for moderate chemical potentials.","feed_headline":"Imaginary-mu trick for QCD boundary fails near critical endpoint","feed_subtitle":"Model test: extrapolating from imaginary to real chemical potential errs by about 150% near the endpoint.","key_machinery":"The object that carries the argument is the reconstruction ansatz for the pseudocritical temperature, $$T_c(\\mu)/T_c = 1 - \\kappa_2(\\mu/T_c)^2 - \\kappa_4(\\mu/T_c)^4,$$ whose coefficients are fitted to the model's imaginary-$\\mu$ data and then continued to real $\\mu$. The comparison metric is the relative error $\\varepsilon_{\\rm rel}=|T_c-T_c^{\\rm(fit)}|/T_c$, and the paper defines the effective convergence radius $\\mu_{\\rm conv}$ as the point where $\\varepsilon_{\\rm rel}$ crosses a 10% threshold. On the model side, the Quark-Meson Lagrangian supplies the chiral phase boundary both at real and imaginary $\\mu$, in mean-field and in FRG with the local-potential approximation.","core_discovery":"The paper's claim is that the standard reconstruction of the chiral phase boundary by continuation from imaginary chemical potential has a finite radius of validity inside the Quark-Meson model: $\\mu_{\\rm conv}\\approx 146$ MeV, essentially the same in the mean-field and FRG-LPA computations. The model's actual critical endpoint lies at more than twice that chemical potential, and at that point the reconstructed boundary disagrees with the true one by a relative error $\\varepsilon_{\\rm rel}\\approx 1.5$, about 150%. The authors conclude that, while the reconstruction is dependable at small and moderate chemical potential, the location of the critical endpoint obtained by this continuation should be treated with caution.","pith_inferences":["An implication the authors leave implicit: because the same 146 MeV radius appears in both mean-field and FRG-LPA, the breakdown is likely set by the polynomial ansatz rather than by the inclusion of fluctuations; a ratio-of-polynomials or resummed fit should move the radius, making that interpretation testable.","If the pattern is generic, then critical-endpoint coordinates quoted from imaginary-$\\mu$ continuation are biased toward small $\\mu$; the size of that bias in any given model could be estimated from the ratio $\\mu_{\\rm conv}/\\mu_{\\rm CEP}$.","A cleaner test would compare reconstruction and direct computation in a regime without the FRG backbending that blurs the endpoint determination, for example with a higher-order derivative expansion; the present paper leaves that as open work."],"forward_implications":["Up to $\\mu \\simeq 146$ MeV, the reconstructed crossover line agrees with the true line within 10% relative error in $T_c$; the reconstruction is therefore a useful tool in the moderate-density regime.","Beyond that radius the extrapolated boundary deviates from the actual one, and the deviation grows to about 150% near the critical endpoint.","Because the same failure appears in both mean-field and FRG-LPA calculations, it is not an artifact of the specific truncation studied.","The coefficients $\\kappa_2,\\kappa_4$ from imaginary-$\\mu$ fits cannot encode the non-analyticity of the critical region; hence any critical-endpoint location extracted from them should be treated as indicative only.","The technique remains a valid benchmark for the crossover region, but improved or hybrid methods are needed for critical-point searches."],"supporting_citations":[{"why":"Defines the analytic-continuation-from-imaginary-mu reconstruction whose reliability the paper tests.","marker":"[20]"},{"why":"Supplies the crossover fit form $T_c/T_c=1-\\kappa_2(\\mu/T_c)^2-\\kappa_4(\\mu/T_c)^4$ used as the reconstruction ansatz.","marker":"[67]"},{"why":"Closest earlier comparison of reconstructed and actual phase boundary at real mu in functional methods; serves as the paper's benchmark.","marker":"[28]"},{"why":"Provides the Quark-Meson model Lagrangian used for the model-side test.","marker":"[25]"},{"why":"Gives the exact flow equation for the effective average action at the core of the FRG-LPA computation.","marker":"[37]"},{"why":"Supplies the optimized regulator used in the bosonic and fermionic flow equations.","marker":"[56]"},{"why":"Earlier effective-model study of imaginary chemical potential whose qualitative behavior the paper confirms.","marker":"[66]"},{"why":"Provides the advection-diffusion formulation of the FRG flow used to integrate the effective potential.","marker":"[5]"}],"fun_headline_variants":["Imaginary-mu reconstruction fails by 150% at QCD endpoint","QCD boundary reconstruction errs by 150% near critical point","Critical endpoint via imaginary mu: off by 150% in model test","Model test: imaginary-mu trick gives 150% error at critical point","Imaginary-to-real mu trick overstated: 150% error at endpoint"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative results transfer to real QCD only if the Quark-Meson model reproduces how QCD's phase boundary actually behaves as chemical potential grows, especially the distance from zero chemical potential to the critical endpoint.","fun_headline_variants_meta":{"raw":{"variants":["Imaginary-mu reconstruction fails by 150% at QCD endpoint","QCD boundary reconstruction errs by 150% near critical point","Critical endpoint via imaginary mu: off by 150% in model test","Model test: imaginary-mu trick gives 150% error at critical point","Imaginary-to-real mu trick overstated: 150% error at endpoint"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000887,"raw_usage":{"total_tokens":3822,"prompt_tokens":929,"completion_tokens":2893,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":2795}},"tokens_in":545,"tokens_out":2893,"duration_ms":17544,"temperature":1.0,"reasoning_tokens":2795,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:25:24.904924+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same reconstruction comparison in a second low-energy effective model whose phase boundary includes deconfinement effects, or replace the polynomial fit with a ratio-of-polynomials resummation: if the effective convergence radius moves far from 146 MeV or the near-endpoint discrepancy moves far from 150%, the numbers are model- or ansatz-dependent; if they persist, the cautionary conclusion is generic.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the optimized regulator used in the bosonic and fermionic flow equations."},{"cited_title":"Role of mesonic fluctuations in the Polyakov loop extended quark-meson model at imaginary chemical potential","cited_arxiv_id":"1108.0735","evidence_quote":"Earlier effective-model study of imaginary chemical potential whose qualitative behavior the paper confirms."},{"cited_title":"Functional Renormalization Group Study of Thermodynamic Geometry Around the Phase Transition of Quantum Chromodynamic","cited_arxiv_id":"2312.00665","evidence_quote":"Provides the advection-diffusion formulation of the FRG flow used to integrate the effective potential."}],"review_version":1}