Necklacelike solitons in optically induced photonic lattices

We report the first observation of stationary necklacelike solitons. Such necklace structures were realized when a high-order vortex beam was launched appropriately into a two-dimensional optically induced photonic lattice. Our theoretical results obtained with continuous and discrete models show that the necklace solitons resulting from a charge-4 vortex have a pi phase difference between adjacent "pearls" and are formed in an octagon shape. Their stability region is identified.

Parametric light bullets supported by quasi-phase-matched quadratically nonlinear crystals

We present a comprehensive analysis of the dynamics of three-dimensional spatiotemporal nonspinning and spinning solitons in quasi-phased-matched (QPM) gratings. By employing an averaging approach based on perturbation theory, we show that the soliton's stability is strongly affected by the QPM-induced third-order nonlinearity (which is always of a mixed type, with opposite signs in front of the corresponding self-phase and cross-phase modulation terms). We study the dependence of the stability of the spatiotemporal soliton (STS) on its energy, spin, the wave-vector mismatch between the fundamental and second harmonics, and the relative strength of the intrinsic quadratic and QPM-induced cubic nonlinearities. In particular, all the spinning solitons are unstable against fragmentation, while zero-spin STS's have their stability regions on the system's parameter space.

Two-dimensional solitons with hidden and explicit vorticity in bimodal cubic-quintic media

We demonstrate that two-dimensional two-component bright solitons of an annular shape, carrying vorticities (m,+/-m) in the components, may be stable in media with the cubic-quintic nonlinearity, including the hidden-vorticity (HV) solitons of the type (m,-m) , whose net vorticity is zero. Stability regions for the vortices of both (m,+/-m) types are identified for m=1 , 2, and 3, by dint of the calculation of stability eigenvalues, and in direct simulations. In addition to the well-known symmetry-breaking (external) instability, which splits the ring soliton into a set of fragments flying away in tangential directions, we report two new scenarios of the development of weak instabilities specific to the HV solitons. One features charge flipping, with the two components exchanging angular momentum and periodically reversing the sign of their spins. The composite soliton does not directly split in this case; therefore, we identify such instability as an intrinsic one. Eventually, the soliton splits, as weak radiation loss drives it across the border of the ordinary strong (external) instability. Another scenario proceeds through separation of the vortex cores in the two components, each individual core moving toward the outer edge of the annular soliton. After expulsion of the cores, there remains a zero-vorticity breather with persistent internal vibrations.

Controlling the motion of dark solitons by means of periodic potentials: application to Bose-Einstein condensates in optical lattices

We demonstrate that the motion of dark solitons (DSs) can be controlled by means of periodic potentials. The mechanism is realized in terms of cigar-shaped Bose-Einstein condensates confined in a harmonic magnetic potential, in the presence of an optical lattice (OL). In the case when the OL period is comparable to the width of the DS, we demonstrate that (a) a moving dark soliton can be captured, switching on the OL, and (b) a stationary DS can be dragged by a moving OL.

Nonlinearity from geometric interactions: a case example

We propose a ladder model wherein dynamical nonlinearity arises from geometry. It includes two strings of particles which are set along rigid rails of a "railroad" and coupled by linear springs. Physical realizations of the model include dust-particle strings in plasma sheaths and chains of microparticles trapped in a strong optical lattice. The transverse couplings between the strings, along with the motion constraint imposed by the rails, generate nonlinearity. It gives rise to robust solitary waves, which are found analytically in the long-wavelength limit, and are obtained in simulations of the full system.

Three-dimensional solitary waves and vortices in a discrete nonlinear Schrödinger lattice

In a benchmark dynamical-lattice model in three dimensions, the discrete nonlinear Schrödinger equation, we find discrete vortex solitons with various values of the topological charge S. Stability regions for the vortices with S=0,1,3 are investigated. The S=2 vortex is unstable and may spontaneously rearranging into a stable one with S=3. In a two-component extension of the model, we find a novel class of stable structures, consisting of vortices in the different components, perpendicularly oriented to each other. Self-localized states of the proposed types can be observed experimentally in Bose-Einstein condensates trapped in optical lattices and in photonic crystals built of microresonators.

Two-dimensional dispersion-managed light bullets in Kerr media

