Spatial Biphoton JPDs in Quantum Optics

Updated 2 January 2026

Spatial Biphoton Joint Probability Distributions represent the spatial correlations of entangled photon pairs, defined as the squared magnitude of the biphoton wavefunction.
They are experimentally reconstructed via spatially resolved coincidence measurements in near- and far-field setups, revealing momentum conservation and entanglement properties.
Analyzing these distributions yields insights into phase-matching, pump beam effects, and supports advances in quantum imaging and high-dimensional communication.

Spatial biphoton joint probability distributions quantify the spatial correlations between photon pairs generated by nonlinear quantum optical processes such as spontaneous parametric down-conversion (SPDC) and spontaneous four-wave mixing (SpFWM). The joint distribution $P(\mathbf{r}_1, \mathbf{r}_2)$ or its momentum-space counterpart $P(\mathbf{k}_1, \mathbf{k}_2)$ encodes the likelihood of detecting one photon at position (or momentum) $\mathbf{r}_1$ and its partner at $\mathbf{r}_2$ , and serves as the central observable for spatial entanglement, entanglement quantification, and quantum imaging protocols.

1. Theoretical Foundation and Mathematical Structure

The generic biphoton wavefunction in the spatial domain is $\Psi(\mathbf{r}_1, \mathbf{r}_2)$ , where $\mathbf{r}_i$ denotes the transverse spatial coordinate of photon $i$ . The spatial biphoton joint probability distribution (JPD) is defined as $P(\mathbf{r}_1,\mathbf{r}_2) = |\Psi(\mathbf{r}_1,\mathbf{r}_2)|^2$ (Zheng et al., 2024). The wavefunction arises from first-order perturbative treatment of the SPDC or SpFWM Hamiltonian under the undepleted pump approximation (Perkins, 2011, Schneeloch et al., 2015, Yuan et al., 2019). For collinear, degenerate type-I SPDC and a Gaussian pump, the wavefunction in the transverse momentum basis is

$\Phi(\mathbf{q}_1,\mathbf{q}_2) = N \ \textrm{sinc}[\alpha|\mathbf{q}_1-\mathbf{q}_2|^2] \exp[-\sigma_p^2|\mathbf{q}_1+\mathbf{q}_2|^2]$

with $\alpha$ determined by phase-matching and crystal parameters (Schneeloch et al., 2015). Under typical experimental conditions, the momentum-space JPD exhibits an elongated ridge of anti-correlation, transforming under Fourier transform in the near-field to a JPD peaked along the diagonal $\mathbf{r}_1 \approx \mathbf{r}_2$ (McFadden et al., 31 Dec 2025, Procopio et al., 2016).

Many SPDC sources admit a "double-Gaussian" approximation,

$\Psi(\boldsymbol{\rho}_1,\boldsymbol{\rho}_2) \propto \exp\left(-\frac{|\boldsymbol{\rho}_1+\boldsymbol{\rho}_2|^2}{4\sigma_+^2} - \frac{|\boldsymbol{\rho}_1-\boldsymbol{\rho}_2|^2}{4\sigma_-^2}\right)$

with the JPD

$P(\boldsymbol{\rho}_1,\boldsymbol{\rho}_2) = \exp\left(-\frac{|\boldsymbol{\rho}_1+\boldsymbol{\rho}_2|^2}{2\sigma_+^2} - \frac{|\boldsymbol{\rho}_1-\boldsymbol{\rho}_2|^2}{2\sigma_-^2}\right)$

where $\sigma_+$ and $\sigma_-$ parameterize the center-of-mass and relative-position widths (Zheng et al., 2024, Schneeloch et al., 2015, Procopio et al., 2016). The covariance matrix and Pearson correlation coefficients extracted from $P$ reveal the strength and nature of spatial entanglement.

2. Joint Probability Distributions in Experiment and Theory

Spatial biphoton JPDs are reconstructed experimentally by spatially resolved coincidence measurements in the near-field (image plane) or far-field (momentum plane). In the thin-crystal limit, the near-field JPD approximates

$P(\mathbf{r}_1, \mathbf{r}_2) \propto |E_p\left(\frac{\mathbf{r}_1 + \mathbf{r}_2}{2}\right)|^2 \, \delta^{(2)}(\mathbf{r}_1 - \mathbf{r}_2)$

where $E_p$ is the pump amplitude, and the $\delta$ is broadened in finite-length crystals to a narrow Gaussian of correlation width $\sigma_x$ governed by longitudinal phase-matching (McFadden et al., 31 Dec 2025). In momentum space,

