## Ultra-high-amplitude Peregrine solitons induced by helicoidal spin-orbit coupling

Cui-Cui Ding, 1 Qin Zhou, 1, 2, ∗ and B. A. Malomed 3, 4

1 Research Group of Nonlinear Optical Science and Technology, Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan 430200, China 2 State Key Laboratory of New Textile Materials and Advanced Processing Technologies, Wuhan Textile University, Wuhan 430200, China 3 Department of Physical Electronics, School of Electrical Engineering, Faculty of Engineering, and the Center for Light-Matter University, Tel Aviv University, Tel Aviv, Israel 4 Instituto de Alta Investigaci´n, o Universidad de Tarapac´, a Casilla 7D, Arica, Chile

In the framework of the model of a spatially non-uniform Bose-Einstein condensate with helicoidal spin-orbit (SO) coupling, we find abnormal Peregrine solitons (PSs) on top of flat and periodic backgrounds, with ultra-high amplitudes. We explore the roles of the SO coupling strength and helicity pitch in the creation of these anomalously tall PSs and find that their amplitude, normalized to the background height, attains indefinitely large values. The investigation of the modulation instability (MI) in the same system demonstrates that these PSs exist in a range of relatively weak MI, maintaining the feasibility of their experimental observation.

Introduction. Rogue waves (RWs), first discovered as extreme events in the ocean [1-6], have been widely studied, due to their unique properties and potential applications, in nonlinear optics [7-10], plasmas [11], BoseEinstein condensates (BECs) [12-14], magnetics [15], financial markets [16], and various other settings [17-26]. A widely recognized RW prototype is provided by the exact Peregrine-soliton (PS) solution of the nonlinear Schr¨dinger o equation (NLSE) [27], whose characteristic features are the threefold peak amplitude and spatiotemporal localization on top of the background field [28]. Several landmark experiments have directly demonstrated this remarkable phenomenon and its ramifications [14, 28-30].

ground height [24, 50, 52, 53]. Nevertheless, recent studies have shown that, under the action of self-steepening effects, the amplitude of fundamental PSs can exceed the threefold limit, reaching up to fivefold the background height [54]. In particular, exceptional PSs, which feature ultra-high peak amplitudes, have been reported too in the vector derivative NLSEs including the self-steepening effect [55].

Spin-orbit (SO) coupling in BECs have drawn much interest since its experimental implementation [31-33], as it offers the realization of the SO-coupling phenomenology in the uniquely clean form [34, 35] and make it possible to create artificial vector gauge potentials [36, 37]. Recently, models of BECs with non-uniform SO coupling have been introduced, as they provide high tunability of this effect, and enhance the role of the intrinsic nonlinearity in the SO-coupled BECs [38-43]. In this context, soliton dynamics in the BEC with non-uniform landscapes of the SO coupling has been investigated [44-47], where, in particular, the helicoidal gauge potential may originate from the light propagation in a helical waveguide array [48]. The propagation of matter-wave solitons in a BEC with a random SO coupling was addressed too [49].

SO-coupled BECs are modeled by systems of two (or several) coupled Gross-Pitaevskii equations (GPEs). In this connection, it is relevant to stress that PSs exist in multi-component NLSE models, such as the famous Manakov system, but, due to the energy transfer between different components, the PS amplitude is no longer fixed, although it still does not exceed the triple back-

∗ qinzhou@whu.edu.cn

In this work, we focus on the following questions: can the fundamental PS with an ultra-high peak amplitude be excited in other ways, besides using higher-order effects, such as self-steepening, and to what extent is it possible to increase the PS amplitude? To answer these questions, we first consider a BEC model with non-uniform helicoidal SO coupling (cf. Ref. [44]), which offers experimental feasibility. We construct its exact fundamental PS solutions on top of flat, alias continuous-wave (CW), and periodic backgrounds. Through the analysis of the PS amplitude, we find that PS with ultra-high peak amplitude, reaching indefinitely large values (as normalized to the background height), can be created with the help of the helicoidal SO coupling.

To explore the PS dynamics under the action of spatially non-uniform gauge potentials, we consider the GPE for the spinor wave function Ψ = (Ψ 1 , Ψ ) 2 T of an effectively one-dimensional two-component BEC, including the helicoidal SO coupling. In the scaled form (with M = /planckover2pi1 = 1, where M is the atomic mass), the GPE is [44, 56, 57]

$$i \frac { \partial \Psi } { \partial t } = \frac { 1 } { 2 } Q ^ { 2 } ( x ) \Psi - ( \Psi ^ { \dagger } \Psi ) \Psi, \quad \quad ( 1 )$$

where the helicoidally molded SO coupling is represented by the generalized momentum operator,

$$Q ( x ) = - i \partial _ { x } + \alpha \sigma \cdot n ( x ). \text{ \quad \ \ } ( 2 )$$

Here α is the SO-coupling strength, which is tunable in the experiment [39-41], σ = ( σ , σ x y , σ z ) is the vector of

the Pauli matrices, and the spatial modulation is represented by vector

$$n ( x ) = ( \cos ( 2 \kappa x ), \sin ( 2 \kappa x ), 0 ), \quad \ \ ( 3 ) \quad ^ { \mathcal { N } }.$$

