1 Introduction
Many approaches to critical phenomena obtain results in continuous space dimension,
although physically relevant dimensions are integer. Most notable is the perturbative renor-
malization group in d= 4−dimensions [1–4]. This is not merely a technical issue: quantities
as functions of real dcan clarify features that are harder to see at discrete values. E.g., one
can follow the topology of the renormalization-group (RG) flow as a function of dimension
and find instances where the universality class changes at non-integer values. This proved
particularly useful for systems with long-range interactions [5–7] or disorder [8–13].
The recent very precise numerical conformal bootstrap [14–16] has been formulated in
continuous dimension [17, 18], in particular for the Ising model in its whole range 4> d ≥
2[19–21]. The interest lies in understanding how the strongly interacting Ising conformal
field theory connects to a free scalar in d= 4 and to the integrable fully-solvable model
in d= 2 [22, 23]. Analytic bootstrap approaches which use the dimension as a tunable
parameter were also developed [24–32]. Initially, the non-unitarity of the theory in non-
integer dimensions [33] was thought to hamper the numerical methods involving positive
quantities. These concerns have been overcome by de facto never observing problems for the
quantities of interest, as explained later.
In this paper, we extend the numerical approach of Ref. [20] using a single correlator,
the SDPB [34] routine for determining the unitarity domain, and the Extremal Functional
Method [35, 36] for solving the bootstrap equations. We obtain improved results for the
scaling dimensions in 4> d ≥3by a denser scanning of the unitary region near the Ising
point, i.e., the kink. The latter gets parametrically sharper as dapproaches 4, allowing for
its better identification. The conformal spectrum in dimensions 4> d ≥2.6has also been
obtained in Ref. [21] via the advanced navigator bootstrap technique [37]. We use these very
precise results in combination with ours to obtain a consistent description of the low-lying
spectrum.
The achieved precision allows us to perform a detailed comparison with state-of-the-art
epsilon expansion in two regimes: for dclose to 4, the series is directly compared to bootstrap
data, using the necessary finer scale for the latter; for intermediate values between 4and
3(included), the divergent perturbative series is resummed using well-established methods
involving the Borel transform [38–41].
The analysis is done on the dimensions of the conformal fields σ, , 0, corresponding to
spin, energy and subleading energy. They determine the critical exponents η, ν, ω. The
precision of our bootstrap data is summarized by the (mostly) d-independent value of the
relative error Err(γ)/γ =O(10−3)for the anomalous dimensions γof the conformal fields σ
and . As the anomalous dimensions are very small for d≈4, the precision for the conformal
dimensions ∆σ,∆is even higher in this region. Regarding the subleading energy, the relative
error Err(∆0)/∆0stays at three digits, as explained later. Some of the structure constants
are determined with a higher O(10−4)accuracy.
We compare our data with recent results of the analytic bootstrap [27–32], Monte Carlo
simulations [42–44] and the non-perturbative RG [45, 46]. We find that the data by all
methods agree very well. This is rather rewarding given the achieved precision. Besides
1