Skip to main content
Sandbox Physics
Time-frequency quantum optics

L2 · Chapter 2

Describing a photon pair with complex amplitudes

Sandbox Physics ·

A pair source determines more than whether a pair appears. It also sets the possible ways to share frequency between the photons. The joint spectral amplitude stores these possibilities with their phases. Heralding and interference start from that amplitude.

Specify the events being retained

Spontaneous parametric down-conversion shares pump energy between two output fields. We call their paths signal and idler. At low gain the chance of producing a pair per pump pulse is small. We first condition on exactly one pair. Most pulses still contain vacuum; two-pair and higher contributions require a photon-number model.

The joint spectral amplitude, or JSA, is F(ωs,ωi)F(\omega_s,\omega_i) after specified spatial and polarization collection. The subscripts label the two arms. Normalization makes the probabilities for all frequency combinations add to one. It does not set the pair rate.

∣Ψ⟩=∬F(ωs,ωi)∣ωs⟩s∣ωi⟩i dωs dωi,∬∣F∣2 dωs dωi=1|\Psi\rangle=\iint F(\omega_s,\omega_i)|\omega_s\rangle_s|\omega_i\rangle_i\,d\omega_s\,d\omega_i,\qquad \iint|F|^2\,d\omega_s\,d\omega_i=1
This is the conditional single-pair state. Absolute count rates also need generation probability, loss, efficiency and repetition rate.

Reading: P. J. Mosley et al. · A. Bhattacharjee et al.

Contributions from the crystal must add coherently

The pump spectrum restricts the sum of the two output frequencies. Pair creation can occur throughout the crystal. Contributions from those positions also need compatible phases to add constructively. Phase matching describes this condition. A finite crystal produces a finite matching width and a phase. Collection optics then select the spatial and polarization amplitudes that reach the experiment.

A scalar propagation approximation organizes these factors as the product below. Here α\alpha is the complex pump envelope, Φ\Phi is the phase-matching amplitude and CC represents the specified collection operation. These are amplitudes, so multiplying intensity pictures alone loses information.

F∝α(ωs+ωi) Φ(ωs,ωi) C(ωs,ωi),Φ∝sinc⁡ ⁣(ΔkL2)eiΔkL/2F\propto\alpha(\omega_s+\omega_i)\,\Phi(\omega_s,\omega_i)\,C(\omega_s,\omega_i),\qquad \Phi\propto\operatorname{sinc}\!\left(\frac{\Delta kL}{2}\right)e^{i\Delta kL/2}
LL is length and Δk\Delta k is the longitudinal wave-vector mismatch including the chosen quasi-phase-matching grating term. The convention is sinc⁡x=sin⁡x/x\operatorname{sinc}x=\sin x/x. The expression assumes the stated propagation and collection approximation.

Reading: P. J. Mosley et al.

Replace the continuum with orthogonal modes

Choose two orthogonal frequency modes in each arm. They might be normalized wave packets in nonoverlapping frequency intervals. Rows label signal modes; columns label idler modes. Each matrix entry is a complex amplitude. Its squared magnitude is the probability of that mode pair.

The coefficient FjkF_{jk} is a projection of the continuous JSA onto mode functions. It is not simply a value at one frequency point. A sampled continuous integral needs quadrature weights absorbed into its discrete coefficients. Our two-mode experiment instead defines orthogonal discrete modes directly, with dimensionless coefficients.

∣Ψ⟩=∑j,k=01Fjk∣j⟩s∣k⟩i,∑j,k∣Fjk∣2=1|\Psi\rangle=\sum_{j,k=0}^{1}F_{jk}|j\rangle_s|k\rangle_i,\qquad \sum_{j,k}|F_{jk}|^2=1
Four complex numbers specify this single-pair model. They expose phase and conditioning without a full crystal-propagation calculation.

Four equally bright entries

Set all four amplitudes to one half. Each mode pair has probability one quarter. Now rotate only the lower-right arrow through an angle θ\theta. Every intensity stays fixed. The next chapter follows this change through to single-photon purity.

Pump, material and collection geometry determine which amplitudes a source can realize. Fixing one specified pair amplitude lets us ask how different herald measurements change the partner state. This separates the physical roles of source design and measurement selection.

F(θ)=12(111eiθ)F(\theta)=\frac12\begin{pmatrix}1&1\\1&e^{i\theta}\end{pmatrix}
Phase is in radians. Multiplying the whole matrix by one common phase changes no prediction. Rotating one entry changes coherent sums between rows.

Try the prediction

Match phase arrows to the intensity cells

  1. Choose the equal-phase preset and check that the four intensities sum to one.
  2. Rotate only the lower-right phase and predict which intensity map changes.
  3. Follow the signal coherence and explain how its state can change with fixed joint intensities.

Exercises and explanations

  1. Must two sources with the same normalized coefficient matrix deliver equal detected pair rates?

    Read the explanation

    No. The matrix describes mode probabilities conditioned on producing and collecting one pair. Repetition rate, pair probability, path loss and detector efficiency affect the number detected per second.

  2. How does multiplying every coefficient by ii differ from multiplying only the lower-right coefficient by ii?

    Read the explanation

    The first adds a common phase to the whole state and leaves all probabilities and density matrices unchanged. The second changes a relative phase that coherent measurements can detect.

Sources and calculation scope

  1. P. J. Mosley et al., Heralded Generation of Ultrafast Single Photons in Pure Quantum States, arXiv:0711.1054v1; Physical Review Letters 100, 133601 (2008)

    Equation (1), Figures 1–4 and the source/collection discussion connect pump, phase matching, heralding and independent-source HOM. Figure 2 explicitly qualifies intensity-based purity by assuming no joint spectral phase variation. No crystal fit or measured points are reproduced.

  2. A. Bhattacharjee et al., Coherent temporal filtering of multimode parametric down-conversion using a quantum pulse gate, arXiv:2608.12544v1 (12 August 2026 preprint)

    Main equations (1)–(5), Figures 1–4 and supplementary sections 1–7 were read for conditional filtering and reconstruction. The bounded reproduction is supplementary equations (12)–(13); (15)–(16) are evaluated analytically. We do not reproduce the measured Figure 3 points, tomography algorithm or full device.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy