# Collision-induced frequency shifts in bright matter-wave solitons and soliton molecules

## Abstract

A recent experiment has detected collision-induced frequency shifts in bright matter-wave solitons for the first time [Nat. Phys. 10, 918 (2014)]. Using a particle model, we derive the frequency shift for two solitons in a harmonic trap, and compare it to the results of the recent experiment and reported theoretical approximation. We find regimes where the frequency shift is much smaller than previously predicted, and propose experiments to test these findings. We also predict that reducing the experimental trap frequency will reveal for the first time soliton molecules or soliton bound states in a cold atoms system. The dynamics of such bound states are found to be both highly phase-dependent, and sensitive to the residual 3D nature of the experiment.

###### pacs:

03.75.Lm,Since the first realisations of bright solitons in Bose-Einstein condensates (BECs) (1); (2); (3), many uses for solitons have been suggested including as surface-probes (4); (5), for Bell-state generation (6) and for interferometry (7); (8); (9). The longstanding interest in bright solitons has recently been revived by an experiment on soliton collisions (10).

Solitons are solitary waves robust with respect to collisions; they survive these collisions unscathed up to a shift in their positions and phases. However, the vital characteristic of bright solitons – the collisional shifts – has not been observed in BEC until recently. The experiment of Rice University (10) of two solitons repeatedly colliding in a harmonic trap, demonstrated these shifts for the first time. The position shifts were revealed by their effect on the solitons’ frequency of oscillation in the trap. The authors of Ref. (10) compared their results with an approximate theoretical curve, which agreed reasonably well with the experiment in the tested regime.

In this Letter we provide an improved prediction of the frequency shift based on a particle model that reproduces the exact position shifts in the untrapped case (11); (12); (13). We identify regimes in which there is a measurable difference between our predictions and the theory of Ref. (10) and propose simple modifications to the experiment to explore these regimes. We also suggest modifications which will produce exotic soliton molecules (bound states) (7). Optical soliton molecules have been produced experimentally (14), but until now such states have not been produced in a cold atoms system. We predict that a reduction in the trap frequency of the experiment (10), along with modifications to control the relative phase, would produce soliton molecules in BEC for the first time. We investigate the 3-dimensional (3D) corrections to the dynamics, which affect both the stability and the mean-field dynamics of the BEC, particularly for the soliton molecules.

Bose-Einstein condensates with tight radial trapping are in most circumstances well-described by the 1-dimensional (1D) Gross-Pitaevskii Equation (GPE) (15):

(1) |

where , is the total atom number, and are the axial and radial trapping angular frequencies, and and are the mass and s-wave scattering length of the atomic species. In experiments, may be tuned to vary the sign and magnitude of the interactions, e.g., negative permits bright soliton solutions (16). We note that the 1D GPE breaks down in regimes where quantum or thermal fluctuations (17) and/or 3D effects become significant (18). Quantum and thermal effects are not expected to play a large role in the situations considered in this paper. However, the slightly 3D nature of the geometry imposes a limit on the interaction strength, beyond which the condensate will collapse (10); (19). This limit is characterised by a critical atom number . The parameter , where is the radial harmonic length of the trapping potential, is also used here and in Ref. (10) to characterise the strength of the interactions, since is proportional to , and is the number of particles per soliton. Adding the quintic nonlinearity to Eq. (1), where , yields a 1D equation which approximates the integrability-broken nature of the full 3D GPE (20). We use the GPE equation modified with this quintic term to verify the regimes of validity of our results based on Eq. (1).

In previous work (11); (12); (13), the following soliton-soliton interaction potential was employed to model solitons as individual particles:

(2) |

where is the dimensionless form of the relative position coordinate, , of the solitons, and are effective soliton mass parameters, which take the value for two solitons of equal size (12); (13). This potential produces the exact position shifts for solitons emerging from collisions within Eq. (1) when the axial potential . The interaction potential was shown numerically also to provide a good description when combined with an external trapping potential in regimes where the solitons were well-separated between collisions, particularly when the solitons collided in-phase (13).

Within the particle model, the oscillation frequency of the solitons’ relative position dynamics is given by

(3) |

Here, are the turning points of the dimensionless relative position coordinate, . For solitons initially at rest (as in the recent experiment (10) and all cases analysed in this paper), is equal to the initial soliton separation. The frequency-shift may then be written . We have found analytical solutions to Eq. (3), in cases where the solitons are initially close or far apart. In the limit , i.e., where the solitons are at all times extremely close, we find:

(4) |

For large initial separation, where and , we can assume that the position-shifts generated by Eq. (2) occur on a timescale much faster than the trap period, and the oscillation period is found by merely adding the time-shifts generated by Eq. (2) during two collisions to the trap period, yielding:

