# Several Localized Waves Induced by Linear Interference between a Nonlinear Plane Wave and Bright Solitons

###### Abstract

We investigate linear interference effects between a nonlinear plane wave and bright solitons, which are admitted by pair-transition coupled two-component Bose-Einstein condensate. We demonstrate the interference effects can induce several localized waves possessing distinctive wave structures, mainly including anti-dark soliton, W-shaped soliton, multi-peak soliton, Kuznetsov-Ma like breather, and multi-peak breather. Especially, the explicit conditions for them are clarified by a phase diagram based on the linear interference properties. Furthermore, the interactions between these localized waves are discussed. The detailed analysis indicate that soliton-soliton interaction induced phase shift brings the collision between these localized waves be inelastic for soliton involving collision, and be elastic for breathers. These characters come from that the profile of solitons depend on relative phase between bright soliton and plane wave, and the profile of breathers do not depend on the relative phase. These results would motivate more discussions on linear interference between other nonlinear waves. Especially, the solitons or breathers obtained here are not related with modulational instability. The underlying reasons are discussed in detail.

###### pacs:

05.45.Yv, 02.30.Ik, 42.65.Tg## I Introduction

Nonlinear wave (NW) has been paid much attention since 1960s Zabusky (); nw (). Many different families of NWs have been found, such as bright soliton BS (), nonlinear plane wave, dark soliton ds (), kink kink (), breather br (), rogue wave rw (), and vector ones of them Kevrekidis (); Lakshmanan (); Bludov (); Zhao (); Zhao1 (); Zhao2 (); Baronio (); Guo (). Many of them have been observed in real experiments or in natural environment rex1 (); rex2 (); rex3 (); Becker (); Engels (). Nowadays, many efforts are still paid to find new families of NWs, uncover the interaction properties of them, and explain the mechanism of them. Furthermore, the applications of them are discussed in many different physical systems. For examples, vector bright soliton was used to generate Bell state in ultra-cold atoms Bells (). Very recently, bright soliton interferometry was proposed in Bose-Einstein condensate (BEC) bsinter1 (); bsinter2 (). The interference between solitons or other NWs can be investigated by their proper superposition nodyzhao (). Since the NWs are admitted by nonlinear systems, the linear superposition form usually fails to describe the interaction between NWs DT1 (); DT2 (); Cheng (). Therefore, it is unusual to investigate dynamical process described by linear superposition between these NWs. In contrast to the nonlinear interference process described by nonlinear superposition nodyzhao (), the interference process described by linear superposition of NWs is seen as a linear interference process. The linear interference effects between them could induce some different wave structures, in contrast to the nonlinear interactions.

We note that it is possible to obtain some linear superposition forms of NWs in a pair-transition coupled two-component BEC system Park (); Tian (); BoTian (); zhaoling (); lingzhao (). This makes it be possible to investigate linear interference between NWs in a real physical system. Moreover, the knowledge of linear interference in wave theory can be used to explain the dynamical process clearly. Therefore, we study on linear interference between NWs in a pair-transition coupled two-component BEC system. Among the well-known NWs, bright soliton and nonlinear plane wave are two of the simplest cases. For simplicity and without losing generality, we would like to discuss the simplest case here to uncover the essential properties of linear interference between NWs.

In this paper, we demonstrate the linear interference effects between bright soliton and nonlinear plane wave can induce several different localized waves possessing distinctive wave structures, mainly including anti-dark soliton, W-shaped soliton, multi-peak soliton, Kuznetsov-Ma (K-M) like breather, and multi-peak breather. It is emphasized that these localized waves are different from anti-dark soliton, W-shaped soliton, multi-peak soliton and K-M breather reported before in other nonlinear systems Ws (); AD (); Liuchong (); K-Mb (). Explicitly, previously reported ones admitted nonlinear superposition form and more complicated expressions. The existence conditions for them are clarified clearly based on the spatial-temporal interference factors (summarized in Fig. 1). The relative phase between bright soliton and plane wave are also discussed in detail for the different four cases. Furthermore, the interactions between these localized waves are discussed. The detailed analysis indicate that soliton-soliton interaction induced phase shift brings the collision between these localized waves can be inelastic for soliton involving collision, and be elastic for breathers. The underlying reasons for these characters are discussed.

