Arxiv:1712.01003V2 [Quant-Ph] 21 Feb 2018 Ρ = a Ρ a Ρ Lator Underlying Each Mode of the Electromagnetic ﬁeld

Nonclassical states of light with a smooth P function

François Damanet1,2, Jonas Kübler3, John Martin2, and Daniel Braun3 1 Department of Physics and SUPA, University of Strathclyde, Glasgow G4 0NG, United Kingdom. 2 Institut de Physique Nucléaire, Atomique et de Spectroscopie, CESAM, Université de Liège, Bâtiment B15, B - 4000 Liège, Belgium. 3 Institut für theoretische Physik, Universität Tübingen, 72076 Tübingen, Germany.

There is a common understanding in quantum optics that nonclassical states of light are states that do not have a positive semidefinite and sufficiently regular Glauber-Sudarshan P function. Al- most all known nonclassical states have P functions that are highly irregular which makes working with them difficult and direct experimental reconstruction impossible. Here we introduce classes of nonclassical states with regular, non-positive-definite P functions. They are constructed by “puncturing” regular smooth positive P functions with negative Dirac-delta peaks, or other sufficiently narrow smooth negative functions. We determine the parameter ranges for which such punctures are possible without losing the positivity of the state, the regimes yielding antibunching of light, and the expressions of the Wigner functions for all investigated punctured states. Finally, we propose some possible experimental realizations of such states.