with κ &lt; 0 and κ &gt; 0 corresponding to the left- and righthanded helicity, respectively [48, 58, 59]. As usual, it is assumed that the inter- and intra-species attractive interactions have equal strengths. Special forms of Eq. (1) include the uniform Rashba-Dresselhaus SO coupling [36] when κ = 0, and the canonical Manakov system [50] when α = 0.

Fundamental PS solutions. Eq. (1) is made gaugeequivalent to the integrable Manakov system,

$$i \mathbf u _ { t } + \frac { 1 } { 2 } \mathbf u _ { x x } + ( \mathbf u ^ { \dagger } \mathbf u ) \mathbf u = 0, \ \mathbf u = ( \mathbf u _ { 1 }, \mathbf u _ { 2 } ) ^ { T }, \ \ ( 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

by means of the transformation [49]

$$\Psi & = \begin{pmatrix} \nu _ { + } e ^ { - i ( k _ { m } + \kappa ) x } & \nu _ { - } e ^ { i ( k _ { m } - \kappa ) x } \\ \nu _ { - } e ^ { - i ( k _ { m } - \kappa ) x } & - \nu _ { + } e ^ { i ( k _ { m } + \kappa ) x } \end{pmatrix} \mathbf u, \quad ( 5 ) & k _ { 2 } ) t \\. \quad. \quad \overbrace { \longrightarrow } \. \quad. \quad. \quad. \quad. \quad. \quad. \end{pmatrix}$$

where k m = √ α 2 + κ 2 is the effective momentum of the lowest-energy states, and

$$\nu _ { + } & = \ s g n ( \alpha ) \sqrt { ( k _ { m } - \kappa ) \, / \, ( 2 k _ { m } ) } \quad \quad ( 6 ) \quad \text{Ir} \\ \nu _ { - } & = \sqrt { ( k _ { m } + \kappa ) \, / \, ( 2 k _ { m } ) } \quad \quad ( 7 ) \quad \text{stra} \\ & \text{com}$$

Below, k m plays a crucial role determining properties of PSs, especially as concerns the amplification of their amplitudes.

$$\nu _ { - } = \quad \sqrt { ( k _ { m } + \kappa ) \, / \, ( 2 k _ { m } ) } \quad ( 7 ) \quad \text{stra} \\ \text{w, } k _ { m } \text{ plays a crucial role determining properties of } \quad \text{PS} \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdots$$

The Manakov system (4) possesses the Lax pair [51] and admits the solution by means of the Darboux dressing method [24]. To begin with, we take the CW seed solution of Manakov system (4), with components

$$u _ { j 0 } = a \exp [ - i ( k _ { j } x - \omega _ { j } t ) ], \ j = 1, 2, \quad ( 8 )$$

which is determined by the amplitude ( a ), wavenumbers ( k j ), and frequencies

$$\omega _ { j } = 2 a ^ { 2 } - k _ { j } ^ { 2 } / 2. \quad \quad \quad ( 9 )$$

Making use of the Manakov system invariance with respect to the rotation of the set of the two components, we choose them in Eq. (8) with equal amplitudes a . Subsequent results demonstrate that the helicoidal SO coupling makes PS heights different in the two components Ψ 1 2 , for the same background amplitudes a , see Eqs. (13) and (15) below.

Utilizing the known PS solutions for Manakov system (4) derived by means of the Darboux transform [24], and substitution (5), we obtain the following exact fundamental PS solutions of the underlying Eq. (1):

$$\begin{array} { l l l } \text{$\text{$\text{$coverage to the main } \mathcal{m}_{1}$} \mathcal{m}_{2}$} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \mathcal{m}_{1} \mathcal{m}_{2} \\ \Psi _ { 1 } & = & a e ^ { - i \kappa x } \left [ \nu _ { + } \left ( 1 - \frac { \mathcal { R } _ { 1 } } { \mathcal { N } _ { 1 } } \right ) e ^ { i \theta _ { 1 } } + \nu _ { - } \left ( 1 - \frac { \mathcal { R } _ { 2 } } { \mathcal { N } _ { 2 } } \right ) e ^ { i \theta _ { 2 } } \right ], & \text{$\text{$the hel}$} \\ \theta _ { 1 } & = & - ( k _ { m } + k _ { 1 } ) x + \omega _ { 1 } t, & \theta _ { 2 } = ( k _ { m } - k _ { 2 } ) x + \omega _ { 2 } t, & ( 1 0 ) & \text{ we first} \end{array}$$

where we define

$$\mu \nu \quad & \text{where we use} \\ ( 3 ) & \quad \mathcal { N } _ { j } \ = \ \left [ ( \theta + \mu t ) ^ { 2 } + \zeta ^ { 2 } t ^ { 2 } + \frac { 4 } { \zeta ^ { 2 } } \right ] \{ [ \delta + ( - 1 ) ^ { j } \mu ] ^ { 2 } + \zeta ^ { 2 } \}, \\ \tilde { \ g h t } \quad & \mathcal { R } _ { j } \ = \ 8 i \{ \zeta ^ { 2 } t - [ \mu + ( - 1 ) ^ { j } \delta ] ( \theta + \mu t ) \} + 1 6, \\ \intertext { i n - 1 } \ \mu \ = \ \pm \frac { \sqrt { 2 } } { 2 } \left [ \sqrt { \delta ^ { 2 } ( 8 a ^ { 2 } + \delta ^ { 2 } ) } - 4 a ^ { 2 } + \delta ^ { 2 } \right ] ^ { 1 / 2 }, \quad ( 1 1 ) \\ [ 3 6 ] \quad & \text{in the case of } | \delta | \geq a, \, \text{with } \delta \equiv k _ { 1 } - k _ { 2 }, \, \text{or}$$

