Dark–Bright Solitons in Coupled Condensates

• Dark–bright solitons are composite nonlinear waves where a localized dark density dip traps a bright pulse via intercomponent repulsion.
• They are modeled by coupled Gross–Pitaevskii and nonlinear Schrödinger equations, enabling analytical, variational, and numerical studies in homogeneous, trapped, and dissipative settings.
• Research reveals rich bifurcation structures, effective potential dynamics, and extensions to optical and Bose–Fermi superfluid systems, highlighting both fundamental insights and practical applications.

Dark-bright solitons are composite nonlinear wave states in coupled defocusing or repulsively interacting media, most prominently in two-component Bose–Einstein condensates, where one component carries a dark soliton—a localized density dip on a nonzero background, typically with a phase jump—while the other carries a bright soliton, a localized hump trapped inside that dip. In this sense they are symbiotic structures: the bright part is supported by the effective potential well created by the dark component, even though the underlying interactions are repulsive. Within the literature surveyed here, dark-bright solitons appear in homogeneous and trapped Gross–Pitaevskii systems, coupled nonlinear Schrödinger equations, spin-orbit-coupled condensates, polariton condensates, photorefractive lattices, microresonators, Bose–Fermi superfluids, and dipolar settings, with both exact and approximate constructions, rich bifurcation structure, and a wide range of dynamical regimes (Yan et al., 2014).

1. Definition, physical mechanism, and canonical structure

In the standard quasi-one-dimensional two-component Gross–Pitaevskii or defocusing vector nonlinear Schrödinger description, a dark-bright soliton consists of a density notch in one field and a localized pulse in the other. The physical binding mechanism is intercomponent repulsion: the density depletion of the dark component acts as an effective trapping well for the bright component, so the bright pulse can persist even though a bright soliton would not exist by itself in the same repulsive medium (Achilleos et al., 2011).

A representative mean-field form in the integrable Manakov case is

$u_d(x,t)=\cos\phi\,\tanh[D(x-x_0(t))]+i\sin\phi,$

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$

with $\phi$ controlling depth and speed, $D$ the inverse width, $\eta$ the bright amplitude, and $x_0(t)$ the center (Katsimiga et al., 2016). Closely related parameter relations recur across the literature, including

$D^2=\cos^2\phi-\eta^2,\qquad \dot x_0=D\tan\phi,$

and a bright-component norm or particle number proportional to $\eta^2/D$ (Achilleos et al., 2011).

The same structural idea extends beyond the equal-coupling limit. In a homogeneous two-component condensate with general interaction coefficients, a stationary dark-bright ansatz of the form

$$u_1=A_1\tanh(bx),\qquad u_2=A_2\,\operatorname{sech}(bx)$$

yields algebraic solvability conditions linking amplitudes and width to $g_{11},g_{12},g_{22}$ and the chemical potentials. In that setting the bright component exists only when

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$0

showing that exact closed-form dark-bright states are not confined to the Manakov point (Yan et al., 2014).

A common misconception is that dark-bright solitons necessarily require attractive self-interaction in the bright component. The surveyed work does not support that view. Rather, in repulsive two-component media the bright structure is typically a trapped filling mode supported by the dark-component depletion, which is why these states are frequently described as symbiotic (Achilleos et al., 2011).

2. Mean-field models and asymptotic reductions

The principal continuum description is a pair of coupled Gross–Pitaevskii equations. One standard dimensionless form is

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$1

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$2

possibly supplemented by impurities, damping, reservoir terms, or unequal nonlinear coefficients depending on the application (Achilleos et al., 2011). In homogeneous settings, the trap is absent; in trapped settings, $u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$3 is the usual harmonic confinement (Yan et al., 2014).

Several reductions clarify the internal mechanics of the state. In treatments of dark-bright soliton pairs, the dark component can be viewed as generating an effective double-well potential for the bright one,

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$4

