Let $\mathbb{F}_2=\{0,1\}$ be the field with two elements, and let $u:\mathbb{F}_2^n\rightarrow \mathbb{F}_2$ be a function.
Is it possible that $$ \sum_{x \in \mathbb{F}_2^n}(-1)^{u(x)+u(x+a)}= 0, $$ for every $a \neq 0$ in $\mathbb{F}_2^n$?
If so, can we characterise the functions satisfying property?
Comment: I treat the sum here as natural (or real) number, not as an element of $\mathbb{F}_2$.
This question arose in the context of this question. (Trying to derive a lower bound on the approximate multiplicativity of a Boolean function).
If $u$ is linear, i.e. $u(x+y)=u(x)+u(y)$, then the sum above is never zero.
Edit:
Jyrki Lahtonen gave a nice example for all even $n$. It turns out that there are no such functions for odd $n$.
I think the following idea I borrowed from the theory of so called bent functions (that Claude Carlet's group is actively researching) gives a construction of an example of such a function $u$. Assuming that $n$ is an even number.
Write $n=2k$, $V=\Bbb{F}_2^k$, and $\langle\ ,\ \rangle:V\times V\to \Bbb{F}_2$ the usual inner product. We identify $\Bbb{F}_2^n$ with $V\times V$ and use $u=\langle\ ,\ \rangle$. Let $(a,b)$ be a non.zero vector of $V\times V$. Let $(x,y)$ be an arbitrary vector. Bilinearity means that $$ \begin{aligned} u((x,y))+u((x,y)+(a,b))&=\langle x,y\rangle+\langle x+a,y+b\rangle\\ &=\langle a,y\rangle+\langle x,b\rangle+\langle a,b\rangle. \end{aligned} $$
Therefore the sum $$ \begin{aligned} \sum_{(x,y)\in\Bbb{F}_2^n}(-1)^{u((x,y))+u((x,y)+(a,b))}&= \sum_{x\in V}\sum_{y\in V}(-1)^{\langle a,y\rangle+\langle x,b\rangle+\langle a,b\rangle}\\ &=(-1)^{\langle a,b\rangle}\sum_{y\in V}(-1)^{\langle a,y\rangle}\sum_{x\in V}(-1)^{\langle x,b\rangle}. \end{aligned} $$ The assumption is that at least one of $a,b$ is non-zero. Therefore at least one of the sums that appear as factors in the final formula vanishes by the usual orthogonality of characters (see below). Hence so does the entire sum.
In other words, if $n=2$, use $u(x_1,x_2)=x_1x_2$. If $n=4$, use $u(x_1,x_2,x_3,x_4)=x_1x_3+x_2x_4$ et cetera.
Orthogonality. Assume that $a\in V$ is non-zero. Then there exists an element $\epsilon\in V$ such that $\langle a,\epsilon\rangle=1$. We fix one such. Therefore $$ \begin{aligned} S(a)&:=\sum_{x\in V}(-1)^{\langle x,a\rangle}\\ &=\sum_{x+\epsilon\in V}(-1)^{\langle x+\epsilon,a\rangle}\\ &=\sum_{x\in V}(-1)^{\langle x+\epsilon,a\rangle}&\text{because $x$ ranges over $V$ while $x+\epsilon$ does}\\ &=(-1)^{\langle\epsilon,a\rangle}\left(\sum_{x\in V}(-1)^{\langle x,a\rangle}\right)\\ &=-S(a). \end{aligned} $$ Hence $S(a)=0$ as claimed.