Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Correlations for subsets of particles in symmetric states: what photons are doing within a beam of light when the rest are ignored

Open Access Open Access

Abstract

Given a state of light, how do its properties change when only some of the constituent photons are observed and the rest are neglected (traced out)? By developing formulas for mode-agnostic removal of photons from a beam, we show how the expectation value of any operator changes when only q photons are inspected from a beam, ignoring the rest. We use this to re-express expectation values of operators in terms of the state obtained by randomly selecting q photons. Remarkably, this only equals the true expectation value for a unique value of q: expressing the operator as a monomial in normally ordered form, q must be equal to the number of photons annihilated by the operator. A useful corollary is that the coefficients of any q-photon state chosen at random from an arbitrary state are exactly the qth-order correlations of the original state; one can inspect the intensity moments to learn what any random photon will be doing and, conversely, one need only look at the n-photon subspace to discern what all of the nth-order correlation functions are. The astute reader will be pleased to find no surprises here, only mathematical justification for intuition. Our results hold for any completely symmetric state of any type of particle with any combination of numbers of particles and can be used wherever bosonic correlations are found.

1. Introduction

Photodetection and photon statistics have been the backbone of quantum optics since its inception [114]. For both single- and multi-mode states, correlation functions are routinely used to characterize the quantum or nonclassical nature of the field [1532], fundamentally requiring simultaneous detection of multiple photons.

Photons, like all bosons, are totally symmetric under particle exchange, so a photodetector should be agnostic as to which photons it perceives. Still, a detector seldom registers all of the photons in a beam, especially since beams tend to possess photon-number uncertainty, making the properties of subsets of photons from a beam crucial to its characterization. How do relevant quantities such as photon correlations change when some number of photons are removed from a beam [answered in Eq. (9)] or when only a certain number of photons from the beam are detected [Eq. (15)]? How many photons must be registered to learn about particular properties of the beam [Eqs. (16) and (18)]? Do the measured properties differ when different numbers of photons are detected [Eq. (17)]? All of these questions have intuitive answers that follow directly from the framework we detail here.

A familiar starting point for our investigation is the correspondence between the Poincaré and Bloch spheres. Classically, a quasi-monochromatic beam of light is represented by the “Stokes vector” pointing somewhere within the unit sphere named after Poincaré, while a qubit state such as a single photon’s polarization state is represented by a “Bloch vector” lying within Bloch’s sphere. Intuition says that a single photon chosen at random from a classical beam of light should have its Bloch vector correspond to the Stokes vector of the entire beam, and this is indeed the case [33]. What happens when more than one photon is selected, especially in the case of quantum polarization where there is more polarization information encoded beyond the Stokes vector [3436]? How does this intuition generalize when there are more than two polarization modes available to the photons, such as spatial or spectral modes? Not only are all of these questions answered here, they also do not depend on whether or not the beam of light has a determinate number of photons to begin with.

The polarization case also exemplifies the crucial role of correlations in quantum optics: the Stokes vector and the beam’s intensity are in one-to-one correspondence with the set of first-order correlations between the two polarization modes. These are the first-order degrees of coherence $g^{(1)}$ that are relevant for any pair of modes and are fundamental to experiments as basic as Young’s double slit [37]. Since the correspondence implies that single photons inspected from a beam are sufficient for revealing $g^{(1)}$ functions, how many photons are required for measuring an arbitrary $g^{(n)}$ function? Could any $q$ photons be inspected to learn the $g^{(1)}$ functions? Our simple yet rigorous framework uniquely answers all of these questions in one fell swoop.

Beyond intuition and mathematical completeness, the formulas we develop are useful for tasks such as quantifying quantumness in optical or other bosonic systems. Knowing how states behave after removing a certain number of photons is crucial to finding optimal quantum rotosensors [3841] and to studies of photon loss in multimode systems such as boson sampling devices [42]. These tools can be added to the quantum optician’s arsenal for a variety of applications.

Nomenclatural preliminaries are required. The removal of photons from a beam of light in different contexts give rise to the terms “loss” and “photon subtraction” that typically describe scenarios different from the present one. Photon subtraction, for example, is the result of acting on a state with an annihilation operator, which has a significant history [4359] including experimental demonstrations that conditionally split off a single photon from a state at a weakly reflective beam splitter [47]. That action may indeed leave the overall state intensity unchanged, or even increased, to appear as though no photons have been “subtracted.” Loss is the result of light coupling linearly with another mode that is later inaccessible; modelled using a beam splitter transformation, it predictably reduces the intensity of input beams yet requires probabilistic combinations to handle situations in which definite numbers of photons have partial loss. Here, in contrast, we deal with beams that have an exact number of photons removed or that remove all photons until an exact number remain, which may be the result of “discarding” or “ignoring” or “neglecting” photons, or “inspecting” or “focusing on” a subset of photons, mathematically described by “tracing out” particles in a first-quantized picture. In the framework of states with a fixed number of photons, our scenario has actually been described as “particle loss” [60,61] and as a “fixed-loss model” [62], terms we avoid due to their conflict with standard descriptions of optical loss. We thus re-emphasize: photons are indeed lost from our states and the resulting number of photons is indeed subtracted from the initial number, yet we refer to removal of photons from a beam in order to distance ourselves from the typical outcomes of loss or photon subtraction.

2. Removing Photons From a Beam

We begin by considering a generic state ${\hat {\rho }}_N$ that contains exactly $N$ photons and asking what happens when one is removed. Considering each photon to potentially belong to one of $d$ modes, such a state is typically expanded in the Fock basis as

$${\hat{\rho}}_N=\sum_{|\mathbf{m}|=|\mathbf{n}|=N}{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\vert{\mathbf{m}}\rangle \langle{\mathbf{n}}\vert.$$
The basis is comprised of states
$$\vert{\mathbf{n}}\rangle=\frac{\hat{a}_1^{{\dagger} n_1}}{\sqrt{n_1!}}\frac{\hat{a}_2^{{\dagger} n_2}}{\sqrt{n_2!}}\cdots \frac{\hat{a}_d^{{\dagger} n_d}}{\sqrt{n_d!}} \vert{\mathrm{vac}\rangle}$$
with $|\mathbf {n}|\equiv \sum _{i=1}^d n_i$, where the creation operators satisfy the usual bosonic commutation relations $[\hat {a}_i,\hat {a}^\dagger _j]=\delta _{ij}$. To discuss removal of a photon, we must rewrite each basis state in a first-quantized picture, with a completely symmetric sum over having $n_i$ of the photons being in state $\vert {i}\rangle$:
$$\vert{\mathbf{n}}\rangle=\sqrt{\frac{1}{\binom{N}{\mathbf{n}}}}\sum_{\mathrm{permutations}}\vert{1}\rangle^{{\otimes} n_1} \otimes\vert{2}\rangle^{{\otimes} n_2}\otimes\cdots\otimes\vert{d}\rangle^{{\otimes} n_d}.$$
All of the photons are in a symmetric state, so we fiducially choose to remove the first one. With this equivalence, we can immediately see that projecting the first photon onto some state $\vert{i}\rangle$ just lowers the index $n_i$ by one in the second-quantized notation:
$$\begin{aligned}(\langle i|\otimes\mathbb{1}^{{\otimes} N-1})\vert{\mathbf{n}\rangle}&=\sqrt{\frac{n_i}{N}}\vert{n_1,n_2,\ldots,n_{i-1},n_i-1,n_{i+1},\ldots,n_d}\rangle\\&=\frac{\hat{a}_i}{\sqrt{N}}\vert{\mathbf{n}\rangle}.\end{aligned}$$
By linearity, any pure state gets projected to $(\langle i|\otimes \mathbb{1}^{\otimes N-1})\vert{\psi }\rangle=\frac {\hat {a}_i}{\sqrt {N}}\vert{\psi }\rangle$. It immediately follows that tracing out the first photon from a state leads to a convex combination of such operations and, since the first photon is equivalent to the rest, we find the generic relation
$$\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)=\frac{1}{N}\sum_{i=1}^d \hat{a}_i{\hat{\rho}}_N \hat{a}_i^\dagger.$$
Removing another photon repeats the procedure with $N\to N-1$, so we employ the notation $\mathop {\mathrm {Tr}} \nolimits _{k}({\hat {\rho }})$ to denote the removal of $k$ photons from a state ${\hat {\rho }}$ (i.e., a state $\mathop {\mathrm {Tr}} \nolimits _{N-q}({\hat {\rho }}_N)$ has $q$ photons remaining and is the result of inspecting only $q$ photons from the beam).

