def METHOD_NAME(self):
    np.random.seed(seed=42)
    num_pairs = 200
    positions = np.random.random((2 * num_pairs, 3))
    self.system.part.add(pos=positions * self.system.box_l,
                         q=num_pairs * [-1, 1])
    self.minimize()