Return the bitphase flip kraus operators

def bitphase_flip_operators(p):
    """
    Return the bitphase flip kraus operators
    """
    k0 = np.sqrt(1 - p) * I
    k1 = np.sqrt(p) * Y
    return k0, k1