I. INTRODUCTION all possible quantum states ρ of the light-field (quantum or classical) can be represented by (1). This is possible Creating and characterizing nonclassical states of light due to the fact that coherent states are overcomplete. are two tasks of substantial interest in a wide range of ap- The price to pay, however, is that in general P (α) is not plications [1, 2], such as sub-shot-noise measurements [3– a smooth function, but a functional. Obviously, indeed, 5], high resolution imaging [6–9], light source and detec- even a coherent state α0 is already represented by a functional, namely a Dirac-delta| i peak, P (α) = δ(α α ). tor calibration [10–12], quantum cryptography [13, 14], − 0 or quantum information processing [15]. Such a level of singularity represents the “boundary of In quantum optics, there is wide-spread agreement that acceptable irregular behavior” of the P function for a the most classical pure states of light are coherent states classical state. Much worse singularities arise for exam- α . These are eigenstates with complex eigenvalue α of ple from simple Fock states n , i.e. energy eigenstates | i of the mode, a†a n = n n , for| i which the P function is the annihilation operator a of the mode of the light-field | i | i under consideration, a α = α α , and are characterized given by the n-th derivative of a delta-function [18]. by the minimal uncertainty| i product| i of the quadratures X = (a + a†)/√2 and Y = (a a†)/(i√2). In the case It would be desirable that one could directly recon- of a mechanical harmonic oscillator− the quadratures cor- struct the P function and in this way show the nonclas- respond simply to phase space coordinates x and p in sicality of states that are genuinely quantum. However, suitable units. Hence, a coherent state resembles most in the case of highly non-regular P functions this is im- closely the notion of a classical phase space point located possible. One way out is to apply filters that do not en- at the mean position of the quadratures. The label α C hance the negativity of the characteristic function, i.e. the determines these mean values as α X α = √2 Re(∈α) Fourier transform of P (α) [19–22]. h | | i and α Y α = √2 Im(α). Coherent states of light, Also for the ease of theoretical work one would like to whichh can| | bei produced by lasers operated far above the have classes of P functions that are regular and smooth, threshold [16], retain the property of minimal uncertainty but represent genuine quantumness due to the fact that p ∆X∆Y = 1/2, where ∆X = X2 X 2 is the stan- they become negative [23]. One class of such states are dard deviation, under the dynamicsh ofi − the h harmonici oscil- single-photon added thermal states [20], i.e. states of the form † , where represents a thermal state arXiv:1712.01003v2 [quant-ph] 21 Feb 2018 ρ = a ρ a ρ lator underlying each mode of the electromagnetic field. N T T Mixing classically coherent states (i.e. choosing randomly and is a normalization factor. The P function of this stateN is given by the product of a linear function of α 2 such states according to a classical probability distri- | | bution) should not increase their quantumness. Hence, and the Gaussian of the thermal state (see Eq. (4.36) (possibly mixed) quantum states given by a density op- in [18]). But as Agarwal writes (p.84 in [18]): “Here we erator ρ as perhaps have the unique case where P (α) exists but is negative. Most known cases of nonclassical P involve P Z functions which do not even exist.“ ρ = P (α) α α d2α (1) | ih | In this paper, we introduce novel whole classes of with a positive P function are considered “classical such genuinely-quantum states with nevertheless smooth states” in quantum optics. While it may appear that (1) P functions (i.e. exhibiting at most singularities in the describes only states that are “diagonal” in a coherent form of Dirac δ functions). The recipe for doing so is state basis, it was shown by Sudarshan [17] that in fact very simple: Start− with a classical quantum state with 2 a smooth positive P function, and then superpose some If the weights are chosen such that the P function narrow, negative peaks. Altogether, one must guarantee, (2) becomes negative for some α and that the result- of course, that the density matrix ρ remains a semidefi- ing state ρ is a positive semidefinite operator with unit nite positive operator. Finding parameter regimes where trace, it follows immediately that the state is nonclassi- this is the case constitutes the main technical difficulty cal. Furthermore, the P function remains “well-behaved” for constructing such states. Obviously, this general (i.e. sufficiently regular), if πi(α) are well-behaved, as we recipe allows for an infinity of possibilities, starting from assumed Pcl(α) smooth. In agreement with the common the choice of the smooth original P function, over the understanding that positive P functions corresponding number and form of the negative peaks superposed, to to simple Dirac-δ peaks are still accepted in the class of their amplitudes and positions. We start with the sim- “well-behaved” functions because they represent coherent plest situation of the rotationally symmetric Gaussian P states, we consider in the following sections Dirac δ(α) function centered at α = 0 of a thermal state, “punc- functions as worst singularities for the punctures πi(α). tured” with a single negative| | delta-peak whose position This is the also the simplest case, as it does not contain and amplitude we can vary. Only the radial position is any further free parameter. relevant here, and we have herewith already three parameters, exemplifying the wide range of tunability of punctured states. B. Conditions for (semi)positive definiteness The paper is organized as follows. In Section II, we introduce the “punctured states”, their general properties, For a given choice of the functions πi(α) for the punc- the conditions for positive definiteness of their density tures, positivity of the state will depend on the position matrix representation and the conditions yielding non- of the αi C and the weights wi. Our task is hence classical states. In Section III, we investigate the partic- to find the∈ allowed parameter ranges for w , α such { i i} ular simple case of a thermal state punctured by a single that ρ > 0, i.e. ρ is a positive semidefinite operator. The delta-peak. In Section IV, we study the case of a squeezed normalization factor ensures that R P (α)d2α = 1, i.e. thermal state punctured by a delta-peak. In Section V, N − we generalize our results by considering a thermal state N ! 1 X punctured by a Gaussian peak of finite width. In each = 1 w . (4) N − i Section, we provide the expression of the Wigner func- i=1 tion, calculate the second-order correlation function and Since does not influence the positivity, we will often identify the regimes of parameters yielding antibunch- skip it.N ing of light, i.e. sub-Poissonian statistics of the photon It is clear that in general the condition for any single counting. In Section VI, we discuss some possibilities of parameter will depend on the values of all the others. experimentally creating some of these states. Finally, in This can be seen immediately from the general neces- Section VII, we give a conclusion of our work. sary and sufficient condition of positive semidefiniteness, 2 ψ ρ ψ > 0 ψ L (C), which for the state (1) withh | | aiP function∀ | i ∈ of H the ≡ form (2) reads II. DEFINITION, GENERAL PROPERTIES, POSITIVITY AND NONCLASSICALITY Z X CONDITIONS P (α) w π (α α ) ψ α 2d2α 0 ψ . cl − i i − i |h | i| > ∀| i ∈ H i (5) A. General punctured states To go beyond this general criterion, we first review briefly a few necessary and sufficient conditions for positivity of Let ρcl represent a classical state, that is a state as a linear operator. defined in (1) with a smooth non-negative P function Necessary and sufficient condition ( ). A neces- Pcl(α) > 0. We define a “punctured ρcl state” as a sary and sufficient condition ( ) forNSC a hermitian ma- state (1) with a P function of the form NSC trix ρ to be semi-positive definite (ρ > 0) is that all its eigenvalues are non-negative. While the eigenvalues of " N # X an infinite dimensional matrix are generally not readily P (α) = P (α) w π (α α ) , (2) N cl − i i − i accessible, the condition nevertheless allows for a numer- i=1 ical approach. Indeed, the density operator (1) can be where is a normalization constant, the are posi- expressed in the Fock state basis n : n N0 using wi {| i ∈ } tive coefficientsN representing the weights of the different the expansion of coherent states punctures i (i = 1,...,N), and π (α) are positive func- i 2 ∞ n − |α| X α tions on the complex plane that determine their shape. α = e 2 n . (6) | i √ | i We take them as normalized according to n=0 n! Z 2 Then, by introducing a cut-off dimension in Hilbert space πi(α) d α = 1 i = 1,...,N. (3) ∀ of sufficiently high excitation nmax, a finite dimensional 3 matrix representation can be obtained and one can check function can be directly measured (see e.g. [26–31]). numerically that when increasing the cut-off, the lowest Since the Wigner function W (α) is the convolution of eigenvalue converges to a value that is sufficiently far the P function with the vacuum state, i.e. from zero for judging the positivity of the state. 2 Z 2 Necessary condition ( ). Since ρ > 0 ψ ρ ψ > 0 W (α) = P (β) e−2|α−β| d2β, (9) for all states ψ , a necessaryNC condition (⇔ h) for| | ρi 0 π | i NC > is that all diagonal elements ρφφ φ ρ φ of the density the existence of negative values for the Wigner function ≡ h | | i matrix are non-negative for a given class of states φ implies the existence of negative values for the func- | i ∈ P where is a subset of , i.e. tion, but the reverse is not true. Thus, the criterion based K ⊂ H K H ρ φ ρ φ 0, φ ρ 0. (7) on negative values of the Wigner function detects fewer φφ ≡ h | | i > ∀ | i ∈ K ⇐ > nonclassical states than the previous criterion. More Sufficient condition ( ). A sufficient condition ( ) specifically, it does not detect states with partly neg- for positivity can be foundSC from Gershgorin’s Circle The-SC ative P function but everywhere positive Wigner func- orem. Consider first the case of finite dimensional matri- tion. This explains why the definitions of punctured ces A. The theorem then says [24]: “Every eigenvalue of a states based on the puncturing of P functions rather than complex square matrix A of coefficients a lies within Wigner functions is more appropriate. Moreover, delta- { ij} at least one of the Gershgorin disks D(aii,Ri) centred at puncturing of a Wigner function would be forbidden as P P aii with radii Ri = j=6 i aij or Ri = j=6 i aji .” An Wigner functions are always bounded (see e.g. p. 73 in interpretation of this theorem| | is that if the off-diagonal| | [32]). elements of A have small norms compared to its diagonal Antibunching. Another standard criterion of nonclas- elements, its eigenvalues cannot be far from the values sicality is based on the second order correlation function of the diagonal elements. A matrix is called (strictly) a†a†aa diagonally dominant if aii > Ri ( aii > Ri) for all i. g(2) = h i, (10) This leads to the following| | sufficient| condition:| a hermi- a†a 2 h i tian (strictly) diagonally dominant matrix with real non- where = Tr ( ρ), or, equivalently, on the Mandel Q negative diagonal entries is positive semidefinite (defi- M parameterh · i [33] · nite). Gershgorin’s Circle Theorem has been generalized to the case of infinite dimensional matrices in [25], where 2 (∆n) n (2) a proof was given for matrices of the space 1. The proof QM = − h i = n g 1 , (11) L n h i − for the space L2 in which our operators live can be done h i analogously and leads to the same result, namely that the where n = a†a. A state is said to be nonclassical if ∞ (2) eigenvalues lie in i=1D(aii,Ri) (just as in the case of a g < 1 (or, equivalently, QM < 0), which corresponds to finite space). Hence,∪ by expressing the density matrix in antibunching (sub-Poissonian statistics of photon count- the Fock state basis, we obtain the sufficient condition ing that cannot be observed for mixtures of coherent ( ) states). SC It is interesting to note the connection between classical- ρ R 0 n ρ 0, (8) nn − n > ∀ | i ⇒ > ity in terms of the P function and in terms of antibunch- P ing, pointed out by Kimble et al. [34], and generalized where ρnn n ρ n and Rn = m=6 n ρnm . ≡ h | | i | | by Vogel to space- and time-dependent correlations [35]: g(2) < 1 can exist if the P function is negative somewhere. C. Conditions for nonclassicality Conversely, if the P function of a state is a valid proba- (2) bility density, then there is no antibunching, i.e. g > 1. There exist various nonclassicality criteria for radiation Proof. This statement can be proved using Jensen’s in- fields. In this work, we shall consider criteria based on equality. We use equation (5’) in the original paper [36]. the negativity of the P function, the negativity of the Let P (α) be a valid probability density defined over the R 2π iφ Wigner function and the existence of antibunching [1]. complex numbers. Then Pr(r) = 0 rP (re ) dφ de- Negativity of the P function. A criterion for a state fines a valid probability density over the real numbers to be nonclassical is that its P function takes negative r [0, ). Using Eq. (1), the second order correlation values, so that it can no longer be interpreted as a prob- function∈ ∞ (10) can be written as ability distribution in phase space. For example, all phys- R 4 2 R ∞ 4 α P (α) d α r Pr(r) dr ical δ-punctured states are nonclassical according to this g(2) = | | = 0 . (12) criterion. Other punctured states, however, could have R 2 2 2 R ∞ 2 2 α P (α) d α 0 r Pr(r) dr a non-negative P function in the whole complex plane, | | so that general punctured states are not necessarily non- Defining the convex functions f(r) = ϕ(r) = r2, we can classical. rewrite the above as Negativity of the Wigner function. Negativity of the R ∞ (2) 0 ϕ(f(r)) Pr(r) dr Wigner function is a widely used criterion for nonclassi- g = R ∞ . (13) cality. This is largely due to the fact that the Wigner ϕ 0 f(r) Pr(r) dr 4

(2) The condition g > 1 is then equivalent to The minimum over γ of the rhs term is attained for γ = n¯+1 α and directly leads to the condition (16). Z ∞ Z ∞ n¯ 1 Suﬃcient condition. In the Fock state basis, the Ger- ϕ(f(r)) Pr(r) dr > ϕ f(r) Pr(r) dr , (14) 0 0 shgorin disks D(ρnn,Rn) have center ρnn n ρ n and radius R given by ≡ h | | i which is true by Jensen’s inequality. n

