Stimulated Parametric Down-Conversion

• StimPDC is the seeded counterpart of parametric down-conversion, using nonlinear χ(2) media to achieve phase-sensitive amplification and coherent idler generation.
• It utilizes three-wave mixing and phase conjugation to enable spatial mode conversion and advanced optical diagnostics, bridging differences between classical DFG and quantum regimes.
• The process facilitates precise quantum process tomography and spatial coding, making it a versatile tool for engineered squeezed states and turbulence mitigation in optical systems.

Stimulated parametric down-conversion (StimPDC) is the seeded counterpart of parametric down-conversion in a $\chi^{(2)}$ medium: the same three-wave mixing process that, in spontaneous PDC, generates signal–idler photons from vacuum fluctuations is instead driven by an injected signal and/or idler field, so that the interaction is observed through phase-sensitive amplification, depletion, mode conversion, and coherent idler generation. In the literature, StimPDC is treated both as seeded PDC and, in the classical high-flux limit, as difference-frequency generation; these descriptions are linked by the same nonlinear interaction kernel that governs spontaneous pair creation, stimulated gain/de-gain, phase conjugation, and spatiotemporal mode transfer (Avenhaus et al., 2014, Xu et al., 2023, Roux, 2021).

1. Process definition and formal descriptions

At the level of energy and momentum conservation, StimPDC obeys the usual three-wave relations

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$

or, in the difference-frequency convention used in several spatial experiments,

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$

with the precise form determined by the experimental labeling of signal and idler and by the phase-matching convention (Xu et al., 2023).

A quantum-optical description treats PDC as a unitary process generated by a three-wave mixing Hamiltonian. In the interaction picture and under the undepleted-pump, weak-coupling approximation, the process is written as

$$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$

where the complex-valued joint spectral amplitude (JSA) $f(\omega_\mathrm{s},\omega_\mathrm{i})$ is the kernel of the process (Avenhaus et al., 2014). In the multimode spatiotemporal description, the same interaction appears as a Bogoliubov transformation with kernel functions $U$ and $V$, so that the output field operators or phase-space variables are linear combinations of input annihilation and creation components across continuous frequency–wavevector space (Roux, 2021).

In the classical seeded limit, the process is described by coupled-wave equations. Under slowly varying envelope, perfect phase matching, and undepleted pump, the signal and idler envelopes satisfy

$$\frac{d\mathcal{E}_i}{dz} = \frac{2 i d_{eff} \omega_i^2}{k_i c^2} \,\mathcal{E}_p\, \mathcal{E}_s^*, \qquad \frac{d\mathcal{E}_s}{dz} = \frac{2 i d_{eff} \omega_s^2}{k_s c^2} \,\mathcal{E}_p\, \mathcal{E}_i^*,$$

or, equivalently, field transformations of Bogoliubov form,

$$E_s^{\text{out}} = \mu\, E_s^{\text{in}} + \nu\, E_i^{\text{in}*}, \qquad E_i^{\text{out}} = \mu\, E_i^{\text{in}} + \nu\, E_s^{\text{in}*},$$

with $|\mu|^2-|\nu|^2=1$ in the idealized lossless model (Xu et al., 2023, Aguilar-Cardoso et al., 3 Nov 2025). These formulae make explicit why StimPDC is simultaneously a parametric amplifier, a mode-converting interaction, and a phase-sensitive measurement channel.

2. Phase conjugation, advanced-wave optics, and spatial-mode conversion

The defining spatial signature of StimPDC is the phase-conjugation property. In the low-gain, thin-crystal limit, the idler transverse field is proportional to the product of the pump beam and the complex conjugate of the seed field,

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$0

so that, for a plane-wave or slowly varying pump, one obtains

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$1

which is the standard phase-conjugate-wave relation (Xu et al., 2023). In momentum space this becomes an anti-correlation relation, and in OAM space it appears as topological-charge inversion.

This phase-conjugation structure underlies Klyshko’s advanced-wave picture in seeded form. With a structured pump, the crystal no longer acts as an ideal flat mirror; instead, the pump angular spectrum plays the role of a reflective transfer function for the advanced wave. StimPDC then becomes a direct experimental proxy for coincidence-based SPDC geometries: the stimulated idler intensity reproduces the advanced-wave prediction that, in the spontaneous case, would otherwise be observed only in fourth-order correlations (Arruda et al., 2018). This is why StimPDC can test optical layouts designed by advanced-wave reasoning before running the corresponding single-photon experiment.