(5) |

where

(6) |

The behaviour of the frequency shift between these two limits can be explained qualitatively by considering the effective potential, illustrated in Fig. 1. The potential has a clear bi-modal character, with a narrower interaction mode within the harmonic potential. In the limit where , the solitons are extremely close with respect to the width of the inter-soliton potential (which becomes large when the solitons’ interaction strength decreases). The dynamics are completely described by small harmonic oscillations, which increase in frequency as the interaction strength increases. When the interactions are increased further, or the solitons’ initial separations are increased, the solitons may still be strongly bound within the interaction potential, but this potential is no longer effectively harmonic, and the frequency increase with interaction strength is less rapid. For even stronger interactions or wider separations, the inter-soliton potential becomes narrow compared with the initial soliton separation, and the solitons effectively ‘escape’ the inter-soliton potential each collision. In this regime, the frequency shift levels off, and starts to decrease as the width of the inter-soliton potential vanishes in comparison with .

We provide for comparison the theoretical approximation derived in Ref. (10) based on the interaction Hamiltonian of a linear superposition of two Gaussian wavepackets:

(7) |

where is the axial harmonic length of the trapping potential. As opposed to the behaviour of the dynamics in the particle model, within Eq. (7), the shift is linear in the interaction strength. The most marked difference between the prediction of the particle model and that of Eq. (7) is that the frequency shift in Eq. (7) diverges to infinity, rather than to zero, as the interaction strength goes to infinity. In the opposite limit the particle model goes to zero quadratically, rather than linearly. However, the important comparison is for experimentally accessible regions (), and we are most interested in regimes where the difference between the predictions is large enough to resolve in experiment. In Ref. (10), the initial soliton separation, m; number of atoms per soliton, ; and radial and axial trap angular frequencies, Hz and Hz. The scattering length was varied between runs such that took values between and . Note that in the regimes of positive , the wavepackets are not strictly solitons, and such regimes are not considered in this paper.

Figure 2 shows the frequency shift from numerical evaluation of Eq. (3), along with the curves for Eqs. (4) and (5) representing the limiting behaviour, and the theoretical approximation provided in Ref. (10). The relevant region of the experimental parameters lies within Fig. 2(a); while Figs. 2(b)-(d) show some other regimes of interest – with different initial soliton separations or trap frequencies than in the experiment. The main deviation between the particle model and Eq. (7) is obvious in the large negative limit, where the frequency in the particle model levels off and then starts to decrease. This discrepancy is particularly apparent in Figs. 2(b) and (d), where the axial trapping frequency is a factor of approximately 10 less than in the recent experiment. However, for the parameters of recent experiment [Fig. 2(a)] the difference between the curves is probably within experimental error; for large initial soliton separation and tight trapping [Fig. 2(c)], the agreement between the curves is extremely close.

Before issuing confident experimental predictions from the above particle model results, there are two reasons to verify the model’s performance. Firstly, the particle model (and also the theoretical approximation of Ref. (10)) assumes that the solitons separate between collisions, and approach at a sufficient speed such that the relative-phase between the solitons plays no effect. It is clear from the results of Ref. (13) that for solitons slow enough that the collision time approaches the order of the trap period, out-of-phase collisions will have smaller frequency-shifts than those predicted by the particle model, whereas in-phase collisions are still well described. In the regimes where the initial separation was m (10), the ansatz used to derive the particle model consists of initially slightly overlapping solitons, even for the narrowest solitons considered in the experiment. This obviously contradicts the assumption that the solitons always separate completely between collisions. Secondly, in some regimes, 3D effects to lead to condensate collapse during collisions, especially for and in-phase solitons (10); (18).

We first evaluate on the phase dependence of the position shifts by integrating the GPE without the quintic term. Fig. 3(a)-(b) shows GPE simulations two solitons for the experimental parameters with for in- and out-of-phase collisions. The in-phase collisions are well described by the particle model, and, as expected, the predicted frequency shift is too large for the out-of-phase collisions.

Interestingly, we find that when the axial trap is weakened to Hz [Fig. 3(c)-(d)], the in-phase solitons form a bound state (soliton molecule) which is surprisingly well-described by the particle model (until it eventually de-phases at longer timescales). The out-of-phase solitons (which seem to repel each other) are badly described by the particle model. Similar solitons with relative phases other than 0 and will have intermediate behaviour. We expect that the theoretical model employed in Ref. (7) would provide a better description of the GPE dynamics of these soliton molecules if it could be extended to incorporate the (weak) external trapping potential.

