# Tomographic measurement of nonclassical radiation states

## I INTRODUCTION

The concept of nonclassical states of light has drawn much attention in quantum optics . The customary definition of nonclassicality is given in terms of the Glauber-Sudarshan $`P`$-function: a nonclassical state does not admit a regular positive $`P`$-function representation, namely, it cannot be written as a statistical mixture of coherent states. Such states produce effects that have no classical analogue. These kinds of states are of fundamental relevance not only for the demonstration of the inadequacy of classical description, but also for applications, e.g., in the realms of information transmission and interferometric measurements .

In this paper we are interested in testing the nonclassicality of a quantum state by means of an operational criterium, which is based on a set of quantities that can be measured experimentally with some given level of confidence, even in the presence of loss, noise, and less-than-unity quantum efficiency. The positivity of the $`P`$-function itself cannot be adopted as a test, since there is no method available to measure it. The $`P`$-function is a Fourier transform on the complex plane of the generating function for the normal-ordered moments; hence, in principle, it could be recovered by measuring all the quadrature components of the field, and subsequently performing an (deconvolved) inverse Radon transform . Currently, there is a well-established quantitative method for such a universal homodyne measurement, and it is usually referred to as quantum homodyne tomography (see Ref.  for a review). However, as proven in Ref. , only the generalized Wigner functions of order $`s<1\eta ^1`$ can be measured, $`\eta `$ being the quantum efficiency of homodyne detection. Hence, through this technique, all functions from $`s=1`$ to $`s=0`$ cannot be recovered, i.e., we cannot obtain the $`P`$-function and all its smoothed convolutions up to the customary Wigner function. For the same reason, the nonclassicality parameter proposed by Lee , namely, the maximum $`s`$-parameter that provides a positive distribution, cannot be experimentally measured.

Among the many manifestations of nonclassical effects, one finds squeezing, antibunching, even-odd oscillations in the photon-number probability, and negativity of the Wigner function . Any of these features alone, however, does not represent the univocal criterium we are looking for. Neither squeezing nor antibunching provides a necessary condition for nonclassicality . The negativity of the Wigner function, which is well exhibited by the Fock states and the Schrödinger-cat-like states, is absent for the squeezed states. As for the oscillations in the photon-number probability, some even-odd oscillations can be simply obtained by using a statistical mixture of coherent states .

Many authors  have adopted the nonpositivity of the phase-averaged $`P`$-function $`F(I)=\frac{1}{2\pi }_0^{2\pi }𝑑\varphi P(I^{1/2}e^{i\varphi })`$ as the definition for a nonclassical state, since $`F(I)<0`$ invalidates Mandel’s semiclassical formula  of photon counting, i.e., it does not allow a classical description in terms of a stochastic intensity. Of course, some states can exhibit a “weak” nonclassicality , namely, a positive $`F(I)`$, but with a non-positive $`P`$-function (a relevant example being a coherent state undergoing Kerr-type self-phase modulation). However, from the point of view of the detection theory, such “weak” nonclassical states still admit a classical description in terms of having the intensity probability $`F(I)>0`$. For this reason, we adopt nonpositivity of $`F(I)`$ as the definition of nonclassicality.

## II SINGLE-MODE NONCLASSICALITY

The authors of Refs.  have recognized relations between $`F(I)`$ and generalized moments of the photon distribution, which, in turn, can be used to test the nonclassicality. The problem is reduced to an infinite set of inequalities that provide both necessary and sufficient conditions for nonclassicality . In terms of the photon-number probability $`p(n)=n|\widehat{\rho }|n`$ of the state with density matrix $`\widehat{\rho }`$, the simplest sufficient condition involves the following three-point relation for $`p(n)`$ 

$`B(n)`$ $``$ $`(n+2)p(n)p(n+2)`$ (2)   
$`(n+1)[p(n+1)]^2<0.`$   
Higher-order sufficient conditions involve five-, seven-, …, $`(2k+1)`$-point relations, always for adjacent values of $`n`$. It is sufficient that just one of these inequalities be satisfied in order to assure the negativity of $`F(I)`$. Notice that for a coherent state $`B(n)=0`$ identically for all $`n`$.

