Let d≥4, let Tr be the d-dimensional torus of side length r, and let Sr be the average field. For s∈R, set β=βc+swr, where
wr={a4(logr)−1/6r−2adr−d/2(d=4),(d>4),hr={b4(logr)1/4r−1bdr−d/4(d=4),(d>4),
with suitably chosen constants ad,bd>0. Let f(s)=∫Rx2dσs, where
dσs=∫Re−x4/4+sx2/2dxe−x4/4+sx2/2dx.
Non-gaussian limit and universal profile conjecture. As r→∞, hr−1Sr converges in distribution under ⟨⋅⟩βc+swrTr to dσs, with convergence of all moments. Moreover,
χTr(βc+swr)∼rdhr2f(s)=bd2f(s){(logr)1/2r2rd/2(d=4),(d>4),
r→∞limgTr(βc+swr)=3−(∫Rx2dσs)2∫Rx4dσs,
and, for x∈Tr,
τβc+swrTr(0,x)−τβc(0,x)∼hr2f(s)=bd2f(s){(logr)1/2r−2r−d/2(d=4),(d>4).
These predictions give a universal non-Gaussian finite-size scaling profile for the average field, susceptibility, plateau, and renormalised coupling constant. They are conjectured by universality from rigorous results for hierarchical ∣φ∣4 models, but remain unproved for the nearest-neighbour Ising model.