def split_xyz(xyz_file: bytes) -> list[bytes]:
    """Split an xyz file into individual conformers."""
    lines = xyz_file.splitlines()
    structures = []
    while True:
        if len(lines) == 0:
            break

        # removed one deck
        natoms = lines.pop(0)
        n = int(natoms.decode())
        comment = lines.pop(0)
        geom = []
        for _ in range(n):
            geom += [lines.pop(0)]

        deck = b"\n".join([natoms, comment] + geom)
        structures += [deck]
    return structures