Identification of nonclassical properties of light with multiplexing layouts

TL;DR

This paper extends a detector-independent method to identify nonclassical light properties using multiplexing layouts, deriving a general theory and applying it to experimental data with superconducting sensors to reveal new nonclassicality criteria.

Contribution

It introduces a comprehensive theoretical framework for nonclassicality detection in multiplexed setups and demonstrates its effectiveness on experimental data, including the novel concept of sub-multinomial light.

Findings

Sub-multinomial light outperforms sub-Poisson and sub-binomial criteria

The theory applies to arbitrary detectors and multiplexing configurations

Experimental data confirms the validity of the new nonclassicality criteria

Abstract

In a recent contribution, we introduced and applied a detector-independent method to uncover nonclassicality. Here, we extend those techniques and give more details on the performed analysis. We derive a general theory of the positive-operator-valued measure that describes multiplexing layouts with arbitrary detectors. From the resulting quantum version of a multinomial statistics, we infer nonclassicality probes based on a matrix of normally ordered moments. We discuss these criteria and apply the theory to our data which are measured with superconducting transition-edge sensors. Our experiment produces heralded multi-photon states from a parametric down-conversion light source. We show that the known notions of sub-Poisson and sub-binomial light can be deduced from our general approach, and we establish the concept of sub-multinomial light, which is shown to outperform the former two concepts of nonclassicality for our data.

Figures

Equations

$$\hat{\rho }=\int d^2\alpha P(\alpha )|\alpha ⟩⟨\alpha |.$$

$$⟨:{\hat{f}}^{†}\hat{f}:⟩\stackrel{cl.}{\ge }0,$$

$${\hat{a}}_{m,\mathrm{out}}=U_{m,1}{\hat{a}}_{1,\mathrm{in}}+\cdots +U_{m,N}{\hat{a}}_{N,\mathrm{in}}.$$

$$|\alpha ⟩\otimes |0{⟩}^{\otimes (N-1)}\stackrel{U(N)}{⟼}|\alpha /\sqrt{N}{⟩}^{\otimes N}.$$

$$p_k(\alpha )=⟨\alpha /\sqrt{N}|:{\hat{\pi }}_k^{\prime }:|\alpha /\sqrt{N}⟩=⟨\alpha |:{\hat{\pi }}_k:|\alpha ⟩,$$

$$p_{k_1}(\alpha )\cdots p_{k_N}(\alpha )=⟨\alpha |:{\hat{\pi }}_{k_1}\cdots {\hat{\pi }}_{k_N}:|\alpha ⟩,$$

$$\begin{array}{rl}{p}_{(k_1,\dots ,k_N)}=& \int d^2\alpha P(\alpha )p_{k_1}(\alpha )\cdots p_{k_N}(\alpha )\\ =& ⟨:{\hat{\pi }}_{k_1}\cdots {\hat{\pi }}_{k_N}:⟩.\end{array}$$

$$\begin{array}{rl}& c_{(N_0,\dots ,N_K)}=\mathrm{tr}[\hat{\rho }{\hat{\Pi }}_{(N_0,\dots ,N_K)}]\\ =& ⟨:\frac{N!}{N_0!\cdots N_K!}{\hat{\pi }}_0^{N_0}\cdots {\hat{\pi }}_K^{N_K}:⟩,\end{array}$$

$$\begin{array}{rl}{c}_{(N_0,\dots ,N_K)}=& \int d^2\alpha P(\alpha )\frac{N!}{N_0!\cdots N_K!}\\ & \times p_0(\alpha {)}^{N_0}\cdots p_K(\alpha {)}^{N_K}.\end{array}$$

$$\hat{f}=\sum _{m_0+\cdots +m_K\le N/2}{f}_{m_0,\dots ,m_K}{\hat{\pi }}_0^{m_0}\cdots {\hat{\pi }}_K^{m_K}.$$

$$\begin{array}{rl}⟨:{\hat{f}}^{†}\hat{f}:⟩=& \sum _{\begin{array}{c}{m}_0+\cdots +m_K\le N/2\\ m_0^{\prime }+\cdots +m_K^{\prime }\le N/2\end{array}}{f}_{m_0,\dots ,m_K}^{\ast }\\ & \times ⟨:{\hat{\pi }}_0^{m_0+m_0^{\prime }}\cdots {\hat{\pi }}_K^{m_K+m_K^{\prime }}:⟩f_{m_0^{\prime },\dots ,m_K^{\prime }}\\ =& {\vec{f}}^{†}M\vec{f},\end{array}$$

$$0\nleqq M.$$

$$\begin{array}{rl}0\nleqq & M^{(2)}={\left(⟨:\Delta {\hat{\pi }}_k\Delta {\hat{\pi }}_{k^{\prime }}:⟩\right)}_{k,k^{\prime }=0,\dots ,K}\\ & =\left(\begin{array}{ccc}⟨:(\Delta {\hat{\pi }}_0{)}^2:⟩& \dots & ⟨:(\Delta {\hat{\pi }}_0)(\Delta {\hat{\pi }}_K):⟩\\ ⋮& \ddots & ⋮\\ ⟨:(\Delta {\hat{\pi }}_0)(\Delta {\hat{\pi }}_K):⟩& \dots & ⟨:(\Delta {\hat{\pi }}_K{)}^2:⟩\end{array}\right).\end{array}$$

$$\begin{array}{rl}& g(z_0,\dots ,z_N)=\overline{z_0^{N_0}\cdots z_K^{N_K}}\\ =& \sum _{N_0+\cdots +N_K=N}{c}_{(N_0,\dots ,N_K)}{z}_0^{N_0}\cdots z_K^{N_K}\\ =& ⟨:{\left(z_0{\hat{\pi }}_0+\cdots +z_K{\hat{\pi }}_K\right)}^N:⟩.\end{array}$$

$${{\partial }_{z_0}^{m_0}\cdots {\partial }_{z_K}^{m_K}g(z_0,\dots ,z_K)|}_{z_0=\cdots =z_K=1}$$

$$=$$

$$⟨:{\hat{\pi }}_k:⟩=\frac{\overline{N_k}}{N}\text{ and }⟨:{\hat{\pi }}_k{\hat{\pi }}_{k^{\prime }}:⟩=\frac{\overline{N_k{N}_{k^{\prime }}}-{\delta }_{k,k^{\prime }}\overline{N_k}}{N(N-1)}$$

$$⟨:\Delta {\hat{\pi }}_k\Delta {\hat{\pi }}_{k^{\prime }}:⟩=\frac{N\overline{\Delta N_k\Delta N_{k^{\prime }}}-\overline{N_k}\left(N{\delta }_{k,k^{\prime }}-\overline{N_{k^{\prime }}}\right)}{N^2(N-1)}.$$

$$0\stackrel{cl.}{\le }$$

