Thanks.
I noticed one thing more, the way I presented above is non coherent since it does not correspond to the average. So I tried to use the closure relationship and it gave that the result of the covariance operator can be either {1,0}, or {-1,0}.
So that studying the operator in fact gives insight to which measurement result it can be (it's in some way more precise than to look at the average).
In fact I noticed a big problem : since the operator is not symmetric, there are left and right eigenproblems. In order to get rid of that, one can symmetrize the operator, leading to
[tex]C=A\otimes B-\frac{1}{2}(A|\Psi\rangle\langle\Psi|B+B|\Psi\rangle\langle\Psi|A)[/tex]
to consider the simplest case, I chose A=B=diag(1,-1). Letting Psi=(a,b,c,d), I got :
[tex]C=\left(\begin{array}{cccc}1-a^2& 0 & 0 & ad\\0&-1+b^2&-bc&0\\0 &-bc&-1+c^2&0\\ad&0&0&1-d^2\end{array}\right)[/tex]
THe eigenproblem is [tex]C|\Psi\rangle=\lambda|\Psi\rangle[/tex] hence :
1) [tex]a-a^3+ad^2=\lambda a[/tex]
2) [tex]-b+b^3-bc^2=\lambda b[/tex]
3) [tex]-b^2c-c+c^3=\lambda c[/tex]
4) [tex]a^2d+d-d^3=\lambda d[/tex]
To solve this problem i supposed a,b,c,d non zero (if they are 0 this solve the equation), then it comes
1+4 implies [tex]\lambda=1[/tex]
2+3 implies [tex]\lambda=-1[/tex]
hence not all a,b,c,d can be 0. If i suppose b=c=0 I get a=d or a=-d
if a=d=0, then b=c or b=-c
Considering equ. 1 and 3 (hence b=d=0) I get
[tex]1-a^2=\lambda[/tex] and [tex]-1+c^2=\lambda[/tex] leading to [tex]a^2+c^2=2[/tex]
and [tex]\lambda=1-a^2\in[-1;1][/tex]
The same comes out for equ. 2 and 4 (hence a=c=0)
Thus the spectrum is composed of 2 discrete eigenvalues 1 and -1 and a continuous spectrum.
We see however that the eigenvectors are not orthogonal.
Writing the closure relationship : [tex]\int\sum_{\lambda}|\lambda\rangle\langle\lambda|=\mathbb{1}[/tex] leads to the fact that the discrete eigenvalues are not considered, because of equation of the form [tex]p^2+q^2+\underbrace{\langle2cos(\theta)^2\rangle}_{=1}=1[/tex] where p and q are the component of the eigenvector corresponding to the value -1.
I get then for the average of the C operator in the singlet state :
[tex]\langle C\rangle=\frac{1}{\pi}\int_0^{2\pi}(1-2cos(\theta)^2)sin(\theta)^2=1/2[/tex]
This is but different than the average computed directly since it simply gives -1.
There is already a sign difference but I don't know if it's a mistake ??