A Reduced Order technique to study bifurcating phenomena: application to the Gross-Pitaevskii equation
We are interested in the steady state solution of the Gross-Pitaevskii equation, which describes the ground state of a quantum system of identical bosons under simplifying hypotheses, and their stability properties. Often referred to as a non-linear Schrödinger equation, the Gross-Pitaevskii equation models certain classes of Bose-Einstein condensates (BECs), a special state of matter formed by bosons at ultra-low temperatures. It is well known that the solutions of the Gross-Pitaevskii equation with a parabolic trap in two dimensions exhibit a rich bifurcating behavior, which includes symmetry-breaking bifurcations and vortex-bearing states. The bifurcation diagram plots the number of bosons in a BEC as a function of the chemical potential. In order to trace a bifurcation diagram, one has to compute multiple solution of a parametrized, non-linear problem. We will start with a one parameter study (the chemical potential) and then consider a two parameter case (chemical potential and the normalized trap strength). In this work, we propose to apply a Reduced Order Modeling (ROM) technique to reduce the demanding computational costs associated with this stability analysis.
READ FULL TEXT