$$={\left(N\overline{\Delta N_k\Delta N_{k^{\prime }}}-\overline{N_k}\left[N{\delta }_{k,k^{\prime }}-\overline{N_{k^{\prime }}}\right]\right)}_{k,k^{\prime }=0,\dots ,K}.$$

$${\vec{f}}^{\mathrm{T}}{M}^{(2)}\vec{f}<0,$$

$$:\hat{\mu }:=f_0:{\hat{\pi }}_0:+\cdots +f_K:{\hat{\pi }}_K:.$$

$$⟨:(\Delta \hat{\mu }{)}^2:⟩<0.$$

$$Q_{\mathrm{multi}}=\underset{\vec{f}:{\vec{f}}^{\mathrm{T}}\vec{f}=1}{\min }{\vec{f}}^{\mathrm{T}}{M}^{(2)}\vec{f}={\vec{f}}_0^{\mathrm{T}}{M}^{(2)}{\vec{f}}_0,$$

$$\begin{array}{rl}& \frac{(N-1){\vec{f}}^{\mathrm{T}}{M}^{(2)}\vec{f}}{⟨:{\hat{\pi }}_0:⟩(1-⟨:{\hat{\pi }}_0:⟩)}=\frac{N\overline{(\Delta B{)}^2}-N\overline{B}+{\overline{B}}^2}{(N-\overline{B})\overline{B}}\\ =& N\frac{\overline{(\Delta B{)}^2}}{\overline{B}(N-\overline{B})}-1=Q_{\mathrm{bin}},\end{array}$$

$$A=\sum _{k=0}^{K}k{N}_k.$$

$$⟨:\hat{\mu }:⟩=\sum _{k=0}^{K}k\frac{\overline{N_k}}{N}=\frac{\overline{A}}{N}.$$

$$⟨:(\Delta \mu {)}^2:⟩={\vec{f}}^{\mathrm{T}}{M}^{(2)}\vec{f}$$

$$=$$

$$-\frac{\left({\sum }_{k=0}^K{k}^2\frac{\overline{N_k}}{N}\right)-{\left({\sum }_{k=0}^{K}k\frac{\overline{N_k}}{N}\right)}^2-\left({\sum }_{k=0}^{K}k\frac{\overline{N_k}}{N}\right)}{N-1}.$$

$$\frac{⟨:(\Delta \mu {)}^2:⟩}{⟨:\mu :⟩}=$$

Full text

I Introduction

The bare existence of photons highlights the particle nature of electromagnetic waves in quantum optics E05 . Therefore, the generation and detection of photon states are crucial for a comprehensive understanding of fundamental concepts in quantum physics; see Refs. BC10 ; CDKMS14 for recent reviews on single photons. Beyond this scientific motivation, the study of nonclassical radiation fields is also of practical importance. For instance, quantum communication protocols rely on the generation and detection of photons GT07 ; S09 . Yet, unwanted attenuation effects—which are always present in realistic scenarios—result in a decrease of the nonclassicality of a produced light field. Conversely, an inappropriate detector model can introduce fake nonclassicality even to a classical radiation field SV11 ; GLLSSMK11 ; SVA12a . For this reason, we seek robust and detector-independent certifiers of nonclassicality Setal2017 .

The basic definition of nonclassicality is that a quantum state of light cannot be described in terms of classical statistical optics. A convenient way to represent general states is given in terms of the Glauber-Sudarshan $P$ function S63 ; G63 . Whenever this distribution cannot be interpreted in terms of classical probability theory, the thereby represented state is a nonclassical one TG86 ; M86 . A number of nonclassicality tests have been proposed; see Ref. MBWLN10 for an overview. Most of them are formulated in terms of matrices of normally ordered moments of physical observables; see, e.g., Ref. SRV05 . For example, the concept of nonclassical sub-Poisson light M79 can be written and even generalized in terms of matrices of higher-order photon-number correlations AT92 . Other matrix-based nonclassicality tests employ the Fourier or Laplace transform of the Glauber-Sudarshan $P$ function RV02 ; SVA16 .

In order to apply such nonclassicality probes, one has to measure the light field under study with a photodetector S07 ; H09 . The photon statistics of the measured state can be inferred if the used detector has been properly characterized. This can be done by a detector tomography LS99 ; AMP04 ; LKKFSL08 ; LFCPSREPW09 ; ZDCJEPW12 —i.e., measuring a comparably large number of well-defined probe states to construct a detection model. Alternatively, one can perform a detector calibration BCDGMMPPP12 ; PHMH12 ; BKSSV17 —i.e., the estimation of parameters of an existing detection model with some reference measurements. Of particular interest are photon-number-resolving detectors of which superconducting transition-edge sensors (TESs) are a successful example LMN08 ; Getal11 ; BCDGLMPRTP12 ; RFZMGDFE12 ; ZCDPLJSPW12 . Independent of the particular realization, photon-number-resolving devices allow for the implementation of quantum tasks, such as state reconstruction LCGS10 ; BGGMPTPOP11 , imaging LBFD08 ; CWB14 , random number generation ATDYRS15 , and the characterization of sources of nonclassical light WDSBY04 ; HPHP05 ; FL13 ; APHAB16 —even in the presence of strong imperfections TKE15 . Moreover, higher-order ALCS10 ; AOB12 ; SBVHBAS15 , spatial BDFL08 ; MMDL12 ; CTFLMA16 , and conditional SBDBJDVW16 quantum correlations have been studied.

So far, we did not distinguish between the detection scheme and the actual detectors. That is, one has to discern the optical manipulation of a signal field and its interaction with a sensor which yields a measurement outcome. Properly designed detection layouts of such a kind render it possible to infer or use properties of quantum light without having a photon-number-resolution capability ZABGGBRP05 ; PADLA10 ; KV16 or they do not require a particular detector model CKS14 ; AGSB16 . For instance, multiplexing layouts with a number of detectors that can only discern between the presence (“on”) or absence (“off”) of absorbed photons can be combined into a photon-number-resolving detection device PTKJ96 ; KB01 ; ASSBW03 ; FJPF03 ; RHHPH03 ; CDSM07 ; SPDBCM07 . Such types of schemes use an optical network to split an incident light field into a number of spatial or temporal modes of equal intensities which are subsequently measured with on/off detectors. The measured statistics is shown to resemble a binomial distribution SVA12a rather than a Poisson statistics, which is obtained for photoelectric detection models KK64 ; see also Refs. I14 ; PZA16 in this context. For such detectors, the positive-operator-valued measure (POVM), which fully describes the detection layout, has been formulated SVA12a ; MSB16 . Recently, the combination of a multiplexing scheme with multiple TESs has been used to significantly increase the maximal number of detectable photons HBLNGS15 .

Based on the binomial character of the statistics of multiplexing layouts with on/off detectors, the notion of sub-binomial light has been introduced SVA12 and experimentally demonstrated BDJDBW13 . It replaces the concept of sub-Poisson light M79 , which applies to photoelectric counting models KK64 , for multiplexing arrangements using on/off detectors. Nonclassical light can be similarly inferred from multiplexing devices with non-identical splitting ratios LFPR16 . In addition, the on-chip realization of optical networks MGHPGSWS15 can be used to produce integrated detectors to verify sub-binomial light HSPGHNVS16 .

