Particle–Particle Entanglement Overview

• Particle–particle entanglement is the quantum correlation between constituent particles or particle-defined subsystems, distinct from mode or spatial entanglement.
• It is characterized using reduced density matrices and concurrence, with its quantification depending on the choice of subsystem decomposition and experimental context.
• The topic spans various regimes including identical particle frameworks, many-body systems, and relativistic field theory, influencing both theoretical analysis and practical measurement techniques.

Particle–particle entanglement denotes entanglement attributed to constituent particles, or to particle-defined subsystems, rather than to spatial regions or field modes alone. In a bipartite pure state $|\psi\rangle \in \mathcal{H}_A\otimes\mathcal{H}_B$, it is diagnosed by the reduced density matrix $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$: if $\rho_A$ is mixed while the global state is pure, $A$ and $B$ are entangled; for mixed states, entanglement is defined by nonseparability, $$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$ (Buniy et al., 2012). In many-body systems, however, the meaning of “particle–particle entanglement” depends on the tensor-product structure one chooses—particle, mode, spatial, or operational—and this dependence becomes especially consequential for identical particles, relativistic fields, and systems subject to superselection constraints (Wasak et al., 2016).

1. Foundational definitions and subsystem structure

A precise discussion of particle–particle entanglement requires distinguishing several subsystem decompositions. One decomposition is the spatial or laboratory bipartition,

$\mathcal H = \mathcal H_A \otimes \mathcal H_B,$

which underlies ordinary bipartite entanglement between separated laboratories. Another is the particle decomposition,

$\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$

where entanglement is defined with respect to constituent particles themselves (Wasak et al., 2016). A many-body state may be separable in one decomposition and entangled in another, so the phrase “particle–particle entanglement” is not interchangeable with “mode entanglement” or “spatial entanglement.”

For multipartite systems, the distinction between direct pairwise entanglement and global many-body entanglement is essential. Two small subsystems $X$ and $Y$ are directly entangled only if their joint reduced density matrix $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$0 is nonseparable. A global pure state can nevertheless be highly entangled even when $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$1 is separable for every small pair. In that case the entanglement is distributed across many degrees of freedom rather than concentrated in pairwise correlations (Buniy et al., 2012).

The broader literature also separates particle–particle entanglement from intraparticle entanglement. In a single massive particle, different degrees of freedom such as spin and spatial mode can define a tensor-product structure,

$\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$2

and one then obtains entanglement between external and internal degrees of freedom rather than between two distinct particles (Aguilar, 2020). This distinction is mathematically straightforward but operationally significant.

2. Identical particles, overlap, and operational entanglement

For identical bosons and fermions, the standard label-based construction obscures which correlations are physically accessible. A label-free particle-based framework avoids fictitious particle labels and defines reduced states through multi-particle probability amplitudes, wedge products, and dot products, while retaining familiar notions such as partial trace, reduced density matrices, von Neumann entropy, concurrence, and Bell tests (Compagno et al., 2018). Within this framework, particle entanglement is defined via reduced states of subsets of particles, even though the particles are not individually labelable.

Spatial overlap plays a central role. In the label-free treatment, independently prepared identical qubits with overlapping spatial wavefunctions can yield Bell-violating spin correlations between two separated regions $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$3 and $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$4. The projected state with one particle in each region takes the form

$\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$5

with concurrence

$\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$6

and the maximal CHSH violation

$\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$7

which exceeds $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$8 whenever the overlap amplitudes are all nonzero (Compagno et al., 2018). This suggests that indistinguishability becomes operationally relevant only when combined with overlap and localized measurements.

A quantitatively sharper result is available for two identical bosons with pseudospins and partially overlapping spatial wavefunctions. After projecting onto the sector with one boson at detector $\rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|$9 and one at detector $\rho_A$0, and tracing out an auxiliary distinguishability degree of freedom $\rho_A$1, the concurrence of the detected two-spin state is

$\rho_A$2

where $\rho_A$3 and $\rho_A$4 are spatial amplitudes and $\rho_A$5 quantifies indistinguishability (Barros et al., 2019). The amount of entanglement is therefore a monotonically increasing function of both spatial overlap and indistinguishability. If either factor vanishes, the particle entanglement vanishes.

Superselection rules sharpen the operational distinction further. When local operations cannot create or detect coherences between different local particle-number sectors, mode entanglement alone is not enough to violate a Bell inequality. In that setting, every particle-separable state becomes effectively $\rho_A$6–$\rho_A$7 separable under the allowed local operations, so no Bell inequality can be violated. The required particle entanglement may arise solely from indistinguishability, as in the identical-boson beamsplitter example analyzed in that framework (Wasak et al., 2016).

