Transfer matrix and lattice dilatation operator for high-quality fixed points in tensor network renormalization group

Tensor network renormalization group maps study critical points of 2D lattice models like the Ising model by finding the fixed point of the RG map. In a prior work [Phys. Rev. X 15, 031047 (2025)]. we showed that by adding a rotation to the RG map, the Newton method could be implemented to find an extremely accurate fixed point. For a particular RG map (Gilt-TNR), we studied the spectrum of the Jacobian of the RG map at the fixed point and found good agreement between the eigenvalues corresponding to relevant and marginal operators and their known exact values. In this companion work, we use two further methods to extract many more scaling dimensions from this Newton method fixed point, and compare the numerical results with the predictions of conformal field theory (CFT). The first method is the well-known transfer matrix (TM), while the second method we refer to as the lattice dilatation operator (LDO). We introduce some extensions of these methods that also provide spins of the CFT operators, modulo an integer. With comparable computing resources, the TM and LDO methods perform equally well. The agreement for the scaling dimensions and spins is excellent up to Δ = 4¹/₈, and reasonably good up to 2 units higher. Some of the eigenvalues of the Jacobian of the RG map can come from perturbations associated with total derivative interactions and so are not universal. In some past studies [Phys. Rev. Res. 3, 023048 (2021); Phys. Rev. E 109, 034111 (2024)] such nonuniversal eigenvalues did not appear in the Jacobian. We explain this surprising result by showing that their RG map has the unusual property that the Jacobian is equivalent to the LDO operator.