In this paper, we derive the quantum-optical click-counting theory for multiplexing layouts which employ arbitrary detectors. Therefore, we formulate nonclassicality tests in terms of normally ordered moments, which are independent of the detector response. This method is then applied to our experiment which produces heralded multi-photon states. Our results are discussed in relation with other notions of nonclassical photon correlations.

In Ref. Setal2017 , we study the same topic as we do in this paper from a classical perspective. There, the treatment of the detector-independent verification of quantum light is performed solely in terms of classical statistical optics. Here, however, we use a complementary quantum-optical perspective on this topic. Beyond that, we also consider higher-order moments of the statistics, present additional features of our measurements, and compare our results with previously known nonclassicality tests as well as simple theoretical models.

This paper is organized as follows. In Sec. II, the theoretical model for our detection layout is elaborated and nonclassicality criteria are derived. The performed experiment is described in Sec. III with special emphasis on the used TESs. An extended analysis of our data, presented in Sec. IV, includes the comparison of different forms of nonclassicality. We summarize and conclude in Sec. V.

II Theory

In this section, we derive the general, theoretical toolbox for describing the multiplexing arrangement with arbitrary detectors and for formulating the corresponding nonclassicality criteria. The measurement layout under study is shown in Fig. 1. Our detection model shows that for any type of employed detector, the measured statistics can be described in the form of a quantum version of a multinomial statistics [Eq. (8)]. This leads to the formulation of nonclassicality criteria in terms of negativities in the normally ordered matrix of moments [Eq. (12)]. Especially, covariance-based criteria are discussed and related to previously known forms of nonclassicality.

II.1 Preliminaries

We apply well-established concepts in quantum optics in this section. Namely, any quantum state of light $\hat{\rho}$ can be written in terms of the Glauber-Sudarshan representation G63 ; S63 ,

[TABLE]

From this diagonal expansion in terms of coherent states $|\alpha\rangle$, one observes that one can formulate the detection theory in terms of coherent states. A subsequent integration over the $P$ function then describes the model for any state. Furthermore, the definition of nonclassicality is also based on this representation. Namely, the state $\hat{\rho}$ is a classical state if and only if $P$ can be interpreted in terms of classical probability theory TG86 ; M86 , i.e., $P(\alpha)\geq 0$. Whenever this cannot be done, $\hat{\rho}$ refers to a nonclassical state.

Moreover, the $P$ function of a state is related to the normal ordering (denoted by ${:}\cdots{:}$) of measurement operators. For a detailed introduction to bosonic operator ordering, we refer to Ref. AW70 . It can be shown in general that any classical state obeys VW06

[TABLE]

for any operator $\hat{f}$. In addition, we may recall that expectation values of normally ordered operators and coherent states can be simply computed by replacing the bosonic annihilation $\hat{a}$ and creation operator $\hat{a}^{\dagger}$ with the coherent amplitude $\alpha$ and its complex conjugate $\alpha^{\ast}$, respectively. A violation of constraint (2) necessarily identifies nonclassicality, which will be also used to formulate our nonclassicality criteria.

II.2 Multiplexing detectors

The optical detection scheme under study, shown in Fig. 1, consists of a balanced multiplexing network which splits a signal into $N$ modes. Those outputs are measured with $N$ identical detectors which can produce $K+1$ outcomes, labeled as $k=0,\ldots,K$. Let us stress that we make a clear distinction between the well-characterized optical multiplexing, the individual and unspecified detectors, and the resulting full detection scheme.

In the multiplexing part, a coherent-state input $|\alpha\rangle$ is distributed over $N$ output modes. Further on, we have vacuum $|\mathrm{vac}\rangle=|0\rangle$ at all other $N-1$ input ports. In general, the $N$ input modes—defined by the bosonic annihilation operators $\hat{a}_{n,\mathrm{in}}$ ($\hat{a}_{1,\mathrm{in}}=\hat{a}$)—are transformed via the unitary $U(N)=(U_{m,n})_{m,n=1}^{N}$ into the output modes

[TABLE]

Taking the balanced splitting into account, it holds that $|U_{m,n}|=1/\sqrt{N}$. Adjusting the phases of the outputs properly, we get the following input-output relation

[TABLE]

Note that a balanced, but lossy network similarly yields $|\tau\alpha\rangle\otimes\cdots\otimes|\tau\alpha\rangle$ for $\tau\leq 1/\sqrt{N}$.

For describing the detector, we do not make any specifications. Nevertheless, we will be able to formulate nonclassicality tests. The probability $p_{k}$ for the $k$th measurement outcome ($0\leq k\leq K$) for any type of detector can be written in terms of the expectation value of the POVM operators ${:}\hat{\pi}^{\prime}_{k}{:}$, $p_{k}=\langle{:}\hat{\pi}^{\prime}_{k}{:}\rangle$. Note that any operator can be written in a normally ordered form AW70 ; VW06 and that the POVM includes all imperfections of the individual detector, such as the quantum efficiency or nonlinear responses. For the coherent states $|\alpha/\sqrt{N}\rangle$, we have

[TABLE]

whereby we also define ${:}\hat{\pi}_{k}{:}$ in terms of ${:}\hat{\pi}^{\prime}_{k}{:}$ through the mapping $\hat{a}\mapsto\hat{a}/\sqrt{N}$.

We find that for a measurement with our $N$ detectors and our coherent output state (4), the probability to measure the outcome $k_{n}$ with the $n$th detector—more rigorously a coincidence $(k_{1},\ldots,k_{N})$ from the $N$ individual detectors—is given by

[TABLE]

where we used the relation $\langle\alpha|{:}\hat{A}{:}|\alpha\rangle\langle\alpha|{:}\hat{B}{:}|\alpha\rangle=\langle\alpha|{:}\hat{A}\hat{B}{:}|\alpha\rangle$ for any two (or more) operators $\hat{A}$ and $\hat{B}$ and Eq. (5). The Glauber-Sudarshan representation (1) allows one to write for any quantum state $\hat{\rho}$

[TABLE]

So far we studied the individual parts, i.e., the optical multiplexing and the $N$ individual detectors, separately. To describe the full detection scheme in Fig. 1, we need some additional combinatorics, which is fully presented in Appendix A. There, the main idea is that one can group the individual detectors into subgroups of $N_{k}$ detectors which deliver the same outcome $k$. Suppose the individual detectors yield the outcomes $(k_{1},\ldots,k_{N})$. Then, $N_{k}$ is the number of individual detectors for which $k_{n}=k$ holds. In other words, $(N_{0},\ldots,N_{K})$ describes the coincidence that $N_{0}$ detectors yield the outcome [math], $N_{1}$ detectors yield the outcome $1$, etc. Note that the total number of detectors is given by $N=N_{0}+\cdots+N_{K}$.