in the case of | δ | ≥ a , with δ ≡ k 1 -k 2 , or

$$\mathcal { N } _ { j } \ & = \ \left [ \theta ^ { 2 } + ( \zeta + \mu ^ { \prime } ) ^ { 2 } t ^ { 2 } + \frac { 4 } { ( \zeta + \mu ^ { \prime } ) ^ { 2 } } \right ] ( 2 a ^ { 2 } + \zeta \mu ^ { \prime } ), \\ \mathcal { R } _ { 1 } \ & = \ 4 i ( 4 a ^ { 2 } - \delta ^ { 2 } + 2 \zeta \mu ^ { \prime } ) t - 4 i ( - 1 ) ^ { j } \delta \theta + 8, \\ 4 ) \quad \mu ^ { \prime } \ & = \ \pm \frac { 1 } { \sqrt { 2 } } \left [ 4 a ^ { 2 } - \delta ^ { 2 } - \sqrt { \delta ^ { 2 } ( 8 a ^ { 2 } + \delta ^ { 2 } ) } \right ] ^ { 1 / 2 }, \quad ( 1 2 ) \\ & \quad \text{in the case of } | \delta | < a. \ In \text{ either case, we set } \theta \equiv 2 x + ( k _ { 1 } + \cdot \end{matrix}$$

in the case of | δ | &lt; a . In either case, we set θ ≡ 2 x +( k 1 + k 2 ) t and ζ ≡ ( 1 / √ 2 ) [ √ δ 2 (8 a 2 + δ 2 ) + 4 a 2 -δ 2 ] 1 / 2 . Using the translational symmetry, we shift the above solutions to the origin, to produce compact expressions for them. Note that these PS solutions are non-singular ones in the entire parameter range.

In addition to the same features which are demonstrated by the conventional PSs, that exist in some multicomponent systems, such as PSs of the bright-dark type, PS doublets, etc ., the helicoidal SO coupling can generate more intricate PS structures, among which the most salient aspect is, as shown below, the possibility of having PSs with uniquely large heights.

The consideration of the exact solution (10) reveals that the PS is generally located on top of a periodic background formed by the superposition of two different CWs. The exact solution for the periodic background is

$$\ e r s & \quad | \Psi _ { 1 } ^ { \text{bg} } | \, = \, a \sqrt { 1 + \frac { \alpha } { k _ { \text{m} } } \cos \left [ ( \delta + 2 k _ { \text{m} } ) x + \frac { k _ { 1 } ^ { 2 } - k _ { 2 } ^ { 2 } } { 2 } t \right ] }, \\ & \quad | \Psi _ { 2 } ^ { \text{bg} } | \, = \, a \sqrt { 1 - \frac { \alpha } { k _ { \text{m} } } \cos \left [ ( \delta + 2 k _ { \text{m} } ) x + \frac { k _ { 1 } ^ { 2 } - k _ { 2 } ^ { 2 } } { 2 } t \right ] } ( 1 3 ) \\ \intertext { w e - \, \begin{array} { c } \iota \, \colon \, \dots \, \dots \, \end{array} \, \frac { \iota \, 2 \, \sum \, \iota \, \iota \, \sigma \, \iota \, \end{array} }$$

It is moving with speed v = ( k 2 2 -k 2 1 ) / [2( δ +2 k m )], where k 1 2 , are the same wavenumbers as in Eq. (8).

Note that, if wavenumbers k 1 2 , and the momentum minimum k m satisfy the following relationship,

$$k _ { 1 } = - k _ { 2 } = - k _ { m }, \text{ \quad \ \ } ( 1 4 )$$

the cos terms vanish in Eq. (13), i.e., the periodic background degenerates into a flat CW. Due to the presence of the helicoidal SO coupling, the constraint (14) is different from similar ones which provide for the flat background in the coupled-NLSE system [24] and the multi-component long-wave-short-wave resonance model [60].

PS on the CW background. To reveal the amplification effect of the helicoidal SO coupling on the PS amplitude, we first address the PS solution on top of the flat CW

background, subject to constraint (14). The respective background amplitude (13) amounts to

$$| \Psi _ { 1 } ^ { c w } | = a \sqrt { \frac { k _ { m } + \alpha } { k _ { m } } }, \ | \Psi _ { 2 } ^ { c w } | = a \sqrt { \frac { k _ { m } - \alpha } { k _ { m } } }. \quad ( 1 5 ) \quad \text{for } \alpha \\ \tau, \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \end{array} \quad ( 1 5 ) \quad \text{for } \alpha \\ \tau, \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdots$$

Under the action of the SO coupling with strength α , the components of the CW background (15) have different heights.

Taking into regard that the center of the PS solution (10) is pinned to the origin, enhancement factor | F j | of component Ψ j is defined as the peak-to-background ratio:

$$\stackrel { \sim } { \sim } F _ { 1 } & \, = \, \frac { \Psi _ { 1 } ( 0, 0 ) } { | \Psi _ { 1 } ^ { c w } | } = \frac { a ( \nu _ { + } f _ { u _ { 1 } } + \nu _ { - } f _ { u _ { 2 } } ) } { | \Psi _ { 1 } ^ { c w } | }, \quad \text{the } \text{$the\,} \\ F _ { 2 } & \, = \, \frac { \Psi _ { 2 } ( 0, 0 ) } { | \Psi _ { 2 } ^ { c w } | } = \frac { a ( \nu _ { - } f _ { u _ { 1 } } - \nu _ { + } f _ { u _ { 2 } } ) } { | \Psi _ { 2 } ^ { c w } | }, \quad ( 1 6 ) \\ \text{where} \, | \Psi _ { 1 } ^ { c w } | \text{ are given in $Eq.$ (15), coefficients $\, \nu_{+}$ are}$$

where | Ψ cw 1 2 , | are given in Eq. (15), coefficients ν ± are same as in Eq. (6), and factors f u 1 and f u 2 are defined, for | δ | ≥ a , as

$$f _ { u _ { 1 } } & = 1 - \frac { 4 \zeta ^ { 2 } } { \zeta ^ { 2 } + ( \delta - \mu ) ^ { 2 } }, \ f _ { u _ { 2 } } = 1 - \frac { 4 \zeta ^ { 2 } } { \zeta ^ { 2 } + ( \delta + \mu ) ^ { 2 } }, \ ( 1 7 ) \quad \frac { \alpha ^ { 2 } } { - 4 } \\ \dots & \ e _ { - - } + \delta - \ e _ { - } + \delta - \ e _ { - }$$

and, for | δ | &lt; a , as

$$f _ { u _ { 1 } } = f _ { u _ { 2 } } = 1 - \frac { 2 ( \zeta + \mu ^ { \prime } ) ^ { 2 } } { 2 a ^ { 2 } + \zeta \mu ^ { \prime } }. \quad \quad ( 1 8 )$$

FIG. 1. ( a , b 1 1 ) An example of the fundamental PS , produced by solution (10) under condition (14), with the exceptionally high peak amplitude of the Ψ 1 component, for α = -1 2, / κ = 2 5. / ( a , b 2 2 ) A generic PS in the Manakov system, for α = 0. ( a , b 3 3 ) The PS with the zero background in Ψ 1 at x = 0, for α = -1 2, / κ = 0. The initial amplitude a = 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-e1df9678-d0dd-44fe-9646-f5440c8cec5d.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-443e641c-9793-41a7-81cb-21f97b8df470.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-5d94ca45-f8c6-4352-801d-7bd1bb40b294.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-9f2afd74-be10-4106-a817-575503eb96f6.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-73b29a71-a595-4494-86a2-61f6b2087be9.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-6984954d-86f1-4ee8-af11-436f00b905b5.jpg)

Characteristic examples of the PSs featuring large enhancement factors are presented in Fig. 1, which includes a PS with nearly fivefold peak amplitude for the component Ψ , with 1 α = -1 / 2 and κ = 2 5, in Fig. / 1( a 1 ). For comparison, two special cases are presented too, viz ., for α = 0 [Figs. 1( a 2 ) and ( b 2 )] and κ = 0 [Figs. 1( a 3 ) and ( b 3 )], which correspond to the Manakov system limit and the uniform SO coupling, respectively. It is observed that, with α = 0, the PS amplitudes are only twice as large as those of the background (in fact, for the Manakov system the peak amplitude cannot exceed three times the background value [24]). In addition, for κ = 0 the PS with zero background in component Ψ 1 or Ψ 2 is produced by Eq. (15), depending on the sign of α . Thus, the helicoidal SO coupling makes it possible to elevate the amplitude of one component to an exceptional level, while suppressing the other component.

FIG. 2. Enhancement factors F 1 and F 2 , as given by Eq. (16) for the PS with the flat background, vs. the SO-coupling strength α for a = 1, κ = 0 4 . , k 1 = -k 2 = -k m . The red and green curves correspond to µ and µ ′ taking signs + or -in Eqs. (11) and (12), respectively. The cyan dashed curves show the CW background values | Ψ cw j | , as given by Eq. (15). The insets exhibit the corresponding PSs at α = -1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-a7a1bb19-abd4-45f2-9973-bd6a78456b9b.jpg)