The outline of this paper is as follows. In Sec. \@slowromancapii@, we present the physical model and linear superposition form between bright solitons and a nonlinear plane wave. In Sec. \@slowromancapiii@, we analyze the distinctive structures of several nonlinear localized waves induced by linear interference effects between one soliton and a plane wave. A phase diagram for these localized waves are presented. In Sec. \@slowromancapiv@, we discuss the interactions between these nonlinear localized waves. The collision between them can be elastic and inelastic. The explicit conditions for elastic or inelastic collision are clarified. The underlying reason for these collision characters is explained well. In Sec. \@slowromancapv@, we explain why soliton or breather obtained here do not involve with modulational instability. Furthermore, we suggest a possible way to clarify which dispersion form is chosen for a weak perturbation, when the perturbations admit more than one dispersion forms. The conclusion and discussion are made in Sec. \@slowromancapvi@.

## Ii The physical model and linear superposition form solution

One-dimensional two-component BEC system with particle transition can be described by the Hamiltonian where is the particle number operator, the symbol represents the Hermite conjugation. and are the intra and external interactions between atoms. and denote single particle and pair particles transition coupling strength separately Fischer1 (); Fischer2 (). In most studies, are set to be zero usually because it was believed that the presence of tunneling makes the systems become non-integrable Baronio (); Zhao (); Lakshmanan (); Bludov (). Recent experimental results in a double-well Bose-Einstein condensate suggested that pair-tunneling can become dominant with strong interaction between atoms pair (); Meyer (). Therefore, we consider that the case for second-order transition is dominant, namely, and . We find integrable CNLS-p can be derived from the Hamiltonian with . It is convenient to set ( correspond to attractive or repulsive interactions between atoms) without losing generality, since there is a trivial scalar transformation for different values.

The corresponding dynamic evolution equation can be derived from the Heisenberg equation for the field operator. Performing the mean field approximation , we can get the following integrable CNLS-p with scale dimensions .

(1) |

where the symbol overbar represents the complex conjugation. What needs mentioning is that the above coupled equations without the last term usually deem as non-integrable CNLS yang-benney (). However, when we add the particles transition term, the non-integrable CNLS become integrable, which was also proven by Painlevé analysis Park ().

Especially, the CNLS-p model can be transformed to two uncoupled NLS equations through a linear transformation Park (); lingzhao (); Lvling (); Kanna (); Xiang (). The solutions of Eq. (1) can be written in the form and , where and are solutions for the scalar NLS (). This is a striking character for the coupled model with particle transition terms, in contrast to the CNLS without particle transition terms Lakshmanan (); Zhao (); Zhao1 (); Zhao2 (); Baronio (); Guo (). The scalar NLS admits many different NWs, such as bright soliton, breather, rogue wave, and nonlinear plane wave. If is one NW, and is another NW, the linear superposition forms can be used to investigate the linear interference between two types of NWs. Here we would like to discuss the simplest case for bright solitons and a nonlinear plane wave. Namely, is a general form of nonlinear plane wave, and is bright soliton solution of the scalar NLS. The bright soliton solution can be one and multi-soliton solutions DT1 (); DT2 (). This makes us possible to investigate linear interference between bright solitons and nonlinear plane wave systemically. Firstly, we discuss the case for one bright soliton and a nonlinear plane wave.

## Iii Nonlinear localized waves induced by linear interference effects

A general linear superposition form of one bright soliton and a nonlinear plane wave can be given as follows:

(2) |

where and denote the plane wave amplitude and soliton amplitude respectively. , and describe the phase evolution of them respectively. The evolution characters are mainly determined by the phase difference between them. The parameter is a constant phase difference between them, which is called as relative phase here. In the following, we discuss properties of the nonlinear waves according to the two parameters and . It is found that there are mainly four cases for interference properties according to the interference period on spatial and temporal direction. The conditions for them are summarized in Fig. 1. It is obtained with variable wave vector and amplitude of the nonlinear plane wave, and fixed soliton’s amplitude and velocity. Similar phase diagram can be obtained in other cases.