The POVM representation $\hat{\Pi}_{(N_{0},\dots,N_{K})}$ for the event $(N_{0},\ldots,N_{K})$ is given in Eq. (A). In combination with Eq. (7), we get for the detection layout in Fig. 1 the click-counting statistics of a state $\hat{\rho}$ as

[TABLE]

which is a normal-ordered version of a multinomial distribution. The click-counting statistics (8) yields the probability that $N_{0}$ times the outcome $k=0$ together with $N_{1}$ times the outcome $k=1$, etc., is recorded with the $N$ individual detectors. Using Eq. (1), we can rewrite the click-counting distribution,

[TABLE]

In this form, we can directly observe that any classical statistics, $P(\alpha)\geq 0$, is a classical average over multinomial probability distributions.

II.3 Higher-order nonclassicality criteria

Our click-counting model (8) describes a multiplexing scheme and applies to arbitrary detectors. One observes that its probability distribution is based on normally ordered expectation values of the form $\langle{:}\hat{\pi}_{0}^{m_{0}}\cdots\hat{\pi}_{K}^{m_{K}}{:}\rangle$. Hence, we can formulate nonclassicality criteria from inequality (2) while expanding

[TABLE]

This operator is chosen such that it solely includes the operators that are actually measured. We can write

[TABLE]

with a vector $\vec{f}=(f_{m_{0},\ldots,m_{K}})_{(m_{0},\ldots,m_{K})}$, using a multi-index notation, and the matrix of normally ordered moments $M$, which is defined in terms of the elements $\langle{:}\hat{\pi}_{0}^{m_{0}+m^{\prime}_{0}}\cdots\hat{\pi}_{K}^{m_{K}+m^{\prime}_{K}}{:}\rangle$. Also note that the order of the moments is bounded by the number of individual detectors, $N\geq m_{0}+\cdots+m_{K}+m^{\prime}_{0}+\cdots+m^{\prime}_{K}$, as the measured statistics (8) only allows for retrieving them.

As the non-negativity of the expression (11) holds for classical states [condition (2)] and for all coefficients $\vec{f}$, we can equivalently write the following: A state is nonclassical if

[TABLE]

Conversely, the matrix of higher-order, normal-ordered moments $M$ is positive semidefinite for classical light. Note, it can be also shown (Appendix A in Ref. LSV15 ) that the matrix of normally ordered moments can be equivalently expressed in a form that is based on central moments, $\langle{:}(\Delta\hat{\pi}_{0})^{m_{0}+m_{0}^{\prime}}\cdots(\Delta\hat{\pi}_{K})^{m_{K}+m_{K}^{\prime}}{:}\rangle$.

For example and while restricting to the second-order submatrix, we get nonclassicality conditions in terms of normal-ordered covariances,

[TABLE]

The relation $\langle{:}\hat{\pi}_{K}{:}\rangle=1-[\langle{:}\hat{\pi}_{0}{:}\rangle+\cdots+\langle{:}\hat{\pi}_{K-1}{:}\rangle]$ of general POVMs implies that the last row of $M^{(2)}$ is linearly dependent on the other ones. This further implies that zero is an eigenvalue of $M^{(2)}$. Hence, we get for any classical state that the minimal eigenvalue of this covariance matrix is necessarily zero.

In order to relate our nonclassicality criteria to the measurement of the click-counting statistics (8), let us consider the generating function, which is given by

[TABLE]

The derivatives of the generating function relate the measured moments with the normally ordered ones,

[TABLE]

for $m_{0}+\cdots+m_{K}\leq N$ and $(x)_{m}=x(x-1)\cdots(x-m+1)=x!/(x-m)!$ being the falling factorial. Having a closer look at the second and third line of Eq. (15), we see that the factorial moments $\overline{\left(N_{0}\right)_{m_{0}}\cdots\left(N_{K}\right)_{m_{K}}}$ can be directly sampled from $c_{(N_{0},\ldots,N_{K})}$. From the last two lines of Eq. (15) follows the relation to the normally ordered moments, which are needed for our nonclassicality tests.

II.4 Second-order criteria

As an example and due to its importance, let us focus on the first- and second-order moments in detail. In addition, our experimental realization implements a single multiplexing step, $N=2$, which yields a restriction to second-order moments [see comment below Eq. (11)]. As a special case of Eq. (15), we obtain

[TABLE]

for $k,k^{\prime}\in\{0,\ldots,K\}$. Hence, our covariances are alternatively represented by

[TABLE]

As the corresponding matrix (13) of normal-ordered moments is nonnegative for classical states, we get

[TABLE]

The violation of this specific constraint for classical states has been experimentally demonstrated for the generated quantum light Setal2017 .

Let us consider other special cases of the general criterion. In particular, let us study the projections that result in a nonclassicality condition

[TABLE]

see also Eqs. (2) and (13). Note that $M^{(2)}$ is a real-valued and symmetric $(K+1)\times(K+1)$ matrix. Thus, it is sufficient to consider real-valued vectors $\vec{f}=(f_{0},\ldots,f_{K})^{\mathrm{T}}$. Further on, let us define the operator

[TABLE]

Then, we can also read condition (19) as

[TABLE]

That is, the fluctuations of the observable ${:}\hat{\mu}{:}$ are below those of any classical light field. In the following, we consider specific choices for $\vec{f}$ to formulate different nonclassicality criteria.

II.4.1 Sub-multinomial light

The minimization of (19) over all normalized vectors yields the minimal eigenvalue $Q_{\mathrm{multi}}$ of $M^{(2)}$. That is

[TABLE]

where $\vec{f}_{0}$ is a normalized eigenvector to the minimal eigenvalue. If we have $M^{(2)}\ngeq 0$, then we necessarily get $Q_{\mathrm{multi}}<0$. For classical states, we get $Q_{\mathrm{multi}}=0$; see the discussion below Eq. (13). As this criterion exploits the maximal negativity from covariances of the multinomial statistics, we refer to a radiation field with $Q_{\mathrm{multi}}<0$ as sub-multinomial light.

II.4.2 Sub-binomial light

We can also consider the vector $\vec{f}=(0,1,\ldots,1)^{\mathrm{T}}$, which yields ${:}\hat{\mu}{:}=\hat{1}-{:}\hat{\pi}_{0}{:}$. Hence, we have effectively reduced our system to a detection with a binary outcome, represented through the POVMs ${:}\hat{\pi}_{0}{:}$ and ${:}\hat{\mu}{:}=\hat{1}-{:}\hat{\pi}_{0}{:}$. Using a proper scaling, we can write

[TABLE]

defining $B=N_{1}+\cdots+N_{K}=N-N_{0}$ and using Eq. (16). The condition $Q_{\mathrm{bin}}<0$ defines the notion of sub-binomial light SVA12 and is found to be a special case of inequality (19).

II.4.3 Sub-Poisson light