To further unveil the specific role of the helicoidal SO coupling in generating PSs with exceptionally high amplitudes, we display the dependence of enhancement factors F j on the SO-coupling strength α and rotation frequency κ in Figs. 2 and 3, respectively. They exhibit an indefinitely large (diverging) enhancement factor for component Ψ 1 or Ψ 2 at | α | →∞ or κ → 0. In particular, the insets to these figures feature the enhancement factor | F 1 | with values close to 10 at α = -1 and κ = 0 4, and . | F 2 | close to 5 at α = 0 6 and . κ = -0 6. . A caveat is that the enhancement factor is diverging when the background amplitude | Ψ CW j | is vanishing, as shown by the cyan dashed curves in Figs. 2 and 3. The absolute values of the PS peak amplitude may be increased by taking values of amplitude a &gt; 1 in Eq. (8) (recall it is currently fixed as a ≡ 1, by means of scaling).

PSs on top of the periodic background. If the constraint (14) does not hold, the above solution (10) produces the PS built on top of the periodic background. Similar to the case of the flat CW background considered above, we

FIG. 3. Enhancement factors F 1 and F 2 , as given by Eq. (16), vs. rotation frequency κ [see Eq. (3)] for a = 1, α = 0 6 . , k 1 = -k 2 = -k m . The red and green curves correspond to µ and µ ′ taking signs + or -in Eqs. (11) and (12), respectively. The cyan dashed curves show the CW background values | Ψ cw j | , as given by Eq. (15). The insets exhibit the corresponding PSs at κ = -0 6. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-1821f8b0-ca3d-4a82-ab15-630efba71097.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-6b1c74db-954b-4dd3-b8fe-55705fa75afb.jpg)

