Stability of nonlinear Dirac solitons under the action of external potential.

The instabilities observed in direct numerical simulations of the Gross-Neveu equation under linear and harmonic potentials are studied. The Lakoba algorithm, based on the method of characteristics, is performed to numerically obtain the two spinor components. We identify non-conservation of energy and charge in simulations with instabilities, and we find that all studied solitons are numerically stable, except the low-frequency solitons oscillating in the harmonic potential over long periods of time. These instabilities, as in the case of the Gross-Neveu equation without potential, can be removed by imposing absorbing boundary conditions. The dynamics of the soliton is in perfect agreement with the prediction obtained using an Ansatz with only two collective coordinates, namely, the position and momentum of the center of mass. We employ the temporal variation of both field energy and momentum to determine the evolution equations satisfied by the collective coordinates. By applying the same methodology, we also demonstrate the spurious character of the reported instabilities in the Alexeeva-Barashenkov-Saxena model under external potentials.

Exact stationary solutions of the parametrically driven and damped nonlinear Dirac equation.

Two exact stationary soliton solutions are found in the parametrically driven and damped nonlinear Dirac equation. The parametric force considered is a complex ac force. The solutions appear when their frequencies are locked to half the frequency of the parametric force, and their phases satisfy certain conditions depending on the force amplitude and on the damping coefficient. Explicit expressions for the charge, the energy, and the momentum of these solutions are provided. Their stability is studied via a variational method using an ansatz with only two collective coordinates. Numerical simulations confirm that one of the solutions is stable, while the other is an unstable saddle point. Consequently, the stabilization of damped Dirac solitons can be achieved via time-periodic parametric excitations.

Soliton ratchet induced by random transitions among symmetric sine-Gordon potentials.

The generation of net soliton motion induced by random transitions among N symmetric phase-shifted sine-Gordon potentials is investigated, in the absence of any external force and without any thermal noise. The phase shifts of the potentials and the damping coefficients depend on a stationary Markov process. Necessary conditions for the existence of transport are obtained by an exhaustive study of the symmetries of the stochastic system and of the soliton velocity. It is shown that transport is generated by unequal transfer rates among the phase-shifted potentials or by unequal friction coefficients or by a properly devised combination of potentials (N>2). Net motion and inversions of the currents, predicted by the symmetry analysis, are observed in simulations as well as in the solutions of a collective coordinate theory. A model with high efficient soliton motion is designed by using multistate phase-shifted potentials and by breaking the symmetries with unequal transfer rates.

Speed-of-light pulses in the massless nonlinear Dirac equation with a potential.

We consider the massless nonlinear Dirac (NLD) equation in 1+1 dimension with scalar-scalar self-interaction g^{2}/2(Ψ[over ¯]Ψ)^{2} in the presence of three external electromagnetic real potentials V(x), a potential barrier, a constant potential, and a potential well. By solving numerically the NLD equation, we find different scenarios depending on initial conditions, namely, propagation of the initial pulse along one direction, splitting of the initial pulse into two pulses traveling in opposite directions, and focusing of two initial pulses followed by a splitting. For all considered cases, the final waves travel with the speed of light and are solutions of the massless linear Dirac equation. During these processes the charge and the energy are conserved, whereas the momentum is conserved when the solutions possess specific symmetries. For the case of the constant potential, we derive exact analytical solutions of the massless NLD equation that are also solutions of the massless linearized Dirac equation. Decay or growth of the initial pulse is also predicted from the evolution of the charge for the case of a non-zero imaginary part of the potential.

Kink ratchet induced by a time-dependent symmetric field potential.