For first-order spatial modes, the phase-conjugation map admits a geometric representation on a spatial Poincaré sphere. When the seed is prepared as

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$2

the idler intensity corresponds to the transformed superposition with

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$3

that is, a specular reflection with respect to the equatorial plane of the sphere (Oliveira et al., 2018). For a Gaussian pump, OAM conservation reduces to $$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$4, and tilted-lens diagnostics reveal the corresponding inversion of topological charge (Oliveira et al., 2018). In the OAM-dichroism formulation, the gain is itself OAM-selective: for a pump with winding number $$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$5, the peak gain occurs at $$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$6, with an OAM bandwidth

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$7

so StimPDC can act as a differential-gain element in OAM space (Lowney et al., 2014).

The same phase-conjugation principle supports image coding and decoding. In orthogonal spatial coding, the seed carries an image with a Fourier-plane phase mask, the transferred idler inherits the conjugated field, and a corrective phase mask rotated by $$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$8 in the idler Fourier plane restores the image. The reconstruction condition is exact only under the plane-wave-pump approximation; if the pump is not a plane wave, the seed and idler momenta are not perfectly anti-correlated and complete decryption becomes impossible (Xu et al., 2023). This point is also central to all-optical turbulence mitigation, where StimPDC is used as a phase-conjugate stage to reverse turbulence-induced spatial distortions in free-space quantum links (Aguilar-Cardoso et al., 3 Nov 2025).

3. Time–frequency process viewpoint and stimulated tomography

In time–frequency language, the core object of PDC and StimPDC is the JSA

$$\omega_p = \omega_s + \omega_i, \qquad \mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i + \mathbf{K}_{\text{PM}},$$9

with a Gaussian pump envelope

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$0

and a phase-matching term

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$1

The modulus of $$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$2 sets which frequency pairs are coupled; its phase determines the temporal correlations and the local phase-sensitive response (Avenhaus et al., 2014).

StimPDC accesses this kernel directly. Seeding the signal and idler with narrowband coherent fields

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$3

the change in intensity of either seed obeys

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$4

where

$$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$5

The oscillation amplitude yields the JSA modulus, while the amplification/de-amplification phase yields the JSA phase (Avenhaus et al., 2014). This is the central reason stimulated emission tomography can reconstruct the full complex JSA rather than only the joint spectral intensity $$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$6.

Experimentally, the method combines ultrafast pulse shaping and linear detection. In the reported PP-KTP waveguide implementation, the seed spectra were shaped by a Fastlite Dazzler with spectral resolution $$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$7 and minimum seed bandwidth $$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$8; intensity oscillations over 100 pulses were Fourier-analyzed to recover $$\omega_i = \omega_p - \omega_s, \qquad \mathbf{k}_p - \mathbf{k}_s = \mathbf{k}_i,$$9, and specific phase scans were used to extract $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$0 (Avenhaus et al., 2014). The result is a time–frequency quantum process tomography of PDC with direct relevance to engineered squeezed states, heralded non-Gaussian states, and factorable photon-pair sources.

A complementary spatiotemporal theory based on Wigner functionals and Bogoliubov kernels shows that the output mean photon number separates into a stimulated term and a spontaneous background term. In this formalism, the transformed seed parameter function contains a signal part $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$1 and an idler part $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$2, while the detector response inherits the full continuous dependence on transverse wavevector, frequency, pump waist, pump bandwidth, and phase matching (Roux, 2021). This makes StimPDC a metrological probe of phase-matching geometry, absolute seed intensity, and multimode gain structure, not merely a source of brighter down-converted light.

4. Spatial-mode tailoring, beam quality, and coherence control

When the pump carries the target spatial mode and the seed is Gaussian, the idler profile at the crystal plane is controlled by the beam-waist ratio

$$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$3

For Gaussian pump and seed, the idler waist satisfies

$$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$4