3. Many-body, continuum, and lattice formulations

In continuum many-body systems, particle entanglement is naturally defined by a particle bipartition rather than a spatial cut. For $\rho_A$8 identical bosons with coordinates $\rho_A$9, one partitions the particles into subsets $A$0 and $A$1 with $A$2, $A$3, traces out the $A$4 particles, and defines Rényi entropies

$A$5

For a contiguous bosonic ground state that is a Bose–Einstein condensate, the non-interacting state has no particle entanglement, although spatial entanglement remains nonzero because of particle-number fluctuations across spatial bipartitions (Herdman et al., 2013). In direct path-integral Monte Carlo simulations of interacting one-dimensional bosons, the single-particle entanglement entropy exhibits logarithmic scaling with $A$6, and the coefficient increases with interaction strength, saturating to unity in the strongly interacting limit (Herdman et al., 2013).

The same work connects particle entanglement to experimentally accessible one-body information. If $A$7 is the condensate fraction, then for the single-particle second Rényi entropy,

$A$8

one has

$A$9

with tighter bounds in the two-mode limit (Herdman et al., 2013). This identifies condensate depletion as a direct constraint on particle bipartition entropies.

In rotating Bose and Fermi gases, particle entanglement is encoded in one- and two-particle reduced density matrices,

$B$0

with entropies $B$1 and $B$2 (Liu et al., 2010). Oscillations of $B$3 as a function of angular momentum $B$4 indicate both particle localization and vortex localization, and an abrupt change in fermionic $B$5 identifies edge reconstruction (Liu et al., 2010). For special bosonic subspaces,

$B$6

the single-particle entanglement scales numerically as

$B$7

which was used to relate these states to fractional quantum Hall trial states (Liu et al., 2010).

A complementary lattice formulation appears in one-dimensional two-particle continuous-time quantum walks governed by Hubbard Hamiltonians. There, particle–particle entanglement is defined from the entropy of the singular values of the time-evolved coefficient matrix: for distinguishable particles,

$B$8

and for indistinguishable bosons,

$B$9

where $$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$0 are eigenvalues of the one-particle reduced density matrix (Mastandrea et al., 1 Jun 2026). The long-time limits of these entanglement measures are typically non-monotonic as the onsite repulsion increases (Mastandrea et al., 1 Jun 2026).

A minimal continuous-variable example is the two-body Hooke’s-law model of a particle coupled to a lattice degree of freedom. There the reduced density matrix is obtained by integrating out the lattice coordinate in a path integral, and the spatial entanglement of the ground state is quantified by the linear entropy

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$1

Increasing the confining potential increases the spatial entanglement in that model (Puttarprom et al., 2013).

4. Generation mechanisms and entanglement witnesses

Particle–particle entanglement need not arise only from direct interactions. A recent photonic experiment demonstrated that independent photons can be entangled by path identity, without direct interaction, prior established entanglement, or Bell-state measurements. In the weak-pumping regime of four type-II SPDC sources, postselection on one photon in each of four paths yields

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$2

so photons $$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$3 and $$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$4 are projected into

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$5

conditioned on the ancillary photons (Wang et al., 2024). Experimentally, the protocol achieved

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$6

for CHSH, a reconstructed Bell-state fidelity

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$7

and concurrence

$$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$8

in the two-ancilla version; in the one-ancilla version the fidelity was $$\rho_{AB}\neq \sum_i w_i\,\rho_A^{(i)}\otimes\rho_B^{(i)}$$9, the witness value was $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$0, and the concurrence was $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$1 (Wang et al., 2024). The mechanism is indistinguishability of emission histories rather than interaction or entanglement swapping.

In mixed fermion–boson many-body systems, a different strand of work characterizes collective particle entanglement through off-diagonal long-range order. The particle–hole reduced density matrix

$\mathcal H = \mathcal H_A \otimes \mathcal H_B,$2

is positive semidefinite, and the largest eigenvalue $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$3 serves as an entanglement witness (Avdic et al., 19 Dec 2025). Values $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$4 indicate that multiple particle–hole pairs occupy the same excitonic mode, and larger $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$5 indicates stronger ODLRO. Sector-resolved witnesses

$\mathcal H = \mathcal H_A \otimes \mathcal H_B,$6

track fermionic and bosonic ODLRO separately, while the interaction