In the following we show that quantum tomography can be used as a powerful tool for performing the nonclassicality test in Eq. (2). For less-than-unity quantum efficiency ($`\eta <1`$), we rely on the concept of a “noisy state” $`\widehat{\varrho }_\eta `$, wherein the effect of quantum efficiency is ascribed to the quantum state itself rather than to the detector. In this model, the effect of quantum efficiency is treated in a Schrödinger-like picture, with the state evolving from $`\widehat{\varrho }`$ to $`\widehat{\varrho }_\eta `$, and with $`\eta `$ playing the role of a time parameter. Such lossy evolution is described by the master equation

$`_t\widehat{\varrho }(t)={\displaystyle \frac{\mathrm{\Gamma }}{2}}\left\{2\widehat{a}\widehat{\varrho }(t)\widehat{a}^{}\widehat{a}^{}\widehat{a}\widehat{\varrho }(t)\widehat{\varrho }(t)\widehat{a}^{}\widehat{a}\right\},`$ (3)   
wherein $`\widehat{\varrho }(t)\widehat{\varrho }_\eta `$ with $`t=\mathrm{ln}\eta /\mathrm{\Gamma }`$.

For the nonclassicality test, reconstruction in terms of the noisy state has many advantages over the true-state reconstruction. In fact, for nonunit quantum efficiency $`\eta <1`$ the tomographic method introduces errors for $`p(n)`$ which are increasingly large versus $`n`$, with the additional limitation that quantum efficiency must be greater than the minimum value $`\eta =0.5`$ . On the other hand, the reconstruction of the noisy-state probabilities $`p_\eta (n)=n|\widehat{\rho }_\eta |n`$ does not suffer such limitations, and even though all quantum features are certainly diminished in the noisy-state description, nevertheless the effect of nonunity quantum efficiency does not change the sign of the $`P`$-function, but only rescales it as follows:

$`P(z)P_\eta (z)={\displaystyle \frac{1}{\eta }}P(z/\eta ^{1/2}).`$ (4)   
Hence, the inequality (2) still represents a sufficient condition for nonclassicality when the original probabilities $`p(n)=n|\widehat{\rho }|n`$ are replaced with the noisy-state probabilities $`p_\eta (n)=n|\widehat{\rho }_\eta |n`$, the latter being given by the Bernoulli convolution

$`p_\eta (n)={\displaystyle \underset{k=n}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{k}{n}}\right)\eta ^n(1\eta )^{kn}p(k).`$ (5)   
Hence, when referred to the noisy-state probabilities $`p_\eta (n)`$, the inequality in Eq. (2) keeps its form and simply rewrites as follows

$`B_\eta (n)`$ $``$ $`(n+2)p_\eta (n)p_\eta (n+2)`$ (7)   
$`(n+1)[p_\eta (n+1)]^2<0.`$   
According to Eq. (7), the quantity $`B_\eta (n)`$ is nonlinear in the density matrix. This means that $`B_\eta (n)`$ cannot be measured by averaging a suitable kernel function over the homodyne data, as for any other observable . Hence, in the evaluation of $`B_\eta (n)`$ one needs to tomographically reconstruct the photon-number probabilities, using the kernel functions

$`K_\eta ^{(n)}(x)`$ $`=`$ $`2\kappa ^2e^{\kappa ^2x^2}{\displaystyle \underset{\nu =0}{\overset{n}{}}}{\displaystyle \frac{()^\nu }{\nu !}}\left({\displaystyle \genfrac{}{}{0pt}{}{n}{n\nu }}\right)(2\nu +1)!\kappa ^{2\nu }`$ (8)   
$`\times `$ $`\text{Re }\{D_{(2\nu +2)}(2i\kappa x)\},`$ (9)   
where $`D_\sigma (z)`$ denotes the parabolic cylinder function and $`\kappa =\sqrt{\eta /(2\eta 1)}`$. The true-state probabilities $`p(n)`$ are obtained by averaging the kernel function in Eq. (9) over the homodyne data. On the other hand, the noisy-state probabilities $`p_\eta (n)`$ are obtained by using the kernel function in Eq. (9) for $`\eta =1`$, namely without recovering the convolution effect of nonunit quantum efficiency. Notice that the expression (9) does not depend on the phase of the quadrature. Hence, the knowledge of the phase of the local oscillator in the homodyne detector is not needed for the tomographic reconstruction, and it can be left fluctuating in a real experiment.

