Schematic of the nonlinear QWZ model. The model has two sublattices (black circles) at each lattice point encircled by the blue ellipse. The green lines represent the linear couplings. We use the notation i ( x , y ) to represent the state variable at each sublattice, where ( x , y ) is the location of the representative point of each lattice point denoted by the red cross.