We propose a scheme for stabilizing spatiotemporal solitons (STSs) in media with cubic self-focusing nonlinearity and "dispersion management," i.e., a layered structure inducing periodically alternating normal and anomalous group-velocity dispersion. We develop a variational approximation for the STS, and verify results by direct simulations. A stability region for the two-dimensional (2D) STS (corresponding to a planar waveguide) is identified. At the borders between this region and that of decay of the solitons, a more sophisticated stable object, in the form of a periodically oscillating bound state of two subpulses, is also found. In the 3D case (bulk medium), all the spatiotemporal pulses spread out or collapse.

Stable vortex solitons supported by competing quadratic and cubic nonlinearities

We address the stability problem for vortex solitons in two-dimensional media combining quadratic and self-defocusing cubic [chi(2):chi(3)- ] nonlinearities. We consider the propagation of spatial beams with intrinsic vorticity S in such bulk optical media. It was earlier found that the S=1 and S=2 solitons can be stable, provided that their power (i.e., transverse size) is large enough, and it was conjectured that all the higher-order vortices with S> or =3 are always unstable. On the other hand, it was recently shown that vortex solitons with S>2 and very large transverse size may be stable in media combining cubic self-focusing and quintic self-defocusing nonlinearities. Here, we demonstrate that the same is true in the chi(2):chi(3)- model, the vortices with S=3 and S=4 being stable in regions occupying, respectively, approximately 3% and 1.5% of their existence domain. The vortex solitons with S>4 are also stable in tiny regions. The results are obtained through computation of stability eigenvalues, and are then checked in direct simulations, with a conclusion that the stable vortices are truly robust ones, easily self-trapping from initial beams with embedded vorticity. The dependence of the stability region on the chi(2) phase-mismatch parameter is specially investigated. We thus conclude that the stability of higher-order two-dimensional vortex solitons in narrow regions is a generic feature of optical media featuring the competition between self-focusing and self-defocusing nonlinearities. A qualitative analytical explanation to this feature is proposed.

Quantum dissociation of a vortex-antivortex pair in a long josephson junction

The thermal and the quantum dissociation of a single vortex-antivortex (VAV) pair in an annular Josephson junction is experimentally observed and theoretically analyzed. In our experiments, the VAV pair is confined in a pinning potential controlled by external magnetic field and bias current. The dissociation of the pinned VAV pair manifests itself in a switching of the Josephson junction from the superconducting to the resistive state. The observed temperature and field dependence of the switching current distribution is in agreement with the analysis. The crossover from the thermal to the macroscopic quantum tunneling mechanism of dissociation occurs at a temperature of about 100 mK. We also predict the specific magnetic field dependence of the oscillatory energy levels of the pinned VAV state.

Temporal solitons in quadratic nonlinear media with opposite group-velocity dispersions at the fundamental and second harmonics

Temporal solitons in quadratic nonlinear media with normal second-harmonic dispersion are studied theoretically. The variational approximation and direct simulations reveal the existence of soliton solutions, and their stability region is identified. Stable solutions are found for large and normal values of the second-harmonic dispersion, and in the presence of large group-velocity mismatch between the fundamental- and second-harmonic fields. The solitons (or solitonlike pulses) are found to have tiny nonlocalized tails in the second-harmonic field, for which an analytic exponential estimate is obtained. The estimate and numerical calculations show that, in the parameter region of experimental relevance, the tails are completely negligible. The results open a way to the experimental observation of quadratic solitons with normal second-harmonic dispersion, and have strong implication to the experimental search for multidimensional "light bullets."

Two-soliton collisions in a near-integrable lattice system

We examine collisions between identical solitons in a weakly perturbed Ablowitz-Ladik (AL) model, augmented by either onsite cubic nonlinearity (which corresponds to the Salerno model, and may be realized as an array of strongly overlapping nonlinear optical waveguides) or a quintic perturbation, or both. Complex dependences of the outcomes of the collisions on the initial phase difference between the solitons and location of the collision point are observed. Large changes of amplitudes and velocities of the colliding solitons are generated by weak perturbations, showing that the elasticity of soliton collisions in the AL model is fragile (for instance, the Salerno's perturbation with the relative strength of 0.08 can give rise to a change of the solitons' amplitudes by a factor exceeding 2). Exact and approximate conservation laws in the perturbed system are examined, with a conclusion that the small perturbations very weakly affect the norm and energy conservation, but completely destroy the conservation of the lattice momentum, which is explained by the absence of the translational symmetry in generic nonintegrable lattice models. Data collected for a very large number of collisions correlate with this conclusion. Asymmetry of the collisions (which is explained by the dependence on the location of the central point of the collision relative to the lattice, and on the phase difference between the solitons) is investigated too, showing that the nonintegrability-induced effects grow almost linearly with the perturbation strength. Different perturbations (cubic and quintic ones) produce virtually identical collision-induced effects, which makes it possible to compensate them, thus finding a special perturbed system with almost elastic soliton collisions.

