Consider a scalar field $\phi(x)$ described by an action $$ S= \frac12 \int d^4x\ \phi (F\phi) $$ where the kernel $F$ is a self-adjoint operator, in the sense that for any two complex-valued functions $\psi_1$ and $\psi_2$ supported on an open neighborhood $U$ of spacetime,
$$ \int d^4x\ [{\psi_1}^*(F\psi_2)-(F\psi_1)^*\psi_2] = 0\,.$$
Now consider a compact region $\Omega \subset U$, so that we can define a two-edged vector operator $\overleftrightarrow{f^\mu}$ corresponding to $F$ in the following manner.
$$ \int_\Omega d^4x\ [{\psi_1}^*(F\psi_2)-(F\psi_1)^*\psi_2] =: \int_{\partial\Omega} d\Sigma_\mu\ {\psi_1}^* \overleftrightarrow{f^\mu} \psi_2 \,,$$
where $d\Sigma_\mu$ is the outward directed area element of $\partial \Omega$.
Using this, we can define an inner product on the space of mode functions $\{u_i\}$ which satisfy the equation of motion for $\phi\,,$ namely that $Fu_i=0$. Given any complete Cauchy hypersurface $\Sigma$ for these equations of motion, define
$$ \langle u_i |u_j \rangle:= -i \int_\Sigma d\Sigma_\mu\ {u_1}^* \overleftrightarrow{f^\mu} u_2 $$
But, if $Fu_i=Fu_j=0$, then don't we quite clearly have that
$$ \langle u_i |u_j \rangle = 0 $$
identically?
For a contrary claim, read B. S. DeWitt, Phys. Rep. 19C, 292 (1975), section 1.1, just before Eq.(7). I am trying to understand what DeWitt really wanted to say and what I am missing.