docs: def remove_atoms(self, indices):
          """Remove the atoms positioned at *indices*. The molecule
          containing the atom is removed as well.
  
          If you have a system of 10 water molecules (and 30 atoms), if
          you remove the atoms at indices 0, 1 and 29 you will remove
          the first and last water molecules.
  
          **Parameters**
  
          indices: np.ndarray((N,), dtype=int)
              Array of integers between 0 and System.n_atoms
  
          """
          mol_indices = self.atom_to_molecule_indices(indices)
          self.copy_from(self.sub(molecule_index=mol_indices))