$P(\mathbf{q}_1, \mathbf{q}_2) \propto |\widetilde{E}_p(\mathbf{q}_1+\mathbf{q}_2)|^2 \, \delta^{(2)}(\mathbf{q}_1+\mathbf{q}_2)$

up to finite divergence and angular bandwidth (McFadden et al., 31 Dec 2025, Perkins, 2011).

For engineered pump profiles or arbitrary pump spatial modes, $E_p$ is decomposed into spatial mode functions, e.g., Elegant Gauss-Hermite or Laguerre-Gaussian modes, and the resulting JPD reflects both the pump structure and the phase-matching function (Perkins, 2011, Dehghan et al., 2024).

In integrated photonic platforms such as 1D quadratic waveguide arrays, the biphoton JPD translates to $P(n_s, n_i) = |\Psi_{n_s,n_i}(L)|^2$ for the amplitude $\Psi_{n_s,n_i}$ to detect the signal in site $n_s$ and the idler in $n_i$ (Gräfe et al., 2012). In ultrathin nonlinear films supporting SpFWM, the absence of longitudinal phase-matching constraints allows a spatial JPD with higher Schmidt number and increased dimensionality (Yuan et al., 2019).

3. Experimental Measurement and Numerical Construction

Experimentally, spatial biphoton JPDs are measured by accumulating spatially resolved coincidence maps on EMCCD, sCMOS, or time-stamping cameras (Reichert et al., 2018, McFadden et al., 31 Dec 2025, Zheng et al., 2024). In photon-counting cameras, the joint probability $P_{ij}$ for photons arriving at pixels $i$ and $j$ is extracted from the measured singles $C_i$ , coincidences $C_{ij}$ , and known noise/background probabilities via formulas such as

$T_{ij} = \frac{C_{ij} - C_i C_j}{a(1-C_i)(1-C_j)}$

with $a$ dependent on source and detector configuration (Reichert et al., 2018). Signal-to-noise ratio optimization requires balancing photon flux, read noise, and background counts according to detector statistics (Reichert et al., 2018).

For sCMOS detectors operating above the photon-count regime, joint distributions are reconstructed from second-order normalized intensity correlations,

$g^{(2)}(\mathbf{r}_1,\mathbf{r}_2) = \frac{\langle \delta I(\mathbf{r}_1)\delta I(\mathbf{r}_2)\rangle}{\bar{I}(\mathbf{r}_1) \bar{I}(\mathbf{r}_2)} \propto P(\mathbf{r}_1, \mathbf{r}_2)$

with background suppression techniques to remove detector artifacts (McFadden et al., 31 Dec 2025).

Phase retrieval and wavefront-sensing protocols (quantum Shack-Hartmann, computational phase retrieval) reconstruct not just the JPD but the full complex biphoton wavefunction modulus and phase, enabling complete spatial state tomography with only spatially resolved second-order measurements (Dehghan et al., 2024, Zheng et al., 2024).

4. Key Physical Features and Parameter Dependence

The shape, scale, and orientation of the spatial JPD reflect underlying physical constraints:

Momentum Conservation and Correlation: The $(\mathbf{q}_1+\mathbf{q}_2)$ dependence enforces momentum anti-correlation; near-field spatial correlation length is set by the pump waist and phase-matching bandwidth (Schneeloch et al., 2015, Procopio et al., 2016).
Crystal Length and Phase-Matching: The width of the phase-matching function (e.g., sinc or Gaussian) sets the conditional width of position- or momentum-correlation, with explicit analytic formulas for the transverse correlation width $\sigma_c$ as functions of crystal length, wavelength, and refractive index (e.g., $\sigma_c = \sqrt{L_z \lambda_p/(18\pi n)}$ for peak-matched Gaussian fits) (Schneeloch et al., 2015).
Pump Beam Structure: Tailored spatial pump modes (e.g., OAM superpositions, aberrated Gaussians) imprint corresponding structures on the JPD, supporting engineered spatial entanglement and multidimensional encodings (Perkins, 2011, Dehghan et al., 2024).
Propagation Dynamics: Free-space evolution transforms the joint amplitude via diffraction, causing amplitude correlations to migrate (decay) and phase correlations to emerge at characteristic propagation distances ( $z = k \sigma_+ \sigma_-$ ) (Zheng et al., 2024, Dehghan et al., 2024).
Integrated Structures: In waveguide arrays, the real-space JPD encodes biphoton quantum walks and their nonclassical correlations, which can be simulated via classical 2D beam propagation (Gräfe et al., 2012).

Covariance analysis of the bivariate Gaussian JPD provides direct access to the entanglement (Schmidt number, mutual information) and the principal axes correspond to the dominant spatial correlation directions—e.g., major axis at 45° in (x, x) (horizontal-horizontal) geometry (Procopio et al., 2016).

