We develop a numerical algorithm for the solution of the Sturm-Liouville differential equation governing the stationary radial oscillations of nonrotating compact stars. Our method is based on the Numerov's method that turns the Sturm-Liouville differential equation in an eigenvalue problem. In our development we provide a strategy to correctly deal with the star boundaries and the interfaces between layers with different mechanical properties. Assuming that the fluctuations obey the same equation of state of the background, we analyze various different stellar models and we precisely determine hundreds of eigenfrequencies and of eigenmodes. If the equation of state does not present an interface discontinuity, the fundamental radial eigenmode becomes unstable exactly at the critical central energy density corresponding to the largest gravitational mass. However, in the presence of an interface discontinuity, there exist stable configurations with a central density exceeding the critical one and with a smaller gravitational mass.

We study the thermodynamic properties of matter at vanishing temperature for non-extreme values of the isospin chemical potential and of the strange quark chemical potential. From the leading order pressure obtained by maximizing the static chiral Lagrangian density we derive a simple expression for the equation of state in the pion condensed phase and in the kaon condensed phase. We find an analytical expression for the maximum of the ratio between the energy density and the Stefan-Boltzmann energy density as well as for the isospin chemical potential at the peak both in good agreement with lattice simulations of quantum chromodynamics. We speculate on the location of the crossover from the Bose-Einstein condensate state to the Bardeen-Cooper-Schrieffer state by a simple analysis of the thermodynamic properties of the system. For $\mu_I \gtrsim 2 m_\pi$ the leading order chiral perturbation theory breaks down; as an example it underestimates the energy density of the system and leads to a wrong asymptotic behavior.

Inhomogeneous phases may appear when a stress is applied to a system and the system can minimize the free energy breaking the rotational invariance. Various examples are known in Nature of this sort, as the paramagnetic to ferromagnetic phase transition, or the fluid/solid phase transition. If the rotational symmetry is broken down to a discrete symmetry, the system is typically named a crystal. We breifly review crystalline color superconductors, which arise in cold quark matter with mismatched Fermi spheres.

Strange stars are one of the possible compact stellar objects that can be formed after a supernova collapse. We consider a model of strange star having an inner core in the color-flavor locked phase surmounted by a crystalline color superconducting layer. These two phases constitute the {\it quarksphere}, which we assume to be the largest and heaviest part of the strange star. The next layer consists of standard nuclear matter forming a ionic crust, hovering on the top of the quarksphere and prevented from falling by a strong dipolar electric field. The dipolar electric field arises because quark matter is confined in the quarksphere by the strong interaction, but electrons can leak outside forming a few hundreds Fermi thick electron layer separating the ionic crust from the underlying quark matter. The ionic matter and the crystalline color superconducting matter constitute two electromagnetically coupled crust layers. We study the torsional oscillations of these two layers. Remarkably, we find that if a fraction larger than $10^{-4}$ of the energy of a Vela-like glitch is conveyed to a torsional oscillation, the ionic crust will likely break. The reason is that the very rigid and heavy crystalline color superconducting crust layer will absorb only a small fraction of the glitch energy, leading to a large amplitude torsional oscillation of the ionic crust.

Strange stars are one of the hypothetical compact stellar objects that can be formed after a supernova explosion. The existence of these objects relies on the absolute stability of strange {\it collapsed} quark matter with respect to standard nuclear matter. We discuss simple models of strange stars with a bare quark matter surface, thus standard nuclear matter is completely absent. In these models an electric dipole layer a few hundreds Fermi thick should exist close to the star surface. Studying the torsional oscillations of the electrically charged layer we estimate the emitted power, finding that it is of the order of $10^{45}$ erg/s, meaning that these objects would be among the brightest compact sources in the heavens. The associated relaxation times are very uncertain, with values ranging between microseconds and minutes, depending on the crust thickness. Although part of the radiated power should be absorbed by the electrosphere surrounding the strange star, a sizable fraction of photons should escape and be detectable.

The crystalline color superconducting phase is believed to be the ground state of deconfined quark matter for sufficiently large values of the strange quark mass. This phase has the remarkable property of being more rigid than any known material. It can therefore sustain large shear stresses, supporting torsional oscillations of large amplitude. The torsional oscillations could lead to observable electromagnetic signals if strange stars have a crystalline color superconducting crust. Indeed, considering a simple model of strange star with a bare quark matter surface, it turns out that a positive charge is localized in a narrow shell about ten Fermi thick beneath the star surface. The electrons needed to neutralize the positive charge of quarks spill in the star exterior forming an electromagnetically bounded atmosphere hundreds of Fermi thick. When a torsional oscillation is excited, for example by a stellar glitch, the positive charge oscillates with typical kHz frequencies, for a crust thickness of about one-tenth of the stellar radius, to hundreds of Hz, for a crust thickness of about nine-tenths of the stellar radius. Higher frequencies, of the order of few GHz, can be reached if the star crust is of the order of few centimeters thick. We estimate the emitted power considering emission by an oscillating magnetic dipole, finding that it can be quite large, of the order of $10^{45}$ erg/s for a thin crust. The associated relaxation times are very uncertain, with values ranging between microseconds and minutes, depending on the crust thickness. The radiated photons will be in part absorbed by the electronic atmosphere, but a sizable fraction of them should be emitted by the star.

