For a Hubbard model with correlation hopping, the hopping term is

$$

H_0=-J\sum_{<i,j>}a_i^{\dagger} F_0(\hat{n}_i-\hat{n}_j)a_j

$$

with

$$

F_0(\hat{n}_i-\hat{n}_j)=J_0(\Omega(\hat{n}_i-\hat{n}_j))

$$

here, @@J_0@@ is first-order Bessel function. SInce @@F_0@@ is non-linear, it is impossible to express @@ F_0 @@ as a product of two local operators on sites i and j, i.e.

$$

F_0(\hat{n}_i-\hat{n}_j) ?= O_i \cdot O_j

$$

Then in this case, how can we perform DMRG as usual using ITensor ?