We consider long-range Bernoulli bond percolation on the $d$-dimensional

hierarchical lattice in which each pair of points $x$ and $y$ are connected by

an edge with probability $1-\exp(-\beta\|x-y\|^{-d-\alpha})$, where

$0<\alpha<d$ is fixed and $\beta \geq 0$ is a parameter. We study the volume of

clusters in this model at its critical point $\beta=\beta_c$, proving precise

estimates on the moments of all orders of the volume of the cluster of the

origin inside a box. We apply these estimates to prove up-to-constants

estimates on the tail of the volume of the cluster of the origin, denoted $K$,

at criticality, namely \[ \mathbb{P}_{\beta_c}(|K|\geq n) \asymp \begin{cases}

n^{-(d-\alpha)/(d+\alpha)} & d < 3\alpha\\ n^{-1/2}(\log n)^{1/4} & d=3\alpha

\\ n^{-1/2} & d>3\alpha. \end{cases} \] In particular, we compute the critical

exponent $\delta$ to be $(d+\alpha)/(d-\alpha)$ when $d$ is below the

upper-critical dimension $d_c=3\alpha$ and establish the precise order of

polylogarithmic corrections to scaling at the upper-critical dimension itself.

Interestingly, we find that these polylogarithmic corrections are not those

predicted to hold for nearest-neighbour percolation on $\mathbb{Z}^6$ by Essam,

Gaunt, and Guttmann (J. Phys. A 1978). Our work also lays the foundations for

the study of the scaling limit of the model: In the high-dimensional case $d

\geq 3\alpha$ we prove that the sized-biased distribution of the volume of the

cluster of the origin inside a box converges under suitable normalization to a

chi-squared random variable, while in the low-dimensional case $d<3\alpha$ we

prove that the suitably normalized decreasing list of cluster sizes in a box is

tight in $\ell^p\setminus \{0\}$ if and only if $p>2d/(d+\alpha)$.