non-local correlated hopping term expressed as MPO

Question

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 ?

Comments

Hi, do you have a question about this?

---

yes, I have edited the question just now.

---

Thanks,  but I don't have a good idea about this. What is Omega? What is F0? Could you provide a link to a reference since I haven't heard of correlation hopping before. Thanks -

---

No problem, you may find details in this paper : PRL 109, 203005, eqn(4)