Scalar conformal primary fields in the Brownian loop soup

The Brownian loop soup is a conformally invariant statistical ensemble of random loops in two dimensions characterized by an intensity $\lambda>0$, with central charge $c=2 \lambda$. Recent progress resulted in an analytic form for the four-point function of a class of scalar conformal primary "layering vertex operators" $\mathcal{O}_{\beta}$ with dimensions $(\Delta, \Delta)$, with $\Delta = \frac{\lambda}{10}(1-\cos\beta)$, that compute certain statistical properties of the model. The Virasoro conformal block expansion of the four-point function revealed the existence of a new set of operators with dimensions $(\Delta+ k/3, \Delta + k'/3)$, for all non-negative integers $k, k'$ satisfying $|k-k'| = 0$ mod 3. In this paper we introduce the edge counting field $\mathcal E(z)$ that counts the number of loop boundaries that pass close to the point $z$. We rigorously prove that the $n$-point functions of $\mathcal E$ are well defined and behave as expected for a conformal primary field with dimensions $(1/3, 1/3)$. We analytically compute the four-point function $\langle \mathcal{O}_{\beta}(z_1) \mathcal{O}_{-\beta}(z_2) \mathcal{E}(z_3) \mathcal{E}(z_4) \rangle$ and analyze its conformal block expansion. The operator product expansions of $\mathcal{E} \times \mathcal{E}$ and $\mathcal{E} \times \mathcal{O}_{\beta}$ produce higher-order edge operators with "charge" $\beta$ and dimensions $(\Delta + k/3, \Delta + k/3)$. Hence, we have explicitly identified all scalar primary operators among the new set mentioned above. We also re-compute the central charge by an independent method based on the operator product expansion and find agreement with previous methods.