n 2 2n n¯ −|α1| α1 As a conclusion, for radiation field states, both the ρnn = w1e | | , (20) negativity of the Wigner function and antibunching im- (¯n + 1)n+1 − n! +∞ ply that the P function is not a valid probability density, 2 n m −|α1| α1 X α1 while the reverse is not necessarily true. Rn = w1e | | | | . (21) √n! m=0 √m! In the following sections, we investigate the cases of m6=n various punctured states. In each case, we apply the conditions of positivity of and nonclassicality of ρ to de- The sufficient condition (8) then implies the following termine the physical nonclassical states. For the sake of condition on w1, simplicity, we restrict ourselves in the remainder of the paper to a single puncture by setting in Eq. (2). 2 ( !n ) N = 1 |α1| n¯ e 1 n¯+1 √ δ,SC w1 < min n! wmax (22) n¯ + 1 f( α ) n∈ α ≡ | 1| N | 1| III. DELTA-PUNCTURED THERMAL STATES where We define a δ-punctured thermal state as a state deter- +∞ m mined through (1) and (2) with X α1 f( α1 ) = | | . (23) 2 | | √m! e−|α| /n¯ m=0 Pcl(α) = , πn¯ (15) δ,SC The bound wmax depends only on the absolute value of π(α) = δ(α), δ,SC α1, not on its phase. Also, we have wmax = 0 only if n¯ = 0. The largest values of w allowed by the bound lie where n¯ = (e~ω/kB T 1)−1 is the average number of 1 in the region where α and n¯ are small, which is shown thermal excitations at− temperature T of the radiation 1 in Fig. 1. Since wδ,SC| >| 0 for all values of n¯ > 0 and α , field with frequency ω. max 1 this is an analytical proof that there always exist physical nonclassical δ-punctured thermal states. δ,NC A. Conditions for (semi)positive definiteness Tightness of wmax . Computations based on the diagonalization of δ-punctured thermal states and the application of suggest that the bound wδ,NC given Necessary condition. A necessary condition for the max in Eq. (16)NSC is tight. This claim is further supported positivity of ρ can be obtained from the condition (7) by the observation that the coherent state n¯+1 α with the set of coherent states. It yields the following | n¯ 1i K from which the bound is derived is an eigenstate of ρ upper bound for the puncture weight w1: δ,NC with w1 = wmax with eigenvalue 0. In Fig. 1, we dis- 2 δ,NC −|α1| /n¯ play both wmax and the value of w1 that cancels the e δ,NC w1 6 wmax . (16) smallest eigenvalue of ρ. In our computations, the den- n¯ + 1 ≡ sity matrix is expressed in the truncated Fock state ba- Proof. Let γ (γ C) denote a coherent state. The sis 0 , 1 ,..., nmax with a sufficiently large cut-off | i ∈ {| i | i | i} density matrix element ργγ γ ρ γ is then given by nmax in order ensure stable numerical results. An ana- ≡ h | | i lytical proof of the tightness for a puncture at α1 = 0 Z 2 1 − |α| −|α−γ|2 2 −|α −γ|2 will be provided in the context of Gaussian punctures ρ = e n¯ d α w e 1 , γγ N πn¯ − 1 (Sec. V A). Note that a delta-puncture at α = 0 with δ,NC (17) maximum weight wmax corresponds to the complete re- 2 where we used γ α 2 = e−|α−γ| . The condition that moval of the ground-state of the oscillator (“the vacuum” |h | i| in quantum-field theory parlance) from the state. We call ργγ must be non-negative for all γ reads such states “vacuum-removed states” for short. 2 R − |α| −|α−γ|2 2 In summary, our results indicate that any δ-punctured 1 e n¯ d α (18) thermal state with satisfying Eq. (16) corresponds to w1 6 −|α −γ|2 γ C. w1 πn¯ e 1 ∀ ∈ a proper physical state. As can be seen in Fig. 1, the After integration, this yields the condition value of the bound of allowed weights w1 decreases as α1 increases, showing that a thermal state cannot be − 1 |γ|2+|α −γ|2 | | e ( n¯+1 ) 1 punctured significantly far from its center. This feature w1 γ C. (19) 6 n¯ + 1 ∀ ∈ is less pronounced as n¯ increases. 5

1 2 ’g2_delta_thermal.dat’ 1 g(2) 0.8 3 0.8 0.6 1.5 δ max

w 0.4 2 0.6 0.2 |α1| 1 0 δ max 0 0.5 1 1.5 2 w 0.4 α1 2 1 | | ’g2_delta_thermal.dat’ 0.5 g(2) 3 0.2 1.5 0 0 0 0.5 1 1.5 2 2 0 n¯ 0 0.5 1 1.5 2 |α1| 1 α1 | | 1 1 0.5 δ,NC 3 = δ,NC FIG. 1. Plot of the upper bounds wmax (Eq. (16), solid w1 wmax g(2) δ,SC δ curves) and wmax (Eq. (22), dashed curves), and of wmax, 2 0 g(2) 3 the w1 which cancels the smallest eigenvalue of the density 0.75 0 1 matrix, as a function of |α1| for diﬀerent n¯ (green squares: 0 0.5 1 1.5 2 2 n¯ = 0.1, blue dots: n¯ = 0.5, orange triangles: n¯ = 1, and 0 n¯ |α1| 0 0.25 0.5 0.75 1 purple diamonds: n¯ = 2). The eigenvalues are computed in 0.5 |α1| 1 a Fock space of dimension nmax = 15. Inset: convergence of δ w (green squares) with increasing cut-oﬀ dimension nmax max 0.250. 0 (nmax = 2, 3, 4 and 15 from top to bottom) for n¯ = 0.1.

0 0 0.25 0.5 0.75 1 w1 B. Conditions for nonclassicality FIG. 2. Density plot of the second order correlation func- Negativity of the P function. All δ-punctured thermal tion g(2) [Eq. (26)] for single δ-punctured thermal states as states are nonclassical due to the presence of the infinite a function of |α1| and n¯ for the maximal puncture weight δ,NC negative δ-peak in the P function. w1 = wmax given by Eq. (16) (top) and as a function of |α1| Negativity of the Wigner function. The Wigner func- and w1 for n¯ = 0.2 (bottom). Only the range where ρ > 0 is shown. The black dashed curve in the top figure, along which tion of δ-punctured thermal states follows from Eq. (9) (2) g = 1, delimits the sets of states displaying antibunching and is given by or not. The inset of the bottom figure represents g(2) as a δ,NC 2 function of α1 for the maximum weight w1 = wmax . − 2|α| ! 2 e 1+2¯n 2 W (α) = N w e−2|α−α1| . (24) π 1 + 2¯n − 1

It takes negative values for large enough puncture weight IV. DELTA-PUNCTURED SQUEEZED THERMAL STATES 2 e−|α1| /n¯ w1 > . (25) 1 + 2¯n To generalize the previous results, we now replace the Such weights are still acceptable as long as they do not thermal state by a squeezed thermal state [38, 39] exceed the (larger) bound (16). Z Antibunching. For δ-punctured thermal states, the 1 2 ρ = e−|α| /n¯ S(r) α α S†(r) d2α, (27) second-order correlation function (10) is given by cl πn¯ | ih |