Robust soliton clusters in media with competing cubic and quintic nonlinearities

Systematic results are reported for dynamics of circular patterns ("necklaces"), composed of fundamental solitons and carrying orbital angular momentum, in the two-dimensional model, which describes the propagation of light beams in bulk media combining self-focusing cubic and self-defocusing quintic nonlinearities. Semianalytical predictions for the existence of quasistable necklace structures are obtained on the basis of an effective interaction potential. Then, direct simulations are run. In the case when the initial pattern is far from an equilibrium size predicted by the potential, it cannot maintain its shape. However, a necklace with the initial shape close to the predicted equilibrium survives very long evolution, featuring persistent oscillations. The quasistable evolution is not essentially disturbed by a large noise component added to the initial configuration. Basic conclusions concerning the necklace dynamics in this model are qualitatively the same as in a recently studied one which combines quadratic and self-defocusing cubic nonlinearities. Thus, we infer that a combination of competing self-focusing and self-defocusing nonlinearities enhances the robustness not only of vortex solitons but also of vorticity-carrying necklace patterns.

Soliton collisions in the discrete nonlinear Schrödinger equation

We report analytical and numerical results for on-site and intersite collisions between solitons in the discrete nonlinear Schrödinger model. A semianalytical variational approximation correctly predicts gross features of the collision, viz., merger or bounce. We systematically examine the dependence of the collision outcome on initial velocity and amplitude of the solitons, as well as on the phase shift between them, and location of the collision point relative to the lattice; in some cases, the dependences are very intricate. In particular, merger of the solitons into a single one, and bounce after multiple collisions are found. Situations with a complicated system of alternating transmission and merger windows are identified too. The merger is often followed by symmetry breaking (SB), when the single soliton moves to the left or to the right, which implies momentum nonconservation. Two different types of the SB are identified, deterministic and spontaneous. The former one is accounted for by the location of the collision point relative to the lattice, and/or the phase shift between the solitons; the momentum generated during the collision due to the phase shift is calculated in an analytical approximation, its dependence on the solitons' velocities comparing well with numerical results. The spontaneous SB is explained by the modulational instability of a quasiflat plateau temporarily formed in the course of the collision.

Standard and embedded solitons in nematic optical fibers

A model for a non-Kerr cylindrical nematic fiber is presented. We use the multiple scales method to show the possibility of constructing different kinds of wave packets of transverse magnetic modes propagating through the fiber. This procedure allows us to generate different hierarchies of nonlinear partial differential equations which describe the propagation of optical pulses along the fiber. We go beyond the usual weakly nonlinear limit of a Kerr medium and derive a complex modified Korteweg-de Vries equation (CM KdV) which governs the dynamics for the amplitude of the wave packet. In this derivation the dispersion, self-focussing, and diffraction in the nematic fiber are taken into account. It is shown that this CM KdV equation has two-parameter families of bright and dark complex solitons. We show analytically that under certain conditions, the bright solitons are actually double-embedded solitons. We explain why these solitons do not radiate at all, even though their wave numbers are contained in the linear spectrum of the system. We study (numerically and analytically) the stability of these solitons. Our results show that these embedded solitons are stable solutions, which is an interesting property since in most systems the embedded solitons are weakly unstable solutions. Finally, we close the paper by making comments on the advantages as well as the limitations of our approach, and on further generalizations of the model and method presented.

Stable spatiotemporal spinning solitons in a bimodal cubic-quintic medium