Inhomogeneous superconductors and inhomogeneous superfluids appear in a variety of contexts including quark matter at extreme densities, fermionic systems of cold atoms, type-II cuprates, and organic superconductors. In the present review the focus is on properties of quark matter at high baryonic density, which may exist in the interior of compact stars. The conditions realized in these stellar objects tend to disfavor standard symmetric BCS pairing and may favor an inhomogeneous color superconducting phase. The properties of inhomogeneous color superconductors are discussed in detail and in particular of crystalline color superconductors. The possible astrophysical signatures associated with the presence of crystalline color superconducting phases within the core of compact stars are also reviewed.

This paper is a brief journey into the amazing realm of crystalline color superconductors. Starting from a qualitative description of superfluids, superconductors and supersolids, we show how inhomogeneous phases may arise when the system is under stress. These basic concepts are then extended to quark matter, in which a richer variety of phases can be realized. Then, the most interesting properties of the crystalline color superconductors are presented. This brief journey ends with a discussion of crystalline color superconductors in compact stars and related astrophysical observables. We aim at providing a pedagogical introduction for nonexpert in the field to a few interesting properties of crystalline color superconductors, without discussing the methods and the technicalities. Thus, the results are presented without a proof. However, we try to give a qualitatively clear description of the main concepts, using standard quantum field theory and analogies with condensed matter systems.

We study the low energy effective action for the collective modes of the color flavor locked phase of QCD. This phase of matter has long been known to be a superfluid because by picking a phase its order parameter breaks the quark-number $U(1)_B$ symmetry spontaneously. We consider the modes describing fluctuations in the magnitude of the condensate, namely the Higgs mode, and in the phase of the condensate, namely the Nambu-Goldstone (or Anderson-Bogoliubov) mode associated with the breaking of $U(1)_B$. By employing as microscopic theory the Nambu-Jona Lasinio model, we reproduce known results for the Lagrangian of the Nambu-Goldstone field to the leading order in the chemical potential and extend such results evaluating corrections due to the gap parameter. Moreover, we determine the interaction terms between the Higgs and the Nambu-Goldstone field. This study paves the way for a more reliable study of various dissipative processes in rotating compact stars with a quark matter core in the color flavor locked phase.

We present some results about a novel damping mechanism of r-mode oscillations in neutron stars due to processes that change the number of protons, neutrons and electrons. Deviations from equilibrium of the number densities of the various species lead to the appearance in the Euler equations of the system of a dissipative mechanism, the so-called rocket effect. The evolution of the r-mode oscillations of a rotating neutron star are influenced by the rocket effect and we present estimates of the corresponding damping timescales. In the description of the system we employ a two-fluid model, with one fluid consisting of all the charged components locked together by the electromagnetic interaction, while the second fluid consists of superfluid neutrons. Both components can oscillate however the rocket effect can only efficiently damp the countermoving r-mode oscillations, with the two fluids oscillating out of phase. In our analysis we include the mutual friction dissipative process between the neutron superfluid and the charged component. We neglect the interaction between the two r-mode oscillations as well as effects related with the crust of the star. Moreover, we use a simplified model of neutron star assuming a uniform mass distribution.

We consider the phonon contribution to the bulk viscosities $\zeta_1, \zeta_2$ and $\zeta_3$ of a cold relativistic superfluid. We assume the low temperature $T$ regime and that the transport properties of the system are dominated by the phonons. We use kinetic theory in the relaxation time approximation and the low energy effective field theory of the corresponding system. The parametric dependence of the bulk viscosity coefficients is fixed once the equation of state is specified, and the phonon dispersion law to cubic order in momentum is known. We first present a general discussion, valid for any superfluid, then we focus on the color-flavor locked superfluid because all the parameters needed in the analysis can be computed in the high density limit of QCD, and also because of the possible astrophysical applications. For the three independent bulk viscosity coefficients we find that they scale with the temperature as $\zeta_i \sim 1/T$, and that in the conformal limit only the third coefficient $\zeta_3$ is non-zero.

