## References

- Dynamics of bright matter wave solitons in a bose–einstein condensate. Int. J. Mod. Phys. B 19, pp. 3415–3473. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Discrete and continuous nonlinear schrödinger systems. Cambridge University Press, Cambridge. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Nonlinear dispersive waves. Cambridge University Press, Cambridge. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Three-dimensional vortex configurations in a rotating bose-einstein condensate. Phys. Rev. A 68, pp. 023603. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Bifurcation analysis of stationary solutions of two-dimensional coupled Gross–Pitaevskii equations using deflated continuation. arXiv preprint arXiv:1912.00023. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Computing stationary solutions of the two-dimensional Gross–Pitaevskii equation with deflated continuation. Commun. Nonlinear Sci. Numer. Simulat. 54, pp. 482–499. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Three-dimensional parallel vortex rings in Bose–Einstein condensates. Phys. Rev. A 70, pp. 033605. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Finite-difference solutions of a non-linear Schrödinger equation. J. Comput. Phys. 44 (2), pp. 277–288. External Links: Link Cited by: Appendix A, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Deflation techniques for finding distinct solutions of nonlinear partial differential equations. SIAM J. Sci. Comput. 37 (4), pp. A2026–A2045. External Links: Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Vortices in a trapped dilute bose-einstein condensate. J. Phys.-Condens. Mat. 13 (12), pp. R135. External Links: Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Rotating trapped bose-einstein condensates. Rev. Mod. Phys. 81, pp. 647–691. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Dark solitons in atomic Bose–Einstein condensates: from theory to experiments. J. Phys. A-Math. Theor. 43 (21), pp. 213001. External Links: ISSN 1751-8121, Link, Document Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Solitons in Optical Communications. Clarendon Press, Oxford. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems. ACM Trans. Math. Softw. 31 (3), pp. 351–362. External Links: Link, Document Cited by: Appendix A, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- The defocusing nonlinear schrodinger equation. Society for Industrial and Applied Mathematics, Philadelphia, PA. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Rogue waves in the ocean. Springer Verlag, Berlin. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Creation and dynamics of knotted vortices. Nat. Phys. 9, pp. 253–258. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Vortex rings and solitary waves in trapped bose–einstein condensates. Eur. Phys. J.-Spec. Top. 147 (1), pp. 133–152. External Links: ISSN 1951-6401, Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Nonlinear physics of plasmas. Springer Verlag, Heidelberg. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Chladni solitons and the onset of the snaking instability for dark solitons in confined superfluids. Phys. Rev. Lett. 113 (25), pp. 255302. External Links: Document, Link Cited by: Appendix A, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Stability and dispersion relations of three-dimensional solitary waves in trapped Bose–Einstein condensates. New J. Phys. 17 (12), pp. 125013. External Links: Link Cited by: Appendix A.
- Stationary vortex clusters in nonrotating bose-einstein condensates. Phys. Rev. A 71, pp. 033626. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Bose–Einstein Condensation in Dilute Gases. Cambridge University Press, Cambridge, United Kingdom. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Firedrake: automating the finite element method by composing abstractions. ACM Trans. Math. Softw. 43 (3), pp. 1–27. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Creating vortex rings and three-dimensional skyrmions in bose-einstein condensates. Phys. Rev. Lett. 86, pp. 3934–3937. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Monopole core instability and alice rings in spinor bose-einstein condensates. Phys. Rev. Lett. 91, pp. 190402. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Practical bifurcation and stability analysis. 3rd edition, Interdisciplinary Applied Mathematics, Vol. 5, Springer. Cited by: Appendix A.
- A Krylov–Schur algorithm for large eigenproblems. SIAM J. Matrix Anal. A. 23 (3), pp. 601–614. External Links: Link, Document Cited by: Appendix A, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Bose–einstein condensation. Oxford University Press, Oxford, United Kingdom. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- The nonlinear Schrödinger equation: self-focusing and wave collapse. Springer, New York. Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Ring dark solitons in three-dimensional bose-einstein condensates. Phys. Rev. A 100, pp. 053621. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.
- Dark spherical shell solitons in three-dimensional bose-einstein condensates: existence, stability, and dynamics. Phys. Rev. A 93, pp. 023630. External Links: Document, Link Cited by: Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation, Deflation-based Identification of Nonlinear Excitations of the 3D Gross–Pitaevskii equation.

