code: def sI(qubit: Qubit, coefficient: complex = 1.0) -> Pauli:
      
      return Pauli.sigma(qubit, 'I', coefficient)
docs: Return the Pauli sigma_I (identity) operator. The qubit is irrelevant,
      but kept as an  argument for consistency