Finally, we study criterion (19) for $\vec{f}=(0,1,\ldots,K)^{\mathrm{T}}$. We have ${:}\hat{\mu}{:}=\sum_{k=0}^{K}k{:}\hat{\pi}_{k}{:}$ and we also define

[TABLE]

Their mean values are related to each other,

[TABLE]

We point out that $\overline{N_{k}}/N$ can be also interpreted as probabilities, being nonnegative $\overline{N_{k}}/N\geq 0$ and normalized $1=\overline{N_{0}}/N+\cdots+\overline{N_{K}}/N$ since $N=N_{0}+\cdots+N_{K}$. Further, we can write the normally ordered variance (21) in the form

[TABLE]

Again, we can use a proper, nonnegative scaling to find

[TABLE]

The parameters $Q_{\mathrm{Pois}}$ and $Q_{\mathrm{Pois}}^{\prime}$, often denoted as the Mandel or $Q$ parameter, relate to the notion of sub-Poisson light M79 . However, we have a difference of two such Mandel parameters in Eq. (27). The second parameter $Q_{\mathrm{Pois}}^{\prime}$ can be considered as a correction, because the statistics of $A$ is only in a rough approximation a Poisson distribution. This is further analyzed in Appendix B.

II.5 Discussion

We derived the click-counting statistics (8) for unspecified POVMs of the individual detectors. This was achieved by using the properties of a well-defined multiplexing scheme. We solely assumed that the $N$ detectors (with $K+1$ possible outcomes) are described by the same POVM. A deviation from this assumption can be treated as a systematic error; see Supplemental Material to Ref. Setal2017 . The full detection scheme was shown to result in a quantum version of multinomial statistics. This also holds true for an infinite, countable ($K=|\mathbb{N}|$) or uncountable ($K=|\mathbb{R}|$) set of outcomes, for which any measurement run can only deliver a finite sub-sample. For coherent light $|\alpha_{0}\rangle$, we get a true multinomial probability distribution; see Eq. (9) for $P(\alpha)=\delta(\alpha-\alpha_{0})$. For a binary outcome, $K+1=2$, we retrieve a binomial distribution SVA12a ; SVA13 , which applies, for example, to avalanche photodiodes in the Geiger mode HSPGHNVS16 ; SBVHBAS15 ; SBDBJDVW16 or superconducting nanowire detectors BKSSV17 ; BKSSV17atm .

Further on, we derived higher-order nonclassicality tests which can be directly sampled from the data obtained from the measurement layout in Fig. 1. Then, we focused on the second-order nonclassicality probes and compared the cases of sub-multinomial [Eq. (22)], sub-binomial [Eq. (23)], and (corrected) sub-Poisson [Eq. (27)] light. The latter notion is related to nonclassicality in terms of photon-number correlation functions (see also Ref. ZM90 ) and is a special case of our general criteria. Additionally, our method can be generalized to multiple multiplexing-detection arrangements to include multimode correlations similar to the approach in Ref. SVA13 .

Recently, another interesting application was reported to characterize spatial properties of a beam profile with multipixel cameras CTFLMA16 . There, the photon-number distribution itself is described in terms of a multinomial statistics, and the Mandel parameter can be used to infer nonclassical light. Here, we show that an balanced multiplexing and any measurement POVM yield a click-counting statistics—describing a different statistical quantity than the photon statistics of a beam profile—in the form of a quantum version of a multinomial distribution leading to higher-order nonclassicality criteria. We also demonstrated that in some special scenarios (Appendix B), a relation between the click statistic and photon statistics can be retrieved which is, however, much more involved in the general case; see also Sec. III.2.

III Experiment

Before applying the derived techniques to our data, we describe the experiment and study some features of our individual detectors in this section. Especially, the response of our detectors is shown to have a nonlinear behavior which underlines the need for our nonclassicality criteria which are applicable to any type of detector. Additional details can be found in Appendix C.

III.1 Setup description and characterization

An outline of our setup is given in Fig. 2. It is divided into a source that produces heralded photon states and a detection stage which represents one multiplexing step. In total, we use three superconducting TESs. For generating correlated photons, we employ a spontaneous parametric down-conversion (PDC) source. Here, we describe the individual parts in some more detail.

III.1.1 The PDC source

Our spontaneous PDC source is a waveguide-written periodically poled potassium titanyl phosphate (PP-KTP) crystal which is $8\,\mathrm{mm}$ long. The type-II spontaneous PDC process is pumped with laser pulses at $775\,\mathrm{nm}$ and a full width at half maximum (FWHM) of $2\,\mathrm{nm}$ at a repetition rate of $75\,\mathrm{kHz}$. The heralding idler mode has a horizontal polarization and it is centered at $1554\,\mathrm{nm}$. The signal mode is vertically polarized and centered at $1546\,\mathrm{nm}$. A PBS spatially separates the output signal and idler pulses. An edge filter discards the pump beam. In addition, the signal and idler are filtered by $3\,\mathrm{nm}$ bandpass filters. This is done in order to filter out the broadband background which is typically generated in dielectric nonlinear waveguides ECMS11 . In general, such PDC sources have been proven to be well-understood and reliable sources of quantum light LCS09 ; KHQBSS13 . Hence, we may focus our attention on the employed detectors.

III.1.2 The TES detectors

We use superconducting TESs as our photon detectors LMN08 . These TESs are micro-calorimeters consisting of $25\,\mathrm{\mu m}\times 25\,\mathrm{\mu m}\times 20\,\mathrm{nm}$ slabs of tungsten located inside an optical cavity with a gold backing mirror designed to maximize absorption at $1500\,\mathrm{nm}$. They are secured within a ceramic ferule as part of a self-aligning mounting system, so that the fiber core is well aligned to the center of the detector MLCVGN11 . The TESs are first cooled below their transition temperature within a dilution refrigerator and then heated back up to their transition temperature by Joule heating caused by a voltage bias, which is self-stabilized via an electro-thermal feedback effect I95 . Within this transition region, the steep resistance curve ensures that the small amount of heat deposited by photon absorption causes a measurable decrease in current flowing through the device. After photon absorption, the heat is then dissipated to the environment via a weak thermal link to the TES substrate.

To read out the signal from this photon absorption process, the current change—produced by photon absorption in the TES—is inductively coupled to a superconducting quantum interference device (SQUID) module where it is amplified, and this signal is subsequently amplified at room temperature. This results in complex time-varying signals of about $5\,\mathrm{\mu s}$ duration. These signals are sent to a digitizer to perform fast analog-to-digital conversion, where the overlap with a reference signal is computed and then binned. This method allows us to process incoming signals at a speed of up to $100\,\mathrm{kHz}$.

Our TESs are installed in a dilution refrigerator operating at a base temperature of about $70\,\mathrm{mK}$ and a cooling power of $400\,\mathrm{\mu W}$ at $100\,\mathrm{mK}$. One of the detectors has a measured detection efficiency of $0.98^{+0.02}_{-0.08}$ HMGHLNNDKW15 . The other two TESs have identical efficiencies within the error of our estimation.

