For $a > 0$ and $n \in \mathbb Z_+$, consider the oscillatory integral

$$\int_{0}^1 f(x) f(ax) \dots f(a^n x) \, dx,$$

where $f$ is an integrable function on $[0, 1]$, which we extend by periodicity to all of $\mathbb R$.

By a Riemann-Lebesgue type argument, one can show that we have

$$\int_{0}^1 f(x) f(ax) \dots f(a^n x) \, dx \to \left ( \int f(x) \, dx \right )^{n+1}$$

as $a \to \infty$. In particular the limit is $0$ for mean zero functions. We can see this as a “zeroth order” result, in the sense that it says the integral is $o(1)$ as $a \to \infty$. 


Assume now additional smoothness on $f$ such as Lipschitz or $C^1$. 

**Question:** Can we obtain sharper asymptotics on the integral for mean zero functions? 
For instance, do we have

$$\int_{0}^1 f(x) f(ax) \dots f(a^n x) \, dx = O(\frac{1}{a^n})$$

as $a \to \infty$, with the implied constants depending only on some Lipschitz or $C^k$ norm of $f$?