Boson-fermion complementarity in a linear interferometer: An identity relating the determinant and permanent of a matrix
Pith reviewed 2026-05-24 04:50 UTC · model grok-4.3
The pith
Bosonic and fermionic transition probabilities in any linear interferometer are linked by a single equation that also connects the squared moduli of the permanent and determinant of a matrix.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The bosonic and fermionic transition probabilities appear together in a single equation which constrains their values, hence expressing a boson-fermion complementarity that is independent of the details of the interferometer; this also provides a heretofore unknown mathematical identity connecting the squared moduli of the permanent and determinant of an arbitrary complex matrix.
What carries the argument
The unifying equation that places bosonic and fermionic transition probabilities together for particles evolving under a unitary single-particle transformation.
If this is right
- For two particles the average of bosonic and fermionic probabilities must coincide with the classical probability for any unitary interferometer.
- The complementarity relation holds for every possible linear interferometer and is therefore independent of the concrete unitary matrix.
- The same equation yields an identity relating the squared modulus of the permanent of an arbitrary complex matrix to the squared modulus of its determinant.
- The identity extends Muir's nineteenth-century result to a new form connecting permanent and determinant.
Where Pith is reading between the lines
- The relation could be used to verify the correct implementation of bosonic or fermionic statistics without precise knowledge of the interferometer unitary.
- The mathematical identity may supply new bounds or computational relations between permanents and determinants that apply outside physics.
- Similar complementarity equations might exist when more than two particles are involved or when the interferometer includes controlled nonlinear elements.
Load-bearing premise
The interferometer is described by a unitary matrix on the single-particle space and the many-particle states are formed by standard symmetrization for bosons or antisymmetrization for fermions with no extra interactions or loss.
What would settle it
Measure the two-particle bosonic and fermionic transition probabilities through any linear interferometer and check whether their average equals the corresponding classical distinguishable-particle probability; any systematic deviation falsifies the claimed relation.
Figures
read the original abstract
Bosonic and fermionic statistics are well known to give rise to antinomic behaviors, most notably boson bunching vs fermion antibunching. Here, we establish a fundamental relation that combines bosonic and fermionic multiparticle interferences in an arbitrary linear interferometer. The bosonic and fermionic transition probabilities appear together in a single equation which constrains their values, hence expressing a boson-fermion complementarity that is independent of the details of the interferometer. For two particles in any interferometer, for example, it implies that the average between the bosonic and fermionic probabilities must coincide with the probability obeyed by classical particles. Crucially, this fundamental relation also provides a heretofore unknown mathematical identity connecting the squared moduli of the permanent and determinant of an arbitrary complex matrix, hence extending an identity by Muir dating from the nineteenth century.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims to establish a relation combining bosonic and fermionic multiparticle transition probabilities in an arbitrary linear interferometer, yielding a complementarity constraint independent of interferometer details. For two particles this implies that the average of the bosonic and fermionic probabilities equals the classical probability. The same construction is asserted to produce a new algebraic identity relating |per(M)|^2 and |det(M)|^2 for an arbitrary complex matrix M, extending Muir's nineteenth-century identity.
Significance. If the central identity holds, the work supplies both a physically motivated constraint on boson-fermion interference and a previously unknown relation between the squared moduli of the permanent and determinant. The derivation from standard symmetrized/antisymmetrized amplitudes in a unitary interferometer, together with the explicit two-particle verification, constitutes a clear strength; the result is falsifiable by direct matrix computation for small dimensions.
minor comments (3)
- [Abstract and §1] The precise algebraic form of the claimed identity (e.g., the coefficient relating |per(M)|^2 and |det(M)|^2) is stated only qualitatively in the abstract and introduction; an explicit equation should appear in the main text before the generalization is asserted.
- [§2] Notation for the submatrix A extracted from the unitary interferometer should be defined once at the outset (including its dimensions relative to the number of modes and particles) to make the passage from the physical setup to the arbitrary-matrix statement transparent.
- [§4] A short appendix or remark verifying the identity by direct expansion for a generic 3×3 complex matrix would strengthen the generalization claim without altering the manuscript's scope.
Simulated Author's Rebuttal
We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No specific major comments appear in the provided report.
Circularity Check
No significant circularity identified
full rationale
The paper derives the claimed identity directly from the standard definitions of bosonic (permanent) and fermionic (determinant) amplitudes under a unitary linear interferometer and symmetrization/antisymmetrization, without any fitted parameters, self-referential definitions, or load-bearing self-citations. The two-particle case reduces to an algebraically verifiable relation for arbitrary complex matrices, and the generalization is presented as following from the same construction. This makes the central result a mathematical consequence of the inputs rather than a reduction to them by construction, with the derivation self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Many-particle states are obtained by symmetrization (bosons) or antisymmetrization (fermions) of products of single-particle amplitudes under a unitary matrix describing the interferometer.
Reference graph
Works this paper leans on
-
[1]
A. Einstein, in Sitzungsberichte der Preussis- chen Akademie der Wissenschaften , Physikalisch- mathematische Klasse, Berlin (1925) p. 314
work page 1925
-
[2]
C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 59, 2044 (1987)
work page 2044
- [3]
-
[4]
H. J. Kimble, M. Dagenais, and L. Mandel, Phys. Rev. Lett. 39, 691 (1977)
work page 1977
-
[5]
Pauli, Zeitschrift f¨ ur Physik31, 765 (1925)
W. Pauli, Zeitschrift f¨ ur Physik31, 765 (1925)
work page 1925
- [6]
- [7]
-
[8]
T. Rom, T. Best, D. van Oosten, U. Schneider, S. F¨ olling, B. Paredes, and I. Bloch, Nature 444, 733 (2006)
work page 2006
- [9]
-
[10]
T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. Monroe, and J. L. O’Brien, Nature 464, 45 (2010)
work page 2010
- [11]
-
[12]
D. Vasyukov, Y. Anahory, L. Embon, D. Halbertal, J. Cuppens, L. Neeman, A. Finkler, Y. Segev, Y. Mya- soedov, M. L. Rappaport, M. E. Huber, and E. Zeldov, Nature Nanotechnology 8, 639 (2013)
work page 2013
-
[13]
S. Aaronson and A. Arkhipov, in Proceedings of the Forty-Third Annual ACM Symposium on Theory of Com- puting, STOC ’11 (Association for Computing Machin- ery, New York, NY, USA, 2011) p. 333–342
work page 2011
-
[14]
M. Tillmann, B. Daki´ c, R. Heilmann, S. Nolte, A. Sza- meit, and P. Walther, Nat. Photonics 7, 540 (2013)
work page 2013
- [15]
-
[16]
M. A. Broome, A. Fedrizzi, S. Rahimi-Keshari, J. Dove, S. Aaronson, T. C. Ralph, and A. G. White, Science 339, 794 (2013)
work page 2013
-
[17]
J. Carolan, J. D. A. Meinecke, P. J. Shadbolt, N. J. Rus- sell, N. Ismail, K. W¨ orhoff, T. Rudolph, M. G. Thomp- son, J. L. O’Brien, J. C. F. Matthews, and A. Laing, Nat. Photonics 8, 621 (2014)
work page 2014
-
[18]
M. Oszmaniec, N. Dangniam, M. E. Morales, and Z. Zim- bor´ as, PRX Quantum3, 020328 (2022)
work page 2022
-
[19]
M. C. Tichy, M. Tiersch, F. Mintert, and A. Buchleitner, New Journal of Physics 14, 093015 (2012)
work page 2012
-
[20]
M. C. Tichy, P. A. Bouvrie, and K. Mølmer, Phys. Rev. Lett. 109, 260403 (2012)
work page 2012
-
[21]
M. C. Tichy and K. Mølmer, Phys. Rev. A 96, 022119 (2017)
work page 2017
- [22]
-
[23]
M. G. Jabbour and N. J. Cerf, Phys. Rev. Research 3, 043065 (2021)
work page 2021
-
[24]
N. J. Cerf and M. G. Jabbour, Proceedings of the National Academy of Sciences 10.1073/pnas.2010827117 (2020)
-
[25]
C. Weedbrook, S. Pirandola, R. Garc´ ıa-Patr´ on, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, Rev. Mod. Phys. 84, 621 (2012)
work page 2012
- [26]
-
[27]
Permanents in linear optical networks
S. Scheel, Permanents in linear optical networks (2004), arXiv:quant-ph/0406127 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[28]
P. A. MacMahon, Combinatory Analysis, Volumes I and II, 3rd ed., Vol. 137 (American Mathematical Society, 2001)
work page 2001
- [29]
-
[30]
Here, ||v||p = P l |vl|p 1/p and ||A||2 = max||v||2 ||Av||2
-
[31]
Muir, Proceedings of the Royal Society of Edinburgh 22, 134–136 (1899)
T. Muir, Proceedings of the Royal Society of Edinburgh 22, 134–136 (1899). 7 Appendix
-
[32]
X k ukiukσi − X k ukiukσi xσi zk #! = X σ∈SN sgn(σ) NY i=1
Generating function for the transition probabilities We define the 2N-variate generating function of the sequence B(i) n (with i, n ∈ NN) to be g(x, z) = Ti,n h B(i) n i (x, z) := X i∈NN X n∈NN B(i) n NY s=1 xis s ! NY r=1 znr r ! , (29) where x, z ∈ {[0, 1)}N. In order to prove Lemma 1, we begin by showing that g(x, z) = 1 det [1N − U †ZU X] . (30) This ...
-
[33]
Recurrence relation for the transition probabilities Proof of Lemma 2. Similarly as in the proof of Lemma 1, here we write M(α, β) to denote the submatrix of M obtained by keeping the rows (columns) whose indices belong to the subset α (β), and we write M(α) for M(α, α). Furthermore, we write R(N) m for Rm and 1(N) α for 1α (see main text) to make the dep...
-
[34]
Consider some L ≥ 2N and ε ≤ 1/||A||2 (L and ε always exist since N is finite)
General relations for matrices Proof of Theorem 2. Consider some L ≥ 2N and ε ≤ 1/||A||2 (L and ε always exist since N is finite). From Lemma 29 in Ref. [13], there exists a unitary matrix V ∈ CL×L that contains εA as a submatrix. Consider such a matrix V constructed from A. Define the vectors p, q ∈ NL +, as follows: if α (β) is the vector corresponding ...
-
[35]
F urther physical applications of Theorem 1 Consider first the trivial case where particles are all sent and received in the first mode (of course, the same is true for any other mode, not necessarily the same for the input and output). Eq. (9) yields B(n) n − B(n−1) n−1 F (1) 1 = 0. (111) By solving this recurrence and using the fact that F (1) 1 = C(1) ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.