We show that the Z_{N} Berry phase (Berry phase quantized into 2π/N) provides a useful tool to characterize symmetry protected topological phases with correlation that can be directly computed through numerics of a relatively small system size. The Z_{N} Berry phase is defined in a N-1-dimensional parameter space of local gauge twists, which we call the "synthetic Brillouin zone," and an appropriate choice of an integration path consistent with the symmetry of the system ensures exact quantization of the Berry phase. We demonstrate the usefulness of the Z_{N} Berry phase by studying two 1D models of bosons, SU(3) and SU(4) Affleck-Kennedy-Lieb-Tasaki models, where topological phase transitions are captured by Z_{3} and Z_{4} Berry phases, respectively. We find that the exact quantization of the Z_{N} Berry phase at the topological transitions arises from a gapless band structure (e.g., Dirac cones or nodal lines) in the synthetic Brillouin zone.
The local
ZN quantized Berry phase for the
SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local
Z2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning a particular bond, and is effective as a tool to characterize and classify the various symmetry-protected topological phases of one-dimensional
SU(N) spin systems. We extend the path-integral quantum Monte Carlo method for the
Z2 Berry phase in order to calculate the
ZN Berry phase numerically. We demonstrate our method by calculating the
Z4 Berry phase for the bond-alternating
SU(4) antiferromagnetic Heisenberg chain, which is represented by the Young diagram of four columns.