code: def concentration(self, pM=False):
          
          concentr = (self.num_particles / NA) / self.box.volume_L
          if pM:
              concentr *= 1e12
          return concentr
docs: Return the concentration (in Moles) of the particles in the box.