2 4 (2) 2¯n w1 α1 †2 2 g = 1 w1 − | | . (26) where S(r) = exp r[a a ]/2 is the squeezing op- − (¯n w α 2)2 − − − 1| 1| erator and r is the squeezing parameter, taken real for simplicity. We remind that the P function of a state ρ For w 0, we recover the known value g(2) = 2 of the 1 → is given by the Fourier transform of its normal-ordered thermal state [37]. In Fig. 2, we show a density plot of † ∗ characteristic function χ(β) = tr[ρ ei(βa +β a)]. For g(2) as a function of n¯ and α for the maximal allowed 1 the squeezed thermal state ρ defined in Eq. (27), the value of the puncture weight| | , i.e. the bound (16). A cl w1 normal-ordered characteristic function reads [38] whole region of parameter space corresponds to punctured states giving rise to antibunching (g(2) < 1). Also, − 2 − 2 g(2) increases beyond 2 for n¯ 0.1, α 0.5. χ(β) = e n¯I βR n¯RβI , (28) ' | 1| ' 6 where βR and βI are the real and imaginary parts of where the modulus squared of the scalar product between β C, and a squeezed coherent state γ, r0 and a coherent state α ∈ is given by [37] | i | i −2r −r n¯R = e n¯ e sinh r, − (29) 0 2 0 − | |2 | |2 ∗ ∗ 0 n¯ = e2r n¯ + er sinh r. γ, r α = sech(r ) e ( α + γ )+(α γ+γ α)sech(r ) I |h | i| 0 2 2 e−th(r ) Re α −γ . (35) The Fourier transform of the characteristic function – and × thus the function – exists only if both the coefficients P The minimum of the rhs term in Eq. (34) is obtained, n¯ and n¯ in front of β2 and β2 are positive, which I R R I after integration over α, for imposes the condition n¯R +n ¯I + 2¯nRn¯I −2|r| −|r| (30) γ = p α1, e n¯ e sinh r > 0 (¯n n¯ + 2¯n n¯ ) (¯n n¯ + 2¯n n¯ ) − | | R − I R I I − R R I (36) n¯ n¯ thereby limiting the range of allowed squeezing parame- r0 = arctanh R − I ters r R for a given n¯. The characteristic function of 2¯nRn¯I a squeezed∈ thermal state being an asymmetric Gaussian distribution, the same goes for its function. There- and leads to the necessary condition (32). P δ,NC fore, a δ-punctured squeezed thermal state is a state de- Tightness of wmax, squ. Computations relying on the termined through (1) and (2) with diagonalization of ρ given by Eq. (1) with Eqs. (2) δ,NC and (31) suggest that the bound wmax, squ given in α2 α2 Eq. (32) is tight. As for the case of thermal states, this − R + I e n¯R n¯I claim is supported by the observation that the squeezed Pcl(α) = , (31) coherent state γ, r0 with γ and r0 given by Eq. (36) π√n¯Rn¯I is an eigenstate| of ρi with eigenvalue 0. Figure 3 shows π(α) = δ(α), δ,NC wmax, squ and the value of w1 that cancels the smallest eigenvalue of ρ. Numerical results converge to the bound where αR and αI are the real and imaginary parts of as the cut-off dimension increases. α C. ∈ 0.6 0.6

A. Conditions for (semi)positive definiteness 0.4 squ , δ 0.4 max 0.2 Necessary condition. Since Pcl(α) is rotationally non- w squ symmetric, it is reasonable to assume that the necessary , 0 0 0.5 1 1.5 2 δ condition (7) will provide us with a tight bound for w1 max w α | 1| only if we consider states φ with non-symmetric 0.2 P functions. By taking | ito ∈ be K the set of squeezed coherent− states [38, 39], weK obtain the following upper bound for the puncture weight w1: 0 α2 α2 − 1R + 1I 0 0.5 1 1.5 2 n¯R n¯I α 2√n¯Rn¯I e δ,NC | 1| w1 6 wmax, squ. (32) n¯R +n ¯I + 2¯nRn¯I ≡ δ,NC FIG. 3. Plot of the upper bound wmax, squ (solid curves) and δ When n¯R =n ¯I =n ¯, this condition reduces to the condi- of wmax, squ, the w1 which cancels the smallest eigenvalue of tion (16), as expected. the density matrix, as a function of |α1| for n¯R = 1, n¯I = 0.1 0 Proof. Let us denote by γ, r a squeezed coherent state (squares) and n¯I = 0.5 (dots). The eigenvalues are computed defined as [38, 39] | i in a Fock space of dimension nmax = 25. Inset: convergence δ of wmax, squ (green squares) with increasing cut-off dimension γ, r0 = D(γ)S(r0) 0 , (33) nmax (nmax = 6, 8, 25 from top to bottom) for n¯I = 0.1. | i | i where D(γ) = exp γa† γ∗a is the displacement operator and S(r0) the squeezing− operator of parameter r0, where the prime symbol distinguishes r0 from r, this lat- B. Conditions for nonclassicality ter defining n¯R and n¯I . The necessary condition (7) with φ = γ, r0 reads | i | i Negativity of the P function. As for δ-punctured ther- R 0 2 2 Pcl(α) γ, r α d α mal states, the presence of a δ peak always enforces non- w |h | i| , (34) 1 6 γ, r0 α 2 classicality. |h | 1i| 7

2 Negativity of the Wigner function. The Wigner func- ’g2_delta_thermal_alpha0.1.dat’ (2) tion of a δ-punctured squeezed thermal state reads g 3 1.5 α2 α2 −2 R + I ! 2 e 1+2nR 1+2nI 2 W (α) = N w e−2|α−α1| , π 1 + 2¯n − 1 2 n¯I 1 (37) and takes negative values if the puncture weight satisﬁes 1 0.5 α2 α2 − R + I e n¯R n¯I 0 w1 > , (38) 0 1 + 2¯n 0 0.5 1 1.5 2 n¯R which is still acceptable as long as it does not exceed the (larger) bound (32). FIG. 4. Density plot of the second order correlation function (2) Antibunching. For δ-punctured squeezed thermal g [Eq. (39)] for single δ-punctured squeezed thermal states states, the second-order correlation function (10) is given as a function of n¯R and n¯I for the maximal puncture√ weight by w1 given in Eq. (32) for α1R = α1I = 0.1/ 2. The black dashed curve delimits the region where g(2) < 1. Note that (2) 1 2 2 4 in the absence of puncture, g > 2 for any squeezed thermal (2) 4 3¯nR + 2¯nRn¯I + 3¯nI w1 α1 state (see text). g = 1 w1 −2 | | . − n¯R+¯nI 2 2 w1 α1 − | | (39) For w1 0, we recover the known value for squeezed thermal→ states [38] A. Vacuum-centered Gaussian punctures

