Since not all particles contribute to all of constraints, J σ σ is a block matrix and can be solved individually to block-unit of the matrix. In other words, J σ σ can be solved individually for each molecule.