Multicolor Solitons in Photonics

• Multicolor solitons are localized nonlinear waves spanning several carrier frequencies that remain bound through nonlinear coupling, dispersion management, and velocity matching.
• They are modeled using coupled equations such as the Lugiato–Lefever framework, three-wave, and NLS-type equations across Kerr microresonators, quadratic media, and integrable systems.
• Their binding mechanisms—including cross-phase modulation, four-wave mixing, and parametric trapping—pave the way for advanced comb generation, synthetic frequency lattices, and broadband photonic integration.

Multicolor solitons are localized nonlinear wave states whose energy is distributed across multiple carrier frequencies, spectral windows, or field components while remaining mutually trapped through nonlinear coupling, dispersion management, and velocity matching. In dissipative Kerr systems, they appear as composite dissipative Kerr solitons with a common repetition rate and distinct carrier or phase offsets; in quadratic media, as mutually trapped fundamental, second-harmonic, and higher-harmonic fields; and in multicomponent integrable models, as bright, dark, or mixed vector states across several components (Silvestri et al., 2024, Shandilya et al., 6 Feb 2026, Qi, 2018, Ramakrishnan et al., 2021). The term therefore denotes a family of bound solitary structures rather than a single canonical solution class.

1. Terminology, scope, and classification

Usage of “multicolor soliton” is system-dependent. In microresonator physics, a multicolor dissipative Kerr soliton is a composite state in which several spectral envelopes travel together with one repetition rate but retain distinct phase-rotation velocities or carrier-envelope offsets. In the language of synthetic-frequency lattices, the field is decomposed into wavepacket components $a_{-},a_{0},a_{+}$ that are nonlinearly coupled by four-wave-mixing Bragg scattering and cross-phase modulation (Moille et al., 2022).

In integrable and near-integrable multicomponent systems, “multicolor” often refers to component-wise bright, dark, or mixed occupancy rather than distinct optical carrier bands. The $(M+1)$-component Yajima–Oikawa system explicitly distinguishes “bright color” and “dark color” in short-wave components, and the multi-component Mel’nikov system uses the same bright–dark taxonomy for its short-wave fields (Chen et al., 2015, Han et al., 2017). In this literature, multicolor solitons include all-bright, all-dark, and mixed bright–dark $N$-soliton families.

A further usage appears in laser and waveguide contexts. “Polychromatic soliton molecules” in a mode-locked laser are coincident solitons with different frequencies but common group velocities, stabilized by an engineered cavity dispersion profile (Lourdesamy et al., 2020). Two-color pulse compounds in waveguides with two anomalous-dispersion windows are likewise described as molecule-like states built from two tightly bound subpulses with a large frequency gap (Melchert et al., 2022). This suggests that the unifying criterion is not the number of frequency bands alone, but the coexistence of localization and nontrivial inter-component binding.

2. Governing equations and reduced descriptions

The dominant microcomb description is the Lugiato–Lefever framework and its multicomponent reductions. A general multiscale theory for multicolor soliton microcombs begins from the non-normalized Lugiato–Lefever equation and derives a slow “meta-envelope” equation,

$$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$

where $\varepsilon_J$ is the nonlinear-enhancement factor and $\sigma_J$ is the pump coupling efficiency (Silvestri et al., 2024). In this description, each dispersion maximum supports one “color,” and the multicolor state is a slow envelope riding on a multi-peak carrier $c(\tau)$.

Dual-pumped Kerr resonators admit a more explicit three-color decomposition. Writing

$$a(\theta,t)=a_0(\theta,t)+a_-(\theta,t)e^{i\varpi_- t}+a_+(\theta,t)e^{-i\varpi_- t},$$

one obtains coupled modified Lugiato–Lefever equations for the main soliton, secondary signal, and idler. In the weak-idler limit, the idler is governed by an approximate equation whose key source term is $-i\gamma a_0^2 a_-^*$, identified as a purely four-wave mixing “Bragg-scattering” drive (Shandilya et al., 6 Feb 2026). A related three-wave formalism describes pump, signal, and idler bands $b_0,b_-,b_+$ in a dual-pumped microresonator with nearly quartic dispersion and yields comb lines of the form

$(M+1)$0

with a single angular group velocity $(M+1)$1 and different angular phase velocities $(M+1)$2 (Menyuk et al., 2024).

Single-pump interband multicolor states in compound resonators are described by coupled Lugiato–Lefever equations for primary soliton $(M+1)$3, secondary soliton $(M+1)$4, and idler $(M+1)$5, including self-phase modulation, cross-phase modulation, and four-wave mixing terms (Ji et al., 24 Jul 2025). Outside resonator settings, the mathematical structure shifts but the multicomponent logic persists. Two-color waveguide molecules are modeled by coupled NLS-type equations for $(M+1)$6 and $(M+1)$7 with self-phase modulation and cross-phase modulation, whereas quasiperiodic quadratic media yield averaged three-wave equations with effective cubic terms such as self-phase modulation, cross-phase modulation, third-harmonic generation, and four-wave mixing (Oreshnikov et al., 2022, Qi, 2018).

