Exploring the free-energy landscape of a rotating superfluid

The equilibrium state of a superfluid in a rotating cylindrical vessel is a vortex crystal-an array of vortex lines, which is stationary in the rotating frame. Experimental realizations of this behavior typically show a sequence of transient states before the free-energy-minimizing configuration is reached. Motivated by these observations, we construct a new method for a systematic exploration of the free-energy landscape via gradient-based optimization of a scalar loss function. Our approach is inspired by the pioneering numerical work of Campbell and Ziff [Phys. Rev. B. 20, 1886 (1979)] and makes use of automatic differentiation, which crucially allows us to include entire solution trajectories in the loss. We first use the method to converge thousands of low free-energy relative equilibria in the unbounded domain for vortex numbers in the range 10≤N≤30, which reveals an extremely dense set of mostly saddle-like solutions. As part of this search, we discover new continuous families of relative equilibria, which are often global minimizers of free energy. These continuous families all consist of crystals arranged in a double-ring configuration, and we assess which state from the family is most likely to be observed experimentally by computing energy-minimizing pathways from nearby local minima-identifying a common entry point into the family. The continuous families become discrete sets of equal-energy solutions when the wall is introduced in the problem. Finally, we develop an approach to compute homoclinic orbits and use it to examine the dynamics in the vicinity of the minimizing state by converging connections for low-energy saddles.