{"id":"1ce6a083-4628-4bbd-8475-02e1b21bf6d5","arxiv_id":"hep-ph/9907240","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A systematic effective-field-theory matching of QCD to potential NRQCD defines the heavy-quark potential cleanly, computes its scale running, and predicts degeneracies among the gluonic excitations of static quarks.","lead":"The authors build a step-by-step effective theory for heavy quark-antiquark pairs, matching QCD first to non-relativistic QCD and then to a theory of potentials plus soft gluons. It gives a clean definition of the quark-antiquark potential and shows where the usual potential-model picture stops being accurate.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The central EFT matching is well-supported; the specific soft spot is the octet IR-freedom claim in §3.4.2, which contradicts refs. [30, 9] and rests on a structural argument rather than the explicit two-loop calculation the authors say is needed.","rationale":"The reader's weakest_assumption — the existence of a perturbative window mv >> μ >> mv², Λ_QCD — is correct and is explicitly assumed in §2.2. The paper is honest about this and about the regimes where it fails (§§4–5). The central EFT matching claim is supported by the explicit leading-log double calculation (Eqs. 33–35) and by the renormalon check in Appendix B. My concern is narrower: the octet IR-freedom claim in §3.4.2 is a specific assertion that contradicts prior literature and is not independently verified, because the two-loop octet potential is unknown. The authors flag this gap themselves, so it is an acknowledged limitation rather than a hidden flaw. It does not undermine the central singlet matching result, but it does mean the paper's claim that the pNRQCD potentials are clean matching coefficients is not fully established for the octet sector. A separate secondary concern is the Π_g anomaly in the lattice comparison of Section 6: the authors explain it away by hypothesizing that the plotted Π_g is a superposition of the real Π_g and the Π′_g, which is plausible but untested. Neither concern changes the verdict on the central claim, which remains well-supported by the explicit singlet matching and the renormalon analysis. The reader's ACCEPT verdict is appropriate; the octet IR-freedom claim should be read as an open item rather than a settled result.","tokens_in":31028,"tokens_out":16200,"duration_ms":287080,"concrete_test":"Compute the two-loop octet potential (v0)_two-loops from the intermediate steps of the Schröder/Peter calculations [16], as the authors themselves suggest. If it contains an IR divergence at O(α_s²), the claim in §3.4.2 fails and the octet potential is not a clean matching coefficient at that order; if it is IR-finite, the claim is confirmed. This is a finite, well-defined calculation that the authors identify as the missing independent test.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.4.2 asserts that the octet matching potential V_o is free of IR divergences at O(α_s²), explicitly contradicting refs. [30, 9]. The evidence offered is Eq. (38), which shows the leading log is O(α_s³), plus the statement that 'up to a color factor the contribution of Eq. (37) is equal, at leading order in α_s, to the one calculated in the singlet case.' But the O(α_s²) piece of V_o is (v0)_two-loops, which the paper states is unknown. The authors call the structural argument a 'proof' but also say the two-loop octet potential 'would provide an independent test to our proof that the octet matching potential is free of IR divergences at two loops.' So the claim that V_o is IR-free at O(α_s²) is asserted, not demonstrated. This matters because if V_o had an O(α_s²) IR divergence, the octet potential would not be a clean matching coefficient at that order, and the paper's broader claim that the pNRQCD potentials are matching coefficients whose scale dependence cancels against ultrasoft gluons would need qualification for the octet sector. The singlet case is checked explicitly by the double calculation of Eqs. (33)–(35); the octet case is not. The authors acknowledge the gap, so this is an unresolved limitation rather than a hidden error, but it is the one place where a stated conclusion outruns the calculation actually performed.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript develops potential NRQCD (pNRQCD), the effective theory obtained by integrating the soft scale mv out of NRQCD, for heavy-quark–antiquark systems with r < 1/Λ_QCD. The singlet and octet fields depend on relative and centre-of-mass coordinates, the potentials Vs, Vo, VA, VB are defined as r-dependent matching coefficients, and ultrasoft gluons are retained so that non-potential (retardation) effects are systematically included. The matching NRQCD → pNRQCD is carried out at leading order in 1/m and up to O(r²) in the multipole expansion, the leading-logarithm running of the potentials and of the singlet normalization factor is computed (Table 1), the cases Λ_QCD ≲ mv² and mv ≫ Λ_QCD ≫ mv² are treated separately (with a pNRQCD′ in the latter), model-independent predictions for the short-distance behaviour and approximate degeneracies of the gluonic excitations between static quarks are derived (Eq. 50, Table 2), and the cancellation of the IR renormalons of the static potential inside the EFT is exhibited (Appendix B).","tokens_in":31304,"tokens_out":5168,"duration_ms":100836,"significance":"If the construction holds, this is an important structural paper: it turns the static potential from an object with ambiguous infrared status into a matching coefficient with a definite scale dependence, and it supplies the calculational framework in which that scale dependence cancels against ultrasoft dynamics in any physical observable. Several features deserve explicit credit. The leading-logarithm result is obtained twice by independent routes — from the pNRQCD correlator (Eq. 33) and from the Wilson loop (Eqs. 34, 35) — and the two agree; this is the strongest internal check in the paper. The renormalon analysis of Appendix B is explicit rather than programmatic: the Borel-non-summable coefficients D_n^(0) and D_n^(2) are computed (Eq. 65), the O(Λ_QCD) and O(Λ_QCD³ r²) ambiguities are identified, and the cancellation against the pNRQCD UV renormalon is verified at next-to-leading order in the multipole expansion (Eqs. 70, 71). The degeneracy predictions of Eq. (50) are falsifiable and are confronted with quenched lattice data (ref. 41) rather than asserted, and the one place where the data do not fit — the Π_g state — is flagged rather than explained away. The paper is also,","major_comments":[{"comment":"The claim that the octet matching potential Vo is free of infrared divergences at relative O(α_s²), which explicitly contradicts refs. [30, 9], is asserted rather than demonstrated. The evidence offered is Eq. (38), which shows only that the leading logarithm is O(α_s³), together with the statement that 'up to a color factor the contribution of Eq. (37) is equal, at leading order in α_s, to the one calculated in the singlet case.' But the O(α_s²) piece of Vo is (v0)_two-loops, which the manuscript states is unknown, and the authors themselves write in the Conclusions that the two-loop octet potential 'would provide an independent test to our proof.' The singlet sector is checked by the double calculation of Eqs. (33)–(35); the octet sector is not. This is a correctness risk rather than an internal inconsistency, and it has a concrete test: (v0)_two-loops should be extractable from the intermediate steps of the two-loop calculation of ref. [16], as the authors note. Either perform that extraction, or soften the wording so that the O(α_s²) IR freedom of Vo is presented as a consequence of the general matching structure subject to a forthcoming check, rather than as an established result that overturns the earlier literature.","section":"§3.4.2, Eq. (38) and the paragraph following it","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The conclusion TE < TB is drawn from the lattice data of ref. [41] as displayed in Fig. 9, but the same paragraph argues that the plotted Π_g is likely 'a superposition of the real Π_g in the doublet and the Π'_g in the triplet.' Since Π_g is generated by an electric-type operator (r × E in Table 2), a misidentified Π_g bears directly on the extraction of the correlation length of the E correlator, which is what Eq. (57) is about. As written, the paper infers a property of the E correlator from a state it simultaneously distrusts. Please either qualify Eq. (57) accordingly, or add an argument that the inference is robust to the doublet/triplet mixing — for example by showing which combination of the measured levels enters the E channel and that the inequality survives for any admissible splitting of the two states.","section":"§6, Eq. (57) and the surrounding discussion of Fig. 9","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"'the widespread believe that Vo suffers from IR divergences' should read 'the widespread belief'.","section":"§3.4.2, first paragraph","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"'it does not matter that the IR behavior is no accurately described' should read 'not accurately described'.","section":"Appendix B, sentence preceding Eq. (70)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The correlator has an unbalanced parenthesis: '⟨H^a(T/2)φ(T/2,−T/2)^{adj ab}H^b(−T/2))⟩' closes a bracket that was never opened.","section":"Eq. (52)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The arguments of the Schwinger strings are garbled in the typesetting, appearing as 'x1 + x22 ; T2' and 'y1 + y22 ; −T2'. These should presumably read (x1+x2)/2 and T/2; please clean up both equations, since the octet matching definition is one of the paper's central objects.","section":"Eqs. (19) and (27)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The exponential structure is unreadable as typeset ('e^{−i(T/2−t)hs} r^{i} e^{−i(t−t')ho} r^{j} e^{−i(t'+T/2)hs}' appears as 'hsri' and 'horje'). Please re-typeset this equation.","section":"Eq. (43)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The factor written as 'r r · E^a' should be made explicit as r^i/r times r·E^a, or the shorthand defined once; as it stands the reader has to reconstruct the index structure.","section":"Eqs. (6) and (7)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that 'αVs is not a short-distance quantity as αs (in the MS scheme), since it depends on μ' is confusing, because α_s in the MS scheme also depends on μ. The intended point is presumably that αVs carries an explicit dependence on the infrared scale 1/r through ln rμ, i.e. that it is a matching coefficient rather than a purely short-distance coupling. Please rephrase.","section":"§4, paragraph following Eq. (41)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The notation for the excited Σ states is inconsistent between the table (Σ^{+ ′g}) and the text (Σ^{+′}_g), and the placement of the prime relative to the g/u label varies. Please adopt one convention throughout.","section":"Table 2 and Eq. (50)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The figure is reproduced from ref. [37], which in turn displays data from ref. [41]. Please confirm that permission has been obtained and state the provenance of the underlying lattice data in the caption.","section":"Fig. 9","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The additive 'constant' inside the curly braces is left unspecified. Either give it or add a sentence explaining why it cannot affect the matching of Vs and Zs.","section":"Eq. (33)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The assertion that ⟨φ(T/2,−T/2)^{adj ab}⟩ 'does not give contribution to the potential at any finite order in perturbation theory' would be easier to accept with one explicit sentence of justification, for instance that it exponentiates into a self-energy-like factor depending only on the endpoints.","section":"§3.2.2, discussion of the adjoint Schwinger line","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is sound and the central matching is checked twice, so my recommendation is minor revision rather than more. The two substantive points are both acknowledged limitations that the authors have left as promissory notes: the two-loop octet potential (which they say would test their own IR-freedom claim) and the Π_g misidentification in the lattice comparison. Both are fixable within the present scope. On novelty relative to the authors' own earlier pNRQCD letter (ref. [11]) and the companion leading-log paper (ref. [12]): the present manuscript is substantially more than a consolidation — the octet matching, the running of VA and VB, the pNRQCD′ construction of §5, the degeneracy analysis of §6 and Appendix B are new — and the 'Note Added' citing ref. [44] is handled appropriately."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Worth your time if you work on quarkonium or precision QCD. The core of the paper is the NRQCD-to-pNRQCD matching at next-to-leading order in the multipole expansion, and it is done properly: the leading-log coefficients for the singlet potential and the singlet normalization factor are obtained twice, once from the Wilson loop (Eqs. 34, 35) and once from the pNRQCD correlator (Eq. 33), and the ln(r mu) pieces agree. That double calculation is the strongest thing in the paper and it is not a trivial check. Appendix B works the renormalon cancellation coefficient by coefficient, with the non-alternating D_n written out and the O(Lambda_QCD) and O(Lambda_QCD^3 r^2) ambiguities matched against the appropriate pNRQCD diagrams. Table 1 gives the one-loop running of V_s, V_o, V_A, V_B and Z_s, which is the piece most people will actually cite. The static-limit degeneracies of Eq. 50 are a genuine prediction, checked against Juge-Kuti-Morningstar lattice data, and the authors flag the Pi_g state that does not fit instead of burying it. That is the right instinct and it is rare enough to be worth noting. On the stress-test concern about §3.4.2: it lands, but smaller than it sounds. The claim that V_o is IR-free at O(alpha_s^2) is argued structurally, not computed, and the O(alpha_s^2) piece (v_0)_two-loops is stated to be unknown. The authors themselves say the two-loop octet potential would provide an independent test of the claim. So this is an acknowledged gap, not a hidden error, and it does not touch the singlet sector where the real evidence sits. The other soft spots are honestly labelled: the massive-gluon model in §5.1.1 is called crude by the authors, and the Lambda_QCD ~ mv regime is deferred to work in preparation. The perturbative window assumed in §2.2 (mv >> mu >> mv^2, Lambda_QCD) is a scope condition, stated plainly, not a flaw. This is a framework paper that will be cited for the running coefficients and the renormalon argument rather than read cover to cover. Send it out.","headline":"The matching and renormalon work is solid and cross-checked; the octet IR-freedom claim in §3.4.2 is the one place a conclusion outruns the calculation, and the authors say so themselves.","tokens_in":32023,"tokens_out":1127,"would_cite":true,"duration_ms":19213,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["12.38.Aw","12.38.Bx","12.39.Hg"],"model":"selfhosted/deepseek-v4.1-flash","headline":"The static potential is redefined as a matching coefficient whose scale dependence cancels against ultrasoft gluons, and the paper checks the cancellation at leading logarithm.","keywords":["potential NRQCD","NRQCD","heavy quarkonium","static potential","effective field theory","multipole expansion","infrared renormalon","gluelump"],"falsifier":"A calculation or lattice measurement showing that the coefficient of $ln(r mu)$ in the singlet potential differs from $C_A^3 alpha_s^4 / 12 pi$, or that the $mu$ dependence of the potential fails to cancel against ultrasoft contributions in a physical observable, would break the central claim; equally, short-distance lattice data on gluonic excitations that show no approach to a common octet-potential shape plus constant, and no $Sigma_g^{+prime} - Pi_g$ or $Sigma_u^+ - Pi_u^{prime} - Delta_u$ degeneracy, would undercut the static-limit predictions.","tokens_in":32314,"feed_emoji":"⚛️","tokens_out":3250,"duration_ms":63598,"temperature":0.7,"pith_summary":"Heavy quarkonium lives at three widely separated scales: the quark mass $m$, the soft scale $mv$ set by the inverse size of the bound state, and the ultrasoft scale $mv^2$ set by the binding energy. The paper's programme is to integrate these out one at a time: QCD to NRQCD at the mass scale, then NRQCD to a new theory, potential NRQCD, at the soft scale. In that second step the potential stops being a fixed function of distance and becomes a matching coefficient with a subtraction point $mu$; the infrared divergences that have long plagued the perturbative static potential are simply the ultraviolet divergences of the effective theory, and they cancel in any physical quantity once ultrasoft gluons are included. The authors verify this explicitly by computing the leading logarithm of the singlet potential and of the singlet normalization factor in two independent ways, from the static Wilson loop and from the effective-theory correlator, and getting the same $ln(r mu)$ coefficients. They then use the same machinery to predict approximate degeneracies among the gluonic excitations of a static quark-antiquark pair and to show how the renormalon ambiguities of the static potential are absorbed inside the effective theory.","feed_headline":"Quarkonium's static potential is redefined as a matching coefficient","feed_subtitle":"Cut off at the soft scale, the potential's infrared divergences match the effective theory's ultraviolet ones.","key_machinery":"The load-bearing object is the pNRQCD Lagrangian at leading order in $1/m$ and at order $r$ in the multipole expansion, built from a colour-singlet and a colour-octet heavy-quark-antiquark field that depend on the relative coordinate $r$ and the centre-of-mass coordinate $R$, coupled to ultrasoft gluons evaluated only at $R$. The potentials $V_s(r)$, $V_o(r)$, $V_A(r)$, $V_B(r)$ are matching coefficients fixed by equating NRQCD and pNRQCD Green functions; at leading order $V_s$ and $V_o$ are read off from the large-$T$ behaviour of the static Wilson loop and of the Wilson loop with colour-matrix insertions, $i/T ln<W> = u_0(r) + i u_1(r) T + ...$. The multipole expansion is what makes the matching systematic: each power of $r$ brings in a definite set of operators, and at order $r^2$ a gluon correlator $<E_i(t) phi(t,t') E_j(t')>$ appears whose $ln(V_o - V_s)$ dependence reproduces the Wilson-loop infrared logarithm. The same expansion, applied to the static limit, gives the classification of gluonic excitations by $O(3) times C$ rather than the full $D_{infty h}$ of NRQCD, and the short-distance formula $V_H(r) = V_o(r) + Lambda_H$.","core_discovery":"The central claim is that potential NRQCD, obtained by integrating the soft scale $mv$ out of NRQCD, is the correct effective field theory at ultrasoft scales, and that the static potential defined within it is a matching coefficient rather than a short-distance observable. Concretely: the three-loop infrared divergence of the potential and the two-loop infrared divergence of the singlet normalization factor computed from the static Wilson loop in NRQCD match the corresponding ultraviolet divergences of pNRQCD, so the $mu$ dependence of the potential cancels against ultrasoft gluons in physical observables. The authors demonstrate this by calculating the leading-logarithm corrections to the singlet potential, $alpha_{V_s}(r,mu) = alpha_s(r) + ... + (C_A^3/12)(alpha_s^3/pi) ln(r mu)$, and to the octet potential, and by obtaining the one-loop running of the singlet-octet and octet-octet matching coefficients $V_A$ and $V_B$. They also show that the widespread belief that the octet potential suffers infrared divergences at relative order $alpha_s^2$ is incorrect: the leading ultrasoft correction to $V_o$ is of relative order $alpha_s^3$.","pith_inferences":["If the same matching logic is carried through for a system where the perturbative window $mv \\gg \\mu \\gg mv^2, \\Lambda_{\\rm QCD}$ does not exist, the effective-theory bookkeeping survives but every coefficient in the running table has to be replaced by a nonperturbative matrix element; the paper flags this for charmonium and bottomonium but does not attempt it.","The prediction that the $\\Pi_g$ curve on the lattice is a superposition of two distinct states is a concrete, checkable claim that a finer short-distance lattice study could settle, and it would be a clean test of the multipole-expansion picture rather than of perturbation theory.","The statement that in pure gluodynamics with a gluelump gap of order $\\Lambda_{\\rm QCD}$ no ultrasoft gluon fields survive in the lower effective theory amounts to saying that quarkonium there obeys an exact potential description; if that held, non-potential effects in pure-glue quarkonium would be an artefact of the truncation rather than a physical effect.","The renormalon-cancellation argument is presented at the level of the leading and next-to-leading multipole orders only, so a natural extension is to ask whether the cancellation persists at the next order in $r$, where the gluon correlator acquires additional tensor structure."],"forward_implications":["The scale dependence of the singlet potential is now known at leading logarithm, which is the missing ingredient for pushing next-to-next-to-leading-order determinations of the bottom and top masses from quarkonium observables to the next logarithmic order.","The octet potential is free of infrared divergences at two loops, so a two-loop calculation of $V_o$ from the intermediate steps of existing singlet calculations would provide an independent check of the matching.","In the short-distance limit the gluonic excitations of a static quark-antiquark pair should fall into approximate multiplets: $Sigma_g^{+\\prime} \\sim \\Pi_g$, $Sigma_g^- \\sim \\Pi_g^{\\prime} \\sim \\Delta_g$, $Sigma_u^- \\sim \\Pi_u$, and $Sigma_u^+ \\sim \\Pi_u^{\\prime} \\sim \\Delta_u$.","Each such static energy should approach the perturbative octet potential plus a nonperturbative constant, which ties the shape of the lattice-measured curves to the correlation lengths of the colour-electric and colour-magnetic gluon correlators, with $T_E < T_B$.","The $O(\\Lambda_{\\rm QCD})$ renormalon of the static potential cancels against the pole-mass self-energy, and the $O(\\Lambda_{\\rm QCD}^3 r^2)$ renormalon cancels against an ultraviolet renormalon of pNRQCD at next-to-leading order in the multipole expansion."],"supporting_citations":[{"why":"Supplies NRQCD itself, the theory from which pNRQCD is obtained by integrating out the soft scale.","marker":"[2]"},{"why":"Introduces pNRQCD, the effective theory this paper develops and matches in detail.","marker":"[11]"},{"why":"Provides the leading-logarithm contribution to the three-loop static potential that the present matching reproduces and extends.","marker":"[12]"},{"why":"Supplies the two-loop static potential, the input denoted $(u_0(r))_{two-loops}$ in the running of $alpha_{V_s}$.","marker":"[16]"},{"why":"Identifies the infrared divergences of the static Wilson loop that the pNRQCD matching is designed to absorb.","marker":"[9]"},{"why":"Provides the quenched lattice data on gluonic excitations between static quarks against which the degeneracy predictions and the $V_o + Lambda_H$ shape are compared.","marker":"[41]"},{"why":"Supplies the Voloshin-Leutwyler local-condensate treatment of nonperturbative corrections in the regime $mv^2 gg Lambda_{QCD}$.","marker":"[20]"},{"why":"Supplies the treatment of the pole-mass renormalon that cancels the $O(Lambda_{QCD})$ ambiguity of the static potential.","marker":"[48]"}],"fun_headline_variants":["In pNRQCD the static potential is a matching coefficient, not an observable","Integrating out the soft scale turns the static potential into a matching coefficient","The potential’s infrared divergence cancels an ultrasoft ultraviolet divergence in pNRQCD","The octet potential’s leading ultrasoft correction is order alpha_s^3, not alpha_s^2","Infrared renormalons in the static potential cancel within potential NRQCD"],"cache_read_input_tokens":0,"weakest_assumption_plain":"Everything computed here assumes there is a scale window in which the soft scale $mv$ is much larger than the matching scale $mu$, which in turn is much larger than $mv^2$ and $Lambda_{QCD}$, so that the matching between NRQCD and pNRQCD can be done order by order in perturbation theory.","fun_headline_variants_meta":{"raw":{"variants":["In pNRQCD the static potential is a matching coefficient, not an observable","Integrating out the soft scale turns the static potential into a matching coefficient","The potential’s infrared divergence cancels an ultrasoft ultraviolet divergence in pNRQCD","The octet potential’s leading ultrasoft correction is order alpha_s^3, not alpha_s^2","Infrared renormalons in the static potential cancel within potential NRQCD"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":3023,"prompt_tokens":884,"completion_tokens":2139,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":2041}},"tokens_in":884,"tokens_out":2139,"duration_ms":40107,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":2041,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T18:26:48.863099+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A calculation or lattice measurement showing that the coefficient of $ln(r mu)$ in the singlet potential differs from $C_A^3 alpha_s^4 / 12 pi$, or that the $mu$ dependence of the potential fails to cancel against ultrasoft contributions in a physical observable, would break the central claim; equally, short-distance lattice data on gluonic excitations that show no approach to a common octet-potential shape plus constant, and no $Sigma_g^{+prime} - Pi_g$ or $Sigma_u^+ - Pi_u^{prime} - Delta_u$ degeneracy, would undercut the static-limit predictions.","supporting_citations":[{"cited_title":"Phys. Lett. B 167 , 437 (1986); G.T. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. D 51 , 1125 (1995), Erratum ibid. D 55 , 5853","cited_arxiv_id":null,"evidence_quote":"Supplies NRQCD itself, the theory from which pNRQCD is obtained by integrating out the soft scale."},{"cited_title":"Nucl. Phys. B (Proc. Suppl.) 64, 428","cited_arxiv_id":null,"evidence_quote":"Introduces pNRQCD, the effective theory this paper develops and matches in detail."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the leading-logarithm contribution to the three-loop static potential that the present matching reproduces and extends."},{"cited_title":"Schr¨ oder, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the two-loop static potential, the input denoted $(u_0(r))_{two-loops}$ in the running of $alpha_{V_s}$."},{"cited_title":"Phys. Rev. D 17 , 2074","cited_arxiv_id":null,"evidence_quote":"Identifies the infrared divergences of the static Wilson loop that the pNRQCD matching is designed to absorb."},{"cited_title":"Gluon excitations of the static-quark potential","cited_arxiv_id":"hep-lat/9809015","evidence_quote":"Provides the quenched lattice data on gluonic excitations between static quarks against which the degeneracy predictions and the $V_o + Lambda_H$ shape are compared."},{"cited_title":"Nucl. Phys. B 154 , 365 (1979); H. Leutwyler, Phys. Lett. B 98 , 447","cited_arxiv_id":null,"evidence_quote":"Supplies the Voloshin-Leutwyler local-condensate treatment of nonperturbative corrections in the regime $mv^2 gg Lambda_{QCD}$."},{"cited_title":"Pineda, PhD thesis, Univ","cited_arxiv_id":null,"evidence_quote":"Supplies the treatment of the pole-mass renormalon that cancels the $O(Lambda_{QCD})$ ambiguity of the static potential."}],"review_version":1}