---
title: Spatial Biphoton JPDs in Quantum Optics
url: https://www.emergentmind.com/topics/spatial-biphoton-joint-probability-distributions
type: topic
---

# Spatial Biphoton JPDs in Quantum Optics

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$ [2406.04973]. The wavefunction arises from first-order perturbative treatment of the SPDC or SpFWM Hamiltonian under the undepleted pump approximation [1102.3261][1502.06996][1903.10936]. 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 [1502.06996]. 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$ [2512.24878][1605.01347].

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 [2406.04973][1502.06996][1605.01347]. 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 [2512.24878]. 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 [2512.24878][1102.3261].

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 [1102.3261][2403.06905].

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$ [1205.4085]. In ultrathin nonlinear films supporting SpFWM, the absence of longitudinal phase-matching constraints allows a spatial JPD with higher Schmidt number and increased dimensionality [1903.10936].

## 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 [1802.00489][2512.24878][2406.04973]. 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 [1802.00489]. Signal-to-noise ratio optimization requires balancing photon flux, read noise, and background counts according to detector statistics [1802.00489].

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 [2512.24878].

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 [2403.06905][2406.04973].

## 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 [1502.06996][1605.01347].
- **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) [1502.06996].
- **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 [1102.3261][2403.06905].
- **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_-$) [2406.04973][2403.06905].
- **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 [1205.4085].

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 [1605.01347].

## 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) [1903.10936].
- **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 [1502.06996].
- **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 [1205.4085].
- **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$ [2512.24878].

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 [2406.04973].

## 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 [2403.06905].
- **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 [2406.04973].
- **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 [2512.24878].
- **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 [1903.10936].

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 [1502.06996][1102.3261][2403.06905].

---

**Citations:**  
[1102.3261]  
[1502.06996]  
[1802.00489]  
[2512.24878]  
[2406.04973]  
[1205.4085]  
[2403.06905]  
[1903.10936]  
[1605.01347]

Source: https://www.emergentmind.com/topics/spatial-biphoton-joint-probability-distributions