NSPT for $O(N)$ non-linear sigma model: the larger $N$ the better

The $O(N)$ non-linear sigma model (NLSM) is an example of field theory on a target space with nontrivial geometry. One interesting feature of NLSM is asymptotic freedom, which makes perturbative calculations interesting. Given the successes in Lattice Gauge Theories, Numerical Stochastic Perturbation Theory (NSPT) is a natural candidate for performing high-order computations also in the case of NLSM. However, in low-dimensional systems NSPT is known to display statistical fluctuations substantially increasing for increasing orders. In this work, we explore how for $O(N)$ NLSM this behaviour is strongly dependent on $N$. As largely expected on general grounds, the larger is $N$, the larger is the order at which a NSPT computation can be effectively performed.