 Solution

Let ${\psi }_0$ be the normalized stationary profile. Multiplying its stationary Gross–Pitaevskii equation by ${\psi }_0^{\ast }$, integrating and integrating the kinetic term by parts yieldsThis is not simply $H/N$: the quartic contribution is twice its contribution to that energy per particle. For an infinite periodic system one uses one cell and the corresponding cell particle number, or equivalently average densities.

At fixed lattice spacing, mean density $n$, coupling and temperature, the commonly used condensate compressibility in an optical lattice isIf compressibility means the thermodynamic volume response instead, the Gibbs-Duhem equation $dp=nd\mu$ at zero temperature gives ${\kappa }_T=[n^2(\partial \mu /\partial n{)}_s{]}^{-1}=\kappa /n$. Stating the convention removes this factor-of-density ambiguity; the lattice-depth derivative has the same sign at fixed $n$.

For the stable repulsive ground-band Bose-Einstein condensate, the intended trend is $\partial \kappa /\partial s<0$ for increasing nonzero lattice depth, with mean density held fixed. Deeper wells localize the field, enhancing the repulsive interaction response. To quantify this, in the weak-interaction band approximation write ${\psi }_0=\sqrt{n}{f}_s$, with the cell average of $|f_s{|}^2$ equal to one. Perturbation in $Un$ givesIn deep wells the localized orbital is approximately a normalized Gaussian of width $a_{\mathrm{ho}}=d/(\pi s^{1/4})$. Neglecting interwell overlap, $I\simeq d/(\sqrt{2\pi }{a}_{\mathrm{ho}})=\sqrt{\pi /2}{s}^{1/4}$, hence $\partial \kappa /\partial s\simeq -\kappa /(4s)<0$. This explicitly establishes the sign in that controlled regime; the repulsive ground-band calculation describes the broader trend.

The sign statement needs the repulsive, stable ground-band interpretation, not merely the literal assumption $U\ne 0$. It also need not be a strictly nonzero derivative at $s=0$: at weak interaction $f_s=1+(s/8)\cos (2\pi x/d)+O(s^2)$, with normalization corrected at second order, gives $I=1+s^2/32+O(s^3)$ and a zero initial slope. A leading density-dominated Thomas–Fermi approximation for a condensate with density positive everywhere gives $\mu =Un+s{E}_R/2$ and no depth dependence of this response at that approximation's order. These limits do not justify a strict universal sign for every unspecified interaction and parameter regime.