The ratchet effect of a sine-Gordon kink is investigated in the absence of any external force while the symmetry of the field potential at every time instant is maintained. The directed motion appears by a time shift of the sine-Gordon potential through a time-dependent additional phase. A symmetry analysis provides the necessary conditions for the existence of net motion. It is also shown analytically, by using a collective coordinate theory, that the novel physical mechanism responsible for the appearance of the ratchet effect is the coupled dynamics of the kink width with the background field. Biharmonic and dichotomic periodic variations of the additional phase of the sine-Gordon potential are considered. The predictions established by the symmetry analysis and the collective coordinate theory are verified by means of numerical simulations. Inversion and maximization of the resulting current as a function of the system parameters are investigated.

Erratum: Nonlinear Dirac equation solitary waves in external fields [Phys. Rev. E 86, 046602 (2012)].

This corrects the article DOI: 10.1103/PhysRevE.86.046602.

Soliton dynamics in optical fibers using the generalized traveling-wave method.

The generalized traveling wave method (GTWM) is applied to the nonlinear Schrödinger (NLS) equation with general perturbations in order to obtain the equations of motion for an ansatz with six collective coordinates, namely the soliton position, the amplitude, the inverse of the soliton width, the velocity, the chirp, and the phase. The advantage of the new ansatz is that it yields three pairs of canonically conjugated coordinates and momenta that all are well-behaved. The new ansatz is applied to model the dynamics of a soliton in a dispersion-shifted optical fiber described by the generalized NLS, including dissipation, higher-order dispersion, Raman scattering, and self-steepening perturbations. It is shown that the GTWM is equivalent to the modified method of moments, which considers the time variation of the norm, the first and the second moment of the norm, the momentum, the first moment of the momentum, and the energy for the perturbed NLS equation.

Nonsinusoidal current and current reversals in a gating ratchet.

In this work, the ratchet dynamics of Brownian particles driven by an external sinusoidal (harmonic) force is investigated. The gating ratchet effect is observed when another harmonic is used to modulate the spatially symmetric potential in which the particles move. For small amplitudes of the harmonics, it is shown that the current (average velocity) of particles exhibits a sinusoidal shape as a function of a precise combination of the phases of both harmonics. By increasing the amplitudes of the harmonics beyond the small-limit regime, departures from the sinusoidal behavior are observed and current reversals can also be induced. These current reversals persist even for the overdamped dynamics of the particles.

General approach for dealing with dynamical systems with spatiotemporal periodicities.

Dynamical systems often contain oscillatory forces or depend on periodic potentials. Time or space periodicity is reflected in the properties of these systems through a dependence on the parameters of their periodic terms. In this paper we provide a general theoretical framework for dealing with these kinds of systems, regardless of whether they are classical or quantum, stochastic or deterministic, dissipative or nondissipative, linear or nonlinear, etc. In particular, we are able to show that simple symmetry considerations determine, to a large extent, how their properties depend functionally on some of the parameters of the periodic terms. For the sake of illustration, we apply this formalism to find the functional dependence of the expectation value of the momentum of a Bose-Einstein condensate, described by the Gross-Pitaewskii equation, when it is exposed to a sawtooth potential whose amplitude is periodically modulated in time. We show that, by using this formalism, a small set of measurements is enough to obtain the functional form for a wide range of parameters. This can be very helpful when characterizing experimentally the response of systems for which performing measurements is costly or difficult.

Soliton stability criterion for generalized nonlinear Schrödinger equations.

A stability criterion for solitons of the driven nonlinear Schrödinger equation (NLSE) has been conjectured. The criterion states that p'(v)<0 is a sufficient condition for instability, while p'(v)>0 is a necessary condition for stability; here, v is the soliton velocity and p=P/N, where P and N are the soliton momentum and norm, respectively. To date, the curve p(v) was calculated approximately by a collective coordinate theory, and the criterion was confirmed by simulations. The goal of this paper is to calculate p(v) exactly for several classes and cases of the generalized NLSE: a soliton moving in a real potential, in particular a time-dependent ramp potential, and a time-dependent confining quadratic potential, where the nonlinearity in the NLSE also has a time-dependent coefficient. Moreover, we investigate a logarithmic and a cubic NLSE with a time-independent quadratic potential well. In the latter case, there is a bisoliton solution that consists of two solitons with asymmetric shapes, forming a bound state in which the shapes and the separation distance oscillate. Finally, we consider a cubic NLSE with parametric driving. In all cases, the p(v) curve is calculated either analytically or numerically, and the stability criterion is confirmed.