Going beyond a fixed $N$, the above relationships can be supplanted by $\langle i\vert{\mathbf {n}}\rangle=\hat {a}_i\frac {1}{\sqrt {\hat {N}}}\vert{\mathbf {n}\rangle}$, where $\hat {N}=\sum _{i=1}^d\hat {a}_i^\dagger \hat {a}_i$ is the total photon-number operator. For an arbitrary state, it then follows that

$$\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}})=\sum_{i=1}^d \hat{a}_i\frac{1}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{1}{\sqrt{\hat{N}}} \hat{a}_i^\dagger{=}\sum_{i=1}^d \frac{1}{\sqrt{\hat{N}+1}}\hat{a}_i{\hat{\rho}} \hat{a}_i^\dagger \frac{1}{\sqrt{\hat{N}+1}},$$
which offers a unique way of removing photons even from superposition states that carry terms of the form $\vert{\mathbf {m}}\rangle \langle{\mathbf {n}}\vert$ with $|\mathbf {m}|\neq |\mathbf {n}|$. This quantity is unchanged via mode transformations and repartitionings of Hilbert space $\hat {a}_i\to \sum _{j=1}^d U_{ij}\hat {a}_j$ for unitary matrices $U$, as expected for a mode-agnostic description of removing photons.

Equation 6 properly describes how to treat a beam of light that has exactly one photon removed; beyond providing a mathematical language for discussing such situations, do the latter arise physically? One can argue that they occur pervasively in nature: when imaging a beam of light with a camera or one’s eye, not all of the photons are always recorded, even if no photons are lost on the way and the detector is perfectly efficient. This occurs, for example, in coincidence detection, when one records events in which two photons from a generic beam of light together impinge on any detector. In that case, one is inspecting a beam that effectively had photons removed according to our formalism. One could consider simulating Eq. 6 using mode-agnostic “photon subtraction” [58] to enact each operation $\rho \to \hat {a}_i\rho \hat {a}_i^\dagger$ and then normalizing the resultant state using a complicated photon-number-dependent operation that acts on all states as $|\psi \rangle \to \frac {1}{\sqrt {\hat {N}+1}}|\psi \rangle$, but this is physically a bit forced and is better considered using the simple photon-subset physical scenario just described.

One can also imagine applying the same formula Eq. 6 to fermionic or anyonic states, but one loses some of the richness upon noticing that a single fermionic mode can only be annihilated once before vanishing.

3. How Correlations Change With Photon Removal

Given an observable quantity, we expect a reasonable relationship between the quantity measured for the entire state and the quantity measured for the state with one or more photons removed. We start with a simple first-order correlation within one mode or between two modes, $\langle \hat {a}_k^\dagger \hat {a}_l\rangle$, and ask how the expectation value changes when taken with respect to ${\hat {\rho }}_N$ versus $\mathop {\mathrm {Tr}} \nolimits _{1}({\hat {\rho }}_N)$. By linearity, these expressions will tell us how any generator of the group SU$(d)$ changes when one $d$-level particle is removed from an $N$-particle symmetric state.

It is easy to compute, using the bosonic commutation relations,

$$\begin{aligned} \langle\hat{a}_k^\dagger\hat{a}_{l}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)} & =\mathop{\mathrm{Tr}} \nolimits(\hat{a}_k^\dagger\hat{a}_{l}\frac{1}{N}\sum_{i=1}^d \hat{a}_i{\hat{\rho}}_N \hat{a}_i^\dagger) \\ & =\frac{1}{N}\sum_{i=1}^d\mathop{\mathrm{Tr}} \nolimits[\hat{a}_k^\dagger(\hat{a}_{l}\hat{a}_i^\dagger{-}\delta_{il}) \hat{a}_i{\hat{\rho}}_N ]\\ & =\frac{1}{N}\mathop{\mathrm{Tr}} \nolimits[\hat{a}_k^\dagger\hat{a}_{l}(\hat{N}-1) {\hat{\rho}}_N ] =\frac{N-1}{N}\langle\hat{a}_k^\dagger\hat{a}_{l}\rangle_{{\hat{\rho}}_N}. \end{aligned}$$
Since this can be repeated for removing more photons, this means that the ratios between all of the first-order correlations for a state ${\hat {\rho }}_N$ are unchanged by removing any number of photons. As well, since the first-order correlations set the value of $\langle \hat {N}\rangle$, the normalization of all of the first-order correlations can be known directly from the total number of photons measured. For the example of the Stokes vector, we see that the classical polarization properties of a beam with a fixed number of photons are unchanged when only a subset of the photons are measured.

There are two immediate directions to generalize this calculation: to more general operators than first-order correlations and to more general states without fixed photon numbers.

A generic operator can always be expressed in normally ordered form [63,64] with components

$$\hat{O}_{\mathbf{k}\mathbf{l}}=\hat{a}_1^{{\dagger} k_1}\hat{a}_2^{{\dagger} k_2}\cdots \hat{a}_d^{{\dagger} k_d}\hat{a}_1^{l_1}\hat{a}_2^{l_2}\cdots \hat{a}_d^{l_d}.$$
Note that each $\hat {a}_i$ acts on mode $i$ such that many of these operators in the product commute with each other. Then, the same calculation as Eq. 7 but using the more general relationship $[\hat {a}_i^{l_i},\hat {a}_i^\dagger ]=l_i \hat {a}_i^{l_i-1}$ dictates that
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)}=\frac{N-|\mathbf{l}|}{N}\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}.$$
For this expression to be nonzero, we also require $|\mathbf {k}|=|\mathbf {l}|$. This reinforces that one requires at least $|\mathbf {l}|$ photons to observe a correlation that has $|\mathbf {l}|$ annihilation or creation operators; after removing one photon at a time until $N=|\mathbf {l}|$, removing one more photon causes these correlations to all vanish from the resulting state.

Convex combinations of the above results suffice to explain how the expectation values of all operators change upon removal of photons from states of the form $\sum _N p_N {\hat {\rho }}_N$ and of operators with $|\mathbf {k}|=|\mathbf {l}|$ for all states. What about expectation values of operators with $|\mathbf {k}|\neq |\mathbf {l}|$ for states with coherences between different photon-number sectors? We simply refrain from the substitution $\hat {N}\to N$ in the derivation of Eq. 7 to find many equivalent expressions:

