Vortices in multilayer stacks of Bose–Einstein condensates with tilted dipoles

In this paper, we consider the formation of vortices in multilayer stacks of Bose–Einstein condensates with tilted dipoles by numerical simulations of the Gross–Pitaevskii equation. Different dependencies of critical rotation frequency (CRF) and optical lattice height, vortex number, and rotation frequency are studied, depending on the direction of the dipole axis and dipole strength. Our results show that the CRF in z = 0 is minimum. When the optical lattice height is gradually increased, the CRF decreases gradually. Reducing of dipole strength in anisotropic dipole–dipole interaction (DDI) favours the formation of vortices, and such decline in isotropic DDI hinders the creation of vortices. The reason for this difference is that the repulsive interaction is favorable and the attractive interaction is disadvantageous for the vortex formation. In addition, we study the first-order correlation function and focus on variation of coherence. For small rotation frequency, the break of coherence occurs earlier in the case of purely repulsive interaction. With an increase in rotation frequency, the coherence concurrently disappears in layer z = 2. Moreover, we also investigate the quenched dynamics, showing that the increase of angular momentum is induced by changing the direction of dipoles and in this process the vortex number remains unchanged.