2 2 We first consider the simple case of a thermal state (2) 3¯nR + 2¯nRn¯I + 3¯nI g = 2 from which we subtract a single Gaussian centered at (¯nR +n ¯I ) α1 = 0. A thermal state with the vacuum component ex- 2 2 2 (40) (2¯n + 1) sinh r cosh r actly removed (case b 0, n¯ = 1, w1 = 1/2) has already = 2 + 2 > 2. → n¯ cosh(2r) + sinh2 r been considered in the literature as a simple example of a nonclassical state [40]. Here we keep a variable width (2) and amplitude for the subtracted Gaussian. Figure 4 shows a density plot of g as a function of n¯R and n¯I for α1 = 0.1 with the maximal allowed value of , i.e. the| bound| (32). Again, a whole region of param- w1 1. Condition for (semi)positive definiteness eter space (delimited by the black dashed curve) corresponds to states displaying antibunching (g(2) < 1), to be contrasted with squeezed thermal states (i.e. without Necessary and sufficient condition. Since both the (2) original and the subtracted state are diagonal in the Fock puncture) for which g > 2. state basis, can be tackled analytically. The P function of the vacuum-centredNSC Gaussian-punctured thermal state reads V. GAUSSIAN-PUNCTURED THERMAL 2 2 − |α| − |α| ! STATES e n¯ e b P (α) = w . (42) N πn¯ − 1 πb We now generalize the results of Sec. III by replacing δ-punctures by single narrow Gaussian punctures. We It takes negative values at α = 0 as soon as w1 > b/n¯. As thus define Gaussian-punctured thermal states as states we will show, only b < n¯ can lead to a positive semidefi- determined through (1) and (2) with nite density operator. Expressing the coherent states in the Fock state basis using Eq. (6) and performing the 2 e−|α| /n¯ Gaussian integrals, we find Pcl(α) = , πn¯ (41) ∞ −| |2 X e α /b ρ = pcl(n) w1 pπ(n) n n , π(α) = , N − | ih | πb n=0 with where b > 0 characterizes the width of the puncture. Note that δ-punctured thermal states correspond to the n¯n bn p (n) = , p (n) = . (43) limit b 0. cl (¯n + 1)n+1 π (b + 1)n+1 → 8

Since ρ is diagonal, the positivity condition reads p (n) attains negative values in some region is obtained if and cl − w1 pπ(n) 0 for all n N0, which is equivalent to only if > ∈ b + 1 nb¯ +n ¯ n b b + 1 G, NSC b < n¯ and < w1 6 . (50) w1 6 inf wmax, vac. (44) n¯ n¯ + 1 n¯ + 1 n∈N0 nb¯ + b ≡ Hence, the puncture weight w must be smaller than We can now distinguish two cases: i) If , then 1 b > n¯ (b + 1)/(¯n + 1) to ensure the physicality of the state, but the fraction in the argument of the infimum is smaller greater than b/n¯ to yield a nonclassical state. Moreover, than 1 and the infimum is obtained in the limit n in the limit b 0, the Gaussian puncture becomes a δ- and evaluates to 0, hence w = 0. This shows that→ → ∞ 1 puncture, and we recover the value given in (16). This we cannot subtract a broader Gaussian than the original provides an analytical proof of the tightness of the bound one. ii) In the other case, where , the minimum b 6 n¯ wδ,NC in the case α = 0. is obtained for n = 0 and evaluates to 1. We can thus max 1 subtract a Gaussian that is tighter than the original one and still obtain a positive semidefinite density operator. B. Arbitrarily centred Gaussian punctures This works as long as We now consider the more general case of an arbitrarily b + 1 G, NSC w1 wmax, vac. centred Gaussian puncture corresponding to a function 6 n¯ + 1 ≡ P of the form

 2 2  − |α| − |α−α1| 2. Conditions for nonclassicality e n¯ e b P (α) =  w1  . (51) N πn¯ − πb Negativity of the P function. In order to obtain a P function which exhibits negative values for some α, we We first obtain an upper bound on the amplitude of the see from Eq. (42) that this requires puncture, and then support the tightness of this bound with numerical results. Finally, we find the regimes of b w > , (45) parameters yielding nonclassicality and antibunching. 1 n¯ showing that in the limit b 0, we recover the fact that 1. Condition for (semi)positive definiteness all δ-punctured states are nonclassical.→ Negativity of the Wigner function. The Wigner func- Necessary condition. Using Eq. (7) with the set of tion of a vacuum-centered Gaussian-punctured thermal K state reads coherent states, we find that we can only subtract a Gaus- sian tighter than the original state (b < n¯) and obtain as  2 2  − 2|α| − 2|α| a necessary condition for positivity 2 e 1+2¯n e 1+2b W (α) =  w1  , (46) 2 N π 1 + 2¯n − 1 + 2b b + 1 − |α1| G NC w e n¯−b w , . (52) 1 6 n¯ + 1 ≡ max and takes negative values if w1 satisfies 2 Proof. Using γ α 2 = e−|α−γ| and Eq. (41), we find |h | i| 1 + 2b for the expectation value of ρ in the coherent state γ w1 > , (47) | i 1 + 2¯n Z 2 |α−α |2 1 − |α| −|α−γ|2 w1 − 1 −|α−γ|2 2 ργγ = e n¯ e b d α. which is a stronger condition than Eq. (45). Note that πn¯ − πb for b 0, Eq. (47) tends to Eq. (25), the bound found (53) for δ-punctured→ thermal states. Evaluation of the Gaussian integrals leads to

Antibunching. From the expression of the second-order 2 |γ−α |2 b + 1 − |γ| + 1 correlation function (10), we have ρ 0 w e n¯+1 b+1 . (54) γγ > ⇔ 1 6 n¯ + 1 2 2 (2) 2¯n w12b That needs to be true for all γ C in order to keep g = (1 w1) − 2 . (48) ∈ − (¯n w1b) the possibility of a positive semideﬁnite density opera- − tor. Minimizing this expression over γ gives us the upper G, NSC G, NC For w1 = wmax, vac, the condition for antibunching reads bound wmax . Since the exponential function is strictly monotonically increasing and the prefactors are positive g(2) < 1 n¯2 + b2 < 1, (49) we can focus on minimizing the exponent: ⇔ γ 2 γ α 2 To summarize, for α1 = 0, a P function that corre- f(γ) | | + | − 1| . sponds to a positive semideﬁnite density operator and ≡ −n¯ + 1 b + 1 9

By applying a rotation in phase-space, the puncture’s 1 1 center can always be brought along the real axis, so that 0.8 we can set α1 = x1 R. We then split γ into its real and 0.8 ∈ 0.6 G imaginary part γ = γR + iγI , leading to max