define the enhancement factor to analyze the effect of the helicoidal SO coupling on the PSs. In Fig. 4 we demonstrate a characteristic example exceeding the threefold contrast between the peak amplitude and periodic background in component Ψ , for 2 k 1 = -k m , k 2 = k m -3 / 2.

FIG. 4. The PS solution (10) built on top of the periodic background. It exceeds the exceeding the threefold enhancement limit, with k 1 = -k m and k 2 = k m -3 2. / The other parameters are a = 1 , α = -1 , κ = 0 4. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-8c65d219-58b9-4218-82c2-d0bec3ecaf33.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-e5a18acf-0f1d-44e0-b3ce-d9d7e46ee2d6.jpg)

Next, we address the modulation instability (MI) of the CWfield, Ψ j 0 = a e j iµt with µ = a 2 1 + a 2 2 , where a 1 and a 2 are the uniform amplitudes and µ is the chemical potential, in the presence of the helicoidal SO-coupled BECs. To this end, we add small perturbations to the CW fields, viz ., Ψ j = Ψ j 0 { 1 + p j exp[ -i ( βx -Ω )] + t q ∗ j exp[ ( i βx -Ω ∗ t )] } , where β and Ω are, respectively, the real and complex parameters, p j and q j being small complex amplitudes. Linearizing the corresponding Eq. (1) with respect to p j and q j , we derive a quartic equation for the perturbation eigenfrequency Ω, which determines the MI gain as γ h = Im(Ω) | | max . In Fig. 5, we display heatmaps for the so found value of γ h in the ( β, α ) and ( β, κ ) parameter planes. The plots reveal that the MI-gain spectra are symmetrically distributed in broad regions of α and κ , which can give rise to the anomalous PS behavior in a broad range of parameters, in comparison to the usual situation underlain by the baseband-MI analysis [61, 62].

For instance, Figs. 5( a 2 ) and ( b 2 ) demonstrate, respectively, that the gain maximum, γ h ≈ 1 91 at . α = -1 and κ = 0 4, corresponds to the enhancement factor . | F 1 | ≈ 10 in Fig. 2, and the maximum γ h ≈ 1 65, at . α = 0 6 and . κ = -0 6, corresponds to . | F 2 | ≈ 5 in Fig. 3. Such relatively small values of the MI gain, corresponding to the ultra-high PS peak amplitudes, suggest that these large amplitude values may be relatively easy to attain in the experiment, as the background will not be vulnerable to the quick destruction by of the MI-driven blowup, hence these PSs are rather robust modes.

FIG. 5. Heatmaps of the MI gain γ h in the ( β, α ) plane ( a 1 ) for κ = 0 4, and in the ( . β, κ ) plane ( b 1 ) for α = 0 6. . Panels ( a 2 ) and ( b 2 ) exhibit, respectively, the gain profile γ h at α = -1 and κ = -0 6, . with the maximum gain marked by black dots. The amplitudes of the underlying CW state are a 1 = a 2 = 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-d24cab70-8db2-417a-89b9-c290a2e32397.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-512bb6b0-05be-45cb-8213-aa98d5c82bd9.jpg)

FIG. 6. The result of the numerical simulations of the fundamental PSs from Figs. 1( a 1 ) and ( b 1 ) under the action of the 2% noise.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00322v1-figure-404f6bfb-fee4-440a-ad80-4ed281c142e2.jpg)

To test the expected robustness of the PSs in the present setting, in Fig. 6 we display results of the numerically simulated evolution of the PSs from Figs. 1( a 1 ) and ( b 1 ) under the action of 2% random disturbances. It is observed that the PSs with the ultra-high amplitudes indeed demonstrate robust propagation.

Conclusion. We have reported the occurrence of abnormal fundamental PSs (Peregrine solitons) with ultrahigh peak amplitudes in the integrable system of GPEs (Gross-Pitaevskii equations) including the helicoidally modulated SO (spin-orbit) coupling, which is a gauge isomer of the Manakov system. The PS solutions are found on top of both the flat and periodic-wave backgrounds. The results demonstrate the existence of RWs

(rogue waves with the ultra-high amplitude in the context of matter waves (BEC)), while previously this was reported in models of nonlinear optics [54, 55]. The helicoidal SO coupling is crucially important for generating this abnormal PSs, and the controllable nature of the SO coupling makes the predicted phenomenology experimentally feasible. The MI (modulational instability) is also studied in the system, demonstrating that the highamplitude PSs readily coexist with moderate MI, thus

