
3
The process to calculate ΛdQCD for given model is then,
ideally, as follows:
•Calculate the IRFP gauge couplings α∗
sand α∗
d.
•Determine Mby running the visible coupling up
from the measured value αs(MZ) = 0.11729 [38] to
αs(M) = α∗
s.
•Evolve the dark coupling down from αd(M) = α∗
d
until the scale ΛdQCD at which it becomes large
enough to trigger confinement. As a rough pertur-
bative condition, the Cornwall-Jackiw-Tomboulis
bound αs> π/4 for Nd= 3 [39] is used to define the
dark confinement scale through αd(ΛdQCD) = π/4.
So, for a given model with a nontrivial IRFP, the dark
confinement scale ΛdQCD can be uniquely determined.
The argument is then that a model with related IRFP
values α∗
sand α∗
dcould lead to compatible confinement
scales for visible and dark QCD of a similar order of
magnitude.
However, this framework relies upon a number of
simplifying assumptions that should be investigated in
greater detail. The first is that the heavy fields decouple
at the mass scale Mwithout any threshold corrections,
which was addressed in Ref. [25]; we summarise their ap-
proach in the next section. The second key assumption
is that the couplings run all the way to their IRFP val-
ues by the decoupling scale Mregardless of the initial
coupling values in the UV. However, this is not true in
general, and in a later section we analyse how the value
of ΛdQCD depends on the UV coupling values αUV
sand
αUV
d.
III. THRESHOLD CORRECTIONS
In MS-like mass-independent renormalization schemes,
the Appelquist-Carazzone decoupling theorem [40] does
not apply in its naive form; heavy fields can continue to
influence coupling constants and β-functions at energy
scales below their masses. In these schemes there also
arises the related issue of large logarithms when working
at energy scales much smaller than the mass Mof the
heavy fields.
To properly account for these problems, the decoupling
is treated explicitly by constructing an effective field the-
ory in which the heavy fields have been integrated out.
Consistency is ensured by matching the full theory onto
the effective field theory at a matching scale µ0, where
µ0∼Mto avoid large logarithms. This procedure leads
to a ‘consistency condition’ relating the coupling con-
stants in the full theory with the coupling constants in
the effective field theory [41,42],
αEFT
s(µ0) = ζ2
c(µ0, αs(µ0))αs(µ0),(4)
αEFT
d(µ0) = ζ2
d(µ0, αs(µ0))αd(µ0).(5)
The decoupling functions ζ2
cand ζ2
dhave been deter-
mined to three- and four-loop order in Refs. [43,44] for
the case of integrating out one heavy quark from QCD
with nfflavours. Since we are working with two-loop β-
functions, we only require the one-loop decoupling func-
tions; in Ref. [25], the results of Ref. [43] were adapted
for integrating out the full selection of heavy field content
at mass M, with the one-loop decoupling functions given
by
ζ2
c(µ, α(µ)) = 1 −αs(µ)
6π˜nclog µ2
M2,(6)
ζ2
d(µ, α(µ)) = 1 −αd(µ)
6π˜ndlog µ2
M2,(7)
where the coefficients ˜ncand ˜ndaccount for the degrees
of freedom of the relevant heavy fields,
˜nc=nfc,h + 3nfj+1
4nsc+ 3nsj,(8)
˜nd=nfd,h + 3nfj+1
4nsd+ 3nsj.(9)
We note that, for a given M, the values of the couplings
at low energies now depend on the choice of decoupling
scale µ0; this renormalization scale dependence is a non-
physical artifact of our perturbative analysis being trun-
cated at finite loop order. As a result, the procedure
of Bai and Schwaller (matching αs(MZ) to its experi-
mental value) no longer uniquely determines the value
of Mwhen we incorporate the threshold corrections to
the coupling constants. Instead, we obtain a relationship
between Mand µ0, and can solve for one of these values
if we specify the other.
In Ref. [25], values of ΛdQCD were determined for
a given model and value of M. By using the consis-
tency condition and again matching the strong coupling
to its measured value at MZ, they solved for the de-
coupling scale µ0, and then ran the dark coupling until
αd(ΛdQCD) = π/4. However, this process does not en-
sure that µ0is on the order of M. We take a different
approach to implementing threshold corrections, as is de-
tailed in the following section.
IV. DEPENDENCE ON INITIAL UV
COUPLINGS
As mentioned earlier, the framework in the previous
sections assumes that the couplings run to their IRFP
values by the decoupling scale µ0regardless of their ini-
tial values in the UV, αUV
sand αUV
d. However, this
is not the case, as illustrated in Fig. 1for the model
(nfc,h , nfd,l , nfd,h , nfj, nsc, nsd, nsj) = (3,3,3,3,3,2,0)
with IRFP (α∗
s, α∗
d) = (0.045,0.077), where we show the
running of a grid of couplings (αs, αd) from the UV scale
ΛUV = 1019 GeV down to a typical infrared scale ΛIR = 1
TeV. While the couplings evolve towards their fixed point
values, they do not precisely reach them, and so the low