I don't know how to prove it for an arbitrary QM system, but in the particular case of anti-correlated spin-1/2 particles, I can prove it.
In that case, we have:
[itex]E(a,b) = - cos(b-a)[/itex]
So the quantity of interest is:
[itex]C(a,b,a',b') = E(a,b) + E(a,b') + E(a', b') - E(a',b)[/itex]
[itex]= -cos(b-a) -cos(b'-a) -cos(b'-a') + cos(b-a')[/itex]
If we let [itex]\alpha = b-a, \beta = b'-a, \gamma = b'-a'[/itex], then we have:
[itex]C(\alpha, \beta, \gamma) = -cos(\alpha) - cos(\beta) -cos(\gamma) + cos(\alpha + \gamma - \beta)[/itex]
For a minimum or maximum, the partial derivatives with respect to [itex]\alpha, \beta, \gamma[/itex] must all be zero. This implies:
- [itex]sin(\alpha) - sin(\alpha + \gamma - \beta) = 0[/itex]
- [itex]sin(\beta) + sin(\alpha + \gamma - \beta) = 0[/itex]
- [itex]sin(\gamma) - sin(\alpha + \gamma - \beta) = 0[/itex]
We can use some trigonometry:
- If [itex]sin(A) = sin(B)[/itex], then either [itex]A=B[/itex] or [itex]B = \pi - A[/itex]
- If [itex]sin(A) = -sin(B)[/itex], then either [itex]A=-B[/itex] or [itex]B = \pi + A[/itex]
So we have 8 possibilities:
- [itex]\alpha = \alpha + \gamma - \beta[/itex] and [itex]\beta = -(\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \alpha + \gamma - \beta[/itex]
- [itex]\alpha = \alpha + \gamma - \beta[/itex] and [itex]\beta = -(\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \pi -(\alpha + \gamma - \beta)[/itex]
- [itex]\alpha = \alpha + \gamma - \beta[/itex] and [itex]\beta = \pi + (\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \alpha + \gamma - \beta[/itex]
- [itex]\alpha = \alpha + \gamma - \beta[/itex] and [itex]\beta = \pi + (\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \pi - (\alpha + \gamma - \beta)[/itex]
- [itex]\alpha = \pi - (\alpha + \gamma - \beta)[/itex] and [itex]\beta = -(\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \alpha + \gamma - \beta[/itex]
- [itex]\alpha = \pi - (\alpha + \gamma - \beta)[/itex] and [itex]\beta = -(\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \pi - (\alpha + \gamma - \beta)[/itex]
- [itex]\alpha = \pi -(\alpha + \gamma - \beta)[/itex] and [itex]\beta = \pi + (\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \alpha + \gamma - \beta[/itex]
- [itex]\alpha = \pi - (\alpha + \gamma - \beta)[/itex] and [itex]\beta = \pi + (\alpha+ \gamma - \beta)[/itex] and [itex]\gamma = \pi - (\alpha + \gamma - \beta)[/itex]
These can be solved to yield the possibilities:
- [itex]\alpha = \beta = \gamma = 0[/itex]
- Impossible
- Impossible
- [itex]\alpha = 0, \beta = \gamma = \pi[/itex]
- Impossible
- [itex]\alpha = \gamma = 0, \beta = -\pi[/itex]
- [itex]\alpha = \beta = \pi, \gamma = 0[/itex]
- [itex]\alpha = \gamma = \frac{3\pi}{4}, \beta = \frac{5\pi}{4}[/itex]
To find out whether these are minima or maxima, we plug them back into the expression for [itex]C(\alpha, \beta, \gamma)[/itex]
- [itex]C(\alpha, \beta, \gamma) = -1-1-1+1 = -2[/itex]
- Impossible
- Impossible
- [itex]C(\alpha, \beta, \gamma) = -1+1+1+1 = +2[/itex]
- Impossible
- [itex]C(\alpha, \beta, \gamma) = -1+1-1-1 = -2[/itex]
- [itex]C(\alpha, \beta, \gamma) = +1+1-1+1 = +2[/itex]
- [itex]C(\alpha, \beta, \gamma) =+\sqrt{2}/2 +\sqrt{2}/2 +\sqrt{2}/2+\sqrt{2}/2 = 2\sqrt{2}[/itex]
So, [itex]2\sqrt{2}[/itex] is the biggest value that you can get with anti-correlated spin-1/2 particles. It's not at all obvious to me why that is the best you can possibly do with any QM system.