The dynamics of a pair of harmonic oscillators represented by three-dimensional fields coupled with a repulsive cubic nonlinearity is investigated through direct simulations of the respective field equations and with the help of the finite-mode Galerkin approximation (GA), which represents the two interacting fields by a superposition of 3+3 harmonic-oscillator p-wave eigenfunctions with orbital and magnetic quantum numbers l=1 and m=1, 0, -1. The system can be implemented in binary Bose-Einstein condensates, demonstrating the potential of the atomic condensates to emulate various complex modes predicted by classical field theories. First, the GA very accurately predicts a broadly degenerate set of the system's ground states in the p-wave manifold, in the form of complexes built of a dipole coaxial with another dipole or vortex, as well as complexes built of mutually orthogonal dipoles. Next, pairs of noncoaxial vortices and/or dipoles, including pairs of mutually perpendicular vortices, develop remarkably stable dynamical regimes, which feature periodic exchange of the angular momentum and periodic switching between dipoles and vortices. For a moderately strong nonlinearity, simulations of the coupled-field equations agree very well with results produced by the GA, demonstrating that the dynamics is accurately spanned by the set of six modes limited to l=1.
The dynamics of a pair of harmonic oscillators represented by three-dimensional fields coupled with a repulsive cubic nonlinearity is investigated through direct simulations of the respective field equations and with the help of the finite-mode Galerkin approximation (GA), which represents the two interacting fields by a superposition of 3+3 harmonic-oscillator p-wave eigenfunctions with orbital and magnetic quantum numbers l=1 and m=1, 0, -1. The system can be implemented in binary Bose-Einstein condensates, demonstrating the potential of the atomic condensates to emulate various complex modes predicted by classical field theories. First, the GA very accurately predicts a broadly degenerate set of the system's ground states in the p-wave manifold, in the form of complexes built of a dipole coaxial with another dipole or vortex, as well as complexes built of mutually orthogonal dipoles. Next, pairs of noncoaxial vortices and/or dipoles, including pairs of mutually perpendicular vortices, develop remarkably stable dynamical regimes, which feature periodic exchange of the angular momentum and periodic switching between dipoles and vortices. For a moderately strong nonlinearity, simulations of the coupled-field equations agree very well with results produced by the GA, demonstrating that the dynamics is accurately spanned by the set of six modes limited to l=1.
Nonlinearity is the driving force for numerous important effects in nature typically showing transitions between different regimes, regular, chaotic or catastrophic behavior. Localized nonlinear modes have been the focus of intense research in areas such as fluid and gas dynamics, photonics, atomic and solid state physics etc. Due to the richness of the behavior of nonlinear systems and due to the severe numerical demands of accurate three-dimensional (3D) numerical simulations presently only little knowledge is available on the dynamics of complex nonlinear modes in 3D. Here, we investigate the dynamics of 3D non-coaxial matter wave vortices that are trapped in a parabolic potential and interact via a repulsive nonlinearity. Our numerical simulations demonstrate the existence of an unexpected and fascinating nonlinear regime that starts immediately when the nonlinearity is switched-on and is characterized by a smooth dynamics representing torque-free precession with nutations. The reported motion is proven to be robust regarding various effects such as the number of particles, dissipation and trap deformations and thus should be observable in suitably designed experiments. Since our theoretical approach, i.e., coupled nonlinear Schrödinger equations, is quite generic, we expect that the obtained novel dynamical behavior should also exist in other nonlinear systems.
Bose-Einstein condensate(BEC) provides a nice stage when the nonlinear Schrödinger equation plays a vital role. We study the dynamics of multi-component repulsive BEC in 2 dimensions with harmonic traps by using the nonlinear Schrödinger (or Gross-Pitaevskii) equation. Firstly we consider a driven two-component BEC with each component trapped in different vertical positions. The appropriate tuning of the oscillation frequency of the magnetic field leads to a striking anti-gravity transport of BEC. This phenomenon is a manifestation of macroscopic non-adiabatic tunneling in a system with two internal(electronic) degrees of freedom. The dynamics splits into a fast complex spatio-temporal oscillation of each condensate wavefunctions together with a slow levitation of the total center of mass. Secondly, we examine the three-component repulsive BEC in 2 dimensions in a harmonic trap in the absence of magnetic field, and construct a model of conservative chaos based on a picture of vortex molecules. We obtain an effective nonlinear dynamics for three vortex cores, which represents three charged particles under the uniform magnetic field with the repulsive inter-particle potential quadratic in the inter-vortex distance rij on short scale and logarithmic in rij on large scale. The vortices here acquire the inertia in marked contrast to the standard theory of point vortices since Onsager. We then explore “the chaos in the three-body problem” in the context of vortices with inertia.
We consider a two-component linearly coupled system with the intrinsic cubic nonlinearity and the harmonic-oscillator (HO) confining potential. The system models binary settings in BEC and optics. In the symmetric system, with the HO trap acting in both components, we consider Josephson oscillations (JO) initiated by an input in the form of the HO’s ground state (GS) or dipole mode (DM), placed in one component. With the increase of the strength of the self-focusing nonlinearity, spontaneous symmetry breaking (SSB) between the components takes place in the dynamical JO state. Under still stronger nonlinearity, the regular JO initiated by the GS input carries over into a chaotic dynamical state. For the DM input, the chaotization happens at smaller powers than for the GS, which is followed by SSB at a slightly stronger nonlinearity. In the system with the defocusing nonlinearity, SSB does not take place, and dynamical chaos occurs in a small area of the parameter space. In the asymmetric half-trapped system, with the HO potential applied to a single component, we first focus on the spectrum of confined binary modes in the linearized system. The spectrum is found analytically in the limits of weak and strong inter-component coupling, and numerically in the general case. Under the action of the coupling, the existence region of the confined modes shrinks for GSs and expands for DMs. In the full nonlinear system, the existence region for confined modes is identified in the numerical form. They are constructed too by means of the Thomas–Fermi approximation, in the case of the defocusing nonlinearity. Lastly, particular (non-generic) exact analytical solutions for confined modes, including vortices, in one- and two-dimensional asymmetric linearized systems are found. They represent bound states in the continuum.
We consider a two-component linearly-coupled system with the intrinsic cubic nonlinearity and the harmonic-oscillator (HO) confining potential. The system models binary settings in BEC and optics. In the symmetric system, with the HO trap acting in both components, we consider Josephson oscillations (JO) initiated by an input in the form of the HO's ground state (GS) or dipole mode (DM), placed in one component. With the increase of the strength of the self-focusing nonlinearity, spontaneous symmetry breaking (SSB) between the components takes place in the dynamical JO state. Under still stronger nonlinearity, the regular JO initiated by the GS input carry over into a chaotic dynamical state. For the DM input, the chaotization happens at smaller powers than for the GS, which is followed by SSB at a slightly stronger nonlinearity. In the system with the defocusing nonlinearity, SSB does not take place, and dynamical chaos occurs in a small area of the parameter space. In the asymmetric half-trapped system, with the HO potential applied to a single component, we first focus on the spectrum of confined binary modes in the linearized system. The spectrum is found analytically in the limits of weak and strong inter-component coupling, and numerically in the general case. Under the action of the coupling, the existence region of the confined modes shrinks for GSs and expands for DMs. In the full nonlinear system, the existence region for confined modes is identified in the numerical form. They are constructed too by means of the Thomas-Fermi approximation, in the case of the defocusing nonlinearity. Lastly, particular (non-generic) exact analytical solutions for confined modes, including vortices, in one- and two-dimensional asymmetric linearized systems are found. They represent bound states in the continuum.
With use of the nonlinear Schr\"odinger (or Gross-Pitaevskii) equation with strong repulsive cubic nonlinearity, dynamics of multicomponent Bose-Einstein condensates (BECs) with a harmonic trap in two dimensions is investigated beyond the Thomas-Fermi regime. In the case when each component has a single vortex, we obtain an effective nonlinear dynamics for vortex cores (particles). The particles here acquire the inertia, in marked contrast to the standard theory of point vortices widely known in the usual hydrodynamics. The effective dynamics is equivalent to that of charged particles under a strong spring force and in the presence of Lorentz force with the uniform magnetic field. The interparticle (vortex-vortex) interaction is singularly repulsive and short ranged with its magnitude decreasing with increasing distance of the center of mass from the trapping center. ``Chaos in the three-body problem'' in the three-vortice system can be seen, which is not expected in the corresponding point vortices without inertia in two dimensions.
We consider an effectively one-dimensional binary Bose-Einstein condensate (BEC) with nonlinear repulsive interactions and linear spin-orbit (SO) and Zeeman-splitting couplings. In the presence of the trapping harmonic-oscillator (HO) potential, we report the existence of even, odd, and asymmetric spatial modes. They feature alternating domains with opposite directions of the pseudospin, i.e., antiferromagnetic structures, which is explained by the interplay of the linear couplings, HO confinement, and repulsive self-interaction. The number of the domains is determined by the strength of the SO coupling. The modes are constructed analytically in the weakly nonlinear system. The dynamical stability of the modes is investigated by means of the Bogoliubov--de Gennes equations and direct simulations. A notable result is that the multi-domain-wall (DW) structures are stable, alternating between odd and even shapes, while the simplest single-DW structure is unstable. Thus, the system features a transition to the complex ground states under the action of the SO coupling. The addition of the Zeeman splitting transforms the odd modes into asymmetric ones via spontaneous symmetry breaking. The results suggest possibilities for switching the binary system between states with opposite (pseudo)magnetization by external fields, and realization of similar stable states and dynamical effects in solid-state and nonlinear-optical settings emulated by the SO-coupled BECs.