3. Binding mechanisms

The principal binding mechanisms are cross-phase modulation, four-wave mixing, group-velocity matching, and, in some systems, phase-matched radiative exchange. In multicolor dissipative Kerr solitons, strong cross-phase modulation locks all wavepacket components to the same repetition rate while each component retains its own phase-rotation velocity; four-wave-mixing Bragg scattering provides the inter-color coupling that links adjacent synthetic-frequency sites (Moille et al., 2022).

A particularly explicit mechanism is parametric binding. In the bright-bright integrated soliton molecule, the idler band may lie in normal dispersion, yet the term $(M+1)$8 replicates the DKS profile in the idler band. When $(M+1)$9 exceeds the maximum $N$0, there is no linear dispersive-wave resonance, so the idler grows only where $N$1 is large. The resulting bright idler is therefore “slave” to the soliton “master,” and the absence of phase matching prevents pedestal growth and forces localization (Shandilya et al., 6 Feb 2026). This directly contradicts the common assumption that a bright companion pulse must inherit the local anomalous-dispersion sign of its center wavelength.

In interband microcomb realizations, the intense primary soliton creates a $N$2 XPM potential in which the secondary field can be trapped. The secondary mode index shift $N$3 aligns the free-spectral range of the secondary pulse with the primary soliton, while the secondary-soliton exponent $N$4 follows from a balance of XPM and dispersion (Ji et al., 24 Jul 2025). In fiber and waveguide molecule models, the existence conditions are written as group-velocity matching, $N$5, together with propagation-constant locking between the two subpulses (Oreshnikov et al., 2022).

Quadratic media add another route. In quasiperiodic quasi-phase-matched superlattices, cascading of $N$6 processes induces effective cubic terms, and the resulting multicolor solitons are stabilized by the combined action of quasi-phase matching, diffraction–nonlinearity balance, and mutually trapped harmonic components (Qi, 2018). Across these platforms, multicolor binding is therefore not reducible to a single mechanism; it can be parametric, XPM-mediated, cascade-induced, or componentwise integrable, depending on the host system.

4. Experimental platforms and observed states

Integrated Kerr microresonators constitute the most developed experimental arena. Dual-pumped devices support multi-color dissipative Kerr solitons, synthetic dispersive waves, and smooth idler envelopes. In a 1 THz repetition rate Si$N$7N$N$8 resonator, four-wave mixing Bragg scattering between wavepackets produced a multi-color soliton spanning over 150 THz and an integrated synthetic frequency lattice, while the electrical spectrum at the drop port remained a single line at $N$9, confirming one repetition rate for all colors (Moille et al., 2022).

A later integrated experiment demonstrated a “Universal Bright-Bright Integrated Soliton Molecule via Parametric Binding,” where a bright idler pulse formed in normal dispersion and remained fundamentally bound to a bright DKS. The three colors locked to the same repetition rate but had fixed carrier-envelope offsets, and optical spectra versus secondary-pump detuning showed a transition from synthetic dispersive waves to a smooth $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$0 envelope in the idler band (Shandilya et al., 6 Feb 2026). Closely related dual-pumped studies used three-wave equations to explain interleaved frequency combs and the experimentally observed soliton-OPO effect, in which the pump frequency comb and signal frequency comb generate an idler frequency comb in a new frequency range (Menyuk et al., 2024).

Single-pump compound resonators provide another route. In a three-coupled silicon-nitride ring resonator with on-chip heaters, a primary soliton was observed first, and above a threshold detuning a secondary soliton appeared at a different frequency together with an idler sideband. The primary and secondary comb lines formed two interleaved sets separated by $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$1, and both soliton trains yielded a single high-SNR RF tone at $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$2 (Ji et al., 24 Jul 2025).

Multicolor pairwise mode locking has also been demonstrated in normal-dispersion coupled-ring resonators. There, pulses cannot exist alone, and instead must phase lock in pairs to form a bright soliton comb. The pulses can form at recurring spectral windows, the pulses in each pair feature different optical spectra, and three-ring systems support three pulses that mode lock through alternating pairwise pulse coupling (Yuan et al., 2023).

In mode-locked fiber lasers, polychromatic soliton molecules were observed by implementing the desired dispersion with a spectral pulse-shaper. The system supported two or more coincident solitons with different frequencies but common group velocities, and spectrogram retrieval showed temporally coincident soliton atoms with piecewise constant phase and $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$3-phase jumps at intensity nodes (Lourdesamy et al., 2020).

5. Stability, collisions, and radiative dynamics

Stability questions depend strongly on model class. In dissipative microresonators, boundedness arises from the interplay of loss, pump detuning, Kerr nonlinearity, and cross-phase modulation. For the bright-bright integrated molecule, small perturbations in pump detunings or powers remain bounded because loss, nonlinearity, and XPM provide restoring forces (Shandilya et al., 6 Feb 2026). In the multiscale microcomb theory, the usual single-color LLE soliton stability condition carries over to the slow envelope $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$4 in the optimal multi-pump case (Silvestri et al., 2024).