5. Quantification and Applications of Spatial Entanglement

The features of $P(\mathbf{r}_1, \mathbf{r}_2)$ directly relate to high-dimensional entanglement quantification:

Schmidt Number: Defined as $K = (\sum_n \lambda_n)^2/(\sum_n \lambda_n^2)$ from the Schmidt decomposition of the spatial amplitude, with $K \gg 1$ for highly entangled and wideband states (e.g., for ultrathin films) (Yuan et al., 2019).
Birth Zone and Mutual Information: The “biphoton birth zone” ( $\sigma_c$ ) quantifies the region where conditional probability is high, with mutual information $I(r_1:r_2) \approx \ln (\sigma_p/\sigma_c)$ for strong entanglement (Schneeloch et al., 2015).
Nonclassicality Tests: Direct violation of Bell-type inequalities can be demonstrated via intensity correlations mapped from the spatial JPD in both quantum and classical 2D photoic lattice simulators (Gräfe et al., 2012).
EPR-Steering Criteria: The product of conditional variances extracted from $P(x_1, x_2)$ and $P(q_1, q_2)$ enables direct checks of spatial EPR entanglement, e.g., $\sigma_{x_-}\sigma_{q_+} < \hbar/2$ (McFadden et al., 31 Dec 2025).

Applications include quantum imaging (ghost imaging, spatially resolved imaging of phase objects), metrology, quantum key distribution with high-dimensional alphabet, and adaptive optics using the phase profile extracted by quantum wavefront-sensing (Zheng et al., 2024).

6. Extensions and Recent Methodological Innovations

Modern experimental methods extend the reconstruction and characterization of spatial biphoton JPDs:

Phase Retrieval Algorithms: Maximum likelihood estimation and genetic algorithms enable reconstruction of both amplitude and phase of arbitrary spatial biphoton states from two-z-plane measurements, with modal decompositions (HyGG, Zernike) quantifying amplitude/phase features (Dehghan et al., 2024).
Quantum Wavefront Sensing: Quantum Shack-Hartmann sensing reconstructs both $P(\rho_1, \rho_2)$ and the partial derivatives of the biphoton phase, allowing full spatial state reconstruction including topological and phase correlations (Zheng et al., 2024).
Mesoscopic and Linear-Mode Imaging: High-flux, non-photon-counting measurements with sCMOS detectors can faithfully capture the JPD and EPR correlations, provided background normalization and detector artifacts are suppressed algorithmically (McFadden et al., 31 Dec 2025).
Ultrathin Nonlinear Media: SpFWM in ultrathin films enables the generation and measurement of spatiotemporally separable, highly multidimensional biphoton JPDs due to the absence of longitudinal phase-matching constraints (Yuan et al., 2019).

These advancements enable efficient, high-dimensional mapping of biphoton JPDs—amplitude and phase—across diverse experimental platforms, removing the need for interferometric stabilization, and supporting applications in quantum communication, imaging, and metrology.

7. Representative Summary Table: Analytical Forms and Physical Parameters

Source/system	JPD analytic form (position or $k$ -space)	Key parameters/regions
Bulk SPDC, Gaussian pump	$P(\mathbf{r}_1, \mathbf{r}_2) \sim \exp\left[-\frac{(\mathbf{r}_1+\mathbf{r}_2)^2}{2\sigma_+^2} - \frac{(\mathbf{r}_1-\mathbf{r}_2)^2}{2\sigma_-^2}\right]$	$\sigma_-$ (correlation width), $\sigma_+$ (pump width)
Thin-crystal limit	$P(\mathbf{r}_1, \mathbf{r}_2) \sim \|E_p\left(\frac{\mathbf{r}_1+\mathbf{r}_2}{2}\right)\|^2\, \delta^{(2)}(\mathbf{r}_1-\mathbf{r}_2)$	Pump waist $w_p$ , phase-matching, crystal thickness
Quadratic waveguide array	$P(n_s, n_i) = \|\Psi_{n_s,n_i}(L)\|^2$	Array size, coupling $C$ , pump distribution $A_n$
Ultrathin SpFWM film	$P(\mathbf{q}_s, \mathbf{q}_i) \sim \exp\left[-r_p^2\|\mathbf{q}_s+\mathbf{q}_i\|^2/4\right]$	No phase-matching filter, large $r_p$ , large Schmidt number

The explicit forms and scaling relationships of $\sigma_+, \sigma_-$ , as well as phase-matching envelopes (e.g., sinc or Gaussian), are system-specific and directly measurable in experiment or reconstructable via coincidence imaging and computational retrieval techniques (Schneeloch et al., 2015, Perkins, 2011, Dehghan et al., 2024).