The core region of a neutron star may feature quark matter in the color-flavor- locked (CFL) phase. The CFL condensate breaks the baryon number symmetry, such that the phenomenon of superfluidity arises. If the core of the star is rotating, vortices will form in the superfluid, carrying the quanta of angular momentum. In a previous study we have solved the question of stability of these vortices, where we found numerical proof of a conjectured instability, according to which superfluid vortices will decay into an arrangement of so-called semi-superfluid fluxtubes. Here we report first results of an extension of our framework that allows us to study multi-vortex dynamics. This will in turn enable us to investigate the structure of semi-superfluid string lattices, which could be relevant to study pinning phenomena at the boundary of the core.

Superfluid vortices in the color-flavor-locked (CFL) phase of dense quark matter are known to be energetically disfavored relative to well-separated triplets of "semi-superfluid" color flux tubes. However, the short-range interaction (metastable versus unstable) has not been established. In this paper we perform numerical calculations using the effective theory of the condensate field, mapping the regions in the parameter space of coupling constants where the vortices are metastable versus unstable. For the case of zero gauge coupling we analytically identify a candidate for the unstable mode, and show that it agrees well with the results of the numerical calculations. We find that in the region of the parameter space that seems likely to correspond to real-world CFL quark matter the vortices are unstable, indicating that if such matter exists in neutron star cores it is very likely to contain semi-superfluid color flux tubes rather than superfluid vortices.