- [1] K. Dysthe, H. E. Krogstad, P. M¨ uller, Oceanic rogue waves, Annual Rev. Fluid Mech. 40 , 287 (2008).
- [2] C. Kharif, E. Pelinovsky, and A. Slunyaev, Rogue Waves in the Ocean (Springer-Verlag: Berlin, 2009).
- [3] A. Chabchoub, N. P. Hoffmann, and N. Akhmediev, Rogue wave observation in a water wave tank, Phys. Rev. Lett. 106 , 204502 (2011).
- [4] A. Chabchoub, N. Hoffmann, M. Onorato, A. Slunyaev, A. Sergeeva, E. Pelinovsky, and N. Akhmediev, Observation of a hierarchy of up to fifth-order rogue waves in a water tank, Phys. Rev. E 86 , 056601 (2012).
- [5] A. Chabchoub, N. Hoffmann, M. Onorato, and N. Akhmediev, Super rogue waves: Observation of a higherorder breather in water waves. Phys. Rev. X 2 , 011015 (2012).
- [6] J. M. Dudley, G. Genty, A. Mussot, A. Chabchoub, F. Dias, Rogue waves and analogies in optics and oceanography, Nature Rev. Phys. 1 , 675 (2019).
- [7] D. R. Solli, C. Ropers, P. Koonath, B. Jalali, Optical rogue waves, nature, 450 , 1054 (2007).
- [8] C. Lecaplain, P. Grelu, J. M. Soto-Crespo, N. Akhmediev, Dissipative rogue waves generated by chaotic pulse bunching in a mode-locked laser, Phys. Rev. Lett. 108 , 233901 (2012).
- [9] W. R. Sun, B. Tian, Y. Jiang, and H. L. Zhen, Optical rogue waves associated with the negative coherent coupling in an isotropic medium, Phys. Rev. E 91 , 023205 (2015).
- [10] S. H. Chen, J. M. Soto-Crespo, F. Baronio, P. Grelu, and D. Mihalache, Rogue-wave bullets in a composite (2+1)D nonlinear medium, Opt. Exp. 24 , 15251-15260 (2016).
- [11] H. Bailung, S. K. Sharma, Y. Nakamura, Observation of Peregrine solitons in a multicomponent plasma with negative ions, Phys. Rev. Lett. 107 , 255005 (2011).
- [12] Y. V. Bludov, V. V. Konotop, N. Akhmediev, Matter rogue waves, Phys. Rev. A 80 , 033610 (2009).
- [13] A. Romero-Ros, G. C. Katsimiga, S. I. Mistakidis, B. Prinari, G. Biondini, P. Schmelcher, P. G. Kevrekidis, Theoretical and numerical evidence for the potential realization of the Peregrine soliton in repulsive two-component Bose-Einstein condensates, Phys. Rev. A 105 , 053306 (2022).
- [14] A. Romero-Ros, G. C. Katsimiga, S. I. Mistakidis, S. Mossman, G. Biondini, P. Schmelcher, P. Engels, P. G. Kevrekidis, Experimental realization of the Peregrine soliton in repulsive two-component Bose-Einstein condensates, Phys. Rev. Lett. 132 , 033402 (2024).
- [15] W. R. Sun, B. Tian, H. L. Zhen, and Y. Sun, Breathers and rogue waves of the fifth-order nonlinear Schr¨dinger o equation in the Heisenberg ferromagnetic spin chain,

preventing a strong background instability and improving chances for the experimental creation of the predicted tall rogue waves.

Acknowledgments This work was supported by the National Natural Science Foundation of China (Grant No. 11975172) and the Natural Science Foundation of Hubei Province of China (Grant No. 2023AFB222). The work of B.A.M. is supported, in part, by the Israel Science Foundation (Grant No. 1695/22).

Nonlin. Dynamics 81 , 725-732 (2015).

- [16] Z. Y. Yan, Financial rogue waves, Comm. Theor. Phys. 54 , 947-949 (2010).
- [17] Z. Y. Yan, Nonautonomous 'rogons' in the inhomogeneous nonlinear Schr¨dinger equation with variable coefo ficients, Phys. Lett. A 374 , 672-679 (2010).
- [18] Z. Yan, V. V. Konotop, and N. Akhmediev, Threedimensional rogue waves in nonstationary parabolic potentials, Phys. Rev. E 82 , 036610 (2010).
- [19] S. W. Xu, J. S. He, and L. H. Wang, The Darboux transformation of the derivative nonlinear Schr¨dinger equao tion, J. Phys. A: Math. Gen. 44 , 305203 (2011).
- [20] Y. S. Tao and J. S. He, Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation, Phys. Rev. E 85 , 026601 (2012).
- [21] J. S. He, H. R. Zhang, L. H. Wang, K. Porsezian, and A. S. Fokas, Generating mechanism for higher-order rogue waves, Phys. Rev. E 87 , 052914 (2013).
- [22] K. W. Chow, H. N. Chan, D. J. Kedziora, and R. H. J. Grimshaw, Rogue wave modes for the long wave-short wave resonance model, J. Phys. Soc. Jpn. 82 , 074001 (2013).
- [23] H. N. Chan, K. W. Chow, D. J. Kedziora, R. H. J. Grimshaw, and E. Ding, Rogue wave modes for a derivative nonlinear Schr¨dinger o model, Phys. Rev. E 89 , 032914 (2014).
- [24] S. Chen, D. Mihalache, Vector rogue waves in the Manakov system: diversity and compossibility, J. Phys. A 48 , 215202 (2015).
- [25] J. G. Rao, K. W. Chow, D. Mihalache, and J. S. He, Completely resonant collision of lumps and line solitons in the Kadomtsev-Petviashvili I equation, Stud. Appl. Math. 127 , 1007-1035 (2021).
- [26] L. Liu, W.-R. Sun, and B. A. Malomed, Formation of rogue waves and modulational instability with zerowavenumber gain in multi-component systems with coherent coupling, Phys. Rev. Lett. 131 , 093801 (2023).
- [27] D. H. Peregrine, Water waves, nonlinear Schr¨dinger o equations and their solutions, J. Aust. Math. Soc. Ser. B 25 , 16 (1983).
- [28] A. Tikan, C. Billet, G. El, A. Tovbis, M. Bertola, T. Sylvestre, F. Gustave, S. Randoux, G. Genty, P. Suret, J. M. Dudley, Universality of the Peregrine soliton in the focusing dynamics of the cubic nonlinear Schr¨dinger o equation, Phys. Rev. Lett. 119 , 033901 (2017).
- [29] B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, G. Genty, N. Akhmediev, J. M. Dudley, The Peregrine soliton in nonlinear fibre optics, Nature Phys. 6 , 790 (2010).
- [30] G. Michel, F. Bonnefoy, G. Ducrozet, G. Prabhudesai,

