I know how to show that if $(f,\mu)$ is weakly mixing then so is $(f^{k},\mu), \forall k \geq 2$. I was trying to proof the converse. I don't even know if this result is true and, at the end, seems to be a real analysis problem:
Fix $A$ and $B$ measurable sets. Define $a_{n}= |\int U^{n}_{f}(\mathbb{1}_{A}) \mathbb{1}_{B}d\mu-\mu(A)\mu(B)|$, where $U_{f}$ is the Koopman operator. By hypothesis, we have
$lim_{n\to \infty} \frac{1}{n}\sum_{j=1}^{n}a_{jk}=0, \forall k\geq 2$.
Clearly, the question is if $lim_{n\to \infty} \frac{1}{n}\sum_{j=1}^{n}a_{j}$ is equal to zero.
Maybe the last assertion is not true. Is there a known counterexample?
Surprisingly the converse does hold. In fact you don't even need $T^k$ to be weakly mixing for all $k\geq 2$. A single one, such as $T^2$, will do. This is because a transformation $S$ being weak mixing is equivalent (See Einsiedler and Ward, theorem 2.36) to the following property:
Now assume this property holds for $S=T^2$ and consider the pair $A, B$ and the pair $A,T^{-1} B$. For the first pair, we get a zero-upper-density $E_0$ such that $$\lim\limits_{n\notin E_0, n\rightarrow \infty} \mu(A\cap T^{-2n}B) = \mu(A)\mu(B). $$ From the second pair, we get a zero-upper-density $E_1$ such that $$\lim\limits_{n\notin E_1, n\rightarrow \infty} \mu(A\cap T^{-2n-1}B) = \mu(A)\mu(B). $$
Together we have a zero-upper-density $E:=E_0\cup E_1$ which satisfies
$$\lim\limits_{n\notin E, n\rightarrow \infty} \mu(A\cap T^{-n}B) = \mu(A)\mu(B). $$
As an aside, this means that @mathworker12's interesting counterexample couldn't possibly arise as the quantities $|\int U^{n}_{f}(\mathbb{1}_{A}) \mathbb{1}_{B}d\mu-\mu(A)\mu(B)|$, for any measure preserving $f$.