We contribute a publicly available set of tables and code to provide Equations of State (EoS) for matter at neutron star densities. Our EoSes are constrained only by input from hadron physics and fundamental principles, without feedback from neutron star observations, and so without relying on General Relativity. They can therefore be used to test General Relativity itself, as well as modified gravity theories, with neutron star observables, without logical circularity. We have adapted state of the art results from NN chiral potentials for the low--density limit, pQCD results for the asymptotically high-density EoS, and use monotonicity and causality as the only restrictions for intermediate densities, for the EoS sets to remain as model-independent as is feasible today.

At very high densities, as for example in the core of a neutron star, matter may appear in the color-flavor locked (CFL) phase, which is a superfluid. This phase features topologically stable vortex solutions, which arise in a spinning superfluid as localized configurations carrying quanta of angular momentum. Despite the topological stability of these vortices they are not the lowest energy state of the system at neutron star densities and decay into triplets of semi-superfluid fluxtubes. In these proceedings we report on the progress of our numerical study in the Ginzburg-Landau approach, where we investigate lattices of semi-superfluid fluxtubes. The fluxtubes are obtained through controlled decay of global vortex configurations in the presence of a gauge field. Understanding the dynamics of semi-superfluid string configurations is important in the context of angular momentum transfer from a quark matter core of a neutron star beyond the core boundary, since not vortex-, but fluxtube pinning seems to be the relevant mechanism in this scenario.

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.