Skip to content

Position-Space Two-Particle States

A two-particle wavefunction is the position representation of a vector in a tensor-product Hilbert space. For two distinguishable particles moving on a line,

H12=L2(R)⊗L2(R)≅L2(R2).\mathcal H_{12} = L^2(\mathbb R)\otimes L^2(\mathbb R) \cong L^2(\mathbb R^2).

The wavefunction

Ψ(x1,x2)=⟨x1,x2∣Ψ⟩\Psi(x_1,x_2) = \langle x_1,x_2\vert\Psi\rangle

is a probability amplitude on configuration space. The variables x1x_1 and x2x_2 label the position outcomes for the two tensor factors. They are not automatically independent random variables; independence is a property of special states.

This page is the bridge from wave mechanics to composite-system language. The general theory of wavefunction probability densities belongs to Wavefunctions and Probability Density; the general tensor-product rule belongs to Tensor Products of Hilbert Spaces. Here the focus is how Ψ(x1,x2)\Psi(x_1,x_2) encodes product states, entanglement, partial traces, and exchange symmetry.

For a normalized two-particle pure state,

∫−∞∞dx1∫−∞∞dx2 ∣Ψ(x1,x2)∣2=1.\int_{-\infty}^{\infty}dx_1 \int_{-\infty}^{\infty}dx_2\, \lvert\Psi(x_1,x_2)\rvert^2 = 1.

The joint probability density for position measurements is

ρ(x1,x2)=∣Ψ(x1,x2)∣2.\rho(x_1,x_2) = \lvert\Psi(x_1,x_2)\rvert^2.

Thus the probability that particle 1 is found in region AA and particle 2 in region BB is

P(A,B)=∫Adx1∫Bdx2 ∣Ψ(x1,x2)∣2.P(A,B) = \int_A dx_1 \int_B dx_2\, \lvert\Psi(x_1,x_2)\rvert^2.

The marginal density for particle 1 is

ρ1(x1)=∫−∞∞dx2 ∣Ψ(x1,x2)∣2,\rho_1(x_1) = \int_{-\infty}^{\infty}dx_2\, \lvert\Psi(x_1,x_2)\rvert^2,

and similarly for particle 2. These marginals describe position statistics only. They do not contain the full reduced quantum state, because phase coherence between different positions is carried by a density-matrix kernel.

In three spatial dimensions, replace xix_i by ri\mathbf r_i and dxidx_i by d3rid^3r_i. The configuration space of two particles is six-dimensional, even when ordinary physical space is three-dimensional.

A pure product state of two distinguishable particles has the form

∣Ψ⟩=∣ψ⟩1⊗∣ϕ⟩2.\lvert\Psi\rangle = \lvert\psi\rangle_1\otimes\lvert\phi\rangle_2.

In position representation this becomes

Ψ(x1,x2)=ψ(x1)ϕ(x2).\Psi(x_1,x_2) = \psi(x_1)\phi(x_2).

If ψ\psi and ϕ\phi are normalized, then Ψ\Psi is normalized:

∫dx1 dx2 ∣ψ(x1)ϕ(x2)∣2=(∫dx1 ∣ψ(x1)∣2)(∫dx2 ∣ϕ(x2)∣2)=1.\int dx_1\,dx_2\, \lvert\psi(x_1)\phi(x_2)\rvert^2 = \left( \int dx_1\,\lvert\psi(x_1)\rvert^2 \right) \left( \int dx_2\,\lvert\phi(x_2)\rvert^2 \right) = 1.

The joint density factorizes:

∣Ψ(x1,x2)∣2=∣ψ(x1)∣2∣ϕ(x2)∣2.\lvert\Psi(x_1,x_2)\rvert^2 = \lvert\psi(x_1)\rvert^2 \lvert\phi(x_2)\rvert^2.

