Finite-temperature critical point of heavy-quark QCD on large lattices

Hiroto Sugawara; Kazuyuki Kanaya; Masakiyo Kitazawa; Ryo Ashikawa; Shinji Ejiri

REVIEW 3 major objections 4 minor 12 references

Finite-temperature critical point of heavy-quark QCD on large lattices

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Using Binder-cumulant scaling on large lattices, this paper locates the finite-temperature critical point of heavy-quark QCD and finds the physical ratio $m_{PS}^{(CP)}/T_c \simeq 17$–$18$ with only mild lattice-spacing dependence.

desk verdict Careful WHOT-QCD proceedings: the new N_t=8 critical point is genuinely new but preliminary, and the continuum estimate rests on a Z(2)-fixed Binder analysis whose contamination-model uncertainty is real though transparently disclosed. read the letter →

arxiv 2412.00598 v1 pith:P3GOK63Z submitted 2024-11-30 hep-lat

Kazuyuki Kanaya , Ryo Ashikawa , Shinji Ejiri , Masakiyo Kitazawa , Hiroto Sugawara This is my paper · ORCID

classification hep-lat

keywords latticeQCDheavy-quarkregioncriticalpointBindercumulantfinite-sizescalinghoppingparameterexpansiondeconfinementtransitionZ(2)universality

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

On the Columbia plot, the map of the QCD transition as a function of quark masses, the heavy-quark corner is expected to be first-order and to end at a critical point. This paper determines that heavy-quark critical point by finite-size scaling of the Binder cumulant of the Polyakov loop, using the hopping parameter expansion to reach spatial volumes as large as $N_s=120$. On $N_t=6$ and $8$ lattices the extracted critical couplings are $\kappa_c=0.08769(7)$ and $0.09024(46)$, giving $m_{PS}^{(CP)}/T_c=18.07(2)$ and $17.2(2)$. The paper argues that this physical ratio changes little between $N_t=6$ and $N_t=8$, so the values are probably close to the continuum limit. If correct, this locates the end of the first-order transition line in physical units to about a percent.

What carries the argument

The load-bearing objects are the Binder cumulant $B_4$ of $\mathrm{Re}\,\hat\Omega$, treated as the magnetization of a $Z(2)$ spin system, and the hopping parameter expansion of the effective quark action. $B_4$ is a normalized fourth moment whose crossing height is universal at a $Z(2)$ critical point, so the critical coupling is found where $B_4$ becomes independent of the spatial aspect ratio $L_T=N_s/N_t$. The hopping expansion makes large spatial extents affordable, and the central mechanism for incorporating terms beyond next-to-leading order is an effective method that replaces truncated high-order terms by shifted low-order couplings, exploiting the strong linear correlation between the orders; convergence is checked by watching effective couplings stabilize as the truncation order $n_L$ is raised and by comparing with exact next-to-leading-order reweighting. Finite-size-scaling fits that include an energy-like contamination term turn the crossing into the quoted critical couplings.

What would settle it

Measure the same crossing ratio on $N_t=6$ lattices with aspect ratio $L_T=20$ or larger and fit it without fixing the universality class; if the fitted crossing height remains more than two standard deviations above the Ising value, the contamination model used here is not removing the extra contribution and the quoted $\kappa_c$ would need to shift.

Watch

Extended reading notes

Core claim

The central claim is that in two-flavor heavy-quark QCD the endpoint of the first-order deconfinement line lies at $\kappa_c=0.08769(7)$ on $N_t=6$ lattices and at $\kappa_c=0.09024(46)$ on $N_t=8$ lattices, after finite-size-scaling analysis of the Binder cumulant on spatial volumes up to aspect ratios $L_T=18$ and $15$, respectively. Using the zero-temperature pseudo-scalar meson mass, these couplings correspond to $m_{PS}^{(CP)}/T_c=18.07(2)$ and $17.2(2)$. The paper finds this ratio to vary by only about five percent between $N_t=6$ and $N_t=8$, and states that this small lattice-spacing dependence suggests the values are not far from the continuum limit. The analysis treats the Polyakov loop as the magnetization of a $Z(2)$ spin system; an apparent excess of the Binder crossing height at $N_t=6$ is accounted for by an added energy-like contamination term, after which the data are consistent with $Z(2)$ scaling.

