Dear all,

I am trying to generalize the TEBD method mentioned here to a case involving fermions. When changing the following part

```
gates = ITensor[]
for j=1:N-1
s1 = s[j]
s2 = s[j+1]
hj = op("Sz",s1) * op("Sz",s2) +
1/2 * op("S+",s1) * op("S-",s2) +
1/2 * op("S-",s1) * op("S+",s2)
Gj = exp(-1.0im * tau/2 * hj)
push!(gates,Gj)
end
```

to include some fermion operators, do we need to manually add JW strings? I am not sure whether these `gates`

already implement JW string like `OpSum`

or not.

Best,

Junsen