Regarding the estimation of statistical errors, they are generally obtained by dividing the set of homodyne data into blocks. However, in the present case, the nonlinear dependence on the photon number probability introduces a systematic error that is vanishingly small for increasingly larger sets of data. Therefore, the estimated value of $`B(n)`$ will be obtained from the full set of data, instead of averaging the mean value of the different statistical blocks.

In Figs. 17 we present some numerical results that are obtained by a Monte-Carlo simulation of a quantum tomography experiment. The nonclassicality criterium is tested either on a Schrödinger-cat state $`|\psi (\alpha )(|\alpha +|\alpha )`$ or on a squeezed state $`|\alpha ,rD(\alpha )S(r)|0`$, wherein $`|\alpha `$, $`D(\alpha )`$, and $`S(r)`$ denote a coherent state with amplitude $`\alpha `$, the displacement operator $`D(\alpha )=e^{\alpha \widehat{a}^{}\overline{\alpha }\widehat{a}}`$, and the squeezing operator $`S(r)=e^{r(\widehat{a}^2\widehat{a}^2)/2}`$, respectively. Figs. 13 show tomographically-obtained values of $`B(n)`$, with the respective error bars superimposed, along with the theoretical values for a Schrödinger-cat state, for a phase-squeezed state ($`r>0`$), and for an amplitude-squeezed state ($`r<0`$), respectively. For the same set of states the results for $`B_\eta (n)`$ \[cf. Eq. (7)\] obtained by tomographic reconstruction of the noisy state are reported in Figs. 46. Let us compare the statistical errors that affect the two measurements, namely, those of $`B(n)`$ and $`B_\eta (n)`$ on the original and the noisy states, respectively. In the first case (Figs. 13) the error increases with $`n`$, whereas in the second (Figs. 46) it remains nearly constant, albeit with less marked oscillations in $`B_\eta (n)`$ than those in $`B(n)`$. Fig. 7 shows tomographically-obtained values of $`B_\eta (n)`$ for the phase-squeezed state (cf. Fig. 5), but for a lower quantum efficiency $`\eta =0.4`$. Notice that, in spite of the low quantum efficiency, the nonclassicality of such a state is still experimentally verifiable, as $`B_\eta (0)<0`$ by more than five standard deviations. In contrast, for coherent states one obtains small statistical fluctuations around zero for all $`n`$.

We remark that the simpler test of checking for antibunching or oscillations in the photon-number probability in the case of the phase-squeezed state considered here (Figs. 2, 5, and 7) would not reveal the nonclassical features of such a state.

## III TWO-MODE NONCLASSICALITY

Quantum homodyne tomography can also be employed to test the nonclassicality of two-mode states. For a two-mode state nonclassicality is defined in terms of nonpositivity of the following phase-averaged two-mode $`P`$-function :

$`F(I_1,I_2,\varphi )={\displaystyle \frac{1}{2\pi }}{\displaystyle _0^{2\pi }}𝑑\varphi _1P(I_1^{1/2}e^{i\varphi _1},I_2^{1/2}e^{i(\varphi _1+\varphi )}).`$ (10)   
In Ref.  it is also shown that a sufficient condition for nonclassicality is

$`C=(\widehat{n}_1\widehat{n}_2)^2(\widehat{n}_1\widehat{n}_2)^2\widehat{n}_1+\widehat{n}_2<0,`$ (11)   
where $`\widehat{n}_1`$ and $`\widehat{n}_2`$ are the photon-number operators of the two modes.

A tomographic test of the inequality in Eq. (11) can be performed by averaging the kernel functions for the operators in the ensemble averages in Eq. (11) over the two-mode homodyne data. For the normal-ordered field operators one can use the Richter formula in Ref. , namely

$`[a^na^m](x,\varphi )=e^{i(mn)\varphi }{\displaystyle \frac{H_{n+m}(\sqrt{2\eta }x)}{\sqrt{(2\eta )^{n+m}}\left(\genfrac{}{}{0pt}{}{n+m}{n}\right)}},`$ (12)   
$`H_n(x)`$ denoting the Hermite polynomial and $`\varphi `$ being the phase of the fiels with respect to the local oscillator of the homodyne detector. Again, as for the kernel function in Eq. (9), the value $`\eta =1`$ is used to reconstruct the ensemble averages of the noisy state $`\widehat{\rho }_\eta `$. Notice that for $`n=m`$ Eq. (12) is independent on the phase $`\varphi `$, and hence no phase knowledge is needed to reconstruct the ensemble averages in Eq. (11). As an example, we consider the twin-beam state at the output of a nondegenerate parametric amplifier

