“…This set is naturally equipped with the structure of a Riemannian manifold, and this allows the use of Riemann optimization algorithms [1,21]. Direct minimization algorithms are preferred for the Gross-Pitaevskii model with magnetic fields [3,18,27,29], for which the Aufbau principle is not satisfied in general. Gradient-type [2,17,49,61,66], Newton-type [5,15,68], and trust-region methods have also been designed to solve (1.1) for larger values of N .…”