For larger initial soliton separations than used in the recent experiment, the particle model agrees very well with the GPE for any value of the relative phase. We illustrate this by simulating the system for an initial separation of four times that of recent experiment: m. For the experimental trap frequency, the frequency shift is too small to detect practically. We instead consider a weaker trap frequency of Hz, resulting in slower solitons with easily observable frequency shifts. The agreement between the model and the GPE is illustrated in Fig. 3(e)-(f). Note that it is this regime in which there is greatest improvement of the particle model over Eq. (7).

We now investigate 3D effects by integrating the 1D GPE with the quintic nonlinearity. We find that for interaction strength , the in-phase collisions are unstable, to collapse, but not the out-of-phase collisions [see Fig. 4(c)-(d)]. This suggests that the in-phase collisions of Fig. 3(a),(c) and (e) would not be realisable. However, slightly reducing the interaction strength to produced qualitatively similar states which do not collapse. We find that the quintic term increases the oscillation frequency by a factor of approximately 1.5 [see figure 4(a)-(b)]. The frequency shifts of the unbound states remain unaffected [see Fig. 4(d)], and it is here that the particle model is most useful.

In conclusion, we demonstrated the existence of regimes described better by the particle model than the theoretical approximation (10), i.e., those of weak axial trap frequency and stronger interactions. We propose the extension of the recent experiment (10) to verify the particle model predictions. We also suggest that soliton molecules can be created if greater control of the soliton phase is possible, e.g., by applying a light-sheet potential to half of the condensate (21); (22); (23). However, such states are highly dependent on relative-phase, and also on 3D effects. In future work we propose extending the formalism of (7) to include the effects of a harmonic trapping potential and 3D effects in order to describe soliton molecules.

We thank D. Schumayer and P. B. Blakie for useful discussions, and P. B. Blakie for assistance with preparation of Fig, 2(a).

### References

- K. E. Strecker, G. B. Partridge, A. G. Truscott, and R. G. Hulet, Nature 417, 150 (2002).
- L. Khaykovich et al., Science 296, 1290 (2002), http://www.sciencemag.org/content/296/5571/1290.full.pdf.
- S. L. Cornish, S. T. Thompson, and C. E. Wieman, Phys. Rev. Lett. 96, 170401 (2006).
- C. Lee and J. Brand, EPL (Europhysics Letters) 73, 321 (2006).
- S. Cornish et al., Physica D: Nonlinear Phenomena 238, 1299 (2009), Nonlinear Phenomena in Degenerate Quantum Gases.
- B. Gertjerenken et al., Phys. Rev. Lett. 111, 100406 (2013).
- U. A. Khawaja and H. T. C. Stoof, New Journal of Physics 13, 085003 (2011).
- A. D. Martin and J. Ruostekoski, New Journal of Physics 14, 043040 (2012).
- J. L. Helm, T. P. Billam, and S. A. Gardiner, Phys. Rev. A 85, 053621 (2012).
- J. H. V. Nguyen, P. Dyke, D. Luo, B. A. Malomed, and R. G. Hulet, Nat Phys 10, 918 (2014).
- R. Scharf and A. R. Bishop, Phys. Rev. A 46, R2973 (1992).
- A. D. Martin, C. S. Adams, and S. A. Gardiner, Phys. Rev. Lett. 98, 020402 (2007).
- A. D. Martin, C. S. Adams, and S. A. Gardiner, Phys. Rev. A 77, 013620 (2008).
- M. Stratmann, T. Pagel, and F. Mitschke, Phys. Rev. Lett. 95, 143902 (2005).
- F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999).
- T. Billam, A. Marchant, S. Cornish, S. Gardiner, and N. Parker, Bright solitary matter waves: Formation, stability and interactions, in Spontaneous Symmetry Breaking, Self-Trapping, and Josephson Oscillations, edited by B. A. Malomed, volume 1 of Progress in Optical Science and Photonics, pp. 403–455, Springer Berlin Heidelberg, 2013.
- N. P. Proukakis and B. Jackson, Journal of Physics B: Atomic, Molecular and Optical Physics 41, 203002 (2008).
- N. G. Parker, A. M. Martin, S. L. Cornish, and C. S. Adams, Journal of Physics B: Atomic, Molecular and Optical Physics 41, 045303 (2008).
- A. Gammal, T. Frederico, and L. Tomio, Phys. Rev. A 64, 055602 (2001).
- L. Khaykovich and B. A. Malomed, Phys. Rev. A 74, 023607 (2006).
- S. Burger et al., Phys. Rev. Lett. 83, 5198 (1999).
- J. Denschlag et al., Science 287, 97 (2000), http://www.sciencemag.org/content/287/5450/97.full.pdf.
- C. Becker et al., Nat Phys 4, 496 (2008).