so the idler approaches the pump beam size for $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$5 and the seed beam size for $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$6 (Aguilar-Cardoso et al., 19 Jul 2025). The idler Rayleigh range is likewise set by both $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$7 and the wavelength ratio, which means that the propagation behavior of the generated mode can be tuned by changing the seed waist without changing the pump mode itself (Aguilar-Cardoso et al., 19 Jul 2025).

The same product rule governs higher-order modes. For Laguerre–Gaussian and Hermite–Gaussian pump modes, the idler reproduces the intended structure with a fidelity

$$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$8

which improves monotonically as the seed becomes wider than the pump (Aguilar-Cardoso et al., 19 Jul 2025). The reported experiments show that increasing $$\hat{U}_\mathrm{PDC} = \exp\!\left[-i\left(\mathcal{B}\int d\omega_\mathrm{s}\, d\omega_\mathrm{i}\, f(\omega_\mathrm{s},\omega_\mathrm{i})\,\hat{a}^\dagger(\omega_\mathrm{s})\,\hat{b}^\dagger(\omega_\mathrm{i}) +\mathrm{h.c.}\right)\right],$$9 improves both the average fidelity and its uniformity across a mode set. For the angular basis

$f(\omega_\mathrm{s},\omega_\mathrm{i})$0

the fidelity depends on the basis dimension $f(\omega_\mathrm{s},\omega_\mathrm{i})$1 and $f(\omega_\mathrm{s},\omega_\mathrm{i})$2 but not on the mode label $f(\omega_\mathrm{s},\omega_\mathrm{i})$3, so all angular modes within a given basis dimension have the same transfer fidelity (Aguilar-Cardoso et al., 19 Jul 2025). This is why the angular basis is singled out as ensuring uniform fidelity across generated modes.

StimPDC also modifies coherence in a characteristic way. If the idler consists of a spontaneous component with mutual coherence $f(\omega_\mathrm{s},\omega_\mathrm{i})$4 and a stimulated component with coherence $f(\omega_\mathrm{s},\omega_\mathrm{i})$5, then

$f(\omega_\mathrm{s},\omega_\mathrm{i})$6

where

$f(\omega_\mathrm{s},\omega_\mathrm{i})$7

For a Gaussian Schell-model spontaneous contribution,

$f(\omega_\mathrm{s},\omega_\mathrm{i})$8

As the seed power increases, $f(\omega_\mathrm{s},\omega_\mathrm{i})$9 increases, the coherence curve acquires a nonzero floor, and the traditional coherence length can diverge for thresholds $U$0 (Santos et al., 2023). This is why the idler field in StimPDC cannot, in general, be described by a single Gaussian Schell-model beam: it is the sum of a partially coherent spontaneous part and a nearly fully coherent stimulated part (Santos et al., 2023).

5. Quantum and multimode regimes beyond bright classical seeding

Although StimPDC is often analyzed in the classical high-flux regime, several works emphasize that the same mechanism extends into explicitly quantum and multimode domains. In a doubly pumped bulk-crystal configuration, two non-collinear pump modes create overlapping phase-matching surfaces; shared modes experience local gain enhancement and appear as bright hot-spots, with a walk-off-controlled transition from three-mode to four-mode coupling. At resonance, the dominant gain scales with the Golden Ratio, and the corresponding quantum interpretation involves tripartite or quadripartite multimode entangled states in simple bulk sources (Jedrkiewicz et al., 2020).

In non-collinear coherent-beam-stimulated PDC, two-mode coherent seeding reshapes the multiphoton number sectors used for interferometric state engineering. For example, in the weak-field regime the reported scheme gives $U$1 for $U$2 at $U$3 and $U$4 for $U$5 at $U$6, while in the strong-field regime with $U$7 it maintains approximately $U$8 fidelity for $U$9 at $V$0 and approximately $V$1 fidelity for $V$2 at $V$3, with much higher usable photon flux than spontaneous-only schemes (Kolkiran, 2011). In this regime, StimPDC is not simply a brighter source; it is a structured quantum–classical resource for selecting specific $V$4-photon sectors.

At the opposite extreme, weak coherent seeding well below one photon per coherence time produces observable three-photon correlations. In the continuous-wave, low-gain degenerate experiment, the seeded output state contains a two-photon spontaneous term and a three-photon stimulated term,

$V$5

leading to

$V$6

