Particle-Hole Asymmetry and Brightening of Solitons in A Strongly Repulsive BEC

We study solitary wave propagation in the condensate of a system of hard-core bosons with nearest-neighbor interactions. For this strongly repulsive system, the evolution equation for the condensate order parameter of the system, obtained using spin coherent state averages is different from the usual Gross-Pitaevskii equation (GPE). The system is found to support two kinds of solitons when there is a particle-hole imbalance: a dark soliton that dies out as the velocity approaches the sound velocity, and a new type of soliton which brightens and persists all the way up to the sound velocity, transforming into a periodic wave train at supersonic speed. Analogous to the GPE soliton, the energy-momentum dispersion for both solitons is characterized by Lieb II modes.

💡 Research Summary

In this paper the authors investigate solitary‑wave propagation in a strongly interacting Bose‑Einstein condensate (BEC) described by a hard‑core boson (HCB) lattice model with nearest‑neighbour (nn) interactions. By mapping the HCB Hamiltonian onto a spin‑½ XXZ model and employing spin‑coherent‑state averages, they derive coupled continuum equations for the condensate order parameter η(r) and the particle‑hole imbalance δ(r)=1−2ρ(r), where ρ is the particle density. Unlike the Gross‑Pitaevskii equation (GPE), the condensate density ρ_c is not equal to the particle density; instead ρ_c=ρ(1−ρ), reflecting an intrinsic particle‑hole duality.

Assuming a uniform background with constant imbalance δ₀, they look for traveling‑wave solutions of the form ρ(z)=ρ₀+f(z) with z=(x−vt)/a. This leads to a nonlinear differential equation (8) for f(z), which can be expressed in terms of an elliptic integral. In a useful approximation the equation admits two analytic soliton branches f₊(z) and f₋(z) (Eq. 9). The f₋ branch corresponds to the familiar dark soliton: its amplitude vanishes as the soliton speed v approaches the Bogoliubov sound speed v_s, reproducing the GPE behaviour. The f₊ branch is a new “brightening” or “persistent” soliton that exists only when the particle‑hole imbalance is non‑zero (ρ₀≠½). Remarkably, its amplitude remains finite even at v=v_s; the soliton does not flatten out but instead becomes broader. For supersonic velocities (v>v_s) the parameter γ²=1−(v/v_s)² becomes negative and the f₊ solution turns into a spatially periodic wave train, indicating a smooth transition from a localized soliton to a nonlinear wave lattice.

The authors compute the canonical momentum P and energy E of both branches by integrating the continuity equation and the Hamiltonian density. The resulting dispersion curves (Fig. 3) display a linear segment near the sound velocity and a saturation at low velocities, exactly the hallmark of Lieb’s type‑II excitations. The lower branch (f₋) behaves like a “hole” excitation with negative effective mass, while the upper branch (f₊) behaves like a “particle” excitation with positive effective mass. Thus the particle‑hole imbalance controls not only the existence of the two soliton families but also their dynamical character.

Numerical integration of Eq. (8) confirms the analytic forms across a wide range of background densities ρ₀ and velocities, with the best agreement away from half‑filling where the particle‑hole symmetry is broken. The nn interaction V merely rescales the soliton width without altering its profile, suggesting that long‑range interactions (e.g., dipolar forces) would not destroy the solitons.

Experimentally, the predictions could be tested in quasi‑one‑dimensional, highly anisotropic traps where hard‑core behaviour is approached (e.g., using strong on‑site repulsion via a Feshbach resonance and a deep optical lattice to enforce the no‑double‑occupancy constraint). By preparing an initial density perturbation with a controlled imbalance (through atom number or external potentials) one could generate either a dark soliton or the persistent brightening soliton and monitor its propagation with phase‑contrast imaging. The transition to a periodic wave train at supersonic speeds would be observable as a train of density ripples moving faster than the sound speed.

In summary, the work extends soliton physics beyond the conventional GPE framework, revealing that hard‑core bosons support a pair of soliton solutions whose properties are dictated by particle‑hole asymmetry. The persistent soliton, which survives up to the sound velocity and morphs into a wave train beyond, provides a concrete example of a nonlinear excitation that directly maps onto Lieb’s type‑II branch. The study opens avenues for exploring non‑GPE dynamics in strongly correlated quantum gases, both theoretically (e.g., via time‑dependent DMRG or quantum Monte‑Carlo) and experimentally, and highlights the rich interplay between many‑body correlations, particle‑hole duality, and nonlinear wave phenomena.