We investigate the formation of stable spatiotemporal three-dimensional (3D) solitons ("light bullets") with internal vorticity ("spin") in a bimodal system described by coupled cubic-quintic nonlinear Schrödinger equations. Two relevant versions of the model, for the linear and circular polarizations, are considered. In the former case, an important ingredient of the model are four-wave-mixing terms, which give rise to a phase-sensitive nonlinear coupling between two polarization components. Thresholds for the formation of both spinning and nonspinning 3D solitons are found. Instability growth rates of perturbation eigenmodes with different azimuthal indices are calculated as functions of the solitons' propagation constant. As a result, stability domains in the model's parameter plane are identified for solitons with the values of the spins of their components s=0 and s=1, while all the solitons with s> or =2 are unstable. The solitons with s=1 are stable only if their energy exceeds a certain critical value, so that, in typical cases, the stability region occupies approximately 25% of their existence domain. Direct simulations of the full system produce results that are in perfect agreement with the linear-stability analysis: stable 3D spinning solitons readily self-trap from initial Gaussian pulses with embedded vorticity, and easily heal themselves if strong perturbations are imposed, while unstable spinning solitons quickly split into a set of separating zero-spin fragments whose number is exactly equal to the azimuthal index of the strongest unstable perturbation eigenmode.

Discrete vector solitons in two-dimensional nonlinear waveguide arrays: solutions, stability, and dynamics

We identify and investigate bimodal (vector) solitons in models of square-lattice arrays of nonlinear optical waveguides. These vector self-localized states are, in fact, self-induced channels in a nonlinear photonic-crystal matrix. Such two-dimensional discrete vector solitons are possible in waveguide arrays in which each element carries two light beams that are either orthogonally polarized or have different carrier wavelengths. Estimates of the physical parameters necessary to support such soliton solutions in waveguide arrays are given. Using Newton relaxation methods, we obtain stationary vector-soliton solutions, and examine their stability through the computation of linearized eigenvalues for small perturbations. Our results may also be applicable to other systems such as two-component Bose-Einstein condensates trapped in a two-dimensional optical lattice.

Coupling fields and underlying space curvature: an augmented Lagrangian approach

We demonstrate a systematic implementation of coupling between a scalar field and the geometry of the space which carries the field. This naturally gives rise to a feedback mechanism between the field and the geometry. We develop a systematic model for the feedback in a general form, inspired by a specific implementation in the context of molecular dynamics (the so-called Rahman-Parrinello molecular dynamics, or RP-MD). We use a generalized Lagrangian that allows for the coupling of the space's metric tensor to the scalar field, and add terms motivated by RP-MD. We present two implementations of the scheme: one in which the metric is only time-dependent (which gives rise to an ordinary differential equation for its temporal evolution), and the other with spatiotemporal dependence (wherein the metric's evolution is governed by a partial differential equation). Numerical results are reported for the (1+1)-dimensional model with a nonlinearity of the sine-Gordon type.

Ring dark solitons and vortex necklaces in Bose-Einstein condensates

We introduce the concept of ring dark solitons in Bose-Einstein condensates. We show that relatively shallow rings are not subject to the snake instability, but a deeper ring splits into a robust ringlike cluster of vortex pairs, which performs oscillations in the radial and azimuthal directions, following the dynamics of the original ring soliton.

Light bullets in quadratic media with normal dispersion at the second harmonic

Stable two- and three-dimensional spatiotemporal solitons (STSs) in second-harmonic-generating media are found in the case of normal dispersion at the second-harmonic (SH). This result, surprising from the theoretical viewpoint, opens the way for experimental realization of STSs. An analytical estimate for the existence of STSs is derived, and full results, including a complete stability diagram, are obtained in numerical form. STSs withstand not only the normal SH dispersion, but also finite walk-off between the harmonics, and readily self-trap from a Gaussian pulse launched at the fundamental frequency.

Discrete nonlinear model with substrate feedback

We consider a prototypical model in which a nonlinear field (continuum or discrete) evolves on a flexible substrate which feeds back to the evolution of the main field. We identify the underlying physics and potential applications of such a model and examine its simplest one-dimensional Hamiltonian form, which turns out to be a modified Frenkel-Kontorova model coupled to an extra linear equation. We find static kink solutions and study their stability, and then examine moving kinks (the continuum limit of the model is studied too). We observe how the substrate effectively renormalizes properties of the kinks. In particular, a nontrivial finding is that branches of stable and unstable kink solutions may be extended beyond a critical point at which an effective intersite coupling vanishes; passing this critical point does not destabilize the kink. Kink-antikink collisions are also studied, demonstrating alternation between merger and transmission cases.

Polychromatic solitons in a quadratic medium