## Appendix A Supplemental Material

Solutions emanating from the 1st excited state. Deflation identified three solutions emanating from the first excited state at . Fig. 6(a) shows a Cartesian state with one cut along the -axis. Fig. 6(b) shows corresponding to a single vortex line with topological charge . The rotations of these solutions along the , , and axes such as the and Cartesian states are also obtained by deflation but are not reported.

Bifurcation diagram of the 2nd excited states. As discussed in the paper, the steady-state solutions to the nonlinear Schrödinger equation are identified by the deflation method at . The branches are then continued backward in down to the linear limit by a standard zero-order continuation method (Seydel, 2010, §4.4.2). A typical example of a relevant solution is shown in Fig. 6(c) and represents a so-called Chladni soliton, previously identified in cylindrical geometry in Mateo and Brand (2014, 2015). We present the bifurcation diagram of the solutions emanating from the second excited state in Fig. 7. Our diagnostic functional is the total number of atoms (or squared norm):

(7) |

The inset panel in the top-left corner of Fig. 7 uses the atom number difference between the branches 2(a) and 2(a)* to illustrate a bifurcation in the diagram. The latter branch bifurcates from the former around and is shown in Fig. 8.

Bogoliubov-de Gennes stability analysis. The stability of a stationary solution to the NLS equation is determined by perturbing it with the following ansatz:

(8) |

where is a small perturbation parameter, is the eigenfrequency, and the corresponding eigenvector. After substituting (8) into the time-dependent NLS equation:

(9) |

we obtain the following complex eigenvalue problem

(10) |

where is the eigenvalue and the matrix elements are given by equationparentequation

(11a) | ||||

(11b) |

We decompose the eigenvector into real and imaginary components as , , and rewrite (10) as

(12) | ||||

where and . The eigenvalues of the matrix on the right-hand side of (12) are (with multiplicity two). Therefore, solving a real eigenvalue problem with the left-hand matrix of (12) yields the same eigenvalues and eigenvectors as the complex eigenvalue problem (10). We use a Krylov–Schur algorithm with a shift-and-invert spectral transformation Stewart (2002), implemented in the SLEPc library Hernandez et al. (2005), to solve the following eigenvalue problem:

(13) |

where the matrices are real and is complex. This problem is discretized with the same piecewise cubic finite element method used to find multiple solutions with deflation. The spectra of the solutions emanating from the 2nd excited states, presented in Figs. 1 and 2 of the paper, are respectively displayed in the different panels of Figs. 10 and 9. It is interesting to observe that some of the states such as the doubly charged vortex line of Fig. 1(e) can be dynamically stable for large values of the chemical potential , while others such as the anisotropic ring of Fig. 1(f) can be stable for sufficiently low .

Dynamics. The NLS equation (9) is integrated in time until using the following perturbed stationary solution as initial state:

(14) |

where is a stationary solution discovered by deflation and is its most unstable eigendirection. The eigenvector is normalized so that

(15) |

and the perturbation parameter is chosen to be . We use a modified Crank-Nicolson method in time Delfour et al. (1981), which preserves to machine precision the square of the norm (i.e., atom number) and the energy of the solutions

(16) |

and piecewise cubic finite elements in space. Given the solution at time , is obtained by solving the following nonlinear PDE:

(17) | ||||

where is the time step used in the numerical simulations. The nonlinear problem is solved with Newton’s method. We present different snapshots of the two vortex rings solution (see Fig. 2(f) of the paper) in Fig. 11. We observe that the original two rings collapse around to form a single vortex line (see also the relevant movie in the Supplemental Material).

Comments

There are no comments yet.