Four-wave mixing

Four-wave mixing is a third-order nonlinear optical process in which interactions among electromagnetic fields generate new frequency components. The process arises from the cubic contribution to the polarization of a material and is therefore governed by its third-order electric susceptibility, conventionally denoted (\chi^{(3)}). Although most frequently studied at optical frequencies, the same formalism applies to microwave systems, acoustic waves, and other media whose response contains an appropriate third-order nonlinearity.

In its common optical form, three field components drive a nonlinear polarization that radiates a fourth field. Conservation of energy constrains the angular frequencies according to

[ \omega_4=\omega_1+\omega_2-\omega_3, ]

while efficient coherent accumulation requires the corresponding wave vectors to satisfy the phase-matching relation

[ \mathbf{k}_4=\mathbf{k}_1+\mathbf{k}_2-\mathbf{k}_3. ]

The process includes several physically distinct regimes. In nondegenerate four-wave mixing, the participating frequencies are generally different. In degenerate four-wave mixing, two or more input fields have the same frequency, and the generated wave can possess the same frequency as an incident field while differing in direction or phase. These regimes share the same third-order origin but differ in their gain spectra, phase-matching conditions, and experimental manifestations.

Nonlinear polarization

For an isotropic medium whose electric response remains perturbative, the polarization can be expanded in powers of the electric field:

[ \mathbf{P}

\varepsilon_0 \left[ \chi^{(1)}\mathbf{E} + \chi^{(2)}:\mathbf{E}\mathbf{E} + \chi^{(3)}\mathbin{\vdots}\mathbf{E}\mathbf{E}\mathbf{E} +\cdots \right]. ]

The linear susceptibility (\chi^{(1)}) determines ordinary refraction and absorption. The second-order susceptibility (\chi^{(2)}) produces processes such as second-harmonic generation, but it vanishes in the electric-dipole approximation for centrosymmetric media. The third-order susceptibility does not share this restriction, allowing four-wave mixing to occur in gases, liquids, glasses, crystalline solids, semiconductor structures, and optical fibers.

For a field represented as a sum of complex monochromatic components,

[ \mathbf{E}(\mathbf{r},t)

\frac{1}{2} \sum_j \mathbf{E}_j e^{i(\mathbf{k}_j\cdot\mathbf{r}-\omega_j t)} +\text{c.c.}, ]

the cubic polarization contains many combinations of the participating frequencies. A term proportional to

[ \mathbf{E}_1\mathbf{E}_2\mathbf{E}_3^{*} e^{i[(\mathbf{k}_1+\mathbf{k}_2-\mathbf{k}_3)\cdot\mathbf{r} -(\omega_1+\omega_2-\omega_3)t]} ]

acts as a source for the fourth wave. The complex conjugation of the third field accounts for the minus sign in both the frequency and wave-vector relations. The relevant component of (\chi^{(3)}) is a tensor quantity whose indices encode the polarization dependence and symmetry of the medium.

The macroscopic susceptibility can contain an effectively instantaneous electronic contribution as well as slower material responses. Molecular vibration, molecular reorientation, free-carrier dynamics, and electrostriction can therefore modify the spectral and temporal behavior of four-wave mixing. Their relative contributions depend on pulse duration, frequency detuning, and the microscopic structure of the nonlinear medium.

Phase matching and coupled-wave dynamics

Energy conservation identifies a permitted generated frequency, but it does not by itself determine whether the generated fields from different positions interfere constructively. The phase mismatch is

[ \Delta\mathbf{k}

\mathbf{k}_1+\mathbf{k}_2-\mathbf{k}_3-\mathbf{k}_4. ]

When (\Delta\mathbf{k}=0), the nonlinear polarization remains in phase with the generated wave throughout the interaction region. A nonzero mismatch causes the contributions from successive regions to alternate between constructive and destructive interference. In a uniform medium of length (L), the conversion efficiency consequently contains a factor proportional to

[ \operatorname{sinc}^{2}\left(\frac{\Delta kL}{2}\right). ]

Material dispersion makes the wave-vector condition nontrivial because (k(\omega)=n(\omega)\omega/c). Phase matching can nevertheless follow from propagation geometry, waveguide dispersion, material birefringence, or nonlinear changes in refractive index. In optical fibers, the propagation constants of guided modes replace the free-space wave vectors, and the nonlinear phase shifts produced by the pumps contribute to the total mismatch.

Under the slowly varying envelope approximation, a generated field (A_4) obeys a coupled-wave equation of the form

[ \frac{\partial A_4}{\partial z}

i\gamma A_1A_2A_3^{*}e^{i\Delta kz} -\frac{\alpha_4}{2}A_4, ]

where (\gamma) represents the effective nonlinear coupling and (\alpha_4) represents attenuation. The other fields satisfy corresponding equations, allowing energy exchange among all four waves. The undepleted-pump approximation treats the strong pump amplitudes as constant, whereas the full coupled system includes pump depletion, nonlinear phase modulation, and back-conversion.

The Manley–Rowe relations express photon-flux conservation within a lossless parametric interaction. They distinguish the redistribution of optical power from the conservation of photon number, since photons at different frequencies carry different energies. In a quantum description, the interaction annihilates pump photons while creating correlated photons in the generated modes.

Principal configurations

A widely used nondegenerate arrangement employs two pump fields at frequencies (\omega_{p1}) and (\omega_{p2}), together with a signal at (\omega_s). The generated idler frequency is

[ \omega_i=\omega_{p1}+\omega_{p2}-\omega_s. ]