Case 1: When and , the temporal spatial structures of solution is presented in Fig. 2(a,b). This case corresponds to the green point (type I soliton) in Fig. 1. A W-shaped soliton appears in the density distribution of the second component , while an anti-dark soliton appears in the density distribution of the first component . Interestingly, the soliton profile depends on the relative phase value. To show how the soliton profile depend on the relative phase, we define relative oscillation to be (where denotes the maximum values of density value at vs the relative phase). This is depicted in Fig. 3(a). For , the W-shaped soliton varies to be an anti-dark soliton. It should be emphasized that the anti-dark soliton and W-shaped soliton here are distinctive from the ones obtained in other systems AD (); Ws (). The superposition forms are different and their spectrums are also different, since the previously reported anti-dark soliton and W-shaped soliton all admit a nonlinear superposition form.

Case 2: When and , densities of both components show multi-peak soliton, such structures are shown in Fig. 2(c,d). This corresponds to the red solid line (type II soliton) in Fig. 1. Similarly, the multi-peak soliton profile also depends on the relative phase . However, the peak change with different in this case is much weak than the change of the case 1 (see Fig.3(b)). The visible peak number is more, and the sensitiveness on relative phase is weaker, vise versa. Especially, the visible peak number is determined by the soliton visible size and the wave vector of plane wave. The spatial distance between humps is , and this value should be much smaller than the soliton visible size to show multi-peak profile. Therefore, when is very small, is relatively larger, multi-peak structure could be observed clearly. We emphasize that the multi-peak soliton here is also distinctive from the ones obtained in Liuchong (), namely, their interference process and profiles are both distinctive.

Case 3: When and , nonlinear localized wave is periodic on temporal direction and the behavior of density distribution is shown in Fig. 4(a,b). This corresponds to the blue dashed line (type I breather) in Fig. 1. We can see that the breathing behavior is analogous to the K-M like breather of scalar NLS K-Mb (). But the amplification rate is much smaller than the K-M breather, since there is no modulational instability gain value in this case. Interestingly, it is seen that there are humps and valleys alternately emerging in each component, and the location of hump in the component just corresponds to the valley in the component . The breathing period is . This indicates that there are particles transition between the two components during the breathing process. The transition process admits a standard Josephson oscillation form, which is distinctive from the nonlinear one reported in zhaoling (). The profile of localized wave varies with time evolution for breather. The spatial-temporal distribution properties do not depend on the relative phase between bright soliton and plane wave anymore.

Case 4: When and , we show the dynamics of the nonlinear excitation in Fig. 4(c,d). This corresponds to regime except the points on the two lines in Fig. 1. The condition makes the localized wave admit multi-peak profile, and makes the localized wave breath with time evolution. To distinct from the type I breather, we call this as multi-peak breather, which corresponds to the type II breather in Fig. 1. Similar to multi-peak soliton, the visible peak number is also determined by the soliton visible size and the wave vector of plane wave. In this case, the spatial-temporal distribution properties also do not depend on the relative phase between bright soliton and plane wave anymore.

Based on the above discussions on soliton types and profiles, one can expect that the interaction between them can be investigated analytically and exactly with the aid of multi-soliton solution of NLS bscollide (). For simplicity, we just discuss the interaction between two of them.

## Iv The interactions between several obtained localized waves

The solution for two bright solitons linearly interfering with a nonlinear plane wave could be presented as follows:

(3) |

where denotes a nonlinear superposition to the bright soliton. , , , , and , , . The parameters and are related soliton’s velocity, and determine peak value of solitons respectively. and determine the initial locations of solitons, and can be used to vary the relative phase between solitons. This solution is a linear superposition of a plane wave and a two-soliton solution, but the two-soliton solution is a nonlinear superposition of two bright solitons. When the related parameters are chosen, the solution will present us the dynamics of two localized waves directly. The localized wave type can be chosen by setting the parameters of bright soliton or plane wave background based on the phase diagram Fig. 1. The collision of arbitrary two of the above localized waves can be observed conveniently. It is found that collision between these localized waves can be inelastic for soliton involving collision, and be elastic for breathers. Nextly, we discuss on the collision between them in details.

The collision between soliton and soliton or breather can be investigated by setting the bright soliton ’s parameters satisfy different conditions in Fig. 1. For examples, we show one case for type I soliton and type II soliton in Fig. 5 (a) and (b) by setting parameters . The type II soliton admit multi-peak here, which is invisible in Fig. 5 (a) and (b). This comes from that the period is much larger than the bright soliton scale. It is shown clearly that solitons’ profiles both vary after colliding process. Similarly, we can investigate collision between soliton and breather. It should be noted that it is not possible to observe collisions between type I soliton and type I breather, since their existence conditions make them admit identical velocity. From the phase diagram, we can observe the collision between type I soliton and type II breather by choosing , , , , , . Fig. 5 (c) and (d) show the collision between them, which correspond to component and component respectively. It is seen that soliton’s profile varies after collision and breather’s profile is kept well.

