It is a well-known result that classical light (which I take here to mean mixtures of coherent states) cannot produce sub-Poissonian photon-counting statistics, with a single beam of coherent light corresponding to a Poissonian photon-counting statistics (as discussed for example here), and other kinds of non-quantum light corresponding to super-Poissonian statistics.
However, I have never seen this fact proven formally. Usually, texts show how some common kinds of classical light, such as thermal light, result in super-Poissonian statistics, and how quantum states can produce sub-Poissonian ones, but they do not tackle the general case.
More specifically, consider a state which is a mixture of coherent states. This corresponds to a photon counting probability $P(n)$ of the form $$P(n)=\sum_\lambda p_\lambda P_\lambda(n),$$ with $\sum_\lambda p_\lambda =1$, and $P_\lambda(n)$ being the Poisson distribution with expected value $\lambda$: $$P_\lambda(n)\equiv e^{-\lambda}\frac{\lambda^n}{n!}.$$ A super-Poissonian distribution is characterised by the property that the variance is greater than the expected value, that is, $\sigma^2\ge\mu$. More precisely, in the considered case this means $$\sum_n(n-\mu)^2P(n)\ge \mu,\quad \mu\equiv\sum_n nP(n).$$ Can this property be shown in full generality, without making reference to specific types of light?
"which I take here to mean mixtures of coherent states" This is quite important clarification of the notion of "classical light". I think in purely classical EM theory, for individual classical electromagnetic waves produced by molecules one can get sub-Poissonian statistics, since reflection and transmission of a polarized wave on a beam splitter is anticorrelated.
– Ján Lalinský Jun 04 '19 at 17:50