III.2 Detector response analysis

Even though we will not use specific detector characteristics for our analysis of nonclassicality, it is nevertheless scientifically interesting to study their response. This will also outline the complex behavior of superconducting detectors. For the time being, we ignore the detection events of the TESs 1 and 2 in Fig. 2 and solely focus on the measurement of the heralding TES.

In Fig. 3, the measurement outcome of those marginal counts is shown. A separation into disjoint energy intervals represents our outcomes $k\in\{0,\ldots,11\}$ (see also Appendix C). The distribution around the peaked structures can be considered as fluctuations of the discrete energy levels (indicated by vertical dark green, solid lines). We observe that the difference between two discrete energies $E_{n}$ is not constant as one would expect from $E_{n+1}-E_{n}=\hbar\omega$, which will be discussed in the next paragraph. In addition, the marginal photon statistics should be given by a geometric distribution for the two-mode squeezed-vacuum state produced by our PDC source; see Appendix D. In the logarithmic scaling in Fig. 3, this would result in a linear function. However, we observe a deviation from such a model; compare light green, dashed and dot-dashed lines in Fig. 3.

This deviation from the expected, linear behavior could have two origins: The source is not producing a two-mode squeezed-vacuum state (affecting the height of the peaks), or the detector, including the SQUID response, is not operating in a linear detection regime (influence on the horizontal axis). To counter the latter, the measured peak energies $E_{n}$—relating to the photon numbers $n$—have been fitted by a quadratic response function $n=aE_{n}^{2}+bE_{n}+c$; see the inset in Fig. 4. As a result of such a calibration, the peaked structure is well described by a linear function in $n$ for the heralding TES as shown in Fig. 4 (top), which is now consistent with the theoretical expectation. The same nonlinear energy transformation also yields a linear $n$ dependence for the TESs 1 and 2 (cf. Fig. 4, bottom). Note that those two detectors only allow for a resolution of $K+1=8$ outcomes and that these two detectors have indeed a very similar response—the depicted linear function is identical for both. In conclusion, it is more likely that the measured nonlinear behavior in Fig. 4 can be assigned to the detectors, and the PDC source is operating according to our expectations.

In summary, we encountered an unexpected, nonlinear behavior of our data. To study this, a nonlinear fit was applied. This allowed us to make some predictions about the detector response in the particular interval of measurement while using known properties of our source. However, a lack of such extra knowledge prevents one from characterizing the detector. In Sec. II, we have formulated nonclassicality tests which are robust against the particular response function of the individual detectors. They are accessible without any prior detector analysis and include the eventuality of nonlinear detector responses and other imperfections, such as quantum efficiency. With this general treatment, we also avoid the time-consuming detector tomography.

IV Application

In this section, we apply the general theory, presented in Sec. II, to our specific experimental arrangement, shown in Fig. 2. In the first step, we perform an analysis to identify nonclassicality which can be related to photon-number-based approaches. In the second step, we also compare the different criteria for sub-multinomial, sub-binomial, and sub-Poisson light for different realizations of our multi-photon states.

IV.1 Heralded multi-photon states

As derived in Appendix B, the connection of the operator (20), for $f_{k}=k$, to the photon-number statistics for the idealized scenario of photoelectric detection POVMs is given by ${:}\hat{\mu}{:}=(\eta/N)\hat{n}$, where $\eta$ is the quantum efficiency of the individual detectors. This also relates—in this ideal case—the quantities

[TABLE]

Recalling ${:}\hat{n}{:}=\hat{n}$, we see that $\langle{:}\hat{\mu}{:}\rangle$ is proportional to the mean photon number in this approximation. Similarly, we can connect $\langle{:}(\Delta\hat{\mu})^{2}{:}\rangle$ to the normally ordered photon-number fluctuations. They are non-negative for classical states and negative for sub-Poisson light [see Eq. (21)].

An ideal PDC source is known to produce two-mode squeezed-vacuum states,

[TABLE]

where $|q|<1$. One mode can be used to produce multi-photon states by conditioning to the $l$th outcome of the heralding detector. Using photoelectric detector POVMs, we get the following mean value and the variances (Appendix D):

[TABLE]

with a transformed squeezing parameter $\tilde{\lambda}=(1-\tilde{\eta})|q|^{2}$ and $\tilde{\eta}$ being the efficiency of the heralding detector. Note, we get the ideal $l$th Fock state, $|l\rangle$, for $\tilde{\lambda}\to 0$.

The experimental result is shown in the top panel of Fig. 5. Using Eq. (16), we directly sampled the mean value and the variance of ${:}\hat{\mu}{:}=0{:}\hat{\pi}_{0}{:}+\cdots+K{:}\hat{\pi}_{K}{:}$ from the measured statistics for a heralding with $l=0,\ldots,5$. In this plot, $l$ increases from left to right relating to the increased mean photon numbers (including attenuations) of the heralded multi-photon states. The idealized theoretical modeling [see Eq. (30)] is shown in the bottom part of Fig. 5. Note, details of the error analysis have been formulated previously in Ref. Setal2017 .

From the variances, we observe no nonclassicality when heralding to the [math]th outcome, which is expected as we condition on vacuum. In contrast, we can infer nonclassicality for the conditioning to higher outcomes of the heralding TES, $\langle{:}(\Delta\mu)^{2}{:}\rangle<0$ for $l>0$. We have a linear relation between the normally ordered mean and variance of ${:}\hat{\mu}{:}$, which is consistent with the theoretical prediction in Eq. (30). In the ideal case, the normal-ordered variance of the photon number for Fock states also decreases linearly with increasing $l$, $\langle l|\hat{n}|l\rangle=l$ and $\langle l|{:}(\Delta\hat{n})^{2}{:}|l\rangle\rangle=-l$. It is also obvious that the errors are quite large for the verification of nonclassicality with this particular test for sub-Poisson light. We will discuss this in more detail in the next subsection.

IV.2 Varying pump power

So far, we have studied measurements for a single pump power of the PDC process. However, the purity of the heralded states depends on the squeezing parameter, which is a function of the pump power. For instance, in the limit of a vanishing squeezing, we have the optimal approximation of the heralded state to a Fock state. However, the rate of the probabilistic generation converges to zero in the same limit (Appendix D). Hence, we have additionally generated multi-photon states for different squeezing levels. The results of our analysis are shown in Fig. 6 and will be discussed in the following.

Suppose we measure the counts $C_{l}$ for the $l$th outcome of the heralding TES. The efficiency of generating this $l$th heralded state reads

[TABLE]

From the model in Appendix D, we expect that

[TABLE]