Load-bearing premise

The central numbers depend on assuming the transition is of the three-dimensional Ising universality class and that a fitted correction term fully explains why the measured crossing ratio is higher than the Ising value; if that correction is incomplete, the critical coupling would move by more than the quoted uncertainty.

Editorial extensions

If this is right

If the quoted $\kappa_c$ values are correct, the first-order deconfinement region in the heavy-quark corner of the Columbia plot ends at a pseudo-scalar meson mass of about $17$–$18$ times $T_c$ in the continuum.
The small $N_t$ dependence of $m_{PS}^{(CP)}/T_c$ between $N_t=6$ and $N_t=8$ means a continuum extrapolation is not yet needed to locate the critical point to about a percent.
The validated effective high-order hopping expansion justifies applying the same method to $N_t\simeq 10$ lattices, where the hopping parameter is even larger and more orders are needed.
The two-flavor result can be translated to $2+1$ flavor QCD through the analytic flavor dependence of the hopping expansion, which the paper states is straightforward.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The near-flatness of $m_{PS}^{(CP)}/T_c$ across $N_t$ suggests a testable prediction: on $N_t=10$ and $12$ lattices the ratio should stay close to $17$–$18$, and a drift beyond a few percent would indicate that the apparent continuum near-independence was accidental.
Because the hopping expansion organizes the fermion determinant as a power series in $\kappa$, the same large-lattice machinery could be applied at finite baryon density, where the heavy-quark region has a milder sign problem; this could map the critical endpoint in the temperature-density plane.
The $Z(3)$ remnant asymmetry visible in the distribution of $\hat\Omega$ implies that finite-volume corrections to $Z(2)$ scaling are dominated by center-symmetry remnants on currently affordable lattices, so future $N_t>8$ studies should check whether even larger $L_T$ are needed before the scaling window opens.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: this is a straightforward, honest proceedings from WHOT-QCD, and the thing to know is that the new number---the N_t=8 critical point, kappa_c = 0.09024(46)---is real new data but officially preliminary, and the continuum extrapolation built on it is a preliminary three-point exercise, not a settled result. The N_t=6 analysis and the method were published earlier in PRD; the proceedings adds little methodologically but does give the first N_t=8 point and the first N_t=4,6,8 continuum estimate of m_PS/T_c.

What the paper does well: the convergence checks for the HPE are concrete and not decorative. They show the relative truncation error as a function of the order, they show the strong linear correlation between high-order terms, and they test the effective resummation against the exact NLO treatment on the N_t=6 phase diagram. The N_t=4 and N_t=6 FSS analyses are already in refereed journals, and the N_t=8 section is explicitly labeled preliminary, with larger errors and fewer volumes. That is decent scientific hygiene for a proceedings.

Where it is softer, and here the stress-test note mostly lands. The central N_t=6 value comes from a fit in which b4 and nu are fixed to Z(2) after the free fit gave b4 = 1.6297(84), more than two sigma above 1.604. An energy-like contamination term is then introduced to absorb the discrepancy. The quoted systematic (+11 over -0) is simply the difference between fits with and without that term; it does not cover the uncertainty in the contamination model itself. If the Z(3) remnant seen in Fig. 5 induces a correction that is not purely energy-like, kappa_c shifts outside the quoted errors. That is a real model-dependence, and the same fixed-Z(2) procedure is used for the N_t=8 point. The continuum statement at 17-18 therefore rests on an assumption, not on data alone. That said, the paper is fully transparent: it shows the free-fit b4, it shows the Z(3) asymmetry in the complex-plane histogram, and it reports the model-dependence as a systematic. This is an under-constrained choice, not a hidden one. My soundness score matches yours at 6.

Who should read this: anyone who needs the heavy-quark critical point as an anchor for the Columbia plot or as a cross-check for finite-density extrapolations. The N_t=8 number and the continuum estimate will be quoted even in this preliminary form, so the paper deserves to be read carefully. As a full paper, I would send it to referees; as a proceedings, this is the right level of detail and candor.

