def segregate(self, diff_mols = True):
    """
    Separate current Mol in a list of Mols of different molecules

    Parameters
    ----------
    diff_mols : bool (optional)
        If all molecules are of the same length, turn this off for a boost in
        speed.

    """
    molecules = []  # list of molecules
    remaining = self.copy()

    natoms_current = 0
    while len(remaining) > 0:
        molecule = remaining.select(0, natoms = natoms_current)
        # diff_mols = True means use slow algorithm
        if diff_mols:
            natoms_current = 0
        else:
            natoms_current = len(molecule)
        molecules.append(molecule)
        for atom in molecule:
            remaining.remove(atom)
    return molecules