$`|\chi (1|\lambda |^2){\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}\lambda ^n|n|n,`$ (13)   
where $`|n|n`$ denotes the joint eigenvector of the number operators of the two modes with equal eigenvalue $`n`$, and the parameter $`\lambda `$ is related to the gain $`G`$ of the amplifier by the relation $`|\lambda |^2=1G^1`$. The theoretical value of $`C`$ for the state in Eq. (13) is $`C=2|\lambda |^2/(1|\lambda |^2)<0`$. A tomographic reconstruction of the twin-beam state in Eq. (13) is particularly facilitated by the self-homodyning scheme, as shown in Ref. . With regard to the effect of quantum efficiency $`\eta <1`$, the same argument still holds as for the single-mode case: one can evaluate $`C_\eta `$ for the twin-beam state that has been degraded by the effect of loss. In this case, the theoretical value of $`C_\eta `$ is simply rescaled to $`C_\eta =2\eta ^2|\lambda |^2/(1|\lambda |^2)`$.

In Fig. 8 we report $`C_\eta `$ vs. $`1\eta `$, $`\eta `$ ranging from 1 to 0.3 in steps of 0.05, for the twin-beam state in Eq. (13) with $`|\lambda |^2=0.5`$, corresponding to the total average photon number equal to 2. The values of $`C_\eta `$ result from a Monte-Carlo simulation of a homodyne tomography experiment with a sample of $`4\times 10^5`$ data, using the theoretical joint homodyne probability of the state $`|\chi `$

$`p_\eta (x_1,x_2,\varphi _1,\varphi _2)={\displaystyle \frac{2\mathrm{exp}\left[\frac{(x_1+x_2)^2}{d_z^2+4\mathrm{\Delta }_\eta ^2}\frac{(x_1x_2)^2}{d_z^2+4\mathrm{\Delta }_\eta ^2}\right]}{\pi \sqrt{(d_z^2+4\mathrm{\Delta }_\eta ^2)(d_z^2+4\mathrm{\Delta }_\eta ^2)}}},`$ (14)   
with

$`z=e^{i(\varphi _1+\varphi _2)}\mathrm{\Lambda },`$ (15)   
$`d_{\pm z}^2={\displaystyle \frac{|1\pm z|^2}{1|z|^2}},`$ (16)   
$`\mathrm{\Delta }_\eta ^2={\displaystyle \frac{1\eta }{4\eta }},`$ (17)   
$`\varphi _1`$ and $`\varphi _2`$ denoting the phases of the two modes relative to the respective local oscillator. Notice that the nonclassicality test in terms of the noisy state gives values of $`C_\eta `$ that are increasingly near the classically positive region for decreasing quantum efficiency $`\eta `$. However, the statistical error remains constant and is sufficiently small to allow recognition of the nonclassicality of the twin-beam state in Eq. (13) up to $`\eta =0.3`$.

## IV CONCLUSIONS

We have shown that quantum homodyne tomography allows one to perform nonclassicality tests for various single- and two-mode radiation states, even when the quantum efficiency of homodyne detection is rather low. The method involves reconstruction of the photon-number probability or of some suitable function of the number operators pertaining to the noisy state, namely, the state degraded by the less-than-unity quantum efficiency. The noisy-state reconstruction is affected by the statistical errors; however, they are sufficiently small that the nonclassicality of the state can be tested even for low values of $`\eta `$. For the cases considered in this paper, we have shown that the nonclassicality of the states can be proven (deviation from classicality by many error bars) with $`10^5`$$`10^7`$ homodyne data. Moreover, since the knowledge of the phase of the local oscillator in the homodyne detector is not needed for the tomographic reconstruction, it can be left fluctuating in a real experiment. Hence, we conclude that the proposed nonclassicality test should be easy to perform experimentally.

## Acknowledgments

This work is supported by the Italian Ministero dell’Universitá e della Ricerca Scientifica e Tecnologica under the program Amplificazione e rivelazione di radiazione quantistica