  We consider ultracold bosons in a 2D square optical lattice described by the
Bose-Hubbard model. In addition, an external time-dependent sinusoidal force is
applied to the system, which shakes the lattice along one of the diagonals. The
effect of the shaking is to renormalize the nearest-neighbor hopping
coefficients, which can be arbitrarily reduced, can vanish, or can even change
sign, depending on the shaking parameter. It is therefore necessary to account
for higher-order hopping terms, which are renormalized differently by the
shaking, and introduce anisotropy into the problem. We show that the
competition between these different hopping terms leads to finite-momentum
condensates, with a momentum that may be tuned via the strength of the shaking.
We calculate the boundaries between the Mott-insulator and the different
superfluid phases, and present the time-of-flight images expected to be
observed experimentally. Our results open up new possibilities for the
realization of bosonic analogs of the FFLO phase describing inhomogeneous
superconductivity.