The efficiency decays exponentially with $l$ and the decay is stronger for smaller squeezing or pump power—i.e., a decreasing $|q|^{2}$. In the left column of Fig. 6, we can observe this behavior. It can be seen in all other parts of Fig. 6 that $\eta_{\mathrm{gen}}$ influences the significance of our results. A smaller $\eta_{\mathrm{gen}}$ value naturally implies a larger error because of a decreased sample size $C_{l}$. This holds for increasing $l$ and for decreasing squeezing.

In the second column in Fig. 6, labeled as “sub-Poisson”, we study the nonclassicality criterion

[TABLE]

$N=2$, and $K=7$, which is related to sub-Poisson light (Sec. II.4.3). The third column in Fig. 6 correspondingly shows “sub-binomial” light (Sec. II.4.2),

[TABLE]

The last column, “sub-multinomial”, depicts the nonclassicality criterion

[TABLE]

where $\vec{f}_{0}$ is a normalized eigenvector to the minimal eigenvalue of $M^{(2)}$ (Sec. II.4.1).

For all notations of nonclassicality under study, the heralding to the [math]th outcome is consistent with our expectation of a classical state, which also confirms that no fake nonclassicality is detected. For instance, applying the Mandel parameter to the data of this [math]th heralded stated without the corrections derived here [Eq. (27)], we would observe a negative value; see also similar discussions in Refs. SVA12a ; SBVHBAS15 . The case of a Poisson or binomial statistics tends to be above zero, whereas the multinomial case is consistent with the value of zero. This expectation has been justified below Eq. (13).

A lot of information on the quantum-optical properties of the generated multi-photon ($l>0$) light fields can be concluded from Fig. 6. Let us mention some of them by focusing on a comparison. We have the trend that the notion of sub-Poisson light has the least significant nonclassicality. This is due to the vector $\vec{f}$ [Eq. (33)], which assigns a higher contribution to the larger outcome numbers. However, those contributions have lower count numbers, which consequently decreases the statistical significance. As depicted in Fig. 6, this effect is not present for sub-binomial light, which is described by a more or less balanced weighting of the different counts; see vector $\vec{f}$ in Eq. (34). Still, this vector is fixed.

The optimal vector is naturally computed by the sub-multinomial criterion in Eq. (35). The quality of the verified nonclassicality is much better than for the other two scenarios of sub-Poisson and sub-binomial light in most of the cases. Let us mention that the normalized eigenvector to the minimal eigenvalue of the sampled matrix $M^{(2)}$ typically, but not necessarily, yields the minimal propagated error. Additionally, a lower squeezing level allows for the heralding of a state which is closer to an ideal Fock state. This results in higher negativities for decreasing squeezing and fixed outcomes $l$ in Fig. 6. However, the heralding efficiency $\eta_{\mathrm{gen}}$ is also reduced, which results in a larger error.

Finally, we may point out that this comparative analysis of sub-Poisson, sub-binomial, and sub-multinomial light from data of a single detection arrangement would not be possible without the technique that has been elaborated in this paper (Sec. II).

V Summary

In summary, we constructed the quantum-optical framework to describe multiplexing schemes that employs arbitrary detectors and to verify nonclassicality of generated multi-photon states. We formulated the theory of such a detection layout together with nonclassicality tests. Further, we set up an experimental realization and applied our technique to the data.

In a first step, the theory was formulated. We proved that the measured click-counting statistics of the scheme under study is always described by a quantum version of the multinomial statistics. In fact, for classical light, this probability distribution can be considered as a mixture of multinomial statistics. This bounds the minimal amount of fluctuations which can be observed for classical radiation fields. More precisely, the matrix of higher-order, normally ordered moments, which can be directly sampled from data, can exhibit negative eigenvalues for nonclassical light. As a particular example, we discussed nonclassicality tests based on the second-order covariance matrix, which led to establishing the concept of sub-multinomial light. Previously studied notions of nonclassicality, i.e, sub-Poisson and sub-binomial light, have been found to be special cases of our general nonclassicality criteria.

In our second part, the experiment was analyzed. Our source produces correlated photon pairs by a parametric-down-conversion process. A heralding to the outcome of a detection of the idler photons with a transition-edge sensor produced multi-photon states in the signal beam. A single multiplexing step was implemented with a subsequent detection by two transition-edge sensors to probe the signal field. The complex function of these detectors was discussed by demonstrating their nonlinear response to the number of incident photons. Consequently, without worrying about this unfavorable feature, we applied our robust nonclassicality criteria to our data. We verified the nonclassical character of the produced quantum light. The criterion of sub-multinomial light was shown to outperform its Poisson and binomial counterparts to the greatest possible extent.

In conclusion, we presented a detailed and more extended study of our approach in Ref. Setal2017 . We formulated the general positive-operator-valued measure and generalized the nonclassicality tests to include higher-order correlations which become more and more accessible with an increasing number of multiplexing steps. In addition, details of our data analysis and a simple theoretical model were considered. Thus, we described a robust detection scheme to verify quantum correlations with unspecified detectors and without introducing fake nonclassicality.

Acknowledgements.

The project leading to this application has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 665148. A. E. is supported by EPSRC EP/K034480/1. J. J. R. is supported by the Netherlands Organization for Scientific Research (NWO). W. S. K is supported by EPSRC EP/M013243/1. S. W. N., A. L., and T. G. are supported by the Quantum Information Science Initiative (QISI). I. A. W. acknowledges an ERC Advanced Grant (MOQUACINO). The authors thank Johan Fopma for technical support. The authors gratefully acknowledge helpful comments by Tim Bartley and Omar Magaña-Loaiza. Contributions of this work by NIST, an agency of the U.S. Government, are not subject to U.S. copyright.

Appendix A Combinatorics and POVM elements

Here, we provide the algebra that is needed to get from Eq. (7) to Eq. (8). More rigorously, we use combinatorial methods to formulate the POVM $\hat{\Pi}_{(N_{0},\ldots,N_{K})}$ in terms of the POVM ${:}\hat{\pi}_{k_{1}}\cdots\hat{\pi}_{k_{N}}{:}$. Say $N_{k}$ is the number of elements of $(k_{1},\ldots,k_{N})$ which take the value $k$. Then, $(N_{0},\ldots,N_{K})$ describes the coincidence that $N_{0}$ detectors yield the outcome [math], $N_{1}$ detectors yield the outcome $1$, etc. One specific and ordered measurement outcome is defined by $(k_{0,1},\ldots,k_{0,N})$, with

[TABLE]