Referee Report

3 major / 4 minor

Summary. The manuscript reports a lattice study of the finite-temperature critical point in the heavy-quark region of two-flavor QCD, extending earlier N_t=4 work to N_t=6 and N_t=8 lattices with spatial aspect ratios up to L_T=18 and L_T=15. To reach these volumes the authors use the hopping parameter expansion combined with an effective shift of the leading-order couplings (eff[LO]/eff[NLO]), with checks of the convergence of the expansion and of the stability of the effective couplings. The Binder cumulant of Re Omega is studied along the transition line, and finite-size scaling fits give kappa_c=0.08769(7)(+11/-0) at N_t=6 and kappa_c=0.09024(46) at N_t=8, corresponding to m_PS^(CP)/T_c=18.07(2)(+11/-0) and 17.2(2), respectively. The paper concludes that the lattice-spacing dependence of this physical ratio is quite small up to N_t=8 and that these values are not far from the continuum limit.

Significance. If the determination is correct, it provides a high-precision anchor for the heavy-quark QCD critical point and demonstrates that the hopping parameter expansion with effective higher-order terms can reach much larger spatial volumes and finer lattices than previously possible. The paper is transparent about the preliminary nature of the N_t=8 result, and the N_t=6 analysis is carefully cross-checked with respect to the definition of the transition line and the convergence of the higher-order terms. The main reservation, which I detail below, is that the central N_t=6 value is obtained by fixing the Binder crossing height b4 to its Z(2) value after the unconstrained fit gives b4=1.6297(84), more than 2 sigma above 1.604. This is not a circularity problem, because b4 is used as an external benchmark rather than a fitted output, but it makes the central value dependent on a contamination model whose systematic uncertainty is not yet quantified.

major comments (3)

[Sec. 4.2, Eq. (8), Fig. 4 (right)] The unconstrained finite-size scaling fit at N_t=6 gives b4=1.6297(84), more than 2 sigma above the Z(2) value 1.604, and the asymmetry associated with the Z(3) remnant is visible even at L_T=15. The central value is then obtained by fixing nu and b4 to their Z(2) values and adding an energy-like contamination term with a fitted amplitude, because the full six-parameter fit is unstable. The quoted systematic error (+11/-0) is only the difference between fits with and without this contamination term; it does not include uncertainty in the form of the contamination or a possible non-Z(2) correction. Since b4 is the observable that identifies the universality class, this is a load-bearing assumption. I ask the authors to demonstrate quantitatively that the contamination model accounts for the 2-sigma discrepancy, for example by reporting the fitted contamination amplitude and its predicted contribution to b4, by varying the form of the contamination term, or by fitting an independent observable with a different contamination pattern.
[Sec. 5, Eqs. (9)-(11), and abstract] The abstract states that the paper attempts a preliminary continuum extrapolation, but no continuum extrapolation is actually shown. Equations (9)-(11) list three values, m_PS/T_c=16.30(3), 18.07(2), and 17.2(2), which are non-monotonic in N_t, and the text only says that the lattice-spacing dependence is small. The difference between the N_t=6 and N_t=8 values is approximately 0.87, which is about four times the quoted error of the N_t=8 point. To support the claim that these values are not far from the continuum limit, an explicit extrapolation fit (for example linear in 1/N_t^2 or including O(a) corrections) with an assessment of systematic uncertainty is needed; otherwise the wording should be limited to a statement about the available lattice spacings.
[Sec. 4.3] The N_t=8 determination, although labeled preliminary, is used in Sec. 5 as evidence for the small lattice-spacing dependence. The quoted kappa_c=0.09024(46) is obtained after fixing nu (and effectively b4) to the Z(2) values, but the text does not state whether the same energy-like contamination correction is applied at N_t=8 or what systematic error is assigned to the contamination model. Because the N_t=8 point carries a large weight in the comparison underlying the continuum statement, this systematic should be quantified before the comparison is used as evidence.

minor comments (4)