Collective coordinates theory for discrete soliton ratchets in the sine-Gordon model.

A collective coordinate theory is developed for soliton ratchets in the damped discrete sine-Gordon model driven by a biharmonic force. An ansatz with two collective coordinates, namely the center and the width of the soliton, is assumed as an approximated solution of the discrete nonlinear equation. The dynamical equations of these two collective coordinates, obtained by means of the generalized travelling wave method, explain the mechanism underlying the soliton ratchet and capture qualitatively all the main features of this phenomenon. The numerical simulation of these equations accounts for the existence of a nonzero depinning threshold, the nonsinusoidal behavior of the average velocity as a function of the relative phase between the harmonics of the driver, the nonmonotonic dependence of the average velocity on the damping, and the existence of nontransporting regimes beyond the depinning threshold. In particular, it provides a good description of the intriguing and complex pattern of subspaces corresponding to different dynamical regimes in parameter space.

Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity.

We consider the nonlinear Dirac equation in 1 + 1 dimension with scalar-scalar self interaction g(2)/κ+1(̅ΨΨ)(κ+1) and with mass m. Using the exact analytic form for rest frame solitary waves of the form Ψ(x,t)=ψ(x)e(-iωt) for arbitrary κ, we discuss the validity of various approaches to understanding stability that were successful for the nonlinear Schrödinger equation. In particular we study the validity of a version of Derrick's theorem and the criterion of Bogolubsky as well as the Vakhitov-Kolokolov criterion, and find that these criteria yield inconsistent results. Therefore, we study the stability by numerical simulations using a recently developed fourth-order operator splitting integration method. For different ranges of κ we map out the stability regimes in ω. We find that all stable nonlinear Dirac solitary waves have a one-hump profile, but not all one-hump waves are stable, while all waves with two humps are unstable. We also find that the time t(c), it takes for the instability to set in, is an exponentially increasing function of ω and t(c) decreases monotonically with increasing κ.

Propulsion efficiency of a dynamic self-assembled helical ribbon.

We study the dynamic self-assembly and propulsion of a ribbon formed from paramagnetic colloids in a dynamic magnetic field. The sedimented ribbon assembles due to time averaged dipolar interactions between the beads. The time dependence of the dipolar interactions together with hydrodynamic interactions cause a twisted ribbon conformation. Domain walls of high twist connect domains of nearly constant orientation and negligible twist and travel through the ribbon. The particular form of the domain walls can be controlled via the frequency and the eccentricity of the modulation. The flux of twist walls-a true ribbon property absent in slender bodies-provides the thrust onto the surrounding liquid that propels this biomimetic flagellum into the opposite direction. The propulsion efficiency increases with frequency and ceases abruptly at a critical frequency where the conformation changes discontinuously to a flat standing ribbon conformation.

Comment on "Ratchet universality in the presence of thermal noise".

A recent paper [P. J. Martínez and R. Chacón, Phys. Rev. E 87, 062114 (2013)] presents numerical simulations on a system exhibiting directed ratchet transport of a driven overdamped Brownian particle subjected to a spatially periodic, symmetric potential. The authors claim that their simulations prove the existence of a universal waveform of the external force that optimally enhances directed transport, hence confirming the validity of a previous conjecture put forth by one of them in the limit of vanishing noise intensity. With minor corrections due to noise, the conjecture holds even in the presence of noise, according to the authors. On the basis of their results the authors claim that all previous theories, which predict a different optimal force waveform, are incorrect. In this Comment we provide sufficient numerical evidence showing that there is no such universal force waveform and that the evidence obtained by the authors otherwise is due to their particular choice of parameters. Our simulations also suggest that previous theories correctly predict the shape of the optimal waveform within their validity regime, namely, when the forcing is weak. On the contrary, the aforementioned conjecture does not hold.