$$\begin{aligned}\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}})}&= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\frac{\hat{N}-|\mathbf{l}|}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{1}{\sqrt{\hat{N}}}\right)= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\frac{1}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{\hat{N}-|\mathbf{k}|}{\sqrt{\hat{N}}}\right)\\&= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\sqrt{\frac{\hat{N}-|\mathbf{l}|}{\hat{N}}}{\hat{\rho}} \sqrt{\frac{\hat{N}-|\mathbf{k}|}{\hat{N}}}\right).\end{aligned}$$
One can choose which one of the formulas to apply, such as $\frac {\hat {N}-|\mathbf {l}|}{\sqrt {\hat {N}}}{\hat {\rho }} \frac {1}{\sqrt {\hat {N}}}$ to the state or $\frac {1}{\sqrt {\hat {N}}}\hat {O}_{\mathbf {k}\mathbf {l}}\frac {\hat {N}-|\mathbf {l}|}{\sqrt {\hat {N}}}$ to the operator, to see how the correlations change upon removal of a photon. To apply this formula, we take every term in ${\hat {\rho }}$ of the form $\vert{\mathbf {m}}\rangle \langle{\mathbf {n}}\vert$ with $|\mathbf {m}|-|\mathbf {l}|=|\mathbf {n}|-|\mathbf {k}|$ that contributed to $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle$ and scale its contribution by $\frac {|\mathbf {m}|-|\mathbf {l}|}{\sqrt {|\mathbf {m}||\mathbf {n}|}}$. Note that the correlations do not simply get multiplied by the same factor, unlike the case for ${\hat {\rho }}_N$, such that the ratios of the correlations change upon removal of photons from a state and thereby the correlations change significantly when one only observes a subset of photons from a state.

The correlations with $|\mathbf {k}|\neq |\mathbf {l}|$ slightly complicated all of the expressions. This seems to imply that measuring a quadrature operator such as $\hat {a}_i+\hat {a}_i^\dagger$ requires more careful calculation than an intensity correlation like $\hat {a}_i^\dagger \hat {a}_j+\hat {a}_j^\dagger \hat {a}_i$. However, it behooves one to recall that homodyne measurements that purport to measure quadrature operators on mode $i$ are actually measuring an intensity correlation between that mode and a field mode $j$, where the latter is taken to be in a coherent state with large amplitude. All routine measurements thus tend to fundamentally be made from operators that annihilate as many photons as they create, such that the simple expressions $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle _{\mathop {\mathrm {Tr}} \nolimits _{1}({\hat {\rho }})}=\sum _N p_N\frac {N-|\mathbf {l}|}{N} \langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle _{{\hat {\rho }}_N}$ can typically be used.

A different sort of correlations was inspected in Refs. [60,61], where states with fixed $N$ were inspected for their entanglement properties when one or more photons are removed or ignored using the same mechanism as the present one. It is likely that various combinations of the correlations $\langle \hat {O}_{\mathbf {kl}}\rangle$ can be used to witness entanglement, but different combinations will be necessary as witnesses for different scenarios and so we leave that investigation for future study.

4. Observing $q$ Photons at a Time

We argue here that there is only one possible integer number of photons $q$ that should be inspected at a time in order to measure any particular correlation $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle$ (see Fig. 1):

$$q=|\mathbf{l}|=|\mathbf{k}|.$$
For operators with $|\mathbf {k}|\neq |\mathbf {l}|$, we simply assume they are homodyne-type operators whose reference modes should be reintroduced until $|\mathbf {k}|= |\mathbf {l}|$. We give an argument for what to do in cases where this is not possible at the end of this section (again, see Fig. 1).

 figure: Fig. 1.

Fig. 1. Schematizing which photons contribute to measuring which operators. A generic superposition state may be comprised of multiple photons (different kets) in multiple modes (different star icons). (a) Inspecting two photons at a time gives rise to all second-order correlations $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle$ with $|\mathbf {k}|=|\mathbf {l}|=2$, with all the relevant pairs of photons circled by the dashed ovals; the pairs must all come from the same ket. (b) Inspecting another set of correlations for the same state, such as those with $|\mathbf {k}|=2$ and $|\mathbf {l}|=1$, requires comparing one photon from one ket and two photons from the ket with one more total photon; we include three representative combinations in blue circles, green squares, and orange pentagons. Generalizations to mixed states and superpositions with different amplitudes follow the intuition from this picture and are discussed in the text.

Download Full Size | PDF

A generic state ${\hat {\rho }}$ can be rewritten in terms of components with fixed and non-fixed photon numbers

$${\hat{\rho}}=\sum_N p_N {\hat{\rho}}_N +\sum_{|\mathbf{m}|\neq |\mathbf{n}|}{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\vert{\mathbf{m}}\rangle \langle{\mathbf{n}}\vert,$$
where each ${\hat {\rho }}_N$ is a normalized density operator with $N$ photons and none of the terms in the latter sum contribute to expectation values of operators with $|\mathbf {k}|=|\mathbf {l}|$. Observing $q$ photons from a state ${\hat {\rho }}_N$ requires computing
$${\hat{\varrho}}(q|{\hat{\rho}}_N)\equiv\mathop{\mathrm{Tr}} \nolimits_{N-q}({\hat{\rho}}_N)$$
by repeatedly using the simple formula Eq. 6 from above. For the more general state ${\hat {\rho }}$, one must weigh each of these contributions by the probability that the photons were taken from that sector,
$$P(q\,\mathrm{from}\,N)\equiv\frac{1}{\mathcal{N}(q)} p_N\binom{N}{q},$$
where $\mathcal {N}(q)=\sum _N p_N\binom {N}{q}$ is a normalization constant equal to the total number of expected $q$-photon events.

It would be nice to express the expectation value of an observable in terms of the $q$ photons that are actually observed. The former should take the form of the value of the observable for $q$ photons that came from a given $N$-photon subspace $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle _{{\hat {\varrho }}(q|{\hat {\rho }}_N)}$, i.e., the expectation value with respect to the state $\mathop {\mathrm {Tr}} \nolimits _{N-q}({\hat {\rho }}_N)$, multiplied by the probability that the photons came from that subspace, all scaled by the total number of times that $q$ photons are actually observed

$$\begin{aligned} \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}} \underset{?}= & \mathcal{N}(q)\sum_N P(q\,\mathrm{from}\,N) \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(q|{\hat{\rho}}_N)}\\ = & \sum_N p_N \binom{N}{q} \frac{q!}{N!}\frac{(N-|\mathbf{l}|)!}{(q-|\mathbf{l}|)!}\langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}\\ = & \sum_N p_N \binom{N-|\mathbf{l}|}{q-|\mathbf{l}|} \langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}. \end{aligned}$$
This equality is correct for arbitrary states if and only if $q=|\mathbf {l}|$. We thus learn that the expectation value of a correlation operator that annihilates $q$ photons from a variety of modes and creates $q$ photons among those modes can uniquely be reinterpreted as the expectation value the operator takes for a randomly selected subset of $q$ photons from the state. We cement this relationship as
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\mathcal{N}(|\mathbf{l}|)\sum_N P(|\mathbf{l}|\,\mathrm{from}\,N) \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(|\mathbf{l}||{\hat{\rho}}_N)}.$$
Immediate consequences include that all intensity moments must be measured one photon at a time. These include the Stokes vector, for example, and are simply scaled by the average number of photons $\mathcal {N}(1)=\langle \hat {N}\rangle _{\hat {\rho }}$. For a correlation function like $g^{(2)}$, two photons must be measured at a time and, similarly, for any $g^{(n)}$, $n$ photons must be inspected simultaneously.

