“…For the SSH4 model, we can choose the unit cell basis {A 1 , B 1 , A 2 , B 2 }, so the total mean chiral displacement operator takes the form Γm = diag(..., 1, −1, 1, −1, 2, −2, 2, −2, ...). At higher dimension, in order to use the mean chiral displacement to measure the winding number, we need to choose an orthogonal and complete basis of a given sub-lattice [77]. In our case, we prepare the initial state at two orthogonal states (1,0,0,0) and (0,0,1,0), and measure the mean chiral displacement C 1 (t) and C 3 (t) respectively.…”