- A. Cazaubiel, F. Copie, A. Tikan, P. Suret, S. Randoux, E. Falcon, Emergence of Peregrine solitons in integrable turbulence of deep water gravity waves, Phys. Rev. Fluids 5 , 082801 (2020).
- [31] Y. J. Lin, K. Jim´ enez-Garc´ ıa, I. B. Spielman, Spinorbit-coupled Bose-Einstein condensates, Nature (London) 471 , 83 (2011).
- [32] V. Galitski, I. B. Spielman, Spin-orbit coupling in quantum gases, Nature (London) 494 , 49 (2013).
- [33] C. Hamner, Y. Zhang, M. A. Khamehchi, M. J. Davis, P. Engels, Spin-Orbit-Coupled Bose-Einstein Condensates in a One-Dimensional Optical Lattice, Phys. Rev. Lett. 114 , 070401 (2015).
- [34] H. Zhai, Degenerate quantum gases with spin-orbit coupling: a review, Rep. Prog. Phys. 78 , 026001 (2015).
- [35] Y. Zhang, M. E. Mossman, T. Busch, P. Engels, and C. Zhang, Properties of spin-orbit-coupled Bose-Einstein condensates, Front. Phys. 11 , 118103 (2016).
- [36] J. Dalibard, F. Gerbier, G. Juzeli¯nas, P. Ohberg, Colu ¨ loquium: Artificial gauge potentials for neutral atoms, Rev. Mod. Phys. 83 , 1523 (2011).
- [37] J. Ruseckas, G. Juzeli¯nas, u P. ¨ hberg, O and M. Fleischhauer, Non-Abelian Gauge Potentials for Ultracold Atoms with Degenerate Dark States, Phys. Rev. Lett. 95 , 010404 (2005).
- [38] J. Struck, C. ¨ lschl¨ger, O a M. Weinberg, P. Hauke, J. Simonet, A. Eckardt, M. Lewenstein, K. Sengstock, P. Windpassinger, Tunable Gauge Potential for Neutral and Spinless Particles in Driven Optical Lattices, Phys. Rev. Lett. 108 , 225304 (2012).
- [39] Y. Zhang, G. Chen, and C. Zhang, Tunable spin-orbit coupling and quantum phase transition in a trapped Bose-Einstein condensate, Sci. Rep. 3 , 1937 (2013).
- [40] K. Jim´ enez-Garc´ ıa, L. J. LeBlanc, R. A. Williams, M. C. Beeler, C. Qu, M. Gong, C. Zhang, I. B. Spielman, Tunable Spin-Orbit Coupling via Strong Driving in Ultracold-Atom Systems, Phys. Rev. Lett. 114 , 125301 (2015).
- [41] X. Luo, L. Wu, J. Chen, Q. Guan, K. Gao, Z.-F. Xu, L. You, R. Wang, Tunable atomic spin-orbit coupling synthesized with a modulating gradient magnetic field, Sci. Rep. 6 , 18983 (2016).
- [42] S. S. S. Hejazi, J. Polo, R. Sachdeva, T. Busch, Symmetry breaking in binary Bose-Einstein condensates in the presence of an inhomogeneous artificial gauge field, Phys. Rev. A 102 , 053309 (2020).
- [43] Y. V. Kartashov, V. V. Konotop, and D. A. Zezyulin, Bose-Einstein condensates with localized spin-orbit coupling: Soliton complexes and spinor dynamics, Phys. Rev. A 90 , 063621 (2014).
- [44] Y. V. Kartashov and V. V. Konotop, Solitons in BoseEinstein Condensates with Helicoidal Spin-Orbit Coupling, Phys. Rev. Lett. 118 , 190401 (2017).
- [45] Y. V. Kartashov, E. Sherman, B. Malomed, and V. Konotop, Stable two-dimensional soliton complexes in BoseEinstein condensates with helicoidal spin-orbit coupling, New J. Phys. 22, 103014 (2020).
- [46] Y. Li, X. Zhang, R. Zhong, Z. Luo, B. Liu, C. Huang, W. Pang, B. A. Malomed, Two-dimensional composite solitons in Bose-Einstein condensates with spatially confined spin-orbit coupling, Comm. Nonlin. Sci. Num. Sim. 73 , 481-489 (2019).
- [47] B. Liu, R. Zhong, Z. Chen, X. Qin, H. Zhong, Y. Li and
- B. A. Malomed, Holding and transferring matter-wave solitons against gravity by spin-orbit-coupling tweezers, New J. Phys. 22 , 043004 (2020).
- [48] M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev, A. Szameit, Photonic Floquet topological insulators, Nature (London) 496 , 196 (2013).
- [49] Y. V. Kartashov, V. V. Konotop, M. Modugno, E. Ya. Sherman, Solitons in Inhomogeneous Gauge Potentials: Integrable and Nonintegrable Dynamics, Phys. Rev. Lett. 122 , 064101 (2019).
- [50] F. Baronio, A. Degasperis, M. Conforti, S. Wabnitz, Solutions of the vector nonlinear Schr¨dinger equations: evo idence for deterministic rogue waves, Phys. Rev. Lett. 109 , 044102 (2012).
- [51] S. V. Manakov, On the theory of two-dimensional stationary self-focusing of electromagnetic waves, Zh. Eksp. Teor. Fiz. 65 , 505-516 (1973) [Sov. Phys. JETP 38 , 248253 (1974)].
- [52] F. Baronio, M. Conforti, A. Degasperis, S. Lombardo, Rogue waves emerging from the resonant interaction of three waves, Phys. Rev. Lett. 111 , 114101 (2013).
- [53] S. Chen, F. Baronio, J. M. Soto-Crespo, Ph. Grelu, M. Conforti, S. Wabnitz, Optical rogue waves in parametric three-wave mixing and coherent stimulated scattering, Phys. Rev. A 92 , 033847 (2015).
- [54] S. Chen, Y. Ye, J. M. Soto-Crespo, P. Grelu, F. Baronio, Peregrine Solitons Beyond the Threefold Limit and Their Two-Soliton Interactions, Phys. Rev. Lett. 121 , 104101 (2018).
- [55] S. Chen, C. Pan, P. Grelu, F, Baronio, N. Akhmediev, Fundamental Peregrine Solitons of Ultrastrong Amplitude Enhancement through Self-Steepening in Vector Nonlinear Systems, Phys. Rev. Lett. 124 , 113901 (2020).
- [56] X. X. Li, R. J. Cheng, A. X. Zhang, J. K. Xue, Modulational instability of Bose-Einstein condensates with helicoidal spin-orbit coupling, Phys. Rev. E 100 , 032220 (2019).
- [57] X. X. Li, R. J. Cheng, J. L. Ma, A. X. Zhang, J. K. Xue, Solitary matter wave in spin-orbit-coupled Bose-Einstein condensates with helicoidal gauge potential, Phys. Rev. E 104 , 034214 (2021).
- [58] S. V. Samsonov, A. D. R. Phelps, V. L. Bratman, G. Burt, G. G. Denisov, A. W. Cross, K. Ronald, W. He, H. Yin, Compression of frequency-modulated pulses using helically corrugated waveguides and its potential for generating multigigawatt rf radiation, Phys. Rev. Lett. 92 , 118301 (2004).
- [59] G. Burt, S. V. Samsonov, K. Ronald, G. G. Denisov, A. R. Young, V. L. Bratman, A. D. R. Phelps, A. W. Cross, I. V. Konoplev, W. He, J. Thomson, C. G. Whyte, Dispersion of helically corrugated waveguides: Analytical, numerical, and experimental study, Phys. Rev. E 70 , 046402 (2004).
- [60] S. Chen, J. M. Soto-Crespo, P. Grelu, Coexisting rogue waves within the (2+1)-component long-wave-shortwave resonance, Phys. Rev. E 90 , 033203 (2014).
- [61] F. Baronio, M. Conforti, A. Degasperis, S. Lombardo, M. Onorato, S. Wabnitz, Vector Rogue Waves and Baseband Modulation Instability in the Defocusing Regime, Phys. Rev. Lett. 113 , 034101 (2014).
- [62] F. Baronio, S. Chen, Ph. Grelu, S. Wabnitz, M. Conforti, Baseband modulation instability as the origin of rogue waves, Phys. Rev. A 91 , 033804 (2015).
Cui-Cui Ding Qin Zhou B. A. Malomed
Ultra-high-amplitude Peregrine solitons induced by helicoidal spin-orbit coupling
In the framework of the model of a spatially non-uniform Bose-Einstein
condensate with helicoidal spin-orbit (SO) coupling, we find abnormal Peregrine
solitons (PSs) on top of flat and periodic backgrounds, with ultra-high
amplitudes. We explore the roles of the SO coupling strength and helicity pitch
in the creation of these anomalously tall PSs and find that their amplitude,
normalized to the background height, attains indefinitely large values. The
investigation of the modulation instability (MI) in the same system
demonstrates that these PSs exist in a range of relatively weak MI, maintaining
the feasibility of their experimental observation.