This is stronger than saying the particles have no position correlation. The wavefunction itself factorizes, including phases. A state can have a factorized probability density but still fail to be a product state if the phase contains nonlocal dependence on both variables.

For example,

Ψ(x1,x2)=ψ(x1)ϕ(x2)eiαx1x2\Psi(x_1,x_2) = \psi(x_1)\phi(x_2)e^{i\alpha x_1x_2}

has the same position density as ψ(x1)ϕ(x2)\psi(x_1)\phi(x_2), but it is not generally a product wavefunction because the phase does not split into a sum of a function of x1x_1 and a function of x2x_2.

A two-particle pure state is entangled across the particle split when its wavefunction cannot be written as one factor depending only on x1x_1 times one factor depending only on x2x_2.

A transparent example uses orthonormal one-particle wavefunctions u0,u1u_0,u_1 for particle 1 and v0,v1v_0,v_1 for particle 2:

Ψ(x1,x2)=12[u0(x1)v0(x2)+u1(x1)v1(x2)].\Psi(x_1,x_2) = \frac{1}{\sqrt2} \left[ u_0(x_1)v_0(x_2) + u_1(x_1)v_1(x_2) \right].

This is the continuous-variable analogue of a two-term Schmidt state. It is entangled because two nonzero matched product terms appear in orthonormal local bases.

More generally, a normalizable bipartite wavefunction may have a Schmidt expansion

Ψ(x1,x2)=∑npn un(x1)vn(x2),\Psi(x_1,x_2) = \sum_{n} \sqrt{p_n}\, u_n(x_1)v_n(x_2),

where

∫dx um∗(x)un(x)=δmn,∫dx vm∗(x)vn(x)=δmn,\int dx\,u_m^*(x)u_n(x) = \delta_{mn}, \qquad \int dx\,v_m^*(x)v_n(x) = \delta_{mn},

and

pn≥0,∑npn=1.p_n\ge0, \qquad \sum_n p_n=1.

The state is product exactly when only one Schmidt probability is nonzero. If two or more are nonzero, the state is entangled.

Continuous variables allow infinite Schmidt rank and continuous-spectrum subtleties. For ordinary square-integrable wavefunctions, the reduced density operator should be trace class before entropies are used. Ideal distributions, such as exact EPR states, must be treated as limiting models rather than as ordinary vectors.

The density operator of a pure two-particle state is

ρ12=∣Ψ⟩⟨Ψ∣.\rho_{12} = \lvert\Psi\rangle\langle\Psi\rvert.

Its position-space kernel is

ρ12(x1,x2;x1′,x2′)=Ψ(x1,x2)Ψ∗(x1′,x2′).\rho_{12}(x_1,x_2;x_1',x_2') = \Psi(x_1,x_2)\Psi^*(x_1',x_2').

Tracing out particle 2 gives the reduced density operator of particle 1. In position representation,

ρ1(x,x′)=∫−∞∞dy Ψ(x,y)Ψ∗(x′,y).\rho_1(x,x') = \int_{-\infty}^{\infty}dy\, \Psi(x,y)\Psi^*(x',y).

The diagonal of this kernel is the marginal position density:

ρ1(x,x)=∫dy ∣Ψ(x,y)∣2.\rho_1(x,x) = \int dy\, \lvert\Psi(x,y)\rvert^2.

The off-diagonal entries carry coherence between different positions of particle 1. They are why the reduced density operator is more informative than the marginal density alone.

For a product wavefunction Ψ(x,y)=ψ(x)ϕ(y)\Psi(x,y)=\psi(x)\phi(y),