After normalization to the accidental background, the predicted peak

$V$7

was experimentally accompanied by a measured value of approximately $V$8, interpreted as time-domain evidence of seed-induced three-photon correlations (Klein et al., 26 May 2026).

A related single-photon-seeded scenario appears in cascaded PDC for OAM conjugation. There, a signal photon from SPDC seeds a second StimPDC stage while the first idler acts as a herald. The stimulated term is constrained by $V$9, and in the cascaded coincidence distribution the phase-conjugation signature appears as a bias toward equal OAM between the heralding idler and the stimulated idler, $$\frac{d\mathcal{E}_i}{dz} = \frac{2 i d_{eff} \omega_i^2}{k_i c^2} \,\mathcal{E}_p\, \mathcal{E}_s^*, \qquad \frac{d\mathcal{E}_s}{dz} = \frac{2 i d_{eff} \omega_s^2}{k_s c^2} \,\mathcal{E}_p\, \mathcal{E}_i^*,$$0 (Luck et al., 1 Sep 2025). These results make clear that the boundary between “classical DFG” and “quantum StimPDC” is a regime boundary, not a change of process.

6. Applications, interpretation, and recurring misconceptions

StimPDC has become a methodological bridge between nonlinear optics and quantum-state engineering. In source characterization, it provides amplitude-sensitive quantum process tomography of PDC, with direct access to the complex JSA needed for two-mode squeezed vacuum, heralded non-Gaussian states, and entangled photon pairs with tailored spectral-temporal properties (Avenhaus et al., 2014). In spatial-mode engineering, it allows controlled generation of idler modes with tunable size, propagation behavior, and fidelity, including basis-uniform transfer in the angular basis (Aguilar-Cardoso et al., 19 Jul 2025). In imaging and coding, the phase-conjugated idler supports orthogonal spatial coding, image transfer across wavelengths, and aberration cancellation in Fourier space (Xu et al., 2023).

The same phase-conjugation mechanism has motivated free-space correction schemes. In turbulence-resilient high-dimensional quantum key distribution, a distorted stimulating signal is converted into a phase-conjugated idler that can reverse atmospheric phase distortions on backward propagation; the reported scheme is explicitly presented as all-optical dynamic correction of spatial-mode distortion, with reduced quantum error rates under strong turbulence (Aguilar-Cardoso et al., 3 Nov 2025). This suggests a broader role for StimPDC as a nonlinear-optical alternative to conventional adaptive optics in reciprocal channels.

Several misconceptions recur in the literature. One is that StimPDC is “just” classical difference-frequency generation. That statement is accurate in the high-flux, undepleted-pump, thin-crystal limit, but incomplete: the same interaction kernel also governs spontaneous PDC, the JSA phase, Bogoliubov mode mixing, and single-photon or weak-coherent-seed phenomena (Avenhaus et al., 2014, Klein et al., 26 May 2026). Another is that phase conjugation is automatically perfect. In fact, exact conjugation requires the specific approximations under which $$\frac{d\mathcal{E}_i}{dz} = \frac{2 i d_{eff} \omega_i^2}{k_i c^2} \,\mathcal{E}_p\, \mathcal{E}_s^*, \qquad \frac{d\mathcal{E}_s}{dz} = \frac{2 i d_{eff} \omega_s^2}{k_s c^2} \,\mathcal{E}_p\, \mathcal{E}_i^*,$$1 is valid—most notably a plane-wave or sufficiently broad pump, appropriate phase matching, and, in turbulence compensation, reciprocity of the propagation channel (Xu et al., 2023, Aguilar-Cardoso et al., 3 Nov 2025). A third is that spatial-mode fidelity is basis-independent. The reported mode-transfer studies show the opposite: LG and HG bases display order-dependent fidelity, whereas the angular basis is distinguished precisely because the fidelity is uniform across its mode labels for fixed basis dimension (Aguilar-Cardoso et al., 19 Jul 2025).

Taken together, these results establish StimPDC as a regime-spanning concept. It is simultaneously the classical seeded limit of PDC, a phase-sensitive amplifier with conjugation symmetry, a diagnostic for the complex structure of spontaneous sources, a platform for spatial and time–frequency mode engineering, and a route into nontrivial quantum-seeded and multimode dynamics.