Using operator covariance to disentangle scaling dimensions in lattice models

In critical lattice models, distance ($r$) dependent correlation functions contain power laws $r^{-2\Delta}$ governed by scaling dimensions $\Delta$ of an underlying continuum field theory. In Monte Carlo simulations, the leading dimensions can be extracted by data fitting, which is difficult when two or more powers contribute significantly. Here a method utilizing covariance between multiple lattice operators is developed where the r dependent eigenvalues of the covariance matrix reflect scaling dimensions of individual field operators. This disentangling is demonstrated explicitly for conformal field theories. The scheme is first tested on the critical point of the 2D Ising model, where the two primary scaling dimensions and their respective two lowest descendant dimensions are extracted. The 3D Ising model is studied next, revealing the two relevant primaries and their lowest descendants to high precision. The 2D tricritical Ising point is studied with the Blume-Capel model. Here the scaling dimensions of all three symmetric primary operators are successfully isolated along with the leading descendants. The eigenvectors are also studied and give useful information on the boundary between the ordered and disordered phases in the neighborhood of the tricritical point. Finally, the crossover from regular to tricritical Ising scaling is investigated on several points on the phase boundary of the Blume-Capel model away from its tricritical point. The scaling of the eigenvalues corresponding to tricritical descendant operators are found to be remarkably stable even far from the tricritical point. The covariance method represents a simple extension of standard analysis of correlation functions and can significantly enhance the utility of Monte Carlo simulations and other computational methods in studies of criticality, in particular conformal critical points.