{'fraction_non_alphanumeric': 0.03723083506778969, 'fraction_numerical': 0.009648116088622614, 'mean_word_length': 4.501926369863014, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 2, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 21, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Linear optics and photodetection achieve near-optimal unambiguous coherent state discrimination

Coherent states of the quantum electromagnetic field, the quantum description of ideal laser light, are prime candidates as information carriers for optical communications. A large body of literature exists on their quantum-limited estimation and discrimination. However, very little is known about the practical realizations of receivers for unambiguous state discrimination (USD) of coherent states. Here we fill this gap and outline a theory of USD with receivers that are allowed to employ: passive multimode linear optics, phase-space displacements, auxiliary vacuum modes, and on-off photon detection. Our results indicate that, in some regimes, these currently-available optical components are typically sufficient to achieve near-optimal unambiguous discrimination of multiple, multimode coherent states.


I. INTRODUCTION
Quantum mechanics places fundamental limits on the distinguishability of non-orthogonal quantum states. This fact underpins applications to quantum information science and technology, most notably in quantum cryptography [1], and constraints the performance of quantum sensors [2], quantum communications [3,4], and of probabilistic algorithms in computation [5]. However, quantum mechanics also provides the tools to identify these limits and approach them through the design of practical receivers [6]. Within the framework of optical quantum information science and technology, coherent states of the electromagnetic field play a prominent role as information carriers in terrestrial and space-based quantum networks [4,7]. This is essentially due to the fact that these states are relatively easy to prepare and control experimentally, and to their robustness to loss. They are in fact the only pure states that have a classical limit, yet they exhibit fundamental quantum properties when sufficiently attenuated [8,9]. This paper focuses on the quantum-limited detection of coherent states. We present a few explicit and practical schemes for unambiguous state discrimination (USD) of coherent states. In particular, instead of considering global bounds, our goal is to design practical receivers that can be realized with linear optics and on-off photodetection. For brevity, we refer to them as linear receivers.
A well-designed quantum receiver combines practicality with high performance, where the latter is quantified through a suitable task-dependent figure of merit. * Electronic address: jsmdrsidhu@gmail.com † Electronic address: cosmo.lupo@poliba.it In ambiguous state discrimination (ASD), the receiver is designed to always provide an output, with the goal of minimising the average error probability. The latter is ultimately limited by the Yuen-Kennedy-Lax (YKL) conditions [20] (which reduces to the Helstrom bound for two states [21]). Alternative figures of merit are preferred for quantum key distribution [22][23][24][25] and quantum digital signatures [26], where an exact state identification improves the transfer of secure information that cannot be forged or repudiated. This framework is known as USD, where the receiver either identifies the state without error, or outputs an inconclusive result, and the goal is to minimize the average probability of such inconclusive result [10,11]. The ultimate bounds to USD can be computed using the theory of Peres and Terno [12] or that of Sun et al. [18] and Bergou et al. [19], which in a sense provide the analogue for the YKL conditions for USD. These bounds can be computed efficiently through semidefinite programming [13], and sometimes allow for an exact analytical solution [18,19]. However, experimentally attaining the bounds, in particular for the discrimination of coherent states, may require high non-linearities that may be un-accessible even with state-of-the-art or nearterm future optical technologies. There is a wide body of literature devoted to establishing the global bound for USD for different families of quantum states [14][15][16][17][18][19]. However, despite these advances, very little is known about practical USD receivers for coherent states, and whether they can achieve the global bounds [24,[27][28][29][30], with most of the available works focusing on phase-shift keying. For classical communications using coherent-state modulation, a joint quantum measurement acting on the received coherentstate code word that performs USD achieves the optimal communication capacity allowed by quantum mechanics [31], known as the Holevo capacity [32]. Furthermore, when acting on a finite-length inner-code comprised of tensor product of coherent states, the USD measurement can even attain a higher channel capacity-Shannon ca-pacity of the super channel induced by the inner code and the receiver-compared to the optimal ASD measurement that minimizes the average probability of error of choosing between the modulated-received inner code words [33][34][35]. For phase-shift encoded signals, the optimal USD measurement in the limit of small photon numbers consists of mode-wise displacement operations followed by photon-number-resolving detectors [24]. For the so-called binary phase-shift-keying (BPSK) alphabet, this scheme is sufficient to attain the optimal USD performance without adaptive pre-detection displacements [36]. Post-selecting on the measurement result can further reduce the error rate for a fixed probability of inconclusive results [37]. For single-mode coherent state constellations with more than two states, there is a substantial gap between the optimal USD performance and non-adaptive receivers [24,27]. An adaptive receiver for quadrature phase-shift-keying (QPSK) has demonstrated an improvement to correct state identification [28] in an effort to close this gap.
Here we outline a general theory of USD of coherent states with a receiver that leverages only readily available technologies, such as multi-mode linear passive optics, phase-space displacement operations, auxiliary vacuum modes, and mode-wise on-off photon detection. We apply our theory to a number of examples of different coherent-state modulations, using as a benchmark the global bounds on USD. The latter are computed explicitly using the theory of Refs. [18,19]. We demonstrate that in some regimes this practical, yet restricted, set of physical operations is typically sufficient to deliver nearoptimal performance. This work establishes a theoretical framework to understand and master the design of receivers to enable near optimal USD of coherent states.

A. Summary of results
Previous works in USD have primarily focused on developing efficient computational methods to determine the global bounds [12,13,19], providing necessary and sufficient conditions that define optimal measurement schemes [14,15,18], or considering specific examples of quantum states to construct measurement operators [16,17,29,38,39].
Here we focus on the design of practical receivers. We introduce a family of receivers based on multimode linear optics and on-off photodetection. A particular strength of these receivers is that they do not require complex adaptive strategies or non-linear optics, which may be challenging to implement. Our work readily addresses the optimal design to minimize the average probability of an inconclusive event. Further, for multimode coherent states, we provide intuitive insights into how the performance of linear receivers depends on the number of modes and on the number of photons detected. This intuition is useful to provide an understanding of how to experimentally realize improvements when using just lin-ear optics.
We test our receiver design on randomly generated multimode coherent states, showing that linear receivers provide near-optimal performance for a range of coherent states. Non-typical, highly-degenerate constellations of coherent states may necessitate coincidence measurements to achieve USD. We also address an alternative figure of merit, the communication capacity, showing that linear receivers achieve near-optimal performance both in the asymptotic regime and for finite block length. Table I provides a high-level summary of the different receiver designs that we introduce in this work, along with their optimal performance bounds, conditions required to saturate these bounds, and how the designs can be adapted to handle general codes. Our results significantly advance results in the literature by clarifying the receiver designs that are both near-optimal for many codes and immediately accessible with current technologies.

B. Outline
We begin in Section II by introducing our notation and the basic theoretical tools to describe our receivers based on linear optics and photodetection. We provide qualitative bounds that introduce key insights into the performance of different detection schemes in Section III. Then, in Section IV we establish a theory to formalize in a quantitative way the problem of USD with limited resources. In Section V, we discuss applications of linear receivers to different constellations of coherent states, including the pulse-position modulation codes, random codes, and non-typical degenerate codes with single and double degeneracy. These applications demonstrate the strengths of linear receivers when benchmarked against the global bound for USD. Finally, conclusions and open questions are provided in Section VI.

II. LINEAR OPTICS IN PHASE SPACE
Before delving into the details of the USD receivers, we first need to introduce our notation and a few basic elements from the toolbox of linear optics [40]. Consider a collection of m bosonic modes with annihilation and creation operators {a j , a † j }, for j = 1, . . . , m, satisfying the canonical commutation relations, [a j , a † j ] = δ jj . These modes may represent a number of physical degrees of freedom, e.g., polarization, transverse wave vector, time of arrival, orbital angular momentum, as long as they are all degenerate in frequency. A coherent state on mode j is denoted as |α j and is characterized by its complex amplitude α. We have  . Each class implements on-off photodetection and differs in either the detection strategy or resources used. Note that class 3 receivers can readily be generalized to code words with arbitrary degeneracy by implementing multimode detection events (greater than two). Note that the performance for each class of our linear optical receivers is defined through an optimisation problem. We developed a numerical optimizer to address each optimization and derive the optimal solutions. Intuition into the performance of each receiver is provided in section III. Note the construction of each receiver class uses linear optics that can be readily implemented, which is in contrast to the globally optimal scenario summarised in the final row. where is the Fock state with k photons on mode j, and |0 is the vacuum state. A multimode coherent state is the direct product of m coherent states: Such a state is uniquely identified by its amplitude vector α = (α 1 , α 2 , . . . , α m ). As |α j | 2 is the mean photon number in mode j, is the total mean photon number in the state. The fundamental mathematical structure in quantum mechanics is the scalar product defined in the Hilbert space. The scalar product between two multimode coherent states |α and |β reads where is the scalar product between the amplitude vectors. The Hilbert-space scalar product is invariant under the action of general unitary transformations in the Hilbert space and is the central mathematical structure underlying the global bounds on USD [12]. Note that this global bound is known to be achieved by some measurements. However, for generic states, we do not expect any particular form for such optimal measurement. Since our focus here is on achieving USD using a particular subset of measurements, which includes linear optics and photodetection, the Hilbert-space scalar product may not be the most useful mathematical tool. Therefore, we consider an alternative notion of scalar product that seems more naturally suited to describe quantum mechanics under a restricted set of allowed measurements. To identify this alternative scalar product, we need to consider in more detail the set of unitary linear optics operations, which preserve the total mean photon number. These operations identify the group of Linear Optical Passive (LOP) unitary transformations [41,42].
LOP unitaries map coherent states into coherent states. The components of the amplitude vector transform as follows, where [U jk ] is a m × m unitary matrix. Multimode LOP transformations can be implemented physically by combining linear optics elements as beam splitters and phase shifters, and are mathematically described as unitary matrices. Given a unitary matrix description of the LOP transformations, there are known, efficient procedures to simulate it as a network of beam splitters and phase shifters [43][44][45]. Note that the scalar product between amplitude vectors in Eq. (6) is invariant under the action of LOP unitaries. Therefore, given a pair of m-mode coherent states |α , |β , with amplitude vectors α = (α 1 , α 2 , . . . , α m ) and β = (β 1 , β 2 , . . . , β m ), we define the phase-space scalar product as The phase-space scalar product will play in our analysis a similar role played by the Hilbert-space scalar product in the theory of USD developed by Peres and Terno [12]. It will guide us in designing our linear receivers for USD of coherent states.
In the rest of the paper, we will consider examples of codes comprising c ≥ 2 m-mode coherent states, identified by c amplitude vectors α 1 , α 2 , . . . , α c . We can arrange these vectors into a rectangular c × m matrix It is well known that pure states can be unambiguously discriminated against if and only if they are linearly independent [14]. In our setting, where we only allow for limited resources, this condition needs to be modified. As a matter of fact, we need to look at the notion of linear independence in phase space, not in the Hilbert space. This alternative notion of linear independence, that we simply call phase-space linear independence, is strictly related to the phase-space scalar product introduced above. Therefore, we will say that the given coherent states are linear independent in phase space if their associated amplitude vectors are linear independent, i.e., if the matrix R in Eq. (9) has maximum rank. Note that this is possible only if c ≤ m. Below we describe a scheme for USD of coherent states that works when the associated R-matrix is full rank. In general, our receivers will also exploit a number of auxiliary modes. We will then introduce m additional optical modes (as we see below, it is sufficient to use m ≤ m auxiliary modes), characterized by the canonical operators {a m+j , a † m+j }, for j = 1 to m . The totality of m+m will be mixed at a LOP unitary over m+m modes, which is represented by a unitary matrix of size m+m . If all modes are initially prepared in a coherent state, with amplitude vector α = (α 1 , . . . , α m , α m+1 , . . . α m+m ), then the output is also a coherent state, with amplitude vector β = (β 1 , . . . , β m , β m+1 , . . . β m+m ), where the ith output amplitude reads where M and N are submatrices of U , with M ij = U ij for j = 1, . . . , m, and N ij = U ij for j = m + 1, . . . , m + m .
For our applications, we are often interested in the case where the m ancillary modes are prepared in the vacuum state, i.e., α m+j = 0. In this setting, we simply have When the output coherent state is measured by modewise photodetection, for any i we will have a certain probability of detecting a photon in the ith mode. This probability is readily computed as where and The photon detection probability in Eq. (12) is the key quantity to characterize our linear receivers, as we discuss in Section IV. Finally, we recall that the displacement operator, denoted as D(γ), maps a coherent state into a coherent state with a displaced amplitude: In the m-mode case, the displacement operator is identified by a complex displacement vector, γ = (γ 1 , γ 2 , . . . , γ m ), such that where α+γ = (α 1 +γ 1 , α 2 +γ 2 , . . . , α m +γ m ) and D(γ) = D(γ 1 ) ⊗ D(γ 2 ) ⊗ · · · D(γ m ).

III. QUALITATIVE ANALYSIS OF LINEAR RECEIVERS FOR USD OF COHERENT STATES
This work addresses the following general question: Given a collection of c, m-mode coherent states |α j , with amplitude vectors α j = (α j 1 , α j 2 , . . . , α j m ) and prior probabilities p j , is it possible to discriminate them unambiguously using linear optics, m auxiliary vacuum modes, and mode-wise photodetection? In this Section, we provide a qualitative answer to this question. These qualitative results will give some insight into the numerical analysis discussed later in the paper.
Consider a set of m-mode coherent states with n mean photon number. If we use m ancillary vacuum modes, then the mean photon number per mode is ν = n/(m + m ). These states are processed through a LOP unitary, which preserves the photon number, and then measured by mode-wise photodetection. On average, the probability of a single photodetection event in a given mode is expected to be about where the factor f depends on the details of the input states and LOP transformation applied. For weak signals, it is unlikely to observe more than one photodetection event. Therefore, we cannot discriminate more than m+m states, i.e., no more than the number of signal plus auxiliary modes, and the receiver yields an inconclusive result every time no photon is detected. This happens with probability Therefore, we expect the probability of an inconclusive event to decrease linearly with n for n 1. This probability can be further decreased if instead of vacuum auxiliary modes we use auxiliary coherent states. In fact, this would increase the mean photon number from n to n > n. This strategy is equivalent to introducing phasespace displacement in the pool of allowed resources [46]. In Section V C, see Fig. 4(a) in particular, we will show that this qualitative behavior is confirmed quantitatively, both with and without phase-space displacement.
For brighter coherent states, or when the first order term in n vanishes, we need to consider joint detection events on pairs of modes. In this case, the number of distinct outcomes is which gives the maximum number of states that can be unambiguously discriminated by observing coincidence detection on two modes. From Eq. (17), and recalling that a multimode coherent state is a product state, see Eq.
(3), we obtain the probability of a coincidence By ignoring joint detection events in more than two modes, the probability of obtaining an inconclusive result is Recalling that ν = n/(m + m ), we obtain In this case, the probability of the inconclusive event decreases quadratically with n and, for m+m large enough, is essentially independent of the number of modes and ancillas. As for the previous case, using coherent states ancillas instead of vacuum ones, may further decrease P 0 , as this is equivalent to replacing n with n > n. This qualitative result will be confirmed by the quantitative analysis of Section V F, see Fig. 7(b) in particular, where we will discuss examples of degenerate codes where the first order term in n vanishes. In principle, a similar qualitative analysis can be applied to USD schemes based on an arbitrary number of joint detection events. We naturally expect that events of joint detection of n d photons will contribute to P 0 with a term proportional to n n d . However, in this work, we only discuss receivers based on single and double detection events. In the following Section, we will focus on photodetection events on at most one mode, and provide quantitative analysis and explicit USD schemes. Later in Section V F we will consider an example of USD based on joint detection on two modes.

IV. EXPLICIT SCHEMES FOR USD OF COHERENT STATES
Below we present three schemes for linear receivers. The first two, in Sections IV A and IV B, exploit a single detection event to achieve USD of coherent states. The third scheme in Section IV C exploits double detection events.

A. Single detection events
We now introduce an explicit scheme for USD of multimode coherent states, based on single photo-detection events. Consider a set of c coherent states over m modes, identified by the amplitude vectors α 1 , α 2 , . . . , α c , where α j = (α j 1 , α j 2 , . . . , α j m ). In order to discriminate these states, we first introduce m auxiliary vacuum modes, then mix them on a LOP transformation identified by the (m + m ) × (m + m ) unitary matrix U , and finally apply mode-wise on-off photodetection. This setup is shown in Fig. 1. As recalled in Section II, the probability that, given the state |α j in input, a photon is detected on mode i in output, is where the vectors M i are determined by the elements of the matrix U as discussed in Section II.
If we want to achieve USD using the information obtained from a single photodetection event, then it is obvious that we need to impose that different input states correspond to photons in different output modes. Without loss of generality, this means that for input state j, only the output mode j has a non-zero probability of containing a photon, and all other output modes are in the Note that the LOP unitaries can be efficiently implemented through a network of beam splitters and phase shifters [43][44][45].
vacuum state. That is, the matrix U needs to be chosen in such a way that Expressed in terms of the vectors M i , this condition reads We note that a non-zero matrix M that solves these equations exists if and only if the amplitude vectors are linearly independent, i.e., the matrix R in Eq. (9) is full rank [54]. The conditions in Eq. (28) ensure that the detection of a photon in output mode j unambiguously identifies the coherent state |α j in input. Otherwise, no conclusion can be drawn in the event that no photon is detected in any of the output modes. This is the inconclusive event. If the input coherent state |α j has probability p j , then the average probability of the inconclusive event is Our goal is to design an explicit receiver that minimizes this probability. This corresponds to performing an optimization on the c × m matrix M , keeping in mind that the latter is by construction of the submatrix of a larger unitary matrix. This latter condition is expressed by the matrix inequality where I is the identity matrix (see Appendix C for a derivation of this condition). In solving the constrained optimization, we parameterize the elements of the matrix M as where k i ≥ 0 are c non-zero coefficients, and v ij are the components of c unit vectors for i = 1, . . . , c. The unit vectors v i are proportional to the vectors M i , therefore the condition (28) becomes In conclusion, using this parameterization, the optimal linear receiver is determined by solving the constrained optimization: The minimal value represents the minimum probability of obtaining an inconclusive event for distinguishing the given set of coherent states using only vacuum auxiliary modes, LOP unitaries, and on-off photodetection. This value can be compared to the global bound obtained when general operations and measurements are allowed, which can be computed analytically or numerically using results already available in literature [19].
Finally, we remark that the optimisation problem (34) is not explicitly dependent on the number of auxiliary modes m . However, once an optimal form for the matrix M is obtained, one needs to find a unitary matrix that extends it. In general, such a matrix exists only if m is chosen sufficiently large. However, one can assume m ≤ m without loss of generality. An explicit construction is given in Appendix C.

B. Phase-space displacement improves USD
A better USD scheme, i.e., one with a lower probability of an inconclusive event, can be obtained by enlarging the set of allowed resources. Here we describe a USD scheme obtained by adding the operation of phase-space displacement. We consider a setup where multimode displacement is first applied to the input modes, followed by the LOP unitary and mode-wise photodetection. This is shown in Fig. 2.
The optimal setup is then obtained by minimizing the probability of the inconclusive event, where now there are m additional complex degrees of freedom, corresponding to the components of the displacement vector γ: The examples of Section V will show that the introduction of phase-space displacement may improve the performance of USD substantially. As we noted above, the constraint can be satisfied if and only if the displaced amplitude vectors β j = α j + γ are linearly independent. The use of displacement operations has the added value of casting a linearly-dependent system with single degeneracy [55] into an independent one. This extends the range of coherent states that can be discriminated unambiguously. To see this, consider a linear dependent set of coherent states with c ≤ m, and rank-deficient R matrix with rank(R) = c − 1. The matrix can be made full rank by adding to each row a linear independent vector γ. The new matrix, has full rank and represents the set of displaced coherent states. Displacement also comes to the rescue when the number of states is larger than the number of modes, in the case c = m + 1 and rank(R) = c − 1. To make the matrix full rank, we first add one auxiliary mode, then displace the m + 1 modes by (γ, γ m+1 ). The new matrix reads and has rank c for suitable choices of (γ, γ m+1 ). An explicit example of this method is presented in Section V G for the discrimination of the binary phase-shift code |α , | − α .

C. Double detection events
When the amplitude vectors are not linearly independent, and the degeneracy of the associated R matrix is higher than 1, we need to consider multiple detection events. The simplest family of doubly degenerate codes is obtained for c = 3 and m = 1, i.e., a code of three coherent states over one mode: |α 1 , |α 2 , |α 3 . To discriminate these states, we introduce an explicit linear receiver design that makes use of two auxiliary vacuum modes. The three modes (one signal and two auxiliary modes) are first mixed in LOP unitary, then they are displaced by γ 1 , γ 2 , γ 3 , and finally each mode is detected as illustrated in Fig. 3. Unlike the receivers of Sections IV A, IV B, here two joint detection events unambiguously determine the input state. Without loss of generality, we require that input state |α 1 yields the vacuum in the first output mode, whereas the other two output modes both have a non-zero probability of photon detection. Similarly, we require that input state |α 2 yields the vacuum state on the second output mode, and state |α 3 yields the vacuum on the third output mode.
To write this condition explicitly, consider the unitary matrix U ki that represents the LOP transformation. Then the amplitude on the output mode k, given the input state |α j is We require that the output amplitude vanishes for j = k, i.e., from which we obtain three constraints: for j = 1, 2, 3. Given the input state |α 1 , the probability of obtaining two joint photodetection events on modes 2 and 3 is where in the last equation we have used condition (41). Note that P (2, 3|1) is the probability of identify the input state |α 1 , as we use double photodetection events to unambiguously discriminate the input states. From this expression, given the prior probability p j for coherent state |α j , we compute the average probability of obtaining an inconclusive event: The optimal performance of linear receivers based on double detection events is then obtained by minimizing this expression for P 0 under the constraint that the matrix U is unitary: This constrained optimization problem can be solved using Lagrange multipliers. Using the same setup, we can generalize this approach to codes with higher degeneracies.

V. EXAMPLES
In this Section, we explore with a few examples the performance of our linear optical USD receivers relative to the global bound, for different coherent states and figures of merits. We illustrate that, despite its simplicity, decoders constructed using only linear components and onoff photodetection can generate near-optimal USD: this is observed for randomly generated coherent states. Phasespace displacement is necessary to achieve near-optimal performances. The situation appears to be different for non-typical coherent states. In fact, our decoders work poorly when the R matrix is degenerate, especially when its rank is smaller than c − 1, where c is the number of states. For such highly degenerate cases, we illustrate a recipe for achieving improved performance based on coincidence events.

A. Discrimination of PPM codes and equivalent codes
A set of orthogonal vectors can always be discriminated against perfectly. Since coherent states are not orthogonal, they can never be perfectly discriminated (though their discrimination can be improved by increasing their distance in phase space). In this Section, we explore a particular family of coherent states that share some formal features with orthogonal states. A pulse-position modulation (PPM) code is a set of c = m coherent states over m modes: such that the matrix R = αI is a full-rank square matrix of dimension m. Note that these coherent states are mutually orthogonal with respect to the phase-space scalar product, with Gram matrix (α i , α j ) = |α| 2 δ ij .
In general, there is a gap between the global bound and the performance of linear receivers. However, for a PPM code, it is easy to show that the gap vanishes as the PPM code can be optimally discriminated against using on-off photodetection only. In fact, to unambiguously discriminate the states in the PPM code, it is sufficient to apply mode-wise photodetection. If a photon is detected on mode j, then we know the input state necessarily was |α j with no ambiguity. The inconclusive event is when no click is recorded, which happens with probability P 0 = e −|α| 2 . It is well known that this inconclusive probability saturates the global bound. To show this, one can apply the results of Bergou et al. in Ref. [19] (reviewed in Appendix B).
The optimality of linear receivers extends to a larger class of codes beyond PPM. First note that, since the phase-space scalar product is invariant under LOP unitaries, it immediately follows that LOP unitaries and onoff photodetection are sufficient to optimally discriminate any set of coherent states that has the same Gram matrix as the PPM code If the Gram matrix is not diagonal, we may try to make it diagonal by applying a phase-space displacement α j → β j = α j + γ, which changes the Gram matrix into Therefore, the code can be optimally discriminated with LOP unitaries, displacements, and photodetection, if there exist γ and τ > 0 such that This is a system of c 2 real equations and 2m + 1 real unknowns (the components of the complex displacement vector γ and τ ). Therefore, in general, we expect this system of equations to have solutions if c 2 ≤ 2m + 1.
As we have seen above, sometimes the use of an auxiliary mode can improve the effectiveness of linear receivers. In fact, adding an auxiliary mode allows us to introduce one additional real degree of freedom, i.e., |γ m+1 | 2 . The new system of equations reads and comprises c 2 equations and 2m+2 unknowns. Therefore, we expect this to generally admit a solution if c 2 ≤ 2m + 2. Note that there is no benefit in adding more than one auxiliary mode.

B. Dual of the PPM code
An example of code that can be reduced to PPM by applying phase space displacement is the following α 1 = |0, α, α . . . , α , where the state α j has the vacuum on mode j and a coherent state of given amplitude α in the other modes. Note that this code is mapped into a PPM code by displacing each mode by −α. Therefore, the linear receiver is optimal and saturates the global bound on the inconclusive event probability, P 0 = e −|α| 2 .

C. Random codes
In this Section, we illustrate the performance of linear receivers on random states. As an example, we consider c = 3 coherent states over m = 3 modes. We sample their amplitudes (α 1 , α 2 , α 3 ) uniformly from a sphere of radius √ n, where n = i |α i | 2 , where for simplicity we restrict to real-valued amplitudes. Figure 4(a) illustrates the minimized inconclusive probability P 0 achieved from our optimized framework with (35) and without (34) displacement. Notice that the dependence of P 0 on the average photon number is compatible with the exponential law obtained in the qualitative analysis of Section III. Without displacement, the linear receivers perform poorly. However, when equipped with displacement the performance improves significantly and nearly matches the global bound. To see this more clearly, Fig. 4(b) illustrates the statistical distributions of P 0 for n = 0.6 with N = 6600 random codes. The distribution of the inconclusive probability for the linear receivers with displacement closely matches the global bound, with small variations at smaller values of P 0 . The global bound has been computed using the method of Ref. [19]. This result reassuringly demonstrates that practical receivers based on linear optics are near-optimal and sufficient to match the global bound, as long as phase-space displacement operations are accessible.

D. Codes with single degeneracy
Random codes are typically non-degenerate. This means the coherent states in the random codes are almost surely linearly independent, i.e., the associated R is full rank. We now consider examples of degenerate codes. As discussed above, our receiver, endowed with phase-space displacement, can map a code with single degeneracy into a non-degenerate one. As an example, consider the code with single-degeneracy [47] α 1 = |α, α , α 2 = |α, −α , α 3 = |−α, α . (54) In Fig. 5(a), we illustrate the optimized inconclusive probabilities for linear receivers with displacement. Since these amplitude vectors are not linearly independent and have single degeneracy, receivers limited to LOP unitaries alone perform poorly, and we must leverage displacement operations. By allowing for phase-space displacement, our receivers match the global bound up to machine precision for small intensities (|α| 2 < ∼ 0.4). The analytic expression of the global bound can be obtained using the method of Ref. [19], which yields The numerical optimality of linear receivers extends to larger signal intensities for alternative figures of merit, as discussed in the next Section.

E. Alternative figures of merit
The gap between the performance of linear receivers and the global bound can be further reduced for alternative figures of merit. Consider for example the communication capacity associated with the given input code words and a given detection strategy. The coherent states in Eq. (54) can be used as code words for a communication protocol where the sender (Alice) uses these coherent states to encode a random variable X that takes values x = 1, 2, 3 with associated probabilities p X (x). The receiver (Bob) decodes this information by applying either the linear receiver or a globally optimal USD receiver. The outcome of the receiver is described by a random variable Y that takes four possible values, y = 0, 1, 2, 3, with probability p Y (y), where y = 0 is the inconclusive event. The maximum asymptotic communication rate achievable in this way is given by the Shannon capacity [48] where is the Shannon entropy of Y , and is the conditional entropy, where p Y (y|x) is the conditional probability of Y = y given X = x.
Since the receivers are unambiguous, P 0 (x) := p Y (y = 0|x) is the probability of inconclusive event for given input, and p(y|x) = δ yx [1 − P 0 (x)] for y = 1, 2, 3. The key quantity that determines the capacity is thus the conditional inconclusive probability P 0 (x). For linear receivers with displacement, this is given by From this we obtain where P 0 = x p X (x)P 0 (x) is the average inconclusive probability. A comparison between our scheme (with LOP and displacement) and the global bound is shown in Fig. 5(b). The two schemes yield nearly equal capacities, the gap being too small to be visualized in the scale of the plot.
In addition to the inconclusive probability and the Shannon capacity, we explore the optimality of our receiver using the finite communication block length rate, F . This rate is the communication rate attainable when both the code length, L, is finite, and there is a block error probability threshold, , imposed on the communication [49]. The normal approximation to the finite block length rate is given by [50] where Q(x) = 1/ √ 2π ∞ x dt exp[−t 2 /2], and V denotes the variance of the information transition probabilities of the channel, and . (63) Figure 6 illustrates the outcome of this maximization with different code lengths L and for difference receivers. We find that linear receivers match the finite-code-length performance generated from the global bounds. Note that in the asymptotic regime of large code block lengths, the finite rate tends towards the channel capacity with = 0, which are illustrated in horizontal lines.  II: Input-output transformation of the coherent state amplitudes through our linear receiver with two-mode photon detection. The receiver is shown in Fig. 7(a).

F. Codes with double degeneracy
To discriminate codes with higher degeneracy we may exploit joint detection events on pair of output modes. Consider the single-mode code with double degeneracy Single-mode codes find important applications in quantum sensing and communications protocols where individual measurements are preferred at each instance. Given the double degeneracy, our displacement-based receiver is not able to unambiguously discriminate signal states in this code. However, we can still discriminate these coherent states with LOP unitaries, vacuum auxiliary modes, displacement, and on-off detection by constructing receivers that take advantage of double detection events as described in Section IV C. The minimal inconclusive event probability of this class of receivers is determined as the solution to optimisation (47). By noting the equivalence between arbitrary auxiliary modes and vacuum auxiliary modes followed by mode-wise displacement operations, the specific receiver that saturates the performance of linear receivers for this code is illustrated in Fig. 7(a). It requires two auxiliary modes: one prepared in the vacuum and the other in a coherent state with amplitude α/ √ 2 (as noted above using coherent state ancillas is equivalent to using vacuum ancillas and phase-space displacements). The modes are then mixed at two 50/50 beam splitters, before the output is measured by mode-wise on-off photodetection. For an input signal amplitude α j chosen from (64), we write the output coherent state amplitude at mode k as ζ j k . This input-output transformation on the coherent state amplitudes is summarised in Table II. By inspection of the table, it is evident that the coherent states can be unambiguously discriminated if two detectors click. For example, a joint detection on modes one and two identifies the input code word |α 1 = |−α without error. From table II, the probability of a joint photodetection event in modes k and l can be determined through For this receiver, an inconclusive event corresponds to the absence of a joint photodetection event. Given a prior distribution p j , the average probability of obtaining an inconclusive event is then given by which for a uniform prior, p j = 1/3, yields This result can be compared with the global bound, which we compute analytically using Ref. [19] (see Appendix D for proof): The performance of the linear receivers is benchmarked with the global bound in Fig. 7(b). We note that the behavior of P 0 for small |α| 2 is compatible with the quadratic law in Eq. (21).

G. Phase-shift keying
A commonly used way to encode information into coherent states is by phase modulation, i.e., phase-shift keying. Previous works have considered practical schemes for USD of phase-shifted coherent states, both with and without feedback [24,27,28]. In this Section we analyze USD of phase-shifted coherent states using our linear receivers, re-obtaining some of the results of van Enk [24].
Binary phase-shift keying (BPSK) is a code with two states on one mode (m = 1, c = 2): It is well known how to optimally discriminate BPSK, here we express this method in the framework of our linear receivers. Indeed, to discriminate these states we can apply the strategy described in Section IV B. Note that these states are not linearly independent in phase-space.
To make them independent, first we add an auxiliary vacuum mode, generating the two-modes states |α, 0 , | − α, 0 . Second, we displace the auxiliary mode by α, which yields two coherent states with linear independent amplitude vectors, namely |α, α and |−α, α . Finally, we note that the Gram matrix of these latter state is diagonal, thus the code is equivalent to the PPM code. Applying the argument of Section V A, we conclude that BPSK code is optimally discriminated with linear receivers. In general, we may consider M -PSK codes with M states over a single mode: for j = 0, . . . , M − 1. Note that for M ≥ 3 the code has multiple degeneracies. This means that it can be discriminated only using multiple detection events. For example, the 3PSK (M = 3) ensemble of coherent states can be discriminated by exploiting double detection events. For this, we make use of two auxiliary vacuum modes and construct linear receivers according to the strategy described in Section IV C. The best performing linear receiver saturates optimisation (47) and is shown in Fig. 8 along with the global bound. Notice again that we recover the quadratic dependency of P 0 with small |α| 2 as predicted by Eq. (21). For USD of 3PSK it is easy to find the analytical expression for the optimized probability of the inconclusive event for linear receivers As discussed in Ref. [24], this probability can be improved by using more ancillary modes. Finally, the global bound, obtained from [19], reads

VI. CONCLUSIONS
While quantum states cannot in general be discriminated against without error, they can be discriminated unambiguously if we allow for a non-zero probability of an inconclusive outcome. The theory of USD allows us to identify the globally optimal measurements and to compute the ultimate bounds on the probability of obtaining an inconclusive outcome, see for example [12,18,19]. However, existing literature has rarely explored USD under realistic experimental constraints and limitations in what measurements are currently feasible, especially for the discrimination of coherent states of the quantum electromagnetic field.
Here we have outlined a theory of USD for multimode coherent states that focuses on practical receivers that can be feasibly realized with current technologies. Our receivers operate under physical resources described entirely through linear optical passive (LOP) unitaries, phase-space displacement operations, auxiliary vacuum modes, and on-off mode-wise photon detection. We have benchmarked the performance of these linear receivers against the global bound and found examples where they are near-optimal. In particular, this happens for random codes. This shows that high-performance USD receivers can be readily realized with currently available technologies, and suggests that, at least for typical, nondegenerate codes, high order non-linearities, feedback, or more advanced quantum technologies may only provide small improvements.
We do not have a complete theoretical explanation of why linear receivers perform so well. However, we think that the reason may be related to the fact that, at least for small amplitudes, the probability of multiple detection events is suppressed. Therefore, using a photon number resolving detector may not be particularly useful unless the code is degenerate. Furthermore, as the photon number is subject to quantum fluctuations in coherent states, it is unlikely that more detailed information on photon statistics can improve USD (whereas it may useful for ASD). Also, we have seen that linear receivers achieve the global bound for the PPM codes, and all codes that can be reduced to them by linear optics. By symmetry, and by an argument of concentration of measure, we may indeed expect that, up to statistical fluctuations, a random code is not too different from a PPM code: this could explain why linear receivers work pretty well on typical codes.
A number of questions remain open. First, an analytical expression for the inconclusive probability, beyond our numerical results, would characterize more clearly the comparison with the global bound, and allow us to apply our analysis to an arbitrary number of modes and code words. In the regime of weak signals, this problem may be more naturally framed within the language of Poisson quantum information [51]. Second, an extension of our theory is required to handle highly degenerate codes, when feedback operations could be useful. This is especially important when the number of states is much larger than the number of optical modes. In this regime, which is of particular interest for quantum communications, USD may be achieved only by exploiting joint photodetection events on multiple modes. Finally, our work might be formulated as a resource theory, similar to Refs. [52,53]. This approach may provide insight and help in comparing different sets of resources, for example by including homodyne or heterodyne detection, photon addition, and subtraction, or some mild non-linear interactions.

Acknowledgments
This work was supported by the EPSRC Quantum Communications Hub, Grant No. EP/T001011/1. SG acknowledges the NSF Center for Quantum Networks, awarded under grant number EEC-1941583. MSB acknowledges the University of Arizona's Information in a Photon course for being introduced to the problem of optimal linear-optic USD receiver designs.

Appendix A: Peres-Terno theory for USD
In this Appendix we review the theory of Peres and Terno [12] for the USD of a set of linearly independent vectors |u j with j = {1, . . . , c}, associated with prior probabilities p j (satisfying c j=1 p j = 1). These vectors span a c-dimensional Hilbert space H c . We define a unique set of (not necessarily normalized) vectors |v j ∈ H c such that We use these vectors to define a POVM with n elements for j = 1, . . . , n. The POVM element corresponding to an inconclusive event is The parameters k j are chosen in such way to ensure A 0 ≥ 0, and I is the identity in H c . For a suitable choice of the parameters k j 's, this POVM allows for unambiguous discrimination. The corresponding probability of an inconclusive outcome is A globally optimal USD measurement corresponds to the one that minimizes P 0 subject to the positivity of A 0 . This is obtained as the solution to the constrained maximization problem: which defines the globally optimal bound on the probability of the inconclusive outcome.

Appendix B: Method of Bergou, Futschik, and Feldman
Consider the c × c matrix with elements Note that the diagonal elements are the probabilities of an inconclusive event conditioned on the vector |u j , and the off-diagonal elements are simply the overlaps between the code words, for i = j.
If the vectors |u j are linearly independent, the condition of non-negativity of the operator A 0 , A 0 ≥ 0, is equivalent to C ≥ 0. Therefore, following the theory of Bergou et al. [19], the optimal average inconclusive probability is obtained by solving the constrained optimisation minimize q 1 , . . . , q j c j=1 p j q j , This problem can be solved analytically or semianalytically. The case of c = 3 is discussed in detail in Ref. [19].
We can immediately apply this method to show that PPM codes are optimally discriminated by mode-wise photodetection. To prove this, recall that the PPM codes are defined such that the matrix R is square, with R = αI and I the c-dimensional identity matrix. The off-diagonal entries of the matrix C are all equal to e −|α| 2 , while the diagonal entries are all equal, i.e., C jj = q. Therefore, the objective function in the minimization (B4) is equal to q. The eigenvalues of C are (c − 1)e −|α| 2 + q (with multiplicity 1) and q − e −|α| 2 (with multiplicity c − 1). The smallest value of q such that the eigenvalues remain non-negative is therefore q = e −|α| 2 , which matches the inconclusive event probability obtained through modewise photodetection in Section V A.
where I is the identity matrix. By multiplying V by U and V we obtain another unitary matrix, which is a unitary extension of M 0 , As M 0 is an extension of M , it follows that U is a unitary extension of M . We conclude that a matrix M can be extended into a unitary, if and only if its singular eigenvalues are not larger than 1. This condition can equivalently be written as Note that with this construction, the unitary matrix has at most size 2m. For our application, this means that we need at most m auxiliary modes to implement the receiver. In particular, the minimum number of auxiliary modes equals the number of singular values of M 0 that are strictly smaller than 1.

(D5)
Note that the first constraint follows from the positivity of q j , the second from the positivity of the failure POVM Π 0 , and the third from the linear independence (LI) of the states. Collectively, these constraints ensure the nonnegativity of C in Eq. (B4).
To analytically solve (D5), note that the failure scaled probabilities satisfyq :=q 1 =q 2 . This symmetry reduces the general three-dimensional feasible region in {q 1 ,q 2 ,q 3 } to a two-dimensional feasible region in {q,q 3 } defined by a cross-section along the planeq 1 =q 2 . The optimization then amounts to finding the minimum solutions to the objective function along the interior, edge, and vertices that defines this two-dimensional feasible constraint region. To see this explicitly, note that the constraints in (D5) constrain the feasible range forq,q 3 that collectively satisfỹ q ∈ [1, e 2|α| 2 ] ,q 3 ∈ [0, e −|α| 2 ] ,q 3 ≥ 2 q + 1 .
A comparison of the two values is illustrated in Fig. 10(a), which shows that P 0 (A) ≤ P 0 (B). Therefore, we discard the point B from the search for the minimised P 0 value. Next, we minimize P 0 over the interior points between A to B alongq 3 = 2(q + 1) −1 : P 0 (q) = 2 3 e −2|α| 2q + e |α| 2 q + 1 , which amounts to a single parameter optimization. This function has a stationary point atq in = e 3|α| 2 /2 − 1 yielding The feasibility of this solution is conditioned on the stationarity point q in residing within the allowed domain between vertices A and B, which holds if and only if This solution is illustrated together with P 0 (A) in Fig. 10(b). The solution P in 0 is plotted only when e |α| 2 ≥ 4. The two solution crosses exactly at |α| 2 = ln(4). Hence, the global minimum for the double-degeneracy (DD CB is characterized by P 0 (A) for e n < 4, and P in 0 for e n ≥ 4: which concludes our proof.