docs: def select_atoms(indices):
      '''Select atoms by their indices.
  
      You can select the first 3 atoms as follows::
  
        select_atoms([0, 1, 2])
      
      Return the current selection dictionary.
  
      '''
      rep = current_representation()
      rep.select({'atoms': Selection(indices, current_system().n_atoms)})
      return rep.selection_state