Interacting Wannier functions in a superlattice: An iterative method

We present a method to compute interacting Wannier functions in an optical superlattice for cold atom gas experiments. We start from a noninteracting Wannier function obtained from the band projected position operator method and perform Trotter method to compute the interacting Wannier function. Then, we use this function to compute the on-site interaction strength in a superlattice, compare the result to a measurement, and find good agreement between simulation and data.