ρ1(x,x′)=ψ(x)ψ∗(x′)∫dy ∣ϕ(y)∣2=ψ(x)ψ∗(x′).\rho_1(x,x') = \psi(x)\psi^*(x') \int dy\, \lvert\phi(y)\rvert^2 = \psi(x)\psi^*(x').

Particle 1 remains in a pure state. For the two-term Schmidt example,

ρ1(x,x′)=12u0(x)u0∗(x′)+12u1(x)u1∗(x′),\rho_1(x,x') = \frac12 u_0(x)u_0^*(x') + \frac12 u_1(x)u_1^*(x'),

so particle 1 is mixed and the entanglement entropy is log⁡2\log2.

For identical particles without spin, physical two-particle wavefunctions must have definite exchange symmetry:

Ψ(x2,x1)=+Ψ(x1,x2)for bosons,\Psi(x_2,x_1) = +\Psi(x_1,x_2) \quad \text{for bosons},

and

Ψ(x2,x1)=−Ψ(x1,x2)for fermions.\Psi(x_2,x_1) = -\Psi(x_1,x_2) \quad \text{for fermions}.

For particles with spin or other internal labels, exchange acts on the complete one-particle label

q=(x,s,…),q=(x,s,\ldots),

not just on position. The symmetry condition is imposed on Ψ(q1,q2)\Psi(q_1,q_2).

The canonical construction is developed in Symmetric and Antisymmetric Wavefunctions. The main caution here is conceptual: the slot labels 11 and 22 used in first-quantized notation are not directly observable particle identities for identical particles. Exchange symmetry by itself should not be confused with operational entanglement between two addressable subsystems.

For identical particles, mode language is often cleaner. A two-boson state with one excitation in each of two orthogonal modes is naturally described in Fock space, while a first-quantized symmetrized wavefunction can make the same state look artificially like a superposition over particle labels. The physical question should specify whether the subsystems are particles, modes, spatial regions, or internal degrees of freedom.

For equal masses, define

R=x1+x22,r=x1−x2.R = \frac{x_1+x_2}{2}, \qquad r = x_1-x_2.

The inverse transformation is

x1=R+r2,x2=R−r2,x_1 = R+\frac r2, \qquad x_2 = R-\frac r2,

and the measure is

dx1 dx2=dR dr.dx_1\,dx_2 = dR\,dr.

Many two-body Hamiltonians separate in RR and rr, so solutions often have the form

Ψ(x1,x2)=Φ(R)χ(r).\Psi(x_1,x_2) = \Phi(R)\chi(r).

This is product in the collective variables R,rR,r, not necessarily product in the particle variables x1,x2x_1,x_2. The coordinate change mixes the two tensor factors:

R=x1+x22,r=x1−x2.R = \frac{x_1+x_2}{2}, \qquad r = x_1-x_2.

Therefore separability in center-of-mass and relative coordinates is a calculational statement, not automatically a statement about entanglement across the particle split. Conversely, an interaction that is simple in relative coordinates can generate correlations between particle positions.

The safe rule is to name the tensor-product structure before using words such as product, separable, or entangled. A change from (x1,x2)(x_1,x_2) to (R,r)(R,r) is not a local change of basis on H1⊗H2\mathcal H_1\otimes\mathcal H_2.

  • Treating a two-variable wavefunction as entangled merely because it depends on two variables.
  • Checking only whether the probability density factorizes and ignoring the phase.
  • Confusing marginal densities with reduced density operators.
  • Using ideal delta-correlated states as if they were normalizable wavefunctions.
  • Treating identical-particle slot labels as observable particle identities.
  • Forgetting spin or internal labels when imposing exchange symmetry.
  • Mistaking center-of-mass separability for particle separability.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • A. Peres, Quantum Theory: Concepts and Methods, Kluwer, 1995.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.
  1. Product normalization. Let Ψ(x1,x2)=ψ(x1)ϕ(x2)\Psi(x_1,x_2)=\psi(x_1)\phi(x_2), with ψ\psi and ϕ\phi normalized. Show that Ψ\Psi is normalized.
Solution

Compute

∫dx1 dx2 ∣Ψ(x1,x2)∣2=∫dx1 ∣ψ(x1)∣2∫dx2 ∣ϕ(x2)∣2.\int dx_1\,dx_2\, \lvert\Psi(x_1,x_2)\rvert^2 = \int dx_1\,\lvert\psi(x_1)\rvert^2 \int dx_2\,\lvert\phi(x_2)\rvert^2.

Each factor is 11, so the product is 11.

  1. Reduced kernel of a product state. For Ψ(x,y)=ψ(x)ϕ(y)\Psi(x,y)=\psi(x)\phi(y), compute ρ1(x,x′)\rho_1(x,x').
Solution

By definition,

ρ1(x,x′)=∫dy Ψ(x,y)Ψ∗(x′,y).\rho_1(x,x') = \int dy\, \Psi(x,y)\Psi^*(x',y).

Substitute the product form:

ρ1(x,x′)=ψ(x)ψ∗(x′)∫dy ∣ϕ(y)∣2=ψ(x)ψ∗(x′).\rho_1(x,x') = \psi(x)\psi^*(x') \int dy\,\lvert\phi(y)\rvert^2 = \psi(x)\psi^*(x').
  1. Two-term Schmidt wavefunction. For
Ψ(x1,x2)=12[u0(x1)v0(x2)+u1(x1)v1(x2)],\Psi(x_1,x_2) = \frac{1}{\sqrt2} \left[ u_0(x_1)v_0(x_2) + u_1(x_1)v_1(x_2) \right],

with orthonormal uu and vv functions, find the nonzero eigenvalues of ρ1\rho_1.

Solution

The wavefunction is already in Schmidt form with probabilities 1/21/2 and 1/21/2. Therefore the reduced density operator has two nonzero eigenvalues:

λ0=λ1=12.\lambda_0 = \lambda_1 = \frac12.

The entanglement entropy is log⁡2\log2.

  1. Exchange symmetry. Let
Ψ±(x1,x2)=12[u(x1)v(x2)±v(x1)u(x2)],\Psi_\pm(x_1,x_2) = \frac{1}{\sqrt2} \left[ u(x_1)v(x_2) \pm v(x_1)u(x_2) \right],

where uu and vv are orthonormal. Show that Ψ+\Psi_+ is symmetric and Ψ−\Psi_- is antisymmetric.

Solution

Exchange the coordinates:

Ψ±(x2,x1)=12[u(x2)v(x1)±v(x2)u(x1)].\Psi_\pm(x_2,x_1) = \frac{1}{\sqrt2} \left[ u(x_2)v(x_1) \pm v(x_2)u(x_1) \right].

Reordering scalar factors gives

Ψ±(x2,x1)=12[v(x1)u(x2)±u(x1)v(x2)].\Psi_\pm(x_2,x_1) = \frac{1}{\sqrt2} \left[ v(x_1)u(x_2) \pm u(x_1)v(x_2) \right].

For the plus sign this equals Ψ+(x1,x2)\Psi_+(x_1,x_2). For the minus sign it equals −Ψ−(x1,x2)-\Psi_-(x_1,x_2).

  1. Center-of-mass caution. Explain why Ψ(x1,x2)=Φ((x1+x2)/2)χ(x1−x2)\Psi(x_1,x_2)=\Phi((x_1+x_2)/2)\chi(x_1-x_2) is not necessarily a product state of particle 1 and particle 2.
Solution

A product state across the particle split must have the form

Ψ(x1,x2)=ψ(x1)ϕ(x2).\Psi(x_1,x_2) = \psi(x_1)\phi(x_2).

The expression Φ((x1+x2)/2)χ(x1−x2)\Phi((x_1+x_2)/2)\chi(x_1-x_2) factors in the collective variables R=(x1+x2)/2R=(x_1+x_2)/2 and r=x1−x2r=x_1-x_2. Those variables each depend on both particle coordinates. Unless the special functions Φ\Phi and χ\chi combine so that all cross-dependence cancels, the result cannot be written as one function of x1x_1 times one function of x2x_2.