As a corollary, we can immediately conclude what the state of a density operator must be when $q$ photons are chosen at random from it. The coefficients of the density operator in this basis must be the values of the correlation functions given by the appropriate observables,

$$\begin{aligned}{\hat{\varrho}}(q|{\hat{\rho}})&\equiv \sum_N P(q\,\mathrm{from}\,N) {\hat{\varrho}}(q|{\hat{\rho}}_N) \\&\quad\propto \sum_{|\mathbf{m}|=|\mathbf{n}|=q}\frac{\langle \hat{O}_{\mathbf{m}\mathbf{n}}\rangle_{{\hat{\rho}}}}{\sqrt{m_1!\cdots m_d!n_1!\cdots n_d!}}\vert{\mathbf{n}}\rangle \langle{\mathbf{m}}\vert.\end{aligned}$$
This is emphatically different from the projection of a state ${\hat {\rho }}$ onto the $q$-photon subspace and allows us to rewrite Eq. 16 as
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\mathcal{N}(|\mathbf{l}|)\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(|\mathbf{l}||{\hat{\rho}})}.$$
The factorial factors are necessary to cancel the values $\langle \mathbf {m}|\hat {O}_{\mathbf {m}\mathbf {n}}|\mathbf {n}\rangle$, while the proportionality constant is the inverse of $\mathcal {N}(q)=\sum _{|\mathbf {m}|=q}\frac {\langle \hat {O}_{\mathbf {m}\mathbf {m}}\rangle _{{\hat {\rho }}}}{m_1!\cdots m_d!}$. We give a more formal proof of this relationship in Supplement 1, in case it is more appealing than our intuitive arguments for some readers.

From this corollary we learn the intuitive equivalence: the expansion coefficients of a randomly chosen $q$ photons from a beam in the Fock space of the $d$ modes tell you what the state’s $q$th-order correlations are and vice versa. Inspecting another number of photons at a time cannot tell you about other correlations and inspecting other correlations cannot tell you how a certain $q$ photons will behave. More forcibly, this means that, for example, should two photons be detected simultaneously among all of the possible modes, their statistics would say nothing about the intensity correlations in the beam from whence they originated.

All of the above should suffice for operators with $|\mathbf {k}|\neq |\mathbf {l}|$ if appropriate reference modes are considered. We finally seek a construction that may work in the absence of such reference modes. There is no unambiguous meaning to tasks such as “removing one photon from $\vert{2,2}\rangle \langle{1,1}\vert$” but we can apply our $\mathop {\mathrm {Tr}} \nolimits _1$ formula and see what ensues. Performing the same calculations as above, we find the unique equivalence