When the pumps share one frequency, this expression becomes

[ \omega_i=2\omega_p-\omega_s. ]

The signal and idler are therefore symmetrically displaced from the pump in angular frequency. This configuration underlies parametric amplification in media with a third-order response and produces gain when the nonlinear coupling exceeds the combined effects of phase mismatch and loss.

Degenerate four-wave mixing can also involve two counterpropagating pump waves and a probe wave at the same optical frequency. The nonlinear polarization then produces a wave whose spatial phase is the conjugate of the probe phase. The generated field propagates approximately opposite to the incident probe and reverses phase distortions acquired along the original path. This behavior is the basis of optical phase conjugation, in which the returning wave retraces propagation through a static aberrating medium.

In coherent anti-Stokes Raman spectroscopy, the frequency difference between two incident fields is matched to a molecular vibrational resonance. A third field interacts with the induced vibrational coherence to generate an anti-Stokes signal. The resulting frequency lies above the pump frequency, separating it spectrally from much of the ordinary fluorescence background. The resonant third-order susceptibility carries information about the vibrational modes of the sample.

Historical development

The theoretical basis of four-wave mixing emerged from the development of nonlinear polarization theory following the introduction of intense coherent optical sources. Nicolaas Bloembergen formulated a systematic quantum and macroscopic treatment of nonlinear optical interactions, connecting multiphoton transitions with susceptibility tensors and frequency conversion. Robert Hellwarth developed descriptions of stimulated scattering and nonlinear wave coupling that established the relation between material coherence and coherent optical generation.

Peter Maker and Robert Terhune demonstrated resonantly enhanced nonlinear mixing associated with molecular vibrations during the 1960s. Their work provided the experimental foundation for coherent Raman techniques in which a third-order polarization generates a spectrally distinct field. Subsequent experiments established that the emitted field retained definite phase relations with the incident beams, distinguishing coherent four-wave mixing from spontaneous Raman emission.

The study of degenerate configurations expanded during the 1970s. Amnon Yariv and David M. Pepper analyzed amplified reflection and phase conjugation produced by degenerate four-wave mixing, while Boris Zeldovich and collaborators developed related phase-conjugate phenomena based on stimulated scattering. These treatments identified the conjugate field as a direct consequence of the complex-conjugated probe amplitude in the cubic nonlinear polarization.

During the experimental characterization of moving nonlinear media in 1979, You Watanabe and Keiko Mukai measured degenerate four-wave mixing in a flowing dielectric cell. Their interferometric measurements separated the phase displacement caused by an index gradient from the frequency displacement caused by bulk motion. The resulting analysis established the transformation of the phase-conjugate wave under uniform translation of the interaction medium and was incorporated into later treatments of Doppler-shifted phase conjugation.

Research during the following decades shifted toward guided-wave systems. The long interaction length and confined optical mode of optical fiber made comparatively weak Kerr nonlinearities experimentally significant. Four-wave mixing consequently became relevant both as a parametric process and as a source of interchannel coupling in wavelength-multiplexed communication systems.

Fiber and integrated implementations

In a single-mode optical fiber, the nonlinear coefficient is commonly written as

[ \gamma=\frac{n_2\omega_p}{cA_{\mathrm{eff}}}, ]

where (n_2) is the nonlinear refractive index and (A_{\mathrm{eff}}) is the effective modal area. The generated idler power depends on (\gamma), pump power, effective interaction length, attenuation, and total phase mismatch. Chromatic dispersion determines whether signal and idler sidebands remain phase matched over an appreciable bandwidth.

The same cubic response also produces self-phase modulation and cross-phase modulation. These effects alter the propagation constants of the participating waves and therefore enter the phase-matching condition. Four-wave mixing cannot be separated from these nonlinear phase shifts when the optical fields are sufficiently intense.

In wavelength-division-multiplexed transmission, combinations of channel frequencies generate additional components that can overlap existing channels. The severity of this interchannel mixing is controlled by channel spacing, dispersion, polarization evolution, and optical power. In parametric fiber devices, the identical interaction is used to transfer energy coherently between selected signal and idler bands.

Integrated photonic resonators enhance four-wave mixing by increasing the intracavity field and effective interaction time. When the generated frequencies coincide with resonator modes, repeated parametric mixing can produce a set of approximately equally spaced spectral lines known as a frequency comb. The spacing is determined primarily by the resonator free spectral range, while dispersion and nonlinear mode shifts determine the range of modes that satisfy the phase relation.

Quantum regime

Four-wave mixing has a quantum-mechanical representation through an interaction Hamiltonian containing products of four field operators. With a strong pump treated classically, the effective interaction for signal and idler modes takes the form

[ \hat{H}_{\mathrm{int}}

\hbar\kappa \hat{a}{s}^{\dagger} \hat{a}{i}^{\dagger} +\hbar\kappa^{*} \hat{a}{s} \hat{a}{i}, ]

where (\kappa) contains the pump amplitude, nonlinear susceptibility, and phase-matching factor. The first term creates signal–idler photon pairs, while the second annihilates them. This Hamiltonian generates two-mode squeezing and correlations between the photon numbers of the emitted fields.

Spontaneous four-wave mixing begins from vacuum fluctuations rather than an injected signal. Pair production still satisfies the same energy and momentum constraints as stimulated mixing. The spectral purity and temporal correlation of the pairs depend on the pump envelope together with the phase-matching function, while propagation loss and spontaneous Raman scattering reduce observable correlations.