w 0.4 2 0.6 0.2 x1 f(γ) = h(γI ) + g(γR) + 0 G b + 1 max 0 0.5 1 1.5 2 w 0.4 α with | 1| 2 1 1 h(γ ) = γ , 0.2 I I b + 1 − n¯ + 1 2 1 1 2γRx1 g(γR) = γR . 0 b + 1 − n¯ + 1 − b + 1 0 0.5 1 1.5 2 α1 Here we have to distinguish two cases (the case b =n ¯ | | is trivial): i) If b > n¯, h(γ ) is an inverted parabola and G, NC I FIG. 5. Plot of the upper bound wmax (solid curve) and can attain arbitrary negative values and the same holds G of wmax, the w1 which cancels the smallest eigenvalue of the for f, thus pushing the exponential in (54) arbitrary close density matrix, as a function of |α1| for b = 1 and n¯ = 1.5 to 0, implying that the above condition can only be ful- (squares) and n¯ = 2 (dots). The eigenvalues are computed in filled for all coherent states if w1 = 0. This means that a Fock space of dimension nmax = 40. Inset: convergence of G we cannot subtract any Gaussian wider than the original wmax (green squares) as a function of the cut-off dimension thermal one and still obtain a valid state. This general- nmax (nmax = 8, 12, 40 from top to bottom). izes what we have already seen in the case of the vacuum centred Gaussian. ii) For b < n¯, the dependence of f which generalizes Eq. (45). Hence, the lower bound of the on the parameters γR, γI splits into the two functions h and g that depend on different independent arguments weight decreases with the distance between the center of the thermal state and the position of the puncture. γR and γI . Thus we can minimize both h and g sepa- Further this is smaller than wG, NC given in (52) since rately. Obviously h is minimized for γI = 0. We find max n¯ > b, implying that for all allowed values of α1, n¯ and that g(γR) is minimized for γR,min = (¯n + 1)x1/(¯n b) and we get the minimal values − b, there is a puncture weight such that the P function is negative somewhere in the complex plane, and the state 2 n¯ + 1 1 is nonclassical. g(γ min) = x , R, − 1 b + 1 n¯ b Negativity of the Wigner function. The Wigner func- 2 − tion of an arbitrarily-centred Gaussian-punctured ther- x1 f(γR,min) = . mal state reads −n¯ b −  2 2  − 2|α| − 2|α−α1| Finally, for the general case α1 C, plugging f(γR,min) 2 e 1+2¯n e 1+2b ∈ W (α) = N  w1  , (56) into Eq. (54) yields the bound (52). π 1 + 2¯n − 1 + 2b G, NC Tightness of wmax . As for the previous bounds obtained from the necessary condition (7), diagonalizing and is negative if w satisfies the density operator of the Gaussian punctured thermal 1 2 state [i.e. Eq. (1) together with Eqs. (2) and (41)] and us- 1 + 2b − |α1| G, NC w1 > e n¯−b , (57) ing the suggest that the bound wmax [Eq. (52)] is 1 + 2¯n tight. ThisNSC result is shown in Fig. 5, where we plotted the analytic bound (lines) and the numerical bound (squares a stronger condition than Eq. (55). Antibunching. The correlation function (10) for the and dots) for different cut-off dimensions nmax as a function of the center of the puncture α . Also, we found state (51) is given by | 1| that the coherent state γ with γ = (¯n + 1)/(¯n b) α1 is 2 4 2 2 | i − (2) 2¯n w1 α1 + 4 α1 b + 2b an eigenstate of ρ with eigenvalue 0 when w1 is given by g = (1 w1) − | | | | . (58) 2 2 the bound (52). − [¯n w1( α1 + b)] − | | Note that in the limit b 0, Eq. (58) tends to Eq. (26) → 2. Conditions for nonclassicality found in Sec. III for a δ-punctured thermal state. We have that for n¯ < 1, the minimum of g(2) is obtained for b 0. For n¯ 0 this minimum tends to 0. We also Negativity of the P function. From Eq. (51), we see have→ that for n¯→ > 1, the minimum of g(2) is obtained for that partly negative P functions can be obtained when a finite value of b. Figure 6 shows Eq. (58) as a function

|α |2 of and for and for the maximal puncture b − 1 n¯ b < n¯ α1 = 0.1 w1 > e n¯−b , (55) G, NC n¯ weight w1 = wmax given in Eq. (52). 10

procedure suitable for cavity-QED or circuit-QED that realizes a vacuum-removed state diagonal in the Fock 2 basis, based on the nondestructive detection of photons ’g2_gauss_thermal.dat’ u 1:2:($2<$1?$3:3) g(2) [41–45] in a mode of the cavity (See Fig. 7). Two-level 3 atoms in the cavity in a superposition + of their two | i 1.5 states g , e , where = ( g e )/√2, experience | i | i |±i | i ± | i a phase shift ϕ = π(n mod 2q) between g and e that 2 q | i | i b 1 depends on the number of photons in the cavity, with an integer q that can be controlled through the interaction 1 time and strength. A projective measurement of an 0.5 atom in the basis (corresponding to measurements results labelled ){|±i} updates our knowledge of the photon 0 number. Since ϕ±is deﬁned modulo 2π, 2q diﬀerent Fock 0 0 0.5 1 1.5 2 states can be distinguished. n¯

FIG. 6. Density plot of the second order correlation function However, we need no information which Fock state is g(2) as a function of n¯ and b < n¯ for single Gaussian-punctured actually realized, but only need to systematically exclude states with α1 = 0.1 for the maximal puncture weight w1 the vacuum state. This can be achieved with an iterative given in Eq. (52). The black dashed curve delimits the region procedure: In the first step one sets 0. If ’ ’ (2) q = 1 = 2 where g < 1. In the top left corner is shown the P function is measured, the state is tagged as part of the ensemble.− for the parameters indicated by a dot in the density plot. If ’+’ is measured, we still need to distinguish the states n with n = 0, 2, 4, 6, 8, ... in order to not truncate |morei states than the vacuum. So we set q = 2 = 21 π(n mod 4) VI. POSSIBLE EXPERIMENTAL yielding a phase shift of ϕ = 2 that allows us to REALIZATIONS OF PUNCTURED STATES distinguish states with n = 2+4k, k N0 photons (mea- ∈ surement result “-”) from states with n = 4k, k N0, ∈ A. Vacuum-removed states including the vacuum (measurement result “+”). This l−1 can be continued with q = 2 , l = 3, 4, ..., lmax, as long as needed to exclude the ambiguity of states with even A δ-punctured thermal state with puncture at α1 = 0 gaus,vac photon numbers as high as wished. and maximum weight of the puncture, b , i.e. a − NSC In general, setting q = 2l 1, with l N, we can vacuum-removed thermal state, definitely distinguish the vacuum from all∈ states with ∞ n n = q+2kq, k N0. If after a finite number of repetitions X n¯ n¯ + 1 1 ∈ 0 0 ρ = n n = ρ 0 0 , lmax all measurements yielded + we need to completely N (1 +n ¯)n+1 | ih | n¯ T − 1 +n ¯ | ih | n=1 discard the remaining state. This obviously causes (59) some errors, since we also exclude possible realizations is probably the example of a punctured state that is other than the vacuum from being further used. The easiest to realize experimentally. Since a thermal state smallest photon number that was not distinguished from has no coherences between Fock states, a projective the vacuum, and was therefore also subtracted in the lmax measurement in the Fock-basis does not destroy the remaining state, is then n0 = 2qmax = 2 , where qmax quantum properties of the state, but only alters its corresponds to the smallest interaction time that has statistical character. Discarding the ground-state when been realized, i.e. the last measurement. The whole lmax it is found leads to an ensemble of states that realizes procedure post-selects states n = 2 k k N0, lmax | i 6 | i∀ ∈ the vacuum-removed thermal state. This generalizes whereas states n = 2 k , k N0 are truncated. to any initial state that is diagonal in the Fock-basis. The procedure| onlyi works| becausei ∈ the thermal state When single-photon sources are available, one can is diagonal in the Fock-basis; all coherences between of course synthesize a given vacuum-removed state different Fock states are lost. In such a case, however, diagonal in the Fock basis from the beginning by mixing the procedure is very effective: If we want to make sure Fock-states with the corresponding statistical weight that states with photon number less than n0 were not without the need for any measurement. This generalizes truncated, we only need to perform lmax = log2(n0) to states where the vacuum is not removed completely different measurements. d e by adjusting the statistical weight of the vacuum-state to the desired value. For a thermal state, occupation for high photon num- The measurement-based procedure faces the problem ber states decays very rapidly for small temperatures, that traditional photon-counting methods destroy the and the overall error can be kept small. To quantify this, photons during the measurement. Here we propose a we can calculate the fidelity of our desired state (59) with 11

