Energy trapped Ising model

In this paper we have considered the 3D Ising model perturbed with the energy operator coupled with a non uniform harmonic potential acting as a trap, showing that this system satisfies the trap-size scaling behavior. Eventually, we have computed the correlators ⟨σ(z)σ(0)⟩, ⟨ϵ(z)ϵ(0)⟩ and ⟨σ(z)ϵ(0)⟩ near the critical point by means of conformal perturbation theory. Combining this result with Monte Carlo simulations, we have been able to estimate the OPE coefficients Cσσϵ, Cϵσσ and Cϵϵϵ, finding a good agreement with the values obtained in [1,2].