increasing iPEPS truncation dimension $D$. Here the single-copy entanglement is calculated based on the one-site and two-site reduced density matrix, obtained by using the iTEBD method.[]{data-label="figXXZ_1"}](xxzspin1sca.eps "fig:"){width="0.7\linewidth"} ![(color online) Single-copy entanglement (a) $S_1(\Delta)$ and (b) $S_2(\Delta)$ as a function of the anisotropy parameter $\Delta$ for the 2D spin-$1$ XXZ model with increasing iPEPS truncation dimension $D$. Here the single-copy entanglement is calculated based on the one-site and two-site reduced density matrix, obtained by using the iTEBD method.[]{data-label="figXXZ_1"}](xxzspin1scb.eps "fig:"){width="0.7\linewidth"}

![(color online) (a) Sublattice components of the local magnetization $m$, (b) amplitude of the local order parameters $M_{xy}$ and $M_z$, (c) total antiferromagnetic order parameter $M_t$ and (d) GE per lattice site ${\mathcal E}(\Delta)$ as a function of anisotropy parameter $\Delta$ for the 2D spin-$1$ XXZ model. Values of the iPEPS truncation dimension $D$ are as indicated, with $D=9$ for (a).[]{data-label="figXXZ_2"}](xxzspin10a.eps "fig:"){width="0.7\linewidth"} ![(color online) (a) Sublattice components of the local magnetization $m$, (b) amplitude of the local order parameters $M_{xy}$ and $M_z$, (c) total antiferromagnetic order parameter $M_t$ and (d) GE per lattice site ${\mathcal E}(\Delta)$ as a function of anisotropy parameter $\Delta$ for the 2D spin-$1$ XXZ model. Values of the iPEPS truncation dimension $D$ are as indicated, with $D=9$ for 