q = 1 q = 2 q = qmax 0 , 2 , ... 0 , 4 , ... 0 , 2qmax , ... | i | i | i | i | i | i

ρT

1 , 3 , ... 2 , 6 , ... qmax , 3qmax , ... | i | i | i | i | i | i

at least one meaurement yielded '-'

ρ0

FIG. 7. Sketch of realization of vacuum-removed thermal state. A thermal state ρT is subjected to multiple QND-measurements that distinguish the vacuum state of other photon number states. If at least one of the measurements yields 0−0, the vacuum has definitely not been realized and the measured state is added to the ensemble. If all measurements yielded 0+0 the state is discarded. Since we only perform a finite number of measurements, the final state ρ0 given in equation (60) can slightly deviate from the desired state ρ given in equation (59). The deviation can be quantified calculating the fidelity (61). respect to the state that was actually constructed with a = a : When a is measured one discards the 0 6 1 0 state, for a1 one continues the experiment and effectively " ∞ lmax # 2 n gets the desired punctured state through post-selection. 0 0 X n¯ lmax lmax ρ = ρT 2 n 2 n . N − (1 +n ¯)2lmax n+1 | ih | Recently a procedure was proposed that realizes the mea- n=0 surement operator A [48]. (60)

Since both density operators are diagonal in Fock basis C. Synthesizing arbitrary states the ﬁdelity is simply evaluated as [46, p. 409]:

∞ The above considerations are only valid for a com- 0 X p 0 1 1 F (ρ, ρ ) = ρnnρnn = 1 l . plete subtraction of the vacuum state. The procedure − n¯ n¯+1 2 max n=0 1 does not translate easily to more general subtractions, n¯ − (61) especially for α1 = 0, and we therefore present another, completely diﬀerent6 approach, based on approximately In the limit of large lmax this reduces to synthesizing our states from scratch as follows: We can

l numerically express the density operator of the desired 2 max − log( n¯+1 )2lmax 1 n¯ e n¯ punctured state up to a chosen dimension N in the Fock F = 1 = 1 , (62) − n¯ n¯ + 1 − n¯ basis, neglecting terms involving states with higher photon number than N: i.e. the ﬁdelity approaches exponentially with lmax . 1 2 N N This shows the theoretical viability of the procedure. 0 X X ρ = c n m . Experimental imperfections, as for example the non per- mn| ih | n=0 m=0 fect linearity of the phase shift in the photon numbers [43], may worsen the result and will have to be eval- This is in general a mixed state and we can express it as uated for a given experimental setup. a sum of pure states: X ρ0 = p ψ ψ i| iih i| B. “Vacuum or not”-measurements i=1 N X The vacuum-state can be removed from a broader class ψ = c(i) n . | ii n | i of states, namely states that do not contain coherences n=0 between the vacuum state and other Fock states. Indeed, We can then create our state by choosing with the right if ρ = 0 in the Fock basis, then ρ = ρ = 0 is implied nn 0n n0 probability one of the pure states and constructing it by the positivity of the state [47]. Then one can remove with the method proposed in [49–52] and demonstrated the vacuum-state with a measurement operator A of the experimentally in [53]. Multiple repetition (creation of an form adequate ensemble of states) then gives us our state up ∞ X to a certain accuracy. In a single run we actually do not A = a 0 0 + a n n . (63) have the complete state, but for further experiments one 0| ih | 1 | ih | n=1 anyway needs typically multiple repetitions for obtaining 12 measurement statistics of any observable. Hence, this lation function can be modiﬁed through an appropriate procedure creates the desired mixed state in the ensemble choice of the punctures and identify the regimes yielding sense. antibunching of light. All states exhibiting antibunching have a negative P function in some region. Finally, we presented possible experimental realizations of punctured states based on vacuum-or-not measurements and VII. CONCLUSION on complete synthesizing.

In this work, we introduced a novel class of nonclassical states with regular non-positive Glauber-Sudarshan ACKNOWLEDGMENTS P function that we call punctured states. These states are obtained from the addition of sufficiently narrow neg- FD acknowledges fruitful discussions with Daniel K. L. ative peaks to smooth positive P functions. We deter- Oi. FD would like to thank the F.R.S.-FNRS for financial mined the regimes of parameters yielding proper phys- support. Computational resources have been provided ical (i.e. positive semidefinite) nonclassical states in the by the Consortium des Équipements de Calcul Intensif case of δ- and Gaussian punctures of thermal or squeezed (CÉCI), funded by the F.R.S.-FNRS under Grant No. thermal states. We showed that their second-order corre- 2.5020.11.

[1] D. V. Strekalov and G. Leuchs, Nonlinear interac- tron. 10, 1112 (1980). tions and nonclassical light, in Advances in Quantum [13] G. Brassard, N. Lütkenhaus, T. Mor, and B. C. Sanders, Photonics: from the First Single-photon and Nonlin- Limitations on practical quantum cryptography, Phys. ear Optical Experiments to Modern Quantum Photon- Rev. Lett. 85, 1330 (2000). ics, edited by Robert W. Boyd, Svetlana G. Lukishova, [14] H. Zheng, D. J. Gauthier, and H. U. Baranger, Decoy- and Victor N. Zadkov (Springer, New York, in press), state quantum key distribution with nonclassical light arXiv:1701.01403. generated in a one-dimensional waveguide, Opt. Lett. 38, [2] A. Pathak and A. Ghatak, Classical light vs. nonclassi- 622 (2013). cal light: characterizations and interesting applications, [15] S. D. Bartlett, B. C. Sanders, Requirement for quantum Journal of Electromagnetic Waves and Applications 32, computation, J. Mod. Opt. 50, 2331 (2003). 229 (2018). [16] R. Chiao and J. Garrison, Quantum optics (Oxford Uni- [3] M. Xiao, L.-A. Wu, and H. J. Kimble, Precision mea- versity Press, 2008). surement beyond the shot-noise limit, Phys. Rev. Lett. [17] E. C. G. Sudarshan, Equivalence of Semiclassical and 59, 278 (1987). Quantum Mechanical Descriptions of Statistical Light [4] L. K. Shalm, R. B. A. Adamson, and A. M. Steinberg, Beams, Phys. Rev. Lett. 10, 277 (1963). Squeezing and over-squeezing of triphotons, Nature (Lon- [18] G. S. Agarwal, Quantum Optics (Cambridge University don) 457, 67 (2009). Press, 2012). [5] T. Ono, J. Sabines-Chesterking, H. Cable, J. L. O’Brien, [19] T. Kiesel and W. Vogel, Nonclassicality filters and quasi- and J. C. F. Matthews, Optical implementation of spin probabilities, Phys. Rev. A 82, 032107 (2010). squeezing, New J. Phys. 19, 053005 (2017). [20] T. Kiesel, W. Vogel, B. Hage, and R. Schnabel, Direct [6] A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. sampling of negative quasiprobabilities of a squeezed state, P. Williams, and J. P. Dowling, Quantum interferometric Phys. Rev. Lett. 107, 113604 (2011). optical lithography - Exploiting entanglement to beat the [21] E. Agudelo, J. Sperling, and W. Vogel, Quasiprobabilities diffraction limit, Phys. Rev. Lett. 85, 2733 (2000). for multipartite quantum correlations of light, Phys. Rev. [7] D. V. Strekalov and J. P. Dowling, Two-Photon Interfer- A 87, 033811 (2013). ometry for High-Resolution Imaging, J. Mod. Opt. 49, [22] B. Kühn and W. Vogel, Visualizing nonclasscal effects in 519 (2002). phase space, Phys. Rev. A 90, 033821 (2014). [8] M. Genovese, Real applications of quantum imaging, J. [23] D. N. Klyshko, Observable signs of nonclassical light, Opt. 18, 073002 (2016). Phys. Lett. A 213, 7 (1996). [9] M. A. Taylor, J. Janousek, V. Daria, J. Knittel, B. Hage, [24] A. Brauer, Limits of the characteristic roots of a matrix, H.-A. Bachor, and Warwick P. Bowen, Subdiffraction- Duke Mathematical Journal 13, 387 (1946). Limited Quantum Imaging within a Living Cell, Phys. [25] P. Shivakumar, J. Williams, and N. Rudraiah, Eigenval- Rev. X 4, 011017 (2014). ues for infinite matrices, Linear Algebra and its Appli- [10] B. Ya. Zel’dovich and D. N. Klyshko, Statistics of field cations 96, 35 (1987). in parametric luminescence, Pis’ma Zh. Eksp. Teor. Fiz. [26] D. T. Smithey, M. Beck, M. G. Raymer, and A. Faridani, 9, 69 (1969) [JETP Lett. 9, 40 (1969)]. Measurement of the Wigner distribution and the density [11] D. C. Burnham and D. L. Weinberg, Observation of matrix of a light mode using optical homodyne tomog- simultaneity in parametric production of optical photon raphy: Application to squeezed states and the vacuum, pairs, Phys. Rev. Lett. 25, 84 (1970). Phys. Rev. Lett. 70, 1244 (1993). [12] D. N. Klyshko, Use of two-photon light for absolute cali- [27] S. Wallentowitz and W. Vogel, Unbalanced homodyning bration of photoelectric detectors, Sov. J. Quantum Elec- for quantum-state measurements, Phys. Rev. A 53, 4528 13