The collision between breather and breather can be investigated by setting the bright soliton ’s parameters satisfy different conditions in Fig. 1. For examples, type II breather and type I breather can be set by choosing , , , , , . Fig. 6 (a) and (b) show the collision between them, which correspond to component and component respectively. However, it is not possible to observe the collision between type I breather and type I breather, because the existence conditions make them admit identical velocity. Type II breather and type II breather can be investigated by choosing , , , , , . Fig. 6 (c) and (d) show the collision between them, which correspond to component and component respectively. It is shown that the collision between breathers are all elastic.

Then, why the soliton’s profile varies after collision? It is well known that there is a phase shift after bright soliton colliding with another one bscollide (); Nguyen (). The phase shift brings relative phase changes between bright soliton and plane wave. Moreover, the soliton type localized waves’ profiles depend on the relative phase. This makes the soliton’s profile change. Therefore, the underlying reason for inelastic collision is soliton type localized waves depend on the relative phase between bright soliton signal and plane wave, and collision between two bright solitons can bring a phase shift on each bright soliton signal. But the breather type localized waves do not depend on the relative phase, the phase shift will not change the spatial-temporal structure of breathers.

## V A discussion on modulational instability branches

It should be noted that the soliton-type perturbation on the plane wave background can be both weak and strong here. We emphasize that the breathing behaviors here do not involve with modulational instability (MI). This can be seen by the above localized wave dynamics for which the perturbation amplitudes are not amplified at all even the soliton-type perturbation is very weak. MI has been shown widely to induce rogue wave or breather on a plane wave background Kibler (); Dudley (); Baronio1 (); Baronio3 (); zhaoling2 (); NMI (). In fact, rogue wave or Akhmediev breather can also exist on the identical plane wave background in this coupled model. This can be proven by the transformation between the coupled model and standard scalar NLSE. Namely, a rogue wave solution of the coupled model can be constructed by one rogue wave solution with the identical plane wave background and one zero solution of the scalar NLSE. This means that there are at least two dispersion relation branches on the plane wave background. We perform standard MI analysis on the plane wave background in the coupled model. We add the perturbation terms on the plane wave background, and (where are small amplitudes of the Fourier modes). It is found that there are two dispersion relation branches, which admit one MI branch and one modulational stability (MS) branch. They can be calculated as , and , which correspond to MI and MS branch respectively. Since the linear stability analysis holds well for weak perturbations, there will be two possible choices for a weak perturbation on the plane wave background. Then, how to understand that the soliton-type perturbation with small amplitude evolves to be type-I or type II soliton or breather which do not involves MI characters? Namely, why the weak soliton-type perturbation choose the MS branch to evolve?

We try to find the essential factors which determine choosing choice for weak perturbations, by calculating the eigenvector of the linearized equations. For MI branch (), the corresponding eigenvector is . must holds since the linear stability analysis involves linearization. It is seen that , for MI branch, namely, the perturbations form added on the backgrounds in the two components admit identical perturbation profile. For an example, we show one case for perturbation amplitude vs perturbation wave vector in Fig. 7 (a). Considering that the rogue wave solution discussed above is constructed by one rogue wave solution with the identical plane wave background and one zero solution of the scalar NLSE, we can understand that the rogue wave can exist on the background and they demonstrate MI characters, since the rogue wave solutions in the two components admit identical perturbation profile on the background and this makes the weak perturbation choose the MI branch to evolve.

For MS branch , the corresponding eigenvector can be expressed as . For , the corresponding eigenvector is . We can see that the perturbation profiles are inverse in the two components for MS branch. As an example, we present the perturbation amplitude and versus the perturbation wave vector for one MS branch in Fig. 7(b). The weak soliton-type perturbations presented above all satisfy the for MS branch. Namely, the type-I or type II soliton or breather obtained here all admit the case that perturbation profiles on the background are inverse in the two components. This makes that the perturbations all choose MS branch to evolve. The breathers discussed above just involves linear interference mechanism and is not related with MI at all.

## Vi Conclusion and discussion

