Strong parameter renormalization from optimum lattice model orbitals

We revisit the old problem of which is the best single particle basis to express a Hubbard-like lattice model. A rigorous variational solution of this problem leads to equations in which the answer depends in a self-consistent manner on the solution of the lattice model itself. Contrary to naive expectations, for arbitrary small interactions, the optimized orbitals differ from the non-interacting ones, leading also to substantial changes in the model parameters as shown analytically and in an explicit numerical solution for a simple double-well one-dimensional case. At strong coupling, we obtain the direct exchange interaction with a very large renormalization with important consequences for the explanation of ferromagnetism with model hamiltonians. Moreover, in the case of two atoms and two fermions we show that the optimization equations are closely related to reduced density matrix functional theory thus establishing an unsuspected correspondence between continuum and lattice approaches.