Nonlinear Dirac equation solitary waves in external fields.

We consider nonlinear Dirac equations (NLDE's) in the 1+1 dimension with scalar-scalar self-interaction g2/κ+1(Ψ[over ¯]Ψ)κ+1 in the presence of various external electromagnetic fields. We find exact solutions for special external fields and we study the behavior of solitary-wave solutions to the NLDE in the presence of a wide variety of fields in a variational approximation depending on collective coordinates which allows the position, width, and phase of these waves to vary in time. We find that in this approximation the position q(t) of the center of the solitary wave obeys the usual behavior of a relativistic point particle in an external field. For time-independent external fields, we find that the energy of the solitary wave is conserved but not the momentum, which becomes a function of time. We postulate that, similarly to the nonlinear Schrödinger equation (NLSE), a sufficient dynamical condition for instability to arise is that dP(t)/dq[over ̇](t)<0. Here P(t) is the momentum of the solitary wave, and q[over ̇] is the velocity of the center of the wave in the collective coordinate approximation. We found for our choices of external potentials that we always have dP(t)/dq[over ̇](t)>0, so, when instabilities do occur, they are due to a different source. We investigate the accuracy of our variational approximation using numerical simulations of the NLDE and find that, when the forcing term is small and we are in a regime where the solitary wave is stable, that the behavior of the solutions of the collective coordinate equations agrees very well with the numerical simulations. We found that the time evolution of the collective coordinates of the solitary wave in our numerical simulations, namely the position of the average charge density and the momentum of the solitary wave, provide good indicators for when the solitary wave first becomes unstable. When these variables stop being smooth functions of time (t), then the solitary wave starts to distort in shape.

Forced nonlinear Schrödinger equation with arbitrary nonlinearity.

We consider the nonlinear Schrödinger equation (NLSE) in 1+1 dimension with scalar-scalar self-interaction g(2)/κ+1(ψ*ψ)(κ+1) in the presence of the external forcing terms of the form re(-i(kx+θ))-δψ. We find new exact solutions for this problem and show that the solitary wave momentum is conserved in a moving frame where v(k)=2k. These new exact solutions reduce to the constant phase solutions of the unforced problem when r→0. In particular we study the behavior of solitary wave solutions in the presence of these external forces in a variational approximation which allows the position, momentum, width, and phase of these waves to vary in time. We show that the stationary solutions of the variational equations include a solution close to the exact one and we study small oscillations around all the stationary solutions. We postulate that the dynamical condition for instability is that dp(t)/dq ̇(t)<0, where p(t) is the normalized canonical momentum p(t)=1/M(t)∂L/∂q ̇, and q ̇(t) is the solitary wave velocity. Here M(t)=∫dxψ*(x,t)ψ(x,t). Stability is also studied using a "phase portrait" of the soliton, where its dynamics is represented by two-dimensional projections of its trajectory in the four-dimensional space of collective coordinates. The criterion for stability of a soliton is that its trajectory is a closed single curve with a positive sense of rotation around a fixed point. We investigate the accuracy of our variational approximation and these criteria using numerical simulations of the NLSE. We find that our criteria work quite well when the magnitude of the forcing term is small compared to the amplitude of the unforced solitary wave. In this regime the variational approximation captures quite well the behavior of the solitary wave.

Refined empirical stability criterion for nonlinear Schrödinger solitons under spatiotemporal forcing.

We investigate the dynamics of traveling oscillating solitons of the cubic nonlinear Schrödinger (NLS) equation under an external spatiotemporal forcing of the form f(x,t)=aexp[iK(t)x]. For the case of time-independent forcing, a stability criterion for these solitons, which is based on a collective coordinate theory, was recently conjectured. We show that the proposed criterion has a limited applicability and present a refined criterion which is generally applicable, as confirmed by direct simulations. This includes more general situations where K(t) is harmonic or biharmonic, with or without a damping term in the NLS equation. The refined criterion states that the soliton will be unstable if the "stability curve" p(v), where p(t) and v(t) are the normalized momentum and the velocity of the soliton, has a section with a negative slope. In the case of a constant K and zero damping, we use the collective coordinate solutions to compute a "phase portrait" of the soliton where its dynamics is represented by two-dimensional projections of its trajectories in the four-dimensional space of collective coordinates. We conjecture, and confirm by simulations, that the soliton is unstable if a section of the resulting closed curve on the portrait has a negative sense of rotation.

Generalized traveling-wave method, variational approach, and modified conserved quantities for the perturbed nonlinear Schrödinger equation.

The generalized traveling wave method (GTWM) is developed for the nonlinear Schrödinger equation (NLSE) with general perturbations in order to obtain the equations of motion for an arbitrary number of collective coordinates. Regardless of the particular ansatz that is used, it is shown that this alternative approach is equivalent to the Lagrangian formalism, but has the advantage that only the Hamiltonian of the unperturbed system is required, instead of the Lagrangian for the perturbed system. As an explicit example, we take 4 collective coordinates, namely the position, velocity, amplitude and phase of the soliton, and show that the GTWM yields the same equations of motion as the perturbation theory based on the Inverse Scattering Transform and as the time variation of the norm, first moment of the norm, momentum, and energy for the perturbed NLSE.

Nonlinear Schrödinger equation with spatiotemporal perturbations.

We investigate the dynamics of solitons of the cubic nonlinear Schrödinger equation (NLSE) with the following perturbations: nonparametric spatiotemporal driving of the form f(x,t)=a exp[iK(t)x], damping, and a linear term which serves to stabilize the driven soliton. Using the time evolution of norm, momentum and energy, or, alternatively, a Lagrangian approach, we develop a collective-coordinate-theory which yields a set of ordinary differential equations (ODEs) for our four collective coordinates. These ODEs are solved analytically and numerically for the case of a constant, spatially periodic force f(x). The soliton position exhibits oscillations around a mean trajectory with constant velocity. This means that the soliton performs, on the average, a unidirectional motion although the spatial average of the force vanishes. The amplitude of the oscillations is much smaller than the period of f(x). In order to find out for which regions the above solutions are stable, we calculate the time evolution of the soliton momentum P(t) and the soliton velocity V(t): This is a parameter representation of a curve P(V) which is visited by the soliton while time evolves. Our conjecture is that the soliton becomes unstable, if this curve has a branch with negative slope. This conjecture is fully confirmed by our simulations for the perturbed NLSE. Moreover, this curve also yields a good estimate for the soliton lifetime: the soliton lives longer, the shorter the branch with negative slope is.

Symmetries shape the current in ratchets induced by a biharmonic driving force.

Equations describing the evolution of particles, solitons, or localized structures, driven by a zero-average, periodic, external force, and invariant under time reversal and a half-period time shift, exhibit a ratchet current when the driving force breaks these symmetries. The biharmonic force f(t)=1 cos(qomegat+phi1)+2 cospomegat+phi2) does it for almost any choice of vphi1 and phi2, provided p and q are two coprime integers such that p+q is odd. It has been widely observed, in experiments in semiconductors, in Josephson junctions, photonic crystals, etc., as well as in simulations, that the ratchet current induced by this force has the shape v proportional, variant1p2q cos(pphi1-qphi2+theta0) for small amplitudes, where theta0 depends on the damping ( theta0=pi/2 if there is no damping, and theta0=0 for overdamped systems). We rigorously prove that this precise shape can be obtained solely from the broken symmetries of the system and is independent of the details of the equation describing the system.