$\mathcal H = \mathcal H_A \otimes \mathcal H_B,$7

transfers coherence between particle types (Avdic et al., 19 Dec 2025). In the mixed Lipkin–Meshkov–Glick model, this framework captures the onset of collective entanglement when order initially localized in one sector becomes shared between both (Avdic et al., 19 Dec 2025).

5. Relativistic fields, particle creation, and non-inertial observers

In quantum field theory, particle–particle entanglement depends explicitly on the choice of physical variables and on how the subsystems are selected. For particle creation in time-dependent backgrounds, the fundamental mode functions obey parametric-oscillator equations, and an initial vacuum evolves into two-mode squeezed states for momentum pairs $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$8. In a particle basis defined by Bogoliubov coefficients $\mathcal H = \mathcal H_A \otimes \mathcal H_B,$9, the entanglement entropy per pair is

$\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$0

whereas in the field-amplitude variables $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$1 the same Gaussian state can be separable (Lin et al., 2010). This shows that particle entanglement and field-mode separability are not contradictory statements; they correspond to different subsystem choices.

For uniformly accelerated observers, the relevant particles are Rindler rather than Minkowski excitations. In a Dirac field, Minkowski particle operators mix Rindler particle and antiparticle operators, so acceleration redistributes entanglement between particle and antiparticle sectors. For fermions, this redistribution plays a key role in the survival of field entanglement in the infinite-acceleration limit (Martin-Martinez et al., 2011). In a charged bosonic field, by contrast, the same kind of redistribution does not prevent entanglement from vanishing in the infinite-acceleration limit (Bruschi et al., 2012). The contrast was traced to the different tensor-product structure induced by fermionic and bosonic statistics.

These relativistic results imply that “particle–particle entanglement” in field theory is observer-dependent and basis-dependent. An inertial particle–particle Bell pair can become, for an accelerated observer, a mixed state with support on particle and antiparticle sectors, and the accessible entanglement then depends on whether the detector is sensitive to both sectors or only one (Martin-Martinez et al., 2011).

6. Cosmological typicality and geometric reformulations

A cosmological application of typical-state arguments leads to a striking separation between global entanglement and pairwise entanglement. If the pre-inflationary region that later expanded into our observable universe was in local thermal equilibrium and in a typical pure state in the constrained Hilbert space, then any small subsystem $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$2 has a reduced state close to the canonical state,

$\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$3

with deviation probability bounded by

$\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$4

when $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$5 (Buniy et al., 2012). The entanglement entropy of any reasonably small region is then near maximal. Yet for two small subsystems $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$6 and $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$7, if

$\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$8

the reduced state $\mathcal H = \bigotimes_{i=1}^N \mathcal H_i,$9 is exponentially likely to be separable (Buniy et al., 2012). In this cosmological setting, a typical particle is entangled with many particles far outside our horizon, while two randomly chosen particles are overwhelmingly unlikely to be directly entangled with each other (Buniy et al., 2012).

A separate conceptual line reformulates two-particle entanglement equations geometrically. In the pilot-wave limit for two entangled particles, the wavefunction is written as

$X$0

and the associated Bohmian equations can be mapped to classical equations on a conformally stretched $X$1-dimensional configuration-space spacetime or, alternatively, to a Finsler geometry (Wasay et al., 2017). In these dual descriptions, the conformal factor is identified with the amplitude $X$2, the Hamilton principal function is matched to $X$3, and the quantum potential is encoded geometrically (Wasay et al., 2017). A related proposal treats the scalar potential in a stationary two-particle Schrödinger equation as equivalent to an isometric entanglement of the position coordinates in a complexified coordinate system with shared complex time,

$X$4

so that the potential term disappears in the transformed coordinates (Ducharme, 2010). These constructions do not introduce new entanglement measures; they recast entangled particle dynamics as geometry on a higher-dimensional or complexified space.

Particle–particle entanglement therefore spans several logically distinct regimes. In few-body quantum mechanics it is a reduced-state property of constituent particles. For identical particles it is an operational resource controlled by overlap, indistinguishability, and superselection. In many-body systems it is a particle-bipartition property that complements spatial and mode entanglement. In field theory it depends on the observer’s particle notion. In cosmology it can be globally overwhelming while remaining pairwise negligible. Across these regimes, the recurring technical lesson is that the existence, magnitude, and operational accessibility of particle–particle entanglement are determined less by the word “particle” alone than by the subsystem structure, the measurement algebra, and the dynamical or kinematical constraints under which the particles are defined.