where the index \(j(i)\)  runs over all the particles that are nearest neighbors of particle \(i\)  (thus the distance is less than 1.6, the first minimum
of the total pair correlation function), \(\vec{r_j}(\gamma )\)  denotes the position of particle \(j\)  after the minimization of the box that is sheared by \(\gamma \) , and \(N_i\)  is the number of nearest neighbors of particle \(i\) . The matrix \({\bf F}_i\)  is chosen such that it minimizes \(D_i^2\)  and it can be calculated based on the spatial coordinates of the corresponding particles [1]}.
