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:
E(a,b) = - cos(b-a)
So the quantity of interest is:
C(a,b,a',b') = E(a,b) + E(a,b') + E(a', b') - E(a',b)
= -cos(b-a) -cos(b'-a) -cos(b'-a') + cos(b-a')
If we let \alpha = b-a, \beta = b'-a, \gamma = b'-a', then we have:
C(\alpha, \beta, \gamma) = -cos(\alpha) - cos(\beta) -cos(\gamma) + cos(\alpha + \gamma - \beta)
For a minimum or maximum, the partial derivatives with respect to \alpha, \beta, \gamma must all be zero. This implies:
- sin(\alpha) - sin(\alpha + \gamma - \beta) = 0
- sin(\beta) + sin(\alpha + \gamma - \beta) = 0
- sin(\gamma) - sin(\alpha + \gamma - \beta) = 0
We can use some trigonometry:
- If sin(A) = sin(B), then either A=B or B = \pi - A
- If sin(A) = -sin(B), then either A=-B or B = \pi + A
So we have 8 possibilities:
- \alpha = \alpha + \gamma - \beta and \beta = -(\alpha+ \gamma - \beta) and \gamma = \alpha + \gamma - \beta
- \alpha = \alpha + \gamma - \beta and \beta = -(\alpha+ \gamma - \beta) and \gamma = \pi -(\alpha + \gamma - \beta)
- \alpha = \alpha + \gamma - \beta and \beta = \pi + (\alpha+ \gamma - \beta) and \gamma = \alpha + \gamma - \beta
- \alpha = \alpha + \gamma - \beta and \beta = \pi + (\alpha+ \gamma - \beta) and \gamma = \pi - (\alpha + \gamma - \beta)
- \alpha = \pi - (\alpha + \gamma - \beta) and \beta = -(\alpha+ \gamma - \beta) and \gamma = \alpha + \gamma - \beta
- \alpha = \pi - (\alpha + \gamma - \beta) and \beta = -(\alpha+ \gamma - \beta) and \gamma = \pi - (\alpha + \gamma - \beta)
- \alpha = \pi -(\alpha + \gamma - \beta) and \beta = \pi + (\alpha+ \gamma - \beta) and \gamma = \alpha + \gamma - \beta
- \alpha = \pi - (\alpha + \gamma - \beta) and \beta = \pi + (\alpha+ \gamma - \beta) and \gamma = \pi - (\alpha + \gamma - \beta)
These can be solved to yield the possibilities:
- \alpha = \beta = \gamma = 0
- Impossible
- Impossible
- \alpha = 0, \beta = \gamma = \pi
- Impossible
- \alpha = \gamma = 0, \beta = -\pi
- \alpha = \beta = \pi, \gamma = 0
- \alpha = \gamma = \frac{3\pi}{4}, \beta = \frac{5\pi}{4}
To find out whether these are minima or maxima, we plug them back into the expression for C(\alpha, \beta, \gamma)
- C(\alpha, \beta, \gamma) = -1-1-1+1 = -2
- Impossible
- Impossible
- C(\alpha, \beta, \gamma) = -1+1+1+1 = +2
- Impossible
- C(\alpha, \beta, \gamma) = -1+1-1-1 = -2
- C(\alpha, \beta, \gamma) = +1+1-1+1 = +2
- C(\alpha, \beta, \gamma) =+\sqrt{2}/2 +\sqrt{2}/2 +\sqrt{2}/2+\sqrt{2}/2 = 2\sqrt{2}
So, 2\sqrt{2} 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.