Minimax and adaptive estimation of the Wigner function in quantum homodyne tomography with noisy data

Abstract: In quantum optics, the quantum state of a light beam is represented through the Wigner function, a density on $\mathbb R^2$ which may take negative values but must respect intrinsic positivity constraints imposed by quantum physics. In the framework of noisy quantum homodyne tomography with efficiency parameter $1/2 < \eta \leq 1$, we study the theoretical performance of a kernel estimator of the Wigner function. We prove that it is minimax efficient, up to a logarithmic factor in the sample size, for the $\mathbb L_\infty$-risk over a class of infinitely differentiable. We compute also the lower bound for the $\mathbb L_2$-risk. We construct adaptive estimator, i.e. which does not depend on the smoothness parameters, and prove that it attains the minimax rates for the corresponding smoothness class functions. Finite sample behaviour of our adaptive procedure are explored through numerical experiments.