A Phase-Fitted Runge-Kutta-Nystr\"om method for the Numerical Solution of Initial Value Problems with Oscillating Solutions

A new Runge-Kutta-Nystr\"om method, with phase-lag of order infinity, for the integration of second-order periodic initial-value problems is developed in this paper. The new method is based on the Dormand and Prince Runge-Kutta-Nystr\"om method of algebraic order four\cite{pa}. Numerical illustrations indicate that the new method is much more efficient than the classical one.Comment: 10 page

Dark soliton scattering in symmetric and asymmetric double potential barriers

The file attached to this record is the author's final peer reviewed version. The Publisher's final version can be found by following the DOI link.Motivated by the recent theoretical study of (bright) soliton diode effects in systems with multiple scatterers, as well as by experimental investigations of soliton-impurity interactions, we consider some prototypical case examples of interactions of dark solitons with a pair of scatterers. In a way fundamentally opposite to the case of bright solitons (but consonant to their “anti-particle character”), we find that dark solitons accelerate as they pass the first barrier and hence cannot be trapped by a second equal-height barrier. A pair of unequal barriers may lead to reflection from the second one, however trapping in the inter-barrier region cannot occur. We also give some examples of dynamical adjusting of the barriers to trap the dark soliton in the inter-barrier region, yet we show that this can only occur over finite time horizons, with the dark soliton always escaping eventually, contrary again to what is potentially the case with bright solitons

Dark soliton scattering in symmetric and asymmetric double potential barriers

Motivated by the recent theoretical study of (bright) soliton diode effects in systems with multiple scatterers, as well as by experimental investigations of soliton-impurity interactions, we consider some prototypical case examples of interactions of dark solitons with a pair of scatterers. In a way fundamentally opposite to the case of bright solitons (but consonant to their “anti-particle character”), we find that dark solitons accelerate as they pass the first barrier and hence cannot be trapped by a second equal-height barrier. A pair of unequal barriers may lead to reflection from the second one, however trapping in the inter-barrier region cannot occur. We also give some examples of dynamical adjusting of the barriers to trap the dark soliton in the inter-barrier region, yet we show that this can only occur over finite time horizons, with the dark soliton always escaping eventually, contrary again to what is potentially the case with bright solitons. © 2017 Elsevier B.V

Dark solitons near potential and nonlinearity steps

We study dark solitons near potential and nonlinearity steps and combinations thereof, forming rectangular barriers. This setting is relevant to the contexts of atomic Bose-Einstein condensates (where such steps can be realized by using proper external fields) and nonlinear optics (for beam propagation near interfaces separating optical media of different refractive indices). We use perturbation theory to develop an equivalent particle theory, describing the matter-wave or optical soliton dynamics as the motion of a particle in an effective potential. This Newtonian dynamical problem provides information for the soliton statics and dynamics, including scenarios of reflection, transmission, or quasitrapping at such steps. The case of multiple such steps and its connection to barrier potentials is additionally touched upon. The range of validity of the analytical approximation and radiation effects are also investigated. Our analytical predictions are found to be in very good agreement with the corresponding numerical results, where appropriate. © 2016 American Physical Society

Spatiotemporal algebraically localized waveforms for a nonlinear Schrödinger model with gain and loss

We consider the asymptotic behavior of the solutions of a nonlinear Schrödinger (NLS) model incorporating linear and nonlinear gain/loss. First, we describe analytically the dynamical regimes (depending on the gain/loss strengths), for finite-time collapse, decay, and global existence of solutions. Then, for all the above parametric regimes, we use direct numerical simulations to study the dynamics corresponding to algebraically decaying initial data. We identify crucial differences between the dynamics of vanishing initial conditions, and those converging to a finite constant background: in the former (latter) case we find strong (weak) collapse or decay, when the gain/loss parameters are selected from the relevant regimes. One of our main results, is that in all the above regimes, non-vanishing initial data transition through spatiotemporal, algebraically decaying waveforms. While the system is nonintegrable, the evolution of these waveforms is reminiscent to the evolution of the Peregrine rogue wave of the integrable NLS limit. The parametric range of gain and loss for which this phenomenology persists is also touched upon.We would like to thank the referees for valuable comments. The authors Z.A.A., G.F., D.J.F., N.I.K., P.G.K., I.G.S. and K.V., acknowledge that this work wasmade possible by NPRP grant # [ 8-764-160 ] from Qatar National Research Fund (a member of Qatar Foundation). The findings achieved herein are solely the responsibility of the authors. D.J.F. and P.G.K. gratefully acknowledge the support of the "Greek Diaspora Fellowship Program" of Stavros Niarchos Foundation. Appendix AScopu

Spatiotemporal algebraically localized waveforms for a nonlinear Schrödinger model with gain and loss

We consider the asymptotic behavior of the solutions of a nonlinear Schrödinger (NLS) model incorporating linear and nonlinear gain/loss. First, we describe analytically the dynamical regimes (depending on the gain/loss strengths), for finite-time collapse, decay, and global existence of solutions. Then, for all the above parametric regimes, we use direct numerical simulations to study the dynamics corresponding to algebraically decaying initial data. We identify crucial differences between the dynamics of vanishing initial conditions, and those converging to a finite constant background: in the former (latter) case we find strong (weak) collapse or decay, when the gain/loss parameters are selected from the relevant regimes. One of our main results, is that in all the above regimes, non-vanishing initial data transition through spatiotemporal, algebraically decaying waveforms. While the system is nonintegrable, the evolution of these waveforms is reminiscent to the evolution of the Peregrine rogue wave of the integrable NLS limit. The parametric range of gain and loss for which this phenomenology persists is also touched upon. © 2017 Elsevier B.V