[Fig. 4 caption] The caption refers to 'Eq. (38)' and 'Eq. (40)', but the manuscript contains no equations with these numbers; these references are presumably to equations in ref. [9] and should be replaced with local equations or a clear citation.
[Sec. 4.3] The phrase 'to be compared with a previous result 1.1135(8)' is unclear: the number 1.1135(8) cannot be kappa_c in the same notation as the rest of the paper, and no conversion or definition is given. Please correct or clarify.
[Sec. 2.1] There are small typographical errors, including 'where M is he Wilson quark kernel' and 'obsereved' in the Fig. 3 caption; these should be corrected.
[Sec. 5] The physical conversion uses zero-temperature pseudo-scalar masses from ref. [11], which is listed as 'unpublished'. Please provide the interpolated mass values used in Eqs. (9)-(11) or cite a published source, so that the conversion is reproducible.

Circularity Check

2 steps flagged · score 2.0 of 10

No significant circularity; the N_t=6 central value is model-conditional (Z(2) crossing height imposed with a one-sided systematic) and the physical-units conversion cites unpublished in-house masses, but the lattice kappa_c values are data-determined and agree with an independent study.

other [Sec. 4.2 (N_t=6 Binder FSS fit); Fig. 4 right panel]
"Because the full six-parameter fit turned out to be unstable, we perform fits fixing the critical exponents and b4 to their Z(2) values. We find that the fits work well with acceptable χ2/dof when LT ≥ 10. We thus conclude that our data obtained on spatially large lattices are consistent with the Z(2) scaling."

The conclusion that the data are 'consistent with the Z(2) scaling' is drawn from a fit that already fixes b4 and the exponents to their Z(2) values, so the consistency test partly assumes what it asserts. The central kappa_c is taken from this constrained fit at the imposed height b4 = 1.604, while the unconstrained fit gives b4 = 1.6297(84), more than 2 sigma above the Z(2) value. The quoted systematic (+11/-0) covers only the difference with the fit disregarding the energy-like contamination, not the possibility of a non-Z2 crossing height.
other [Sec. 5, Eqs. (9)-(11) and ref. [11]]
"Using previous results of the pseudo-scalar meson mass mPS at T = 0 [5, 11], we express the results of the critical point in terms of physical observables."

Ref. [11] is an unpublished in-house WHOT-QCD result used as the zero-temperature input converting kappa_c into mPS^(CP)/Tc in Eqs. (9)-(11). The continuum-limit suggestion ('the lattice-spacing dependence is quite small up to Nt = 8 in this combination') depends partly on this unverified self-citation. However, the lattice kappa_c values are independent of [11] and the conversion also uses external ref. [5], so this is a data-provenance weakness rather than a load-bearing circular argument.

full rationale

The paper's core derivation is a Binder-cumulant finite-size-scaling determination of the heavy-quark QCD critical point from new N_t=6 and 8 simulations, extending the earlier N_t=4 study. No equation reduces the output (kappa_c, mPS/Tc) to an input by construction: the FSS fits locate kappa_c from the shape of data-fitted B4 curves, and the Z(2) universal value 1.604 is an external benchmark, not a fitted output. The hopping-parameter-expansion method with effective high-order incorporation (eff[LO]/eff[NLO]) is cited from the same group's prior work [8,9], but the convergence is checked numerically in this paper (Fig. 1 truncation stability; Fig. 3 agreement between eff[LO] and eff[NLO] lines), so the self-citation is backed by in-paper evidence and does not smuggle in the target result. The N_t=6 tension (free-fit b4 = 1.6297(84), >2 sigma above Z(2)) is reported transparently, and the (+11/-0) systematic is honest about fit-model dependence; that is a correctness/model-selection risk, not circularity, since the alternative was shown explicitly. Agreement with the independent study [5] (kappa_c = 0.0877(9) at N_t=6) provides an external benchmark confirmation. The two minor caveats are the mildly self-referential 'consistent with Z(2) scaling' inference after imposing Z(2) values and the unpublished in-house zero-T mass input [11] used only for the physical-units conversion. Overall the central lattice claim has independent content, so the circularity burden is small.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The extraction uses a controlled HPE plus an empirical linear-correlation rescaling, a Z(2) scaling assumption, and a fitted energy-like contamination term. These are the prices paid to reach large volumes on fine lattices. No new physical entities are introduced; the effective shifts and contamination amplitude are data-determined scheme parameters.