which results in a given $(N_{0},\ldots,N_{K})$, where the total number of detectors is $N=N_{0}+\cdots+N_{K}$. This specific example can be used to represent all similar outcomes as we will show now. The $(k_{1},\ldots,k_{N})$ for the same combination $(N_{0},\ldots,N_{K})$ can be obtained from $(k_{0,\sigma(1)},\ldots,k_{0,\sigma(N)})$ via a permutation $\sigma\in\mathcal{S}_{N}$ of the elements. Here $\mathcal{S}_{N}$ denotes the permutation group of $N$ elements which has a cardinality of $N!$. Note that all permutations $\sigma$ which exchange identical outcomes result in the same tuple. This means for the outcome defined in Eq. (40) that $(k_{0,\sigma(1)},\ldots,k_{0,\sigma(N)})=(k_{0,1},\ldots,k_{0,N})$ for any permutation of the form $\sigma\in\mathcal{S}_{N_{0}}\times\cdots\times\mathcal{S}_{N_{K}}$. Therefore, the POVM element for a given $(N_{0},\ldots,N_{K})$ can be obtained by summing over all permutations $\sigma\in\mathcal{S}_{N}$ of the POVMs of individual outcomes ${:}\hat{\pi}_{k_{0,1}}\cdots\hat{\pi}_{k_{0,N}}{:}$ [Eq. (7)] while correcting for the $N_{0}!\cdots N_{K}!$ multi-counts. More rigorously, we can write

[TABLE]

where relations of the form ${:}\hat{A}\hat{B}\hat{A}{:}={:}\hat{A}^{2}\hat{B}{:}$ have been used.

Appendix B Corrected Mandel parameter

For the nonclassicality test in Sec. II.4.3, we could assume a detector which can discriminate $K=\infty$ measurement outcomes, which are related to measurement operators of a Poisson form, ${:}\hat{\pi}_{k}^{\prime}{:}={:}\hat{\Gamma}^{k}e^{-\hat{\Gamma}}{:}/k!$ KK64 , where $\hat{\Gamma}=\eta\hat{n}$ is an example of a linear detector response function ($\eta$ quantum efficiency). Using the definition (5), we get ${:}\hat{\pi}_{k}{:}={:}(\hat{\Gamma}/N)^{k}e^{-\hat{\Gamma}/N}{:}/k!$, where the denominator $N$ accounts for the splitting into $N$ modes SVA13 . This idealized model yields $\langle{:}\hat{\mu}{:}\rangle=\langle{:}(\hat{\Gamma}/N){:}\rangle$ and

[TABLE]

Hence, we have $Q_{\mathrm{Pois}}=\langle{:}(\Delta\hat{\Gamma})^{2}{:}\rangle/\langle{:}\hat{\Gamma}{:}\rangle=NQ_{\mathrm{Pois}}^{\prime}$ and

[TABLE]

Thus, we have shown that for photoelectric detection models, we retrieve the notation of sub-Poisson light, $Q_{\mathrm{Pois}}<0$, from the general form (27), which includes a correction term.

Appendix C Binning and measured coincidences

The data in Fig. 4 (Sec. III) are grouped in disjoint intervals around the peaks, representing the photon numbers. They define the outcomes $k=0,\ldots,K$. Because we are free in the choice of the intervals, we studied different scenarios and found that the given one is optimal from the information-theoretic perspective. On the one hand, if the current intervals are divided into smaller ones, we distribute the data of one photon number among several outcomes. This produces redundant information about this photon number. On the other hand, we have a loss of information about the individual photon numbers if the interval stretches over multiple photon numbers. This explains our binning as shown in Fig. 4.

An example of a measured coincidence statistics for outcomes $(k_{1},k_{2})$ is shown in Fig. 7. There, we consider a state which is produced by the simplest conditioning to the [math]th outcome of the heralding TES.

Based on this plot, let us briefly explain how these coincidences for $(k_{1},k_{2})$ result in the statistics $c_{(N_{0},\dots,N_{K})}$ for $(N_{0},\ldots,N_{K})$ and $K=7$. The counts on the diagonal, $k_{1}=k_{2}=k$, of the plot yield $c_{(N_{0},\dots,N_{K})}$ for $N_{k}=2$ and $N_{k^{\prime}}=0$ for $k^{\prime}\neq k$. For example, the highest counts are recorded for $(k_{1},k_{2})=(0,0)$ in Fig. 7 which gives $c_{(2,0,\ldots,0)}$ when normalized to all counts. Off-diagonal combinations, $k_{1}\neq k_{2}$, result in $c_{(N_{0},\dots,N_{K})}$ for $N_{k_{1}}=N_{k_{2}}=1$ and $N_{k}=0$ otherwise. For example, the normalized sum of the counts for $(k_{1},k_{2})\in\{(0,1),(1,0)\}$ yields $c_{(1,1,0,\ldots,0)}$. As we have $N=2$ TESs in our multiplexing scheme and $N_{0}+\cdots+N_{K}=N$, the cases $k_{1}=k_{2}$ and $k_{1}\neq k_{2}$ already define the full distribution $c_{(N_{0},\dots,N_{K})}$.

The asymmetry in the counting statistics between the two detectors results in a small systematic error $\lesssim 1\%$. One should keep in mind that the counts are plotted in a logarithmic scale. For all other measurements of heralded multi-photon states, this error is in the same order Setal2017 .

Appendix D Simplified theoretical model

Let us analytically compute the quantities which are used for the simplified description of physical system under study. The PDC source produces a two-mode squeezed-vacuum state (29), where the first mode is the signal and the second mode is the idler or herald. In our idealized model, the heralding detector is supposed to be a photon-number-resolving detector with a quantum efficiency $\tilde{\eta}$. A multiplexing and a subsequent measurement with $N$ photon-number-resolving detectors ($K=\infty$) are employed for the click counting. Each of the photon-number-resolving detector’s POVM elements is described by

[TABLE]

In addition, we will make use of the relations ${:}e^{y\hat{n}}{:}=(1+y)^{\hat{n}}$ (cf., e.g., Ref. SVA14 ) and

[TABLE]

For this model, we can conclude that the two-mode generating function for the considered two-mode squeezed-vacuum state reads

[TABLE]

where $z\in[0,1]$ relates the heralding mode and the components of $\vec{x}\in[0,1]^{N}$ (recall that $\|\vec{x}\|=\sum_{n}x_{n}$) to the outcomes of the $N$ detectors in the multiplexing scheme. From this generating function, we directly deduce the different properties that are used in this paper for comparing the measurement with our model. The needed derivatives are

[TABLE]

where $k=\|\vec{k}\|$, $x^{\prime}=1-\eta+\eta\|\vec{x}\|/N$, and $z^{\prime}=1-\tilde{\eta}+\tilde{\eta}z$. It is also worth mentioning that the case $N=1$ yields the result for photon-number-resolving detection without multiplexing.

The marginal statistics of the heralding detector reads

[TABLE]

The marginal statistics of the $n$th detector is

[TABLE]

In addition, the case of no multiplexing ($N=1$ and $x\cong\vec{x}$) yields for the $l$th heralded state the following first and second normally ordered photon numbers:

[TABLE]

with $\tilde{\lambda}=(1-\tilde{\eta})|q|^{2}$. The corresponding photon distribution (i.e., for $\eta=1$) of the $l$th multi-photon state reads

[TABLE]

For $\tilde{\lambda}\to 0$, we have $\tilde{p}_{k|l}=\delta_{k,l}$, which is the photon statistics of the $l$th Fock state.