Dissipative processes acting in rotating neutron stars are essential in preventing the growth of the r-mode instability. We estimate the damping time of r-modes of an hypothetical compact quark star made up by color flavor locked quark matter at a temperature $T \lesssim 0.01$ MeV. The dissipation that we consider is due to the the mutual friction force between the normal and the superfluid component arising from the elastic scattering of phonons with quantized vortices. This process is the dominant one for temperatures $T \lesssim 0.01$ MeV where the mean free path of phonons due to their self-interactions is larger than the radius of the star and they can be described as an ideal bosonic gas. We find that r-modes oscillations are efficiently damped by this mechanism for pulsars rotating at frequencies of the order of 1 Hz at most. Our analysis rules out the possibility that cold pulsars rotating at higher frequencies are entirely made up by color flavor locked quark matter.

We write a covariant transport equation for the phonon excitations of a relativistic superfluid valid at small temperatures. The hydrodynamical equations for this system are derived from the effective field theory associated to the superfluid phonons. We describe how to construct the kinetic theory for the phonon quasiparticles using a relativistic generalization of the analogue model of gravity developed by Unruh. This gravity analogy relies on the equivalence between the action of a phonon field moving in a superfluid background with that of a boson propagating in a given curved space-time. Exploiting this analogy we obtain continuity equations for the phonon current, entropy and energy-momentum tensor in a covariant form, valid in any reference frame. Our aim is to shed light on some aspects of transport phenomena of relativistic superfluidity. In particular, we are interested in the low temperature regime of the color flavor locked phase, which is a color superconducting and superfluid phase of high density QCD that may be realized in the core of neutron stars.

We present a study of the shear modulus of the crystalline color superconducting phase of quark matter, showing that this phase of dense, but not asymptotically dense, quark matter responds to shear stress as a very rigid solid. This phase is characterized by a gap parameter $\Delta$ that is periodically modulated in space and therefore spontaneously breaks translational invariance. We derive the effective action for the phonon fields that describe space- and time-dependent fluctuations of the crystal structure formed by $\Delta$, and obtain the shear modulus from the coefficients of the spatial derivative terms. Within a Ginzburg-Landau approximation, we find shear moduli which are 20 to 1000 times larger than those of neutron star crusts. This phase of matter is thus more rigid than any known material in the universe, but at the same time the crystalline color superconducting phase is also superfluid. These properties raise the possibility that the presence of this phase within neutron stars may have distinct implications for their phenomenology. For example, (some) pulsar glitches may originate in crystalline superconducting neutron star cores.

We calculate the shear modulus of crystalline color superconducting quark matter, showing that this phase of dense, but not asymptotically dense, three-flavor quark matter responds to shear stress like a very rigid solid. To evaluate the shear modulus, we derive the low energy effective Lagrangian that describes the phonons that originate from the spontaneous breaking of translation invariance by the spatial modulation of the gap parameter $\Delta$. These massless bosons describe space- and time-dependent fluctuations of the crystal structure and are analogous to the phonons in ordinary crystals. The coefficients of the spatial derivative terms of the phonon effective Lagrangian are related to the elastic moduli of the crystal; the coefficients that encode the linear response of the crystal to a shearing stress define the shear modulus. We analyze the two particular crystal structures which are energetically favored over a wide range of densities, in each case evaluating the phonon effective action and the shear modulus up to order $\Delta^2$ in a Ginzburg-Landau expansion, finding shear moduli which are 20 to 1000 times larger than those of neutron star crusts. The crystalline color superconducting phase has long been known to be a superfluid -- by picking a phase its order parameter breaks the quark-number $U(1)_B$ symmetry spontaneously. Our results demonstrate that this superfluid phase of matter is at the same time a rigid solid. We close with a rough estimate of the pinning force on the rotational vortices which would be formed embedded within this rigid superfluid upon rotation. Our results raise the possibility that (some) pulsar glitches could originate within a quark matter core deep within a neutron star.

We calculate the specific heats of quasi-particles of two-flavor QCD in its crystalline phases for low temperature. We show that for the different crystalline structures considered here there are gapless modes contributing linearly in temperature to the specific heat. We evaluate also the phonon contributions which are cubic in temperature. These features might be relevant for compact stars with an inner shell in a color superconducting crystalline phase.

The gluon self-energies and dispersion laws in the color superconducting phase of QCD with two massless flavors are calculated using the effective theory near the Fermi surface. These quantities are calculated at zero temperature for all the eight gluons, those of the remaining SU(2) color group and those corresponding to the broken generators. The construction of the effective interaction is completed with the one loop calculation of the three- and four-point gluon interactions.

We present an effective field theory for high density, low temperature QCD in the crystalline colour superconductive phase (LOFF phase). This interesting phase of QCD is characterized by a gap parameter with a crystalline pattern, breaking traslational and rotational invariance, and could have astrophysical applications. In the effective theory the fermions have a Majorana mass, which, besides colour, breaks translation and rotation symmetries. Fermions couple to the three phonons arising from the breaking of rotation and translation invariance. Integrating out the fermions leads eventually to an effective lagrangian in terms of the phonon fields only, which satisfies an anisotropic dispersion relation.