$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\sum_{\mathbf{m}\mathbf{n}}\sqrt{{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\binom{|\mathbf{m}|}{|\mathbf{l}|}}\sqrt{{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\binom{|\mathbf{n}|}{|\mathbf{k}|}} \langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{|\mathbf{m}|-|\mathbf{l}|}(\vert{\mathbf{m}}\rangle \langle{\mathbf{n}\vert})}.$$
The expectation value for the overall state is the same as the expectation value in the reduced state that has $|\mathbf {k}|$ photons remaining in the bra and $|\mathbf {l}|$ photons remaining in the ket, weighed by a sort-of probability amplitude that the photons came from the $|\mathbf {m}|$- and $|\mathbf {n}|$-photon subspaces. No other terms can contribute to this expectation value (notice that $|\mathbf {n}|-|\mathbf {k}|$ must equal $|\mathbf {m}|-|\mathbf {l}|$) and this is the unique weight distribution that allows for a reinterpretation of expectation values for all states. We can also use this to construct some sector of a density operator that has “$q_1$ photons in its bras and $q_2$ photons in its kets” whose coefficients would be the corresponding correlations $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle _{\hat {\rho }}$, should we be willing to slightly abuse notation.

5. Agreement With Intuition From Photodetection

The origins of the correlation functions $\langle \hat {a}_i^{\dagger l_i}\hat {a}_i^{l_i}\rangle$ stem from a photodetection model where $l_i$ photons are absorbed by the detector and thereby annihilated from mode $i$. All of the remaining correlation functions to the same order can be obtained via a similar argument among a variety of modes $\hat {a}_i\to \sum _{j=1}^d U_{ij}\hat {a}_j$, so the most basic argument of standard photodetection is that $|\mathbf {l}|$ photons must be absorbed by a detector to measure a correlation of order $|\mathbf {l}|$. In this work, we have shown a separate, equivalent, mathematically rigorous method for motivating the same physical picture: the state of a beam of light that would be obtained by randomly selecting $|\mathbf {l}|$ photons to be sent to a detector is exactly the state that exclusively contains the information about the $|\mathbf {l}|$th-order correlations. Photodetection removes $|\mathbf {l}|$ photons to inspect them; our work shows that the $|\mathbf {l}|$ photons exactly carry the appropriate correlation information. As heralded in the abstract, the astute quantum optician will be pleased to find no surprises here, only justification for intuition.

6. Incorporating and Comparing With Standard Optical Loss

How does one deal with the situation in which $q$ photons are observed simultaneously, but some number of photons $r$ were lost prior to the detection? Should one use the state ${\hat {\varrho }}(q+r|{\hat {\rho }})$ or ${\hat {\varrho }}(q|{\hat {\rho }})$ for predicting what values the correlations will take?

There is actually no ambiguity here: one simply uses the state ${\hat {\rho }}$ that arises from the original state having lost $r$ photons and inputs it into the formula for ${\hat {\varrho }}(q|{\hat {\rho }})$. If one uses the state before it lost the $r$ photons, one must be able to have access to the loss modes. Then one would simply increase $d$ and detect $q$ photons among all of the modes; given that loss modes tend to be lost, this method will seldom work and one should instead treat the state as having lost $r$ photons prior to computing ${\hat {\varrho }}(q|{\hat {\rho }})$.

The other question is how this work relates to standard loss channels, which enact the input–output relations

$$\hat{a}_i\underset{\eta_i}{\to}\sqrt{\eta_i}\hat{a}_i+\sqrt{1-\eta_i^2}\hat{b}_i$$
and then trace over the “vacuum modes” annihilated by $\hat {b}_i$ that begin in their vacuum states. We show in Supplement 2 that having equal loss $\eta _i=\eta$ in all $d$ modes is equivalent to a state undergoing the transformation
$$\begin{aligned} \hat{\rho}_N\underset{\eta}{\to} & \sum_k\binom{N}{k}\left(1-\eta\right)^{k} \eta^{N-k} \mathop{\mathrm{Tr}} \nolimits_{N-k}( \hat{\rho}_N);\\ \hat{\rho}\underset{\eta}{\to} & \sum_k\frac{\left(\frac{1-\eta}{\eta}\right)^{k}}{k!} \underbrace{\mathop{\mathrm{Tr}} \nolimits_1(\sqrt{\hat{N}}\cdots \mathop{\mathrm{Tr}} \nolimits_1(\sqrt{\hat{N}}}_{k\,\mathrm{times}} \sqrt{\eta}^{\hat{N}}\hat{\rho} \sqrt{\eta}^{\hat{N}} \underbrace{\sqrt{\hat{N}})\cdots\sqrt{\hat{N}})}_{k\,\mathrm{times}}. \end{aligned}$$
The first transformation of Eq. 21 matches the one found in Refs. [62,65] in the contexts of boson sampling and quantum metrology. This provides an alternate formula for loss channels and directly explains why equal loss on $d$ modes commutes with all linear optical networks that enact $\hat {a}_i\to \sum _j U_{ij}\hat {a}_j$ (because the $\mathop {\mathrm {Tr}} \nolimits _1$ and $\hat {N}$ operations commute with such networks).

7. Example Applications

7.1 Purity of Reduced State

One indicator of quantum polarization properties is how mixed a pure state ${\hat {\rho }}_N$ is after removing a certain number of photons [66]. With our formula from Eq. 6, we can express the purity of a reduced pure state as

$$\scalebox{0.93}{$\begin{aligned} \mathop{\mathrm{Tr}} \nolimits\{[\mathop{\mathrm{Tr}} \nolimits_{N-q}(\vert{\psi_N}\rangle \langle{\psi_N}\vert)]^2\}& =\frac{q!^2}{N!^2}\sum_{i_1\cdots i_{N-q},j_1\cdots j_{N-q}}|\langle \hat{a}^\dagger_{i_1}\cdots \hat{a}^\dagger_{i_{N-q}} \hat{a}_{j_1}\cdots\hat{a}_{j_{N-q}}\rangle|^2\\ & =\frac{q!^2}{N!^2}\sum_{\mathbf{k}\mathbf{l}}\binom{N-q}{\mathbf{k}}\binom{N-q}{\mathbf{l}}|\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle|^2. \end{aligned}$}$$
This is just a sum of the ($N-q$)th-order moments of the state, as the multinomial coefficients vanish unless $|\mathbf {k}|=|\mathbf {l}|=N-q$. For example, in the two-mode case of polarization, the purity of a state after tracing out a single photon is determined by the Stokes parameters because it exclusively depends on $\langle \hat {a}^\dagger _i\hat {a}_j\rangle$; the purity after tracing out two photons depends on the covariances of the Stokes parameters; etc. We can explicitly calculate for the Stokes parameters that the purity of the state after tracing out a single photon is proportional to $S_0^2+S_1^2+S_2^2+S_3^2$, such that it is maximal for spin-coherent states and minimal for pure states whose Stokes vector vanishes.

7.2 Photon-Number Projection

Our $q$-photon state ${\hat {\varrho }}(q|{\hat {\rho }})$ from Eq. 17 is not the same as the projection of a state ${\hat {\rho }}$ onto the $q$-photon subspace; randomly choosing $q$ photons is different from projecting onto $q$ photons. How do these two ideas compare?

Due to the relation $\vert{0}\rangle \langle{0}\vert=:\exp (-\hat {a}^\dagger \vphantom {a}\hat {a}):$ that uses the normal ordering operation $::$, we can express a projector onto a certain $m$-photon state as

$$\vert{m}\rangle \langle{m}\vert=\frac{1}{m!}\hat{a}^{{\dagger} m}\sum_{n=0}^\infty ({-}1)^n\frac{\hat{a}^{{\dagger} n}\hat{a}^n}{n!}\hat{a}^m=\sum_{n=m}^\infty\frac{({-}1)^{n-m}}{m!(n-m)!}\hat{O}_{nn}.$$
Evaluating the expectation value of the projector for a state will thus require all states ${\hat {\varrho }}(q|{\hat {\rho }})$ with $q\geq m$. Intuitively, this means that there is a contribution to the projection from all components of the state from which $m$ photons could have arisen: we start with a contribution from the $m$-photon random state, adjust for contributions that may have arisen from states that truly had $m+1$ photons out of which $m$ were randomly chosen by considering the $m+1$-photon random state, and so on.

This result holds for multimode states as well. To project onto the $m$-photon subspace of a multimode state, one needs to account for all of the components of the state from which $m$-photons may be randomly selected using our above prescription.

8. Conclusion

There is a unique unambiguous method for describing how a beam of light behaves when $k$ photons are removed from it or when $q$ photons are selected from it at random. We have provided a number of simple formulas for these computations, with the result that the $q$th-order correlation functions of a $d$-mode state are in one-to-one correspondence with the expansion coefficients in the Fock basis for $q$ photons selected at random from the state. This closes the loop between the mathematics and intuition of why singles counts, doubles counts, and more are required for observing particular correlations; a detector effectively sees a $q$-photon state chosen at random from the original state. All of our work has been derived in the language of quantum optics, but the results should find equal application in the study of any totally symmetric $d$-level system across quantum information theory. Any other construction of ignoring particles from symmetric states would be disingenuous.

Funding

Natural Sciences and Engineering Research Council of Canada Postdoctoral Fellowship.

Acknowledgments

The NRC headquarters is located on the traditional unceded territory of the Algonquin Anishinaabe and Mohawk peoples.

Disclosures

The authors declare no conflicts of interest.

Data availability

No data were generated or analyzed in the presented research.

Supplemental document

See Supplement 1 for supporting content.

References

1. R. J. Glauber, “Coherent and incoherent states of the radiation field,” Phys. Rev. 131(6), 2766–2788 (1963). [CrossRef]  

2. R. J. Glauber, “The quantum theory of optical coherence,” Phys. Rev. 130(6), 2529–2539 (1963). [CrossRef]  

3. R. J. Glauber, “Photon correlations,” Phys. Rev. Lett. 10(3), 84–86 (1963). [CrossRef]  

4. E. C. G. Sudarshan, “Equivalence of semiclassical and quantum mechanical descriptions of statistical light beams,” Phys. Rev. Lett. 10(7), 277–279 (1963). [CrossRef]  

5. C. L. Mehta and E. C. G. Sudarshan, “Relation between quantum and semiclassical description of optical coherence,” Phys. Rev. 138(1B), B274–B280 (1965). [CrossRef]  

6. U. M. Titulaer and R. J. Glauber, “Correlation functions for coherent fields,” Phys. Rev. 140(3B), B676–B682 (1965). [CrossRef]  

7. K. E. Cahill and R. J. Glauber, “Ordered expansions in Boson amplitude operators,” Phys. Rev. 177(5), 1857–1881 (1969). [CrossRef]  

8. G. S. Agarwal and E. Wolf, “Calculus for functions of noncommuting operators and general phase-space methods in quantum mechanics. I. Mapping theorems and ordering of functions of noncommuting operators,” Phys. Rev. D 2(10), 2161–2186 (1970). [CrossRef]  

9. H. J. Carmichael and D. F. Walls, “Proposal for the measurement of the resonant Stark effect by photon correlation techniques,” J. Phys. B: At. Mol. Phys. 9(4), L43–L46 (1976). [CrossRef]  

10. R. Loudon, “Photon bunching and antibunching,” Phys. Bull. 27(1), 21–23 (1976). [CrossRef]  

11. H. J. Kimble, M. Dagenais, and L. Mandel, “Photon antibunching in resonance fluorescence,” Phys. Rev. Lett. 39(11), 691–695 (1977). [CrossRef]  

12. G. Leuchs, Photon Statistics, Antibunching and Squeezed States (Springer, 1986), pp. 329–360.

13. P. Grangier, G. Roger, and A. Aspect, “Experimental evidence for a photon anticorrelation effect on a beam splitter: a new light on single-photon interferences,” Europhys. Lett. 1(4), 173–179 (1986). [CrossRef]  

14. G. Rempe, F. Schmidt-Kaler, and H. Walther, “Observation of sub-Poissonian photon statistics in a micromaser,” Phys. Rev. Lett. 64(23), 2783–2786 (1990). [CrossRef]  

15. L. Mandel, “Sub-Poissonian photon statistics in resonance fluorescence,” Opt. Lett. 4(7), 205–207 (1979). [CrossRef]  

16. L. Mandel, “Non-classical states of the electromagnetic field,” Phys. Scr. T12, 34–42 (1986). [CrossRef]  

17. G. S. Agarwal and K. Tara, “Nonclassical properties of states generated by the excitations on a coherent state,” Phys. Rev. A 43(1), 492–497 (1991). [CrossRef]  

18. D. Klyshko, “Observable signs of nonclassical light,” Phys. Lett. A 213(1-2), 7–15 (1996). [CrossRef]  

19. C. T. Lee, “Simple criterion for nonclassical two-mode states,” J. Opt. Soc. Am. B 15(3), 1187–1191 (1998). [CrossRef]  

20. T. Richter and W. Vogel, “Nonclassicality of quantum states: a hierarchy of observable conditions,” Phys. Rev. Lett. 89(28), 283601 (2002). [CrossRef]  

21. E. V. Shchukin and W. Vogel, “Nonclassical moments and their measurement,” Phys. Rev. A 72(4), 043808 (2005). [CrossRef]  

22. D. Achilles, C. Silberhorn, and I. A. Walmsley, “Direct, loss-tolerant characterization of nonclassical photon statistics,” Phys. Rev. Lett. 97(4), 043602 (2006). [CrossRef]  

23. A. Zavatta, V. Parigi, and M. Bellini, “Experimental nonclassicality of single-photon-added thermal light states,” Phys. Rev. A 75(5), 052106 (2007). [CrossRef]  

24. A. Miranowicz, M. Bartkowiak, X. Wang, Y.-x. Liu, and F. Nori, “Testing nonclassicality in multimode fields: a unified derivation of classical inequalities,” Phys. Rev. A 82(1), 013824 (2010). [CrossRef]  

25. S. Ryl, J. Sperling, E. Agudelo, M. Mraz, S. Köhnke, B. Hage, and W. Vogel, “Unified nonclassicality criteria,” Phys. Rev. A 92(1), 011801 (2015). [CrossRef]  

26. J. Sahota and N. Quesada, “Quantum correlations in optical metrology: Heisenberg-limited phase estimation without mode entanglement,” Phys. Rev. A 91(1), 013808 (2015). [CrossRef]  

27. S. Hartmann, F. Friedrich, A. Molitor, M. Reichert, W. Elsäßer, and R. Walser, “Tailored quantum statistics from broadband states of light,” New J. Phys. 17(4), 043039 (2015). [CrossRef]  

28. J. Peřina, V. Michálek, and O. c. v. Haderka, “Higher-order sub-Poissonian-like nonclassical fields: theoretical and experimental comparison,” Phys. Rev. A 96(3), 033852 (2017). [CrossRef]  

29. V. Sudhir, R. Schilling, S. A. Fedorov, H. Schütz, D. J. Wilson, and T. J. Kippenberg, “Quantum correlations of light from a room-temperature mechanical oscillator,” Phys. Rev. X 7(3), 031055 (2017). [CrossRef]  

30. N. Korolkova and G. Leuchs, “Quantum correlations in separable multi-mode states and in classically entangled light,” Rep. Prog. Phys. 82(5), 056001 (2019). [CrossRef]  

31. P. Malpani, N. Alam, K. Thapliyal, A. Pathak, V. Narayanan, and S. Banerjee, “Lower- and higher-order nonclassical properties of photon added and subtracted displaced Fock states,” Ann. Phys. 531(2), 1800318 (2019). [CrossRef]  

32. L. S. Madsen, F. Laudenbach, M. F. Askarani, F. Rortais, T. Vincent, J. F. F. Bulmer, F. M. Miatto, L. Neuhaus, L. G. Helt, M. J. Collins, A. E. Lita, T. Gerrits, S. W. Nam, V. D. Vaidya, M. Menotti, I. Dhand, Z. Vernon, N. Quesada, and J. Lavoie, “Quantum computational advantage with a programmable photonic processor,” Nature 606(7912), 75–81 (2022). [CrossRef]  

33. A. Z. Goldberg, “Chapter Three - Quantum polarimetry,” Progress in Optics (Elsevier, 2022), pp. 185–274. [CrossRef]  

34. D. M. Klyshko, “Polarization of light: Fourth-order effects and polarization-squeezed states,” J. Exp. Theor. Phys. 84(6), 1065–1079 (1997). [CrossRef]  

35. A. Luis, “Polarization in quantum optics,” Prog. Opt. 61, 283–331 (2016). [CrossRef]  

36. A. Z. Goldberg, P. de la Hoz, G. Björk, A. B. Klimov, M. Grassl, G. Leuchs, and L. L. Sánchez-Soto, “Quantum concepts in optical polarization,” Adv. Opt. Photonics 13(1), 1–73 (2021). [CrossRef]  

37. T. Young, “I. The Bakerian Lecture. Experiments and calculations relative to physical optics,” Philos. Trans. R. Soc. London 94, 1–16 (1804). [CrossRef]  

38. L. Arnaud and N. J. Cerf, “Exploring pure quantum states with maximally mixed reductions,” Phys. Rev. A 87(1), 012319 (2013). [CrossRef]  

39. D. Baguette, T. Bastin, and J. Martin, “Multiqubit symmetric states with maximally mixed one-qubit reductions,” Phys. Rev. A 90(3), 032314 (2014). [CrossRef]  

40. O. Giraud, D. Braun, D. Baguette, T. Bastin, and J. Martin, “Tensor representation of spin states,” Phys. Rev. Lett. 114(8), 080401 (2015). [CrossRef]  

41. D. Baguette and J. Martin, “Anticoherence measures for pure spin states,” Phys. Rev. A 96(3), 032304 (2017). [CrossRef]  

42. S. Aaronson and D. J. Brod, “BosonSampling with lost photons,” Phys. Rev. A 93(1), 012335 (2016). [CrossRef]  

43. M. Ueda, N. Imoto, and T. Ogawa, “Quantum theory for continuous photodetection processes,” Phys. Rev. A 41(7), 3891–3904 (1990). [CrossRef]  

44. M. Dakna, T. Anhut, T. Opatrný, L. Knöll, and D.-G. Welsch, “Generating Schrödinger-cat-like states by means of conditional measurements on a beam splitter,” Phys. Rev. A 55(4), 3184–3194 (1997). [CrossRef]  

45. T. Opatrný, G. Kurizki, and D.-G. Welsch, “Improvement on teleportation of continuous variables by photon subtraction via conditional measurement,” Phys. Rev. A 61(3), 032302 (2000). [CrossRef]  

46. S. S. Mizrahi and V. V. Dodonov, “Creating quanta with an ‘annihilation’ operator,” J. Phys. A: Math. Gen. 35(41), 8847–8857 (2002). [CrossRef]  

47. J. Wenger, R. Tualle-Brouri, and P. Grangier, “Non-Gaussian statistics from individual pulses of squeezed light,” Phys. Rev. Lett. 92(15), 153601 (2004). [CrossRef]  

48. A. Ourjoumtsev, R. Tualle-Brouri, J. Laurat, and P. Grangier, “Generating optical Schrödinger kittens for quantum information processing,” Science 312(5770), 83–86 (2006). [CrossRef]  

49. V. Parigi, A. Zavatta, M. Kim, and M. Bellini, “Probing quantum commutation rules by addition and subtraction of single photons to/from a light field,” Science 317(5846), 1890–1893 (2007). [CrossRef]  

50. A. Zavatta, V. Parigi, M. S. Kim, and M. Bellini, “Subtracting photons from arbitrary light fields: experimental test of coherent state invariance by single-photon annihilation,” New J. Phys. 10(12), 123006 (2008). [CrossRef]  

51. H. Takahashi, K. Wakui, S. Suzuki, M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, “Generation of large-amplitude coherent-state superposition via ancilla-assisted photon subtraction,” Phys. Rev. Lett. 101(23), 233605 (2008). [CrossRef]  

52. T. Gerrits, S. Glancy, T. S. Clement, B. Calkins, A. E. Lita, A. J. Miller, A. L. Migdall, S. W. Nam, R. P. Mirin, and E. Knill, “Generation of optical coherent-state superpositions by number-resolved photon subtraction from the squeezed vacuum,” Phys. Rev. A 82(3), 031802 (2010). [CrossRef]  

53. C. Navarrete-Benlloch, R. García-Patrón, J. H. Shapiro, and N. J. Cerf, “Enhancing quantum entanglement by photon addition and subtraction,” Phys. Rev. A 86(1), 012328 (2012). [CrossRef]  

54. T. J. Bartley, P. J. D. Crowley, A. Datta, J. Nunn, L. Zhang, and I. Walmsley, “Strategies for enhancing quantum entanglement by local photon subtraction,” Phys. Rev. A 87(2), 022313 (2013). [CrossRef]  

55. L. Fan and M. S. Zubairy, “Quantum illumination using non-Gaussian states generated by photon subtraction and photon addition,” Phys. Rev. A 98(1), 012319 (2018). [CrossRef]  

56. N. Stiesdal, H. Busche, K. Kleinbeck, J. Kumlin, M. G. Hansen, H. P. Büchler, and S. Hofferberth, “Controlled multi-photon subtraction with cascaded Rydberg superatoms as single-photon absorbers,” Nat. Commun. 12(1), 4328 (2021). [CrossRef]  

57. K. Takase, J.-I. Yoshikawa, W. Asavanant, M. Endo, and A. Furusawa, “Generation of optical Schrödinger cat states by generalized photon subtraction,” Phys. Rev. A 103(1), 013710 (2021). [CrossRef]  

58. M. F. Melalkia, L. Brunel, S. Tanzilli, J. Etesse, and V. D’Auria, “Theoretical framework for photon subtraction with non–mode-selective resources,” Phys. Rev. A 105(1), 013720 (2022). [CrossRef]  

59. C. M. Nunn, J. D. Franson, and T. B. Pittman, “Modifying quantum optical states by zero-photon subtraction,” Phys. Rev. A 105(3), 033702 (2022). [CrossRef]  

60. A. Neven, J. Martin, and T. Bastin, “Entanglement robustness against particle loss in multiqubit systems,” Phys. Rev. A 98(6), 062335 (2018). [CrossRef]  

61. G. M. Quinta, R. André, A. Burchardt, and K. Życzkowski, “Cut-resistant links and multipartite entanglement resistant to particle loss,” Phys. Rev. A 100(6), 062329 (2019). [CrossRef]  

62. M. Oszmaniec and D. J. Brod, “Classical simulation of photonic linear optics with lost particles,” New J. Phys. 20(9), 092002 (2018). [CrossRef]  

63. V. V. Mikhailov, “Ordering of some boson operator functions,” J. Phys. A: Math. Gen. 16(16), 3817–3827 (1983). [CrossRef]  

64. P. Blasiak, K. Penson, and A. Solomon, “The general boson normal ordering problem,” Phys. Lett. A 309(3-4), 198–205 (2003). [CrossRef]  

65. M. Oszmaniec, R. Augusiak, C. Gogolin, J. Kołodyński, A. Acín, and M. Lewenstein, “Random bosonic states for robust quantum metrology,” Phys. Rev. X 6(4), 041044 (2016). [CrossRef]  

66. M. Rudziński, A. Burchardt, and K. Zyczkowski, “Orthonormal bases of extreme spin coherence,” (2023).

Supplementary Material (1)

NameDescription
Supplement 1       Details of some derivations

Data availability

No data were generated or analyzed in the presented research.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (1)

Fig. 1.
Fig. 1. Schematizing which photons contribute to measuring which operators. A generic superposition state may be comprised of multiple photons (different kets) in multiple modes (different star icons). (a) Inspecting two photons at a time gives rise to all second-order correlations $\langle \hat {O}_{\mathbf {k}\mathbf {l}}\rangle$ with $|\mathbf {k}|=|\mathbf {l}|=2$, with all the relevant pairs of photons circled by the dashed ovals; the pairs must all come from the same ket. (b) Inspecting another set of correlations for the same state, such as those with $|\mathbf {k}|=2$ and $|\mathbf {l}|=1$, requires comparing one photon from one ket and two photons from the ket with one more total photon; we include three representative combinations in blue circles, green squares, and orange pentagons. Generalizations to mixed states and superpositions with different amplitudes follow the intuition from this picture and are discussed in the text.

Equations (23)

Equations on this page are rendered with MathJax. Learn more.

$${\hat{\rho}}_N=\sum_{|\mathbf{m}|=|\mathbf{n}|=N}{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\vert{\mathbf{m}}\rangle \langle{\mathbf{n}}\vert.$$
$$\vert{\mathbf{n}}\rangle=\frac{\hat{a}_1^{{\dagger} n_1}}{\sqrt{n_1!}}\frac{\hat{a}_2^{{\dagger} n_2}}{\sqrt{n_2!}}\cdots \frac{\hat{a}_d^{{\dagger} n_d}}{\sqrt{n_d!}} \vert{\mathrm{vac}\rangle}$$
$$\vert{\mathbf{n}}\rangle=\sqrt{\frac{1}{\binom{N}{\mathbf{n}}}}\sum_{\mathrm{permutations}}\vert{1}\rangle^{{\otimes} n_1} \otimes\vert{2}\rangle^{{\otimes} n_2}\otimes\cdots\otimes\vert{d}\rangle^{{\otimes} n_d}.$$
$$\begin{aligned}(\langle i|\otimes\mathbb{1}^{{\otimes} N-1})\vert{\mathbf{n}\rangle}&=\sqrt{\frac{n_i}{N}}\vert{n_1,n_2,\ldots,n_{i-1},n_i-1,n_{i+1},\ldots,n_d}\rangle\\&=\frac{\hat{a}_i}{\sqrt{N}}\vert{\mathbf{n}\rangle}.\end{aligned}$$
$$\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)=\frac{1}{N}\sum_{i=1}^d \hat{a}_i{\hat{\rho}}_N \hat{a}_i^\dagger.$$
$$\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}})=\sum_{i=1}^d \hat{a}_i\frac{1}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{1}{\sqrt{\hat{N}}} \hat{a}_i^\dagger{=}\sum_{i=1}^d \frac{1}{\sqrt{\hat{N}+1}}\hat{a}_i{\hat{\rho}} \hat{a}_i^\dagger \frac{1}{\sqrt{\hat{N}+1}},$$
$$\begin{aligned} \langle\hat{a}_k^\dagger\hat{a}_{l}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)} & =\mathop{\mathrm{Tr}} \nolimits(\hat{a}_k^\dagger\hat{a}_{l}\frac{1}{N}\sum_{i=1}^d \hat{a}_i{\hat{\rho}}_N \hat{a}_i^\dagger) \\ & =\frac{1}{N}\sum_{i=1}^d\mathop{\mathrm{Tr}} \nolimits[\hat{a}_k^\dagger(\hat{a}_{l}\hat{a}_i^\dagger{-}\delta_{il}) \hat{a}_i{\hat{\rho}}_N ]\\ & =\frac{1}{N}\mathop{\mathrm{Tr}} \nolimits[\hat{a}_k^\dagger\hat{a}_{l}(\hat{N}-1) {\hat{\rho}}_N ] =\frac{N-1}{N}\langle\hat{a}_k^\dagger\hat{a}_{l}\rangle_{{\hat{\rho}}_N}. \end{aligned}$$
$$\hat{O}_{\mathbf{k}\mathbf{l}}=\hat{a}_1^{{\dagger} k_1}\hat{a}_2^{{\dagger} k_2}\cdots \hat{a}_d^{{\dagger} k_d}\hat{a}_1^{l_1}\hat{a}_2^{l_2}\cdots \hat{a}_d^{l_d}.$$
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}}_N)}=\frac{N-|\mathbf{l}|}{N}\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}.$$
$$\begin{aligned}\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{1}({\hat{\rho}})}&= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\frac{\hat{N}-|\mathbf{l}|}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{1}{\sqrt{\hat{N}}}\right)= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\frac{1}{\sqrt{\hat{N}}}{\hat{\rho}} \frac{\hat{N}-|\mathbf{k}|}{\sqrt{\hat{N}}}\right)\\&= \mathop{\mathrm{Tr}} \nolimits\left(\hat{O}_{\mathbf{k}\mathbf{l}}\sqrt{\frac{\hat{N}-|\mathbf{l}|}{\hat{N}}}{\hat{\rho}} \sqrt{\frac{\hat{N}-|\mathbf{k}|}{\hat{N}}}\right).\end{aligned}$$
$$q=|\mathbf{l}|=|\mathbf{k}|.$$
$${\hat{\rho}}=\sum_N p_N {\hat{\rho}}_N +\sum_{|\mathbf{m}|\neq |\mathbf{n}|}{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\vert{\mathbf{m}}\rangle \langle{\mathbf{n}}\vert,$$
$${\hat{\varrho}}(q|{\hat{\rho}}_N)\equiv\mathop{\mathrm{Tr}} \nolimits_{N-q}({\hat{\rho}}_N)$$
$$P(q\,\mathrm{from}\,N)\equiv\frac{1}{\mathcal{N}(q)} p_N\binom{N}{q},$$
$$\begin{aligned} \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}} \underset{?}= & \mathcal{N}(q)\sum_N P(q\,\mathrm{from}\,N) \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(q|{\hat{\rho}}_N)}\\ = & \sum_N p_N \binom{N}{q} \frac{q!}{N!}\frac{(N-|\mathbf{l}|)!}{(q-|\mathbf{l}|)!}\langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}\\ = & \sum_N p_N \binom{N-|\mathbf{l}|}{q-|\mathbf{l}|} \langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\rho}}_N}. \end{aligned}$$
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\mathcal{N}(|\mathbf{l}|)\sum_N P(|\mathbf{l}|\,\mathrm{from}\,N) \langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(|\mathbf{l}||{\hat{\rho}}_N)}.$$
$$\begin{aligned}{\hat{\varrho}}(q|{\hat{\rho}})&\equiv \sum_N P(q\,\mathrm{from}\,N) {\hat{\varrho}}(q|{\hat{\rho}}_N) \\&\quad\propto \sum_{|\mathbf{m}|=|\mathbf{n}|=q}\frac{\langle \hat{O}_{\mathbf{m}\mathbf{n}}\rangle_{{\hat{\rho}}}}{\sqrt{m_1!\cdots m_d!n_1!\cdots n_d!}}\vert{\mathbf{n}}\rangle \langle{\mathbf{m}}\vert.\end{aligned}$$
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\mathcal{N}(|\mathbf{l}|)\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{{\hat{\varrho}}(|\mathbf{l}||{\hat{\rho}})}.$$
$$\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\hat{\rho}}=\sum_{\mathbf{m}\mathbf{n}}\sqrt{{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\binom{|\mathbf{m}|}{|\mathbf{l}|}}\sqrt{{\hat{\rho}}_{\mathbf{m}\mathbf{n}}\binom{|\mathbf{n}|}{|\mathbf{k}|}} \langle\hat{O}_{\mathbf{k}\mathbf{l}}\rangle_{\mathop{\mathrm{Tr}} \nolimits_{|\mathbf{m}|-|\mathbf{l}|}(\vert{\mathbf{m}}\rangle \langle{\mathbf{n}\vert})}.$$
$$\hat{a}_i\underset{\eta_i}{\to}\sqrt{\eta_i}\hat{a}_i+\sqrt{1-\eta_i^2}\hat{b}_i$$
$$\begin{aligned} \hat{\rho}_N\underset{\eta}{\to} & \sum_k\binom{N}{k}\left(1-\eta\right)^{k} \eta^{N-k} \mathop{\mathrm{Tr}} \nolimits_{N-k}( \hat{\rho}_N);\\ \hat{\rho}\underset{\eta}{\to} & \sum_k\frac{\left(\frac{1-\eta}{\eta}\right)^{k}}{k!} \underbrace{\mathop{\mathrm{Tr}} \nolimits_1(\sqrt{\hat{N}}\cdots \mathop{\mathrm{Tr}} \nolimits_1(\sqrt{\hat{N}}}_{k\,\mathrm{times}} \sqrt{\eta}^{\hat{N}}\hat{\rho} \sqrt{\eta}^{\hat{N}} \underbrace{\sqrt{\hat{N}})\cdots\sqrt{\hat{N}})}_{k\,\mathrm{times}}. \end{aligned}$$
$$\scalebox{0.93}{$\begin{aligned} \mathop{\mathrm{Tr}} \nolimits\{[\mathop{\mathrm{Tr}} \nolimits_{N-q}(\vert{\psi_N}\rangle \langle{\psi_N}\vert)]^2\}& =\frac{q!^2}{N!^2}\sum_{i_1\cdots i_{N-q},j_1\cdots j_{N-q}}|\langle \hat{a}^\dagger_{i_1}\cdots \hat{a}^\dagger_{i_{N-q}} \hat{a}_{j_1}\cdots\hat{a}_{j_{N-q}}\rangle|^2\\ & =\frac{q!^2}{N!^2}\sum_{\mathbf{k}\mathbf{l}}\binom{N-q}{\mathbf{k}}\binom{N-q}{\mathbf{l}}|\langle \hat{O}_{\mathbf{k}\mathbf{l}}\rangle|^2. \end{aligned}$}$$
$$\vert{m}\rangle \langle{m}\vert=\frac{1}{m!}\hat{a}^{{\dagger} m}\sum_{n=0}^\infty ({-}1)^n\frac{\hat{a}^{{\dagger} n}\hat{a}^n}{n!}\hat{a}^m=\sum_{n=m}^\infty\frac{({-}1)^{n-m}}{m!(n-m)!}\hat{O}_{nn}.$$
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.