We introduce the simplest model to describe parametric interactions in a quadratically nonlinear optical medium with the fundamental harmonic containing two components with (slightly) different carrier frequencies [which is a direct analog of wavelength-division multiplexed models, well known in media with cubic nonlinearity]. The model takes a closed form with three different second-harmonic components, and it is formulated in the spatial domain. We demonstrate that the model supports both polychromatic solitons (PCSs), with all the components present in them, and two types of mutually orthogonal simple solitons, both types being stable in a broad parametric region. An essential peculiarity of PCS is that its power is much smaller than that of a simple (usual) soliton (taken at the same values of control parameters), which may be an advantage for experimental generation of PCSs. Collisions between the orthogonal simple solitons are simulated in detail, leading to the conclusion that the collisions are strongly inelastic, converting the simple solitons into polychromatic ones, and generating one or two additional PCSs. A collision velocity at which the inelastic effects are strongest is identified, and it is demonstrated that the collision may be used as a basis to design a simple all-optical XOR logic gate.

Stable three-dimensional spinning optical solitons supported by competing quadratic and cubic nonlinearities

We show that the quadratic interaction of fundamental and second harmonics in a bulk dispersive medium, combined with self-defocusing cubic nonlinearity, gives rise to completely localized spatiotemporal solitons (vortex tori) with vorticity s=1. There is no threshold necessary for the existence of these solitons. They are found to be stable if their energy exceeds a certain critical value, so that the stability domain occupies about 10% of the existence region of the solitons. On the contrary to spatial vortex solitons in the same model, the spatiotemporal ones with s=2 are never stable. These results might open the way for experimental observation of spinning three-dimensional solitons in optical media.

Solitons in triangular and honeycomb dynamical lattices with the cubic nonlinearity

We study the existence and stability of localized states in the discrete nonlinear Schrödinger equation on two-dimensional nonsquare lattices. The model includes both the nearest-neighbor and long-range interactions. For the fundamental strongly localized soliton, the results depend on the coordination number, i.e., on the particular type of lattice. The long-range interactions additionally destabilize the discrete soliton, or make it more stable, if the sign of the interaction is, respectively, the same as or opposite to the sign of the short-range interaction. We also explore more complicated solutions, such as twisted localized modes and solutions carrying multiple topological charge (vortices) that are specific to the triangular and honeycomb lattices. In the cases when such vortices are unstable, direct simulations demonstrate that they typically turn into zero-vorticity fundamental solitons.

Targeted transfer of solitons in continua and lattices

We propose a robust mechanism of targeted energy transfer along a line, as well as on a surface, in the form of transport of coherent solitary-wave structures, driven by a moving, spatially localized external ac field ("arm") in a lossy medium. The efficiency and robustness of the mechanism are demonstrated analytically and numerically in terms of the nonlinear Schrödinger (NLS) equation, and broad regions of stable operation are identified in the model's parameter space. Direct simulations show that the driving arm can manipulate solitons equally well in a lattice NLS model. A salient feature, which is revealed by simulations and explained analytically, is a resonant character of the operation of the driving arm in the lattice medium, both integer and fractional resonances being identified. Numerical experiments also demonstrate that the same solitary-wave-transport mechanism works well in two-dimensional lattice media.

Multichannel pulse dynamics in a stabilized Ginzburg-Landau system

We study the stability and interactions of chirped solitary pulses in a system of nonlinearly coupled cubic Ginzburg-Landau (CGL) equations with a group-velocity mismatch between them, where each CGL equation is stabilized by linearly coupling it to an additional linear dissipative equation. In the context of nonlinear fiber optics, the model describes transmission and collisions of pulses at different wavelengths in a dual-core fiber, in which the active core is furnished with bandwidth-limited gain, while the other, passive (lossy) one is necessary for stabilization of the solitary pulses. Complete and incomplete collisions of pulses in two channels in the cases of anomalous and normal dispersion in the active core are analyzed by means of perturbation theory and direct numerical simulations. It is demonstrated that the model may readily support fully stable pulses whose collisions are quasielastic, provided that the group-velocity difference between the two channels exceeds a critical value. In the case of quasielastic collisions, the temporal shift of pulses, predicted by the analytical approach, is in semiquantitative agreement with direct numerical results in the case of anomalous dispersion (in the opposite case, the perturbation theory does not apply). We also consider a simultaneous collision between pulses in three channels, concluding that this collision remains quasielastic, and the pulses remain completely stable. Thus, the model may be a starting point for the design of a stabilized wavelength-division-multiplexed transmission system.