In quasiperiodic quadratic media, stability is developed more formally. The averaged three-wave equations possess the Manley–Rowe power integral, a Hamiltonian, and a Lyapunov functional, from which orbital stability is established. A virial identity rules out finite-distance collapse, and a multiscale analysis yields the Vakhitov–Kolokolov marginal-stability condition $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$5 (Qi, 2018).

Radiative dynamics are especially rich in two-color molecules. Cherenkov radiation obeys resonance conditions of the form $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$6, and additional four-wave-mixing-mediated radiation occurs at frequencies satisfying $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$7 (Oreshnikov et al., 2022). If the molecule oscillates, the resonance spectrum acquires sidebands indexed by the oscillation harmonic $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$8, producing a Kushi-comb-like multi-frequency radiation pattern (Melchert et al., 2022). These works also show that external dispersive waves can excite internal oscillations and that the scattering process can be used to probe the internal mode of a two-color soliton molecule.

In multicomponent integrable systems, the emphasis shifts from dissipation to collision structure and parameter constraints. Nondegenerate fundamental bright solitons of the $$\frac{\partial f}{\partial T} = \frac1{t_R}\Bigl[ -(\alpha + i\,\delta)\,f -\frac{i\,L\,\bar\beta_2}{2}\,\frac{\partial^2 f}{\partial X^2} +i\,\gamma\,L\,\varepsilon_J\,|f|^2f +\sigma_J\,\sqrt{\theta}\,E_0 \Bigr],$$9-coupled Manakov-type CNLS system exhibit up to $\varepsilon_J$0 humps when all $\varepsilon_J$1 are distinct, and numerical Crank–Nicolson simulations for the $\varepsilon_J$2-CNLS case showed that the triple-hump profile remains stable against perturbations of $\varepsilon_J$3 and $\varepsilon_J$4 white noise (Ramakrishnan et al., 2021). In mixed Yajima–Oikawa and Mel’nikov systems, inelastic collision can occur only among short-wave components that carry bright solitons in at least two components, whereas dark short-wave parts and the long-wave bright soliton undergo elastic collision with phase or position shifts (Chen et al., 2015, Han et al., 2017). A plausible implication is that “multicolor” does not imply a universal collision law; the outcome is encoded in the component architecture of the model.

6. Applications, extensions, and related directions

The application space is broad because multicolor solitons combine localization with large spectral separation. In integrated Kerr combs, the bright-bright parametric molecule enables on-chip spectral extension of Kerr combs into the visible without complex dispersion engineering or low-Q waveguide design, and the idler can address atomic transitions such as Rb and Cs while the DKS remains in telecom bands (Shandilya et al., 6 Feb 2026). Dual-pumped three-wave models further identify interleaved combs with shared repetition rate and distinct offset frequencies, a configuration relevant to low-noise microwave generation and potentially to a chip-scale clockwork (Menyuk et al., 2024).

Synthetic-dimension and topological applications form a second cluster. The integrated dispersive multi-color soliton realizes an all-optical synthetic frequency lattice compatible with octave-spanning microcombs and has been proposed as a route toward Harper–Hofstadter and SSH models, quantum emulation, and graph-based optical processing (Moille et al., 2022). In Floquet topological insulators made of helical $\varepsilon_J$5 waveguide arrays, multicolor edge solitons bifurcate from topological edge states and, away from resonance, are described by a period-averaged nonlinear Schrödinger equation with an effective cubic coefficient $\varepsilon_J$6 whose sign and magnitude depend on the overall phase mismatch (Ivanov et al., 2021).

A third cluster concerns communications, spectroscopy, and waveform generation. Multicolor interband soliton combs can generate THz-rate intensity modulation from interference between two femtosecond pulses with THz-scale spectral separation, with proposed conversion to a THz-band comb via photoconductive or optical-rectification techniques (Ji et al., 24 Jul 2025). Pairwise mode-locked normal-dispersion microresonators suggest all-optical soliton buffers and memories over multiple storage rings, as well as platforms for quantum combs and topological photonics (Yuan et al., 2023). In mode-locked lasers, the dispersion-shaped cavity supporting polychromatic soliton molecules offers optical analogies to mutual trapping and spectral tunneling in quantum mechanics (Lourdesamy et al., 2020).

Beyond resonators and lasers, the concept extends to supercontinuum physics and spatially modulated quadratic media. Simultaneous bright and dark solitons in photonic-crystal fibers can seed a two-octave-spanning continuum, where weak trapped radiation produced by bright–dark collisions defines the short- and long-wavelength limits of the continuum (Milián et al., 2016). In two-dimensional $\varepsilon_J$7 media with triangular singular modulation, stable multi-soliton sets, asymmetric states, and vortex rings provide a spatial counterpart to multicolor pattern formation (Lutsky et al., 2018).

Taken together, these developments indicate that multicolor solitons occupy a boundary region between solitary-wave theory, comb physics, synthetic dimensions, and topological photonics. This suggests that future work will continue to focus on universal reduction methods, phase- and group-matching design rules, and the controlled conversion between bound multicolor states, interleaved combs, and radiative sideband spectra.