free parameters (3)

effective Polyakov-loop coupling shift in eff[NLO] method = not tabulated; stable for truncation order n_L > 14 at N_t=6
Section 3.2: high-order PLT terms are replaced by rescaled LO/NLO terms; the rescaling slope is determined from the simulated configurations.
effective gauge-coupling shift in eff[LO]/eff[NLO] method = not tabulated
Section 3.2: Wilson-loop high-order contributions are absorbed into a shifted beta; the shift is read off from the data.
energy-like operator contamination amplitude in B4 FSS fit = not quoted
Section 4.2: central kappa_c at N_t=6 is taken from a fit with Z(2) exponents fixed and an additional contamination term with fitted amplitude.

assumptions (6)

domain assumption The hopping parameter expansion of the quark determinant converges sufficiently well at kappa_c for N_t = 6 and 8, and the truncation after effective resummation is accurate.
Invoked in Secs. 2.2 and 3.2; at N_t=8 'we need NNLO and higher orders to achieve a good accuracy', and the paper relies on the effective method to supply them.
domain assumption The deconfinement critical point in the heavy-quark region is in the Z(2) (Ising) universality class, with Binder cumulant crossing height b4 = 1.604 and critical exponent nu = 0.630.
Used in Secs. 4.2 and 4.3, where fits fix b4 and nu to the Z(2) values; the data are consistent only after adding an energy-like contamination term.
ad hoc to paper Re(Omega) (the real part of the Polyakov loop) behaves as the Z(2) magnetization up to a small contamination from an energy-like operator.
Introduced in Sec. 4.2 when the unconstrained Binder cumulant fit gives b4 more than 2 sigma above the Z(2) value; the contamination model is not derived from first principles.
domain assumption For large spatial volumes, spatial winding contributions to the hopping parameter expansion are negligible.
Stated in Sec. 2.1: 'we assume that the spatial extent is sufficiently large such that the effects of spatial windings are negligible.'
domain assumption Finite-size scaling is dominated by the leading singularity; the transition/crossover line can be defined by the minimum of B4, and the three definitions converge for large L_T.
Assumed in Sec. 4 and checked for three definitions; the B4 FSS determines the CP.
domain assumption Reweighting brings the NLO effects and continuous parameter variation into the measurement without residual overlap problems.
Sec. 3.1 states the overlap problem is 'largely resolved' by including LO in the generation; this is a standard reweighting assumption.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Finite-temperature critical point of heavy-quark QCD on large lattices." pith.science (2026). https://pith.science/paper/P3GOK63Z

@misc{pith2026241200598,
  author       = {Pith},
  title        = {Pith review of: Finite-temperature critical point of heavy-quark QCD on large lattices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/P3GOK63Z}},
  note         = {Machine review of arXiv:2412.00598}
}

abstract

We study the finite-temperature critical point of QCD in the heavy-quark region by a scaling study of the Binder cumulant on large lattices. Extending our previous study at $N_t=4$, we perform simulations on $N_t=6$ and 8 lattices with spatial volumes up to the aspect ratio $LT=N_s/N_t=18$ and 15 ($N_s=108$ and 120), respectively, to determine the critical point in the thermodynamic limit with a high precision. To enable simulations with large spatial volumes, we adopt the hopping parameter expansion combined with a method to effectively incorporate high order terms of the expansion. The reliability of the method is confirmed by examining the effect of high order terms. Using the results of the critical point at $N_t=4$, 6, and 8, we also attempt a preliminary continuum extrapolation of the critical point in physical units.

Figures

Figures reproduced from arXiv: 2412.00598 by the authors.

