code: def collective_dr_squared( self ):
          
          return sum( np.square( sum( [ atom.dr for atom in self.atoms ] ) ) )
docs: Squared sum of total displacements for these atoms.
  
          Args:
              None
  
          Returns:
              (Float): The square of the summed total displacements for these atoms.