so that the bright field is described in a two-mode basis of the lowest symmetric and antisymmetric eigenstates of this induced well (Karamatskos et al., 2014). This connects dark-bright soliton pairs to Bosonic Josephson Junction phenomenology, including plasma oscillations, $u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$5-oscillations, and self-trapping in the bright population imbalance and phase-difference variables (Karamatskos et al., 2014).

In spin-orbit-coupled condensates with vanishing Raman coupling, a multiscale expansion reduces the system to the Mel'nikov model, from which approximate traveling dark-bright solitons can be reconstructed. In that reduction one imposes boundary conditions

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$6

and obtains a KdV-type equation coupled to a stationary Schrödinger equation for the bright component (D'Ambroise et al., 2017). This formalism demonstrates that the existence of traveling dark-bright structures can persist even when the original model is not itself integrable.

A more recent variational treatment emphasizes that equal spatial scales are not fundamental. A family of different-width dark-bright solitons is built from the ansatz

$u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$7

where $u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$8 and $u_b(x,t)=\eta\,\sech[D(x-x_0(t))]\,e^{i[kx+\theta(t)+(\mu-1)t]},$9 are allowed to differ. The associated effective wells

$\phi$0

need not coincide, in sharp contrast to the equal-width case (Mao et al., 2024).

3. Existence theory: integrable, nonintegrable, and variational regimes

The Manakov limit,

$\phi$1

occupies a special place because it is integrable and supports exact single- and multi-dark-bright soliton solutions (Katsimiga et al., 2018). Much of the subsequent literature studies what survives when this equality is broken, either by varying interspecies or intraspecies interactions, by introducing unequal dispersion coefficients, or by adding trapping and dissipation.

For single solitons in homogeneous condensates with general interaction coefficients, exact analytical solutions exist only in restricted parameter windows, but numerically the existence region is far wider (Yan et al., 2014). In that work the analytical single-dark-bright solution is corroborated numerically, while lattice generalizations are also found well beyond the analytically tractable regime. A related message appears in the different-width theory: the allowed nonlinear-parameter region is much broader than for previously known equal-width exact solutions, and the equal-width family is recovered only in the special limit $\phi$2 (Mao et al., 2024).

Rigorous existence theory has recently been developed for a defocusing one-dimensional Gross–Pitaevskii system with asymptotic conditions

$\phi$3

Using constrained minimization of an energy functional at fixed bright mass and modified dark momentum, one proves existence of traveling-wave dark-bright solitons with $\phi$4, $\phi$5, even and radially nonincreasing modulus defects, and subsonic speed

$\phi$6

under the hypotheses denoted $\phi$7 or $\phi$8 in the theorem (López-Martínez, 26 Aug 2025). This result places dark-bright solitons in a fully nonperturbative variational framework beyond the traditional integrable setting.

Another route to nonintegrable existence arises from unequal dispersion. In a two-component defocusing nonlinear Schrödinger system with $\phi$9 and $D$0, the dark background

$D$1

induces a Pöschl–Teller well for the bright component. The bound states are indexed by $D$2, and the $D$3-th one exists for

$D$4

The $D$5 branch is the usual dark-bright soliton, while higher $D$6 yield excited dark-in-bright bound states (Charalampidis et al., 2014).

These results collectively suggest that the integrable Manakov soliton is only one member of a much broader class. That inference is explicit in the numerical and variational literature, where dark-bright states persist beyond exact solvability, unequal width can arise self-consistently, and constrained variational arguments now establish existence in a broad defocusing regime (Mao et al., 2024).

4. Bound states, lattices, and bifurcation structure

Dark-bright solitons support higher-order organizations, including periodic lattices and bound pairs. In homogeneous condensates, periodic arrays can be written in Jacobi elliptic form. Two canonical examples are the sn-cn lattice,

$D$7

with out-of-phase bright neighbors, and the sn-dn lattice,

$D$8

with in-phase bright neighbors. Their inter-soliton spacing is

$D$9

and both reduce to the single-dark-bright soliton as $\eta$0 (Yan et al., 2014). Numerically, these lattice states are generally unstable, but the instability growth rates are smallest near the miscibility-immiscibility threshold (Yan et al., 2014).

For two-soliton bound states, relative phase in the bright component is decisive. In a homogeneous repulsive condensate beyond the integrable limit, a variational interaction theory decomposes the total energy as

$\eta$1

where the cross-component term $\eta$2 is of the same asymptotic order as $\eta$3 and cannot be neglected. Out-of-phase bright components, $\eta$4, support a genuine equilibrium molecule, whereas in-phase equilibria predicted by the ansatz are not supported by direct numerical fixed-point computations or dynamics (Katsimiga et al., 2016).

The stability of these pair states is organized by symmetry-breaking bifurcations. In the homogeneous beyond-integrable problem, increasing $\eta$5 drives an exponential instability of the out-of-phase bound state near

$\eta$6

associated with a subcritical pitchfork bifurcation and an asymmetric redistribution of bright mass that can end in a dark / dark-bright configuration (Katsimiga et al., 2016). In a trapped system, by contrast, the out-of-phase branch of a dark-bright pair undergoes a symmetry-breaking pitchfork bifurcation at

$\eta$7

for $\eta$8, and two new stable asymmetric branches emerge (Karamatskos et al., 2014). In that setting the effective double-well picture captures bright-atom tunneling between the two dark-soliton wells and recovers the standard Bosonic Josephson Junction phase portrait in $\eta$9 (Karamatskos et al., 2014).

A different bifurcation scenario appears when one varies a single intra-species coupling through the Manakov limit. Antisymmetric and asymmetric dark-bright pairs are degenerate at $x_0(t)$0, but away from that point they separate in a transcritical bifurcation with symmetry. The antisymmetric branch is stable for $x_0(t)$1 and unstable for $x_0(t)$2, while the asymmetric branch exhibits the opposite exchange of stability (Katsimiga et al., 2018).

A frequent misconception is that out-of-phase bright components are always the stable bound configuration. The literature does not sustain so general a claim. Stability depends on setting and control parameter: in homogeneous nonintegrable systems the out-of-phase pair is the physically relevant molecule but can destabilize through a subcritical pitchfork, whereas in trapped systems the out-of-phase branch can lose stability and generate stable asymmetric states through a supercritical pitchfork-like scenario (Katsimiga et al., 2016).

5. Dynamics in traps, at interfaces, and under external perturbations

In a harmonic trap, a single dark-bright soliton behaves as a nonlinear quasiparticle with an anomalous mode governing its oscillation. Perturbation theory for a trapped condensate yields an effective law

$x_0(t)$3

which reduces to an oscillatory motion for $x_0(t)$4 (Yan et al., 2014). Related formulations write the center dynamics as

$x_0(t)$5

or, at finite temperature,

$x_0(t)$6

the latter describing an anti-damped oscillator in dissipative Gross–Pitaevskii theory (Achilleos et al., 2011). In that framework the bright filling reduces the effective anti-damping and partially stabilizes the dark soliton against thermal dissipation, thereby extending the lifetime relative to a bare dark soliton (Achilleos et al., 2011).

External inhomogeneities can strongly reshape transport. In the presence of delta-like impurities acting only on the bright component, an attractive impurity produces an effective barrier for the composite soliton, while a repulsive impurity produces an effective well. This sign reversal originates from the repulsive back-action of the bright atoms on the dark component and is confirmed by perturbation theory, Bogoliubov–de Gennes analysis, and direct simulation (Achilleos et al., 2011).

Domain walls between immiscible condensates provide another nonlinear scattering environment. In a one-dimensional immiscible mixture satisfying

$x_0(t)$7

a dark soliton incident on the interface can be transmitted, reflected, or converted into a dark-bright soliton by capturing atoms from the opposite component in the overlap region of the domain wall (Arazo et al., 2021). The resulting dynamically generated dark-bright soliton appears even outside the regime where stationary analytical dark-bright solutions are known, and its subsequent trap motion is harmonic-like. A semi-analytical frequency formula,

$x_0(t)$8

agrees well with numerical data for approximately

$x_0(t)$9

(Arazo et al., 2021).

Collisions reveal a sharp distinction between topological and non-topological sectors. In numerical Hong–Ou–Mandel-type experiments with repulsive condensates, pure dark solitons interacting with a barrier show no significant asymmetries in the outgoing waves, whereas dark-bright solitons can exhibit pronounced asymmetries in the bright component, including cases where one outgoing bright packet carries almost all the bright mass (Sun et al., 2016). This difference is attributed to the topological nature of the dark notch and the non-topological character of the bright filling.

6. Generalizations, dissipative realizations, and beyond-mean-field dynamics

Dark-bright solitons persist well beyond conservative binary Gross–Pitaevskii theory. In spinor polariton condensates under nonresonant pumping, the open-dissipative dynamics is described by coupled condensate equations and a reservoir rate equation. Hamiltonian perturbation theory in the fast-reservoir limit yields a Newton-like law

$D^2=\cos^2\phi-\eta^2,\qquad \dot x_0=D\tan\phi,$0

and a lifetime scale

$D^2=\cos^2\phi-\eta^2,\qquad \dot x_0=D\tan\phi,$1

In that setting the dark-bright soliton relaxes by blending with the background at finite time, but remains long-lived enough for experimental observation, even in the presence of Langevin noise (Xu et al., 2018).

In passive optical microresonators, mutually trapped dark-bright cavity solitons arise under bichromatic pumping of two mode families with opposite dispersion signs and nearly matching free spectral ranges. The governing model consists of two coupled generalized Lugiato–Lefever equations with cross-phase modulation, group-velocity mismatch, and higher-order dispersion. One field forms a bright dissipative Kerr soliton in the anomalous-dispersion band, while the other forms a dark pulse in the normal-dispersion band through Kerr-induced cross-phase modulation with the bright pulse (Zhang et al., 2021). This establishes a dissipative, cavity-based realization distinct from conservative atomic systems.

The concept also extends to mixed-statistics superfluids. In a Bose–Fermi superfluid mixture at unitarity, a fermionic dark soliton can be filled by a broad bosonic component in the partially separated phase, yielding a dark-bright soliton that preserves full spatial coherence. In the fully separated phase, by contrast, the bosons completely fill the central region and the left and right fermionic domains become independent, so the coherent dark-bright soliton no longer persists in the same sense (Tylutki et al., 2016).

Higher-dimensional analogues reveal stabilization by filling. In quasi-two-dimensional two-component condensates, a dark-bright ring soliton consists of a circular dark depletion in one component and a bright filling in the other. The effective radial dynamics admits an equilibrium radius and breathing oscillations, while the bright population $D^2=\cos^2\phi-\eta^2,\qquad \dot x_0=D\tan\phi,$2 shifts symmetry-breaking bifurcations to larger $D^2=\cos^2\phi-\eta^2,\qquad \dot x_0=D\tan\phi,$3, thereby extending the stability window relative to the one-component ring dark soliton (Stockhofe et al., 2011). The same paper emphasizes that stabilization by a filling component is not limited to radially symmetric structures and also occurs in cross-like dark-bright configurations.

Beyond mean field, the composite picture changes qualitatively. Multi-orbital many-body simulations using ML-MCTDHB show that a mean-field dark-bright soliton in a trap can fragment and split into a fast and a slow daughter solitary wave, while collisions can produce multiple entangled dark-bright fragments, dark-antidark states, and domain-wall-bright complexes (Katsimiga et al., 2016). This directly contradicts the expectation that the parent dark-bright soliton remains an intact elementary object under all weakly perturbed evolutions.

Taken together, these developments show that dark-bright solitons are not a narrow integrable curiosity. They form a broad class of multicomponent coherent structures whose existence, stability, and dynamics depend sensitively on coupling asymmetry, trapping, dissipation, dimensionality, and many-body correlations, while retaining a unifying core mechanism: localization of a bright component inside the self-induced defect of a dark one (López-Martínez, 26 Aug 2025).