In this paper, we show that linear interference effects between a nonlinear plane wave and bright soliton, can induce anti-dark soliton, W-shaped soliton, K-M like breather, multi-peak soliton, and multi-peak breather. Their profile properties are discussed in detail. It is shown that soliton type localized waves’ profile depends on relative phase between bright soliton and plane wave, but the breathers do not. Furthermore, the interactions between these localized waves are discussed. The detailed analysis indicate that soliton interaction induced phase shift brings the collision between these localized waves can be inelastic and elastic. The underlying reason for these characters are discussed. Additionally, considering the soliton profile depends on the relative phase between plane wave and bright soliton, we expect that the inelastic collision property could be used to measure the phase shift during the collision between bright solitons. Similar studies can be extended to linear interference between other nonlinear waves in other coupled systems. Especially, the soliton or breather obtained here is proven to be not related with MI. A possible way is suggested to clarify which MI branch is chosen for a weak perturbation, when the perturbations admit more than one MI branches.

## Acknowledgments

This work is supported by National Natural Science Foundation of China (Contact No. 11775176), and Shaanxi Province Science association of colleges and universities (Contact No. 20160216).

## References

- (1) N.J. Zabusky, M.D. Kruskal, Phys. Rev. Lett. 15, 240 (1965).
- (2) R. Carretero-Gonz¨¢lez , D. J. Frantzeskakis, P. G. Kevrekidis, Nonlinearity 21, R139 (2008).
- (3) Khaykovich L, Schreck F, Ferrari G, et al., Science 296, 1290-1293 (2002).
- (4) Burger S, Bongs K, Dettmer S, et al., Phys. Rev. Lett. 83, 5198 (1999).
- (5) P. O. Fedichev, A. E. Muryshev, and G. V. Shlyapnikov, Phys. Rev. A 60, 3220 (1999).
- (6) M. Matuszewski, E. Infeld, B.A. Malomed, et al., Phys. Rev. Lett. 95, 050403 (2005).
- (7) Yu. V. Bludov, V. V. Konotop, and N. Akhmediev, Phys. Rev. A 80, 033610 (2009).
- (8) P. G. Kevrekidis, D. Frantzeskakis, and R. Carretero-Gonzalez, Emergent Nonlinear Phenomena in Bose-Einstein Condensates: Theory and Experiment (Springer, Berlin Heidelberg, 2009).
- (9) T. Kanna and M. Lakshmanan, Phys. Rev. Lett. 86, 5043-5046(2001); M. Vijayajayanthi, T. Kanna, and M. Lakshmanan, Phys. Rev. A 77, 013820 (2008).
- (10) Y.V. Bludov, V.V. Konotop, and N. Akhmediev, Eur. Phys. J. Special Topics 185, 169 (2010).
- (11) L.C. Zhao, J. Liu, J. Opt. Soc. Am. B 29, 3119-3127 (2012).
- (12) L.C. Zhao, J. Liu, Phys. Rev. E 87, 013201 (2013).
- (13) F. Baronio, A. Degasperis, M. Conforti, and S. Wabnitz, Phys. Rev. Lett. 109, 044102 (2012).
- (14) L.C. Zhao, G.G. Xin, Z.Y. Yang, Phys. Rev. E 90, 022918 (2014).
- (15) B.L. Guo and L.M. Ling, Chin. Phys. Lett. 28, 110202 (2011).
- (16) Kibler, B., Fatome, J., Finot, C., Millot, G., et al. Nat. Phys. 6, 790 (2010).
- (17) Bailung, H., Sharma, S.K., Nakamura, Y., Phys. Rev. Lett. 107, 255005 (2011).
- (18) Chabchoub, A., Hoffmann, N.P., Akhmediev, N.,Phys. Rev. Lett. 106, 204502 (2011).
- (19) C. Becker, S. Stellmer, P.S. Panahi, S. Dorscher, M. Baumert, Eva-Maria Richter, J. Kronjager, K. Bongs, K. Sengstock, Nature Phys. 4, 496-501 (2008).
- (20) C. Hamner, J. J. Chang, and P. Engels, Phys. Rev. Lett. 106, 065302 (2011); M. A. Hoefer, J. J. Chang, C. Hamner, and P. Engels, Phys. Rev. A 84, 041605(R) (2011).
- (21) B. Gertjerenken, T.P. Billam, C.L. Blackley, et al., Phys. Rev. Lett. 111, 100406 (2013).
- (22) McDonald, G.D., Kuhn, C.C.N., Hardman, K.S., Bennetts, S., Everitt, P.J., Altin, P.A., Debs, J.E., Close, J.D., Robins, N.P., Phys. Rev. Lett. 113, 013002 (2014).
- (23) Polo, J., Ahufinger, V., Phys. Rev. A 88, 053628 (2013).
- (24) L.C. Zhao, L. Ling, Z.Y. Yang, J. Liu, Nonlinear Dyn. 83, 659-665 (2016).
- (25) V.B. Matveev and M.A. Salle, Darboux Transformation and Solitons (Springer-Verlag, Berlin, 1991).
- (26) B.L. Guo, L.M. Ling, Q. P. Liu , Phys. Rev. E 85, 026607 (2012).
- (27) X.-P. Cheng, S. Y. Lou, C. Chen, and X.-Y. Tang, Phys. Rev. E 89, 043202 (2014).
- (28) Q-Han Park and H.J. Shin, Phys. Rev. E 59, 2373 (1999).
- (29) X. Lü, and B. Tian, Phys. Rev. E 85, 026117 (2012).
- (30) W.R. Sun, B. Tian, Y. Jiang, and H.L. Zhen, Phys. Rev. E 91, 023205 (2015).
- (31) L. Ling and L.C. Zhao, Phys. Rev. E 92, 022924 (2015).
- (32) L.C. Zhao, L. Ling, Z.Y. Yang, J. Liu, Commun. Nonlinear Sci. Numer. Simulat. 23, 21-27 (2015).
- (33) L.C. Zhao, S.C. Li, L. Ling, Phys. Rev. E 89, 023210 (2014).
- (34) L. Duan, L.C. Zhao, W.H. Xu, et al., Phys. Rev. E 95, 042212 (2017).
- (35) C. Liu, Z.Y. Yang, L.C. Zhao, et al., Phys. Rev. E 94, 042221 (2016).
- (36) E. Kuznetsov, Sov. Phys. Dokl. 22, 507-508 (1977); Y.C. Ma, Stud. Appl. Math. 60, 43-58 (1979).
- (37) P. Bader, U.R. Fischer, Phys. Rev. Lett. 103, 060402 (2009).
- (38) U.R. Fischer, K.S. Lee, B. Xiong, Phys. Rev. A 84, 011604 (2011).
- (39) S. Flling, S. Trotzky, P. Cheinet, et al., Nature 448, 06112 (2007).
- (40) S. Zöllner, H.D. Meyer, and P. Schmelcher, Phys. Rev. Lett. 100, 040401 (2008).
- (41) J. Yang and D. J. Benney, Stud. Appl. Math. 96, 111 (1996).
- (42) R. Babu Mareeswaran, T. Kanna, Phys. Lett. A 380, 3244-3252 (2016).
- (43) Xing Lü, and Liming Ling, CHAOS 25, 123103 (2015).
- (44) Rui Xiang, Liming Ling, Xing Lü, App. Math. Lett. 68, 163-170 (2017).
- (45) N.G. Parker, A.M. Martin, S.L. Cornish, C.S. Adams, J. Phys. B: At. Mol. Opt. Phys. 41, 045303 (2008).
- (46) J.H.V. Nguyen, P. Dyke, D. Luo, B.A. Malomed, R.G. Hulet, Nature Phys. 10, 918-922 (2014).
- (47) B. Kibler, J. Fatome, C. Finot, G. Millot, et al., Nature Phys. 6,790 (2010).
- (48) J. M. Dudley, G. Genty, F. Dias, B. Kibler, and N. Akhmediev, Opt. Express 17, 21497 (2009).
- (49) F. Baronio, M. Conforti, A. Degasperis, S. Lombardo, M. Onorato, and S. Wabnitz, Phys. Rev. Lett. 113, 034101 (2014).
- (50) F. Baronio, S. Chen, P. Grelu, S. Wabnitz and M. Conforti, Phys. Rev. A 91, 033804 (2015).
- (51) L.C. Zhao, L. Ling, J. Opt. Soc. Am. B 33, 850-856 (2016).
- (52) G. Biondini, and D. Mantzavinos, Phys. Rev. Lett. 116, 043902 (2016).