if not op:
        return spin_operators
    elif op == "x":
        return sx
    elif op == "y":