**Figure 1.** Left: Relative deviation from the exact value due to truncation of the HPE for the PLT contribution in the effective action at 𝑁𝑡 = 6 [9]. Right: Scatter plot of 𝑛th order PLT terms of the HPE it vs. the LO Polyakov-loop, observed on an 𝑁𝑡 = 6 lattice near the CP [9]. 2.1 Hopping parameter expansion In the heavy-quark region 𝜅 ≪ 1, we may adopt the hopping parameter expansion (HPE): 𝑆eff = 𝑆g − 𝑁f𝑁 3 𝑠 𝑁𝑡 ∑︁∞ 𝑛=1 [… view at source ↗

**Figure 2.** Simulation of heavy-quark QCD incorporating LO and NLO terms of the HPE [7]. From the figure we find that, around the CP at 𝑁𝑡 = 4, 𝜅𝑐 = 0.0603(4) [7] for 𝑁f = 2, the effective action in the LO approximation (red curve) may have an error of 10% maximally, while the NLO (blue curve) is guaranteed to be fairly accurate. Around 𝜅𝑐 = 0.08769(7) (+11 −0 ) [9] at 𝑁𝑡 = 6, the LO and NLO approximations may have error of 25 … view at source ↗

**Figure 3.** Left: Phase diagram of two-flavor QCD in the (𝛽, 𝜅) plane at 𝑁𝑡 = 6 [9]. The full and dotted lines represent the first-order transition line and the crossover line, while the symbols at the end of the first-order lines represent the CP. The green, blue, and red lines show the result of the NLO analysis discussed in Sec. 3.1, the eff[LO] method, and the eff[NLO] method discussed in Sec. 3.2, respectively. Right: Effe… view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: Binder cumulant of the Polyakov loop along the transition line, obtained with various spatial volumes. Left: Results at 𝑁𝑡 = 4 [7]. Right: Results at 𝑁𝑡 = 6 [9]. The symbols with both vertical and horizontal error bars are the results of conventional FSS fits identifyi…

**Figure 5.** Figure 5: Left: Two-dimensional histogram of Ωˆ on the complex plane at the critical point at 𝑁𝑡 = 6 and 𝐿𝑇 = 6 [9]. Right: Preliminary results of the Binder cumulant obtained at 𝑁𝑡 = 8 [10]. The symbol with both vertical and horizontal error bars is the crossing point estimated…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 4 canonical work pages

[1]

Y.Kuramashi, Y.Nakamura, H.Ohno, andS.Takeda, Natureofthephasetransitionforfinite temperature𝑁f = 3QCDwithnonperturbatively 𝑂(𝑎) improvedWilsonfermionsat 𝑁𝑡 = 12, Phys. Rev. D101(2020) 054509 [hep-lat/2001.04398]

arXiv 2020
[2]

O.Philipsen, LatticeconstraintsontheQCDchiralphasetransitionatfinitetemperatureand baryon density, Symmetry13(2021) 2079 [hep-lat/2111.03590]

arXiv 2021
[3]

L. Dini, P. Hegde, F. Karsch, A. Lahiri, C. Schmidt, and S. Sharma,The chiral phase transition in three-flavor QCD from lattice QCD, Phys. Rev. D 105 (2022) 034510 [hep-lat/2111.12599]

arXiv 2022
[4]

Binder,Finite size scaling analysis of Ising model block distribution functions, Z

K. Binder,Finite size scaling analysis of Ising model block distribution functions, Z. Phys. B 43 (1981) 119

work page 1981
[5]

F.Cuteri,O.Philipsen,A.Schön,andA.Sciarra, DeconfinementcriticalpointoflatticeQCD with𝑁f = 2 Wilson fermions, Phys. Rev. D103 (2021) 014513 [hep-lat/2009.14033]

arXiv 2021
[6]

Ejiri, S

S. Ejiri, S. Itagaki, R. Iwami, K. Kanaya, M. Kitazawa, A. Kiyohara, M. Shirogane, and T. Umeda, End point of the first-order phase transition of QCD in the heavy quark region by reweightingfromquenchedQCD,Phys.Rev.D 101(2020)054505[ hep-lat/1912.10500]

arXiv 2020
[7]

A.Kiyohara,M.Kitazawa,S.Ejiri,andK.Kanaya, Finite-sizescalingaroundthecriticalpoint intheheavyquarkregionofQCD,Phys.Rev.D 104(2021)114509[ hep-lat/2108.00118]

arXiv 2021
[8]

Wakabayashi, S

N. Wakabayashi, S. Ejiri, K. Kanaya, and M. Kitazawa, Scope and convergence of the hopping parameter expansionin finite-temperature quantum chromodynamics with heavy quarks around the critical point, Prog. Theor. Exp. Phys. 2022 (2022) 033B05 [hep-lat/2112.06340]

arXiv 2022

Show all 12 references

[9]

R.Ashikawa,M.Kitazawa,S ˙Ejiri,andK.Kanaya, High-precisionanalysisofthecriticalpoint in heavy-quark QCD at𝑁𝑡 = 6, Phys. Rev. D110 (2024) 074508 [hep-lat/2407.09156]

2024 arXiv
[10]

Sugawara et al., for the WHOT-QCD Collaboration, in preparation

H. Sugawara et al., for the WHOT-QCD Collaboration, in preparation
[11]

Itagaki et al., for the WHOT-QCD Collaboration, unpublished

S. Itagaki et al., for the WHOT-QCD Collaboration, unpublished
[12]

Ser.523 (2014) 012046;http://bridge.kek.jp/Lattice-code/

S.Uedaetal., DevelopmentofanobjectorientedlatticeQCDcode’Bridge++’,J.Phys.Conf. Ser.523 (2014) 012046;http://bridge.kek.jp/Lattice-code/. 9

2014

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [1]

Y.Kuramashi, Y.Nakamura, H.Ohno, andS.Takeda, Natureofthephasetransitionforfinite temperature𝑁f = 3QCDwithnonperturbatively 𝑂(𝑎) improvedWilsonfermionsat 𝑁𝑡 = 12, Phys. Rev. D101(2020) 054509 [hep-lat/2001.04398]

arXiv 2020

[2] [2]

O.Philipsen, LatticeconstraintsontheQCDchiralphasetransitionatfinitetemperatureand baryon density, Symmetry13(2021) 2079 [hep-lat/2111.03590]

arXiv 2021

[3] [3]

L. Dini, P. Hegde, F. Karsch, A. Lahiri, C. Schmidt, and S. Sharma,The chiral phase transition in three-flavor QCD from lattice QCD, Phys. Rev. D 105 (2022) 034510 [hep-lat/2111.12599]

arXiv 2022

[4] [4]

Binder,Finite size scaling analysis of Ising model block distribution functions, Z

K. Binder,Finite size scaling analysis of Ising model block distribution functions, Z. Phys. B 43 (1981) 119

work page 1981

[5] [5]

F.Cuteri,O.Philipsen,A.Schön,andA.Sciarra, DeconfinementcriticalpointoflatticeQCD with𝑁f = 2 Wilson fermions, Phys. Rev. D103 (2021) 014513 [hep-lat/2009.14033]

arXiv 2021

[6] [6]

Ejiri, S

S. Ejiri, S. Itagaki, R. Iwami, K. Kanaya, M. Kitazawa, A. Kiyohara, M. Shirogane, and T. Umeda, End point of the first-order phase transition of QCD in the heavy quark region by reweightingfromquenchedQCD,Phys.Rev.D 101(2020)054505[ hep-lat/1912.10500]

arXiv 2020

[7] [7]

A.Kiyohara,M.Kitazawa,S.Ejiri,andK.Kanaya, Finite-sizescalingaroundthecriticalpoint intheheavyquarkregionofQCD,Phys.Rev.D 104(2021)114509[ hep-lat/2108.00118]

arXiv 2021

[8] [8]

Wakabayashi, S

N. Wakabayashi, S. Ejiri, K. Kanaya, and M. Kitazawa, Scope and convergence of the hopping parameter expansionin finite-temperature quantum chromodynamics with heavy quarks around the critical point, Prog. Theor. Exp. Phys. 2022 (2022) 033B05 [hep-lat/2112.06340]

arXiv 2022