(1996). [43] C. Guerlin, J. Bernu, S. Deleglise, C. Sayrin, S. Gleyzes, [28] K. Banaszek and K. Wódkiewicz, Direct probing of quan- S. Kuhr, M. Brune, J.-M. Raimond, and S. Haroche, Pro- tum phase space by photon counting, Phys. Rev. Lett. 76, gressive field-state collapse and quantum non-demolition 4344 (1996). photon counting, Nature (London) 448, 889 (2007). [29] L. G. Lutterbach and L. Davidovich, Method for direct [44] Y. F. Xiao, S. K. Özdemir, V. Gaddam, C. H. Dong, N. measurement of the Wigner function in cavity QED and Imoto and L. Yang, Quantum nondemolition measure- ion traps, Phys. Rev. Lett. 78, 2547 (1997). ment of photon number via optical Kerr effect in an ultra- [30] K. Banaszek, C. Radzewicz, K. Wódkiewicz, and J. S. high-Q microtoroid cavity, Opt. Exp. 16, 21462 (2008). Krasiński, Direct measurement of the Wigner function [45] F. Helmer, M. Mariantoni, E. Solano, and F. Marquardt, by photon counting, Phys. Rev. A 60, 674 (1999). Quantum nondemolition photon detection in circuit QED [31] P. Bertet, A. Auffeves, P. Maioli, S. Osnaghi, T. Meunier, and the quantum Zeno effect, Phys. Rev. A 79, 052115 M. Brune, J. M. Raimond, and S. Haroche, Direct Mea- (2009). surement of the Wigner Function of a One-Photon Fock [46] M. A. Nielsen and I. L. Chuang, Quantum computation State in a Cavity, Phys. Rev. Lett. 89, 200402 (2002). and quantum information (Cambridge University Press, [32] W. P. Schleich, “Quantum Optics in Phase space”, Wiley- 2010). VCH Verlag Berlin GmbH, Berlin (2001). [47] R. A. Horn and C. R. Johnson, Matrix Analysis (Cam- [33] L. Mandel, Sub-Poissonian photon statistics in resonance bridge University Press, New York, 1985). fluorescence, Opt. Lett. 4, 205 (1979). [48] D. K. L. Oi, V. Potoček, and J. Jeffers, Nondemolition [34] H. J. Kimble, M. Dagenais, and L. Mandel, Photon An- Measurement of the Vacuum State or its Complement, tibunching in Resonance Fluorescence, Phys. Rev. Lett. Phys. Rev. Lett. 110, 210504 (2013). 39, 691 (1977). [49] M. Hofheinz, E. M. Weig, M. Ansmann, R. C. Bialczak, [35] W. Vogel, Nonclassical correlation properties of radiation E. Lucero, M. Neeley, A. D. O’Connell, H. Wang, J. M. fields, Phys. Rev. Lett. 100, 013605 (2008). Martinis, and A. N. Cleland, Generation of Fock states [36] J. L. W. V. Jensen, Sur les fonctions convexes et les in- in a superconducting quantum circuit, Nature (London) égalités entre les valeurs moyennes, Acta Math. 30, 175 454, 310 (2008). (1906). [50] A. Ben-Kish, B. DeMarco, V. Meyer, M. Rowe, J. Brit- [37] M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- ton, W. M. Itano, B. M. Jelenković, C. Langer, D. bridge University Press, 2001). Leibfried, T. Rosenband, and D. J. Wineland, Experi- [38] M. S. Kim, F. A. M. de Oliveira, and P. L. Knight, Prop- mental Demonstration of a Technique to Generate Ar- erties of squeezed number states and squeezed thermal bitrary Quantum Superposition States of a Harmonically states, Phys. Rev. A 40, 2494 (1989). Bound Spin-1/2 Particle, Phys. Rev. Lett. 90, 037902 [39] S. M. Barnett and P. M. Radmore, Methods in Theoret- (2003). ical Quantum Optics (Oxford University Press, Oxford, [51] C. K. Law and J. H. Eberly, Arbitrary Control of a Quan- 2002). tum Electromagnetic Field, Phys. Rev. Lett. 76, 1055 [40] L. Diósi, Comment on "Nonclassical States: An Observ- (1996). able Criterion", Phys. Rev. Lett. 85, 2841 (2000). [52] Y.-X. Liu, L.-F. Wei, and F. Nori, Generation of non- [41] G. Gabrielse, S. Peil, B. Odom, and B. D’Urso, QND ob- classical photon states using a superconducting qubit in servation of quantum jumps between Fock states: a one- a quantum electrodynamic microcavity, Europhys. Lett. electron cyclotron oscillator at 70 mK to 4.2 K, in Sum- 67, 941 (2004) maries of Papers Presented at the Quantum Electronics [53] M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, E. and Laser Science Conference, 1999. QELS’99. Technical Lucero, M. Neeley, A. O’Connell, D. Sank, J. Wenner, J. Digest. (IEEE,1992) pp. 19-20. M. Martinis, and A. N. Cleland, Synthesizing arbitrary [42] A. Reiserer, S. Ritter, and G. Rempe, Non-destructive quantum states in a superconducting resonator, Nature detection of an optical photon, Science 342, 1349 (2013). (London) 459, 546 (2009).