The transport coefficients of quark matter at non-zero chemical potential and temperature are computed within the two-flavor Nambu--Jona-Lasinio model. We apply the Kubo formalism to obtain the thermal ($\kappa$) and electrical ($\sigma$) conductivities as well as an update of the shear viscosity ($\eta$) by evaluating the corresponding equilibrium two-point correlation functions to leading order in the $1/N_c$ expansion. The Dirac structure of the self-energies and spectral functions is taken into account as these are evaluated from the meson-exchange Fock diagrams for on-mass-shell quarks. We find that the thermal and electrical conductivities are decreasing functions of temperature and density above the Mott temperature $T_{\rm M}$ of dissolution of mesons into quarks, the main contributions being generated by the temporal and vector components of the spectral functions. The coefficients show a universal dependence on the ratio $T/T_{\rm M}$ for different densities, i.e., the results differ by a chemical-potential dependent constant. We also show that the Wiedemann-Franz law for the ratio $\sigma/\kappa$ does not hold. The ratio $\eta/s $, where $s$ is the entropy density, is of order of unity (or larger) close to the Mott temperature and, as the temperature increases, approaches the AdS/CFT bound $1/4\pi$. It is also conjectured that the ratio $\kappa T/c_V $, with $c_V$ being the specific heat, is bounded from below by $1/18$.

In this paper, we discuss the damping of density oscillations in dense nuclear matter in the temperature range relevant to neutron star mergers. This damping is due to bulk viscosity arising from the weak interaction ``Urca'' processes of neutron decay and electron capture. The nuclear matter is modelled in the relativistic density functional approach. The bulk viscosity reaches a resonant maximum close to the neutrino trapping temperature, then drops rapidly as temperature rises into the range where neutrinos are trapped in neutron stars. We investigate the bulk viscous dissipation timescales in a post-merger object and identify regimes where these timescales are as short as the characteristic timescale $\sim$10 ms, and, therefore, might affect the evolution of the post-merger object. Our analysis indicates that bulk viscous damping would be important at not too high temperatures of the order of a few MeV and densities up to a few times saturation density.

We compute the mass, radius and tidal deformability of stars containing phase transitions from hadronic to quark phase(s). These quantities are computed for three types of hadronic envelopes: purely nuclear, hyperonic, and $\Delta$-resonance--hyperon admixed matter. We consider either a single first-order phase transition to a quark phase with a maximally stiff equation of state (EOS) or two sequential first-order phase transitions mimicking a transition from hadronic to a quark matter phase followed by a second phase transition to another quark phase. We explore the parameter space which produces low-mass twin and triplet configurations where equal-mass stars have substantially different radii and tidal deformabilities. We demonstrate that while for purely hadronic stiff EOS the obtained maximum mass is inconsistent with the upper limit on this quantity placed by GW170817, the inclusion of the hyperonic and $\Delta$-resonance degrees of freedom, as well as the deconfinement phase transition at sufficiently low density, produces a configuration of stars consistent with this limit. The obtained hybrid star configurations are in the mass range relevant for the interpretation of the GW170817 event. We compare our results for the tidal deformability with the limits inferred from GW170817 showing that the onset of non-nucleonic phases, such as $\Delta$-resonance--hyperon admixed phase and/or the quark phase(s), are favored by this data if the nuclear EOS is stiff. Also, we show that low-mass twins and especially triplets proliferate the number of combinations of possible types of stars that can undergo a merger event, the maximal number being six in the case of triplets. The prospects for uncovering the first-order phase transition(s) to and in quark matter via measurements of tidal deformabilities in merger events are discussed.

We study bulk viscosity arising from weak current Urca processes in dense baryonic matter at and beyond nuclear saturation density. We consider the temperature regime where neutrinos are trapped and therefore have nonzero chemical potential. We model the nuclear matter in a relativistic density functional approach, taking into account the trapped neutrino component. We find that the resonant maximum of the bulk viscosity would occur at or below the neutrino trapping temperature, so in the neutrino trapped regime the bulk viscosity decreases with temperature as $T^{-2}$, this decrease being interrupted by a drop to zero at a special temperature where the proton fraction becomes density-independent and the material scale-invariant. The bulk viscosity is larger for matter with lower lepton fraction, i.e., larger isospin asymmetry. We find that bulk viscosity in the neutrino-trapped regime is smaller by several orders of magnitude than in the neutrino-transparent regime, which implies that bulk viscosity in neutrino-trapped matter is probably not strong enough to affect significantly the evolution of neutron star mergers. This also implies weak damping of gravitational waves emitted by the oscillations of the postmerger remnant in the high-temperature, neutrino-trapped phase of evolution.

Nuclear matter and finite nuclei exhibit the property of superfluidity by forming Cooper pairs. We review the microscopic theories and methods that are being employed to understand the basic properties of superfluid nuclear systems, with emphasis on the spatially extended matter encountered in neutron stars, supernova envelopes, and nuclear collisions. Our survey of quantum many-body methods includes techniques that employ Green functions, correlated basis functions, and Monte Carlo sampling of quantum states. With respect to empirical realizations of nucleonic and hadronic superfluids, this review is focused on progress that has been made toward quantitative understanding of their properties at the level of microscopic theories of pairing, with emphasis on the condensates that exist under conditions prevailing in neutron-star interiors. These include singlet $S$-wave pairing of neutrons in the inner crust, and, in the quantum fluid interior, singlet-$S$ proton pairing and triplet coupled $P$-$F$-wave neutron pairing. Additionally, calculations of weak-interaction rates in neutron-star superfluids within the Green function formalism are examined in detail. We close with a discussion of quantum vortex states in nuclear systems and their dynamics in neutron-star superfluid interiors.

A set of hadronic equations of state (EoSs) derived from relativistic density functional theory and constrained by terrestrial experiments, astrophysical observations, in particular by the GW170817 event, and chiral effective field theory ($\chi$EFT) of neutron matter is used to explore the sensitivity of the EoS parameterization on the few nuclear matter characteristics defined at the saturation density. We find that the gross properties of compact stars are most sensitive to the isoscalar skewness coefficient $Q_{\text{sat}}$ and the isovector slope coefficient $L_{\text{sym}}$ around saturation density, since the higher order coefficients, such as $K_{\text{sym}}$, are fixed by our model. More specifically, (i) among these $Q_{\text{sat}}$ is the dominant parameter controlling both the maximum mass and the radii of compact stars while $L_{\rm sym}$ is constrained somewhat by $\chi$EFT of neutron matter; (ii) massive enough ($M\sim 2.0~M_{\odot})$ compact stars featuring both hyperons and $\Delta$ resonances can be obtained if the value of $Q_{\text{sat}}$ is large enough; (iii) the emergence of $\Delta$'s reduces the radius of a canonical mass ($M\sim 1.4~M_{\odot})$ compact star thus easing the tension between the predictions of the relativistic density functionals and the inferences from the X-ray observation of nearby isolated neutron stars.

The thermal evolution of hypernuclear compact stars is studied for stellar models constructed on the basis of covariant density functional theory in Hartree and Hartree-Fock approximation. Parametrizations of both types are consistent with the astrophysical mass constraints on compact stars and available hypernuclear data. We discuss the differences of these density functionals and highlight the effects they have on the composition and on the cooling of hypernuclear stars. It is shown that hypernuclear stars computed with density functional models that have a low symmetry energy slope, $L$, are fairly consistent with the cooling data of observed compact stars. The class of stellar models based on larger $L$ values gives rise to the direct Urca process at low densities, which leads to significantly faster cooling. We conjecture high-density pairing for protons and $\Lambda$'s in the $P$-wave channel and provide simple scaling arguments to obtain these gaps. As a consequence the most massive stellar models with masses $1.8 \le M/M_{\odot} \le2$ experience slower cooling by hyperonic dUrca processes which involve $\Lambda$'s and protons.

The effects of $\Delta$ isobars on the equation of state of dense matter and structure of compact stars (CSs) are explored within the covariant density functional theory and confronted with the data on tidal deformability (TD) extracted from the GW170817 event. We show that the presence of $\Delta$ isobars substantially softens the tension between the predictions of the hypernuclear density functionals and the inference from the observations of relatively small radius and small TD of canonical mass CSs. The TDs deduced from GW170817 are compatible with the existence of hypernuclear CSs containing an admixture of $\Delta$ isobars. We thus argue that the GW170817 event is consistent with a merger of a binary CS system having both strangeness (hyperons) and $\Delta$ isobars in the stellar core.

We study the axion cooling of neutron stars within the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) model, which allows for tree level coupling of electrons to the axion {and locks the Peccei-Quinn charges of fermions via an angle parameter}. This extends our previous study [Phys. Rev. D 93, 065044 (2016)] limited to hadronic models of axions. We explore the two-dimensional space of axion parameters within the DFSZ model by comparing the theoretical cooling models with the surface temperatures of a few stars with measured surface temperatures. It is found that axions masses $m_a\ge 0.06$ to 0.12 eV can be excluded by x-ray observations of thermal emission of neutron stars (in particular by those of Cas A), the precise limiting value depending on the angle parameter of the DFSZ model. It is also found that axion emission by electron bremsstrahlung in neutron star crusts is negligible except for the special case where neutron Peccei-Quinn charge is small enough, so that the coupling of neutrons to axions can be neglected.

We examine the range of rest-mass densities, temperatures and magnetic fields involved in simulations of binary neutron-star mergers and identify the conditions under which the ideal-magnetohydrodynamics approximation breaks down and hence the magnetic-field decay should be accounted for. We use recent calculations of the conductivities of warm correlated plasma in envelopes of compact stars and find that the magnetic-field decay timescales are much larger than the characteristic timescales of the merger process for lengthscales down to a meter. Because these are smaller than the currently available resolution in numerical simulations, the ideal-magnetohydrodynamics approximation is effectively valid for all realistic simulations. At the same time, we find that the Hall effect can be important at low densities and low temperatures, where it can induce a non-dissipative rearrangement of the magnetic field. Finally, we mark the region in temperature and density where the hydrodynamic description breaks down.

The hypernuclear matter is studied within the relativistic Hartree-Fock theory employing several parametrizations of the hypernuclear density functional with density-dependent couplings. The equations of state and compositions of hypernuclear matter are determined for each parametrization and compact stars are constructed by solving their structure equations in spherical symmetry. We quantify the softening effect of Fock terms on the equation of state, as well as discuss the impact of tensor interactions, which are absent in the Hartree theories. Starting from models of density functionals which are fixed in the nuclear sector to the nuclear phenomenology, we vary the couplings in the hyperonic sector around the central values which are fitted to the hyperon potentials in nuclear matter. We use the SU(6) spin-flavor and SU(3) flavor symmetric quark models to relate the hyperonic couplings to the nucleonic ones. We find, consistent with previous Hartree studies, that for the SU(6) model the maximal masses of compact stars are below the two-solar mass limit. In the SU(3) model we find sufficiently massive compact stars with cores composed predominantly of $\Lambda$ and $\Xi$ hyperons and a low fraction of leptons (mostly electrons). The parameter space of the SU(3) model is identified where simultaneously hypernuclear compact stars obey the astrophysical limits on pulsar masses and the empirical hypernuclear potentials in nuclear matter are reproduced.

The $\Delta$-isobar degrees of freedom are included in the covariant density functional (CDF) theory to study the equation of state (EoS) and composition of dense matter in compact stars. In addition to $\Delta$'s we include the full octet of baryons, which allows us to study the interplay between the onset of delta isobars and hyperonic degrees of freedom. Using both the Hartree and Hartree-Fock approximation we find that $\Delta$'s appear already at densities slightly above the saturation density of nuclear matter for a wide range of the meson-$\Delta$ coupling constants. This delays the appearance of hyperons and significantly affects the gross properties of compact stars. Specifically, $\Delta$'s soften the EoS at low densities but stiffen it at high densities. This softening reduces the radius of a canonical $1.4 M_\odot$ star by up to 2~km for a reasonably attractive $\Delta$ potential in matter, while the stiffening results in larger maximum masses of compact stars. We conclude that the hypernuclear CDF parametrizations that satisfy the 2$M_\odot$ maximum mass constraint remain valid when $\Delta$ isobars are included, with the important consequence that the resulting stellar radii are shifted toward lower values, which is in agreement with the analysis of neutron star radii.

Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, that give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star--neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.

We study the two-flavor color superconductivity of low-temperature quark matter in the vicinity of chiral phase transition in the quark-meson model where the interactions between quarks are generated by pion and sigma exchanges. Starting from the Nambu-Gor'kov propagator in real-time formulation we obtain finite temperature (real axis) Eliashberg-type equations for the quark self-energies (gap functions) in terms of the in-medium spectral function of mesons. Exact numerical solutions of the coupled nonlinear integral equations for the real and imaginary parts of the gap function are obtained in the zero temperature limit using a model input spectral function. We find that these components of the gap display a complicated structure with the real part being strongly suppressed above $2\Delta_0$, where $\Delta_0$ is its on-shell value. We find $\Delta_0\simeq 40$ MeV close to the chiral phase transition.

We study the thermal evolution of hypernuclear compact stars constructed from covariant density functional theory of hypernuclear matter and parameterizations which produce sequences of stars containing two-solar-mass objects. For the input in the simulations, we solve the Bardeen-Cooper-Schrieffer gap equations in the hyperonic sector and obtain the gaps in the spectra of $\Lambda$, $\Xi^0$ and $\Xi^-$ hyperons. For the models with masses $M/M_{\odot} \ge 1.5$ the neutrino cooling is dominated by hyperonic direct Urca processes in general. In the low-mass stars the $(\Lambda p)$ plus leptons channel is the dominant direct Urca process, whereas for more massive stars the purely hyperonic channels $(\Sigma^-\Lambda)$ and $(\Xi^-\Lambda)$ are dominant. Hyperonic pairing strongly suppresses the processes on $\Xi^-$s and to a lesser degree on $\Lambda$s. We find that intermediate-mass $1.5 \le M/M_{\odot} \le 1.8$ models have surface temperatures which lie within the range inferred from thermally emitting neutron stars, if the hyperonic pairing is taken into account. Most massive models with $M/M_{\odot} \simeq 2$ may cool very fast via the direct Urca process through the $(\Lambda p)$ channel because they develop inner cores where the $S$-wave pairing of $\Lambda$s and proton is absent.

The phase diagram of isospin-asymmetrical nuclear matter may feature a number of unconventional phases, which include the translationally and rotationally symmetric, but isospin-asymmetrical BCS condensate, the current-carrying Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase, and the heterogeneous phase-separated phase. Because the Cooper pairs of the condensate carry a single unit of charge, these phases are charged superconductors and respond to electromagnetic gauge fields by either forming domains (type-I superconductivity) or quantum vortices (type-II superconductivity). We evaluate the Ginzburg-Landau (GL) parameter across the phase diagram and find that the unconventional phases of isospin-asymmetrical nuclear matter are good type-II superconductors and should form Abrikosov vortices with twice the quantum of magnetic flux. We also find that the LOFF phase at the boundary of the transition to the type-I state, with the GL parameter being close to the critical value $1/\sqrt{2}$.

Compact stars may contain quark matter in their interiors at densities exceeding several times the nuclear saturation density. We explore models of such compact stars where there are two first-order phase transitions: the first from nuclear matter to a quark-matter phase, followed at higher density by another first-order transition to a different quark matter phase [e.g., from the two-flavor color superconducting (2SC) to the color-flavor-locked (CFL) phase). We show that this can give rise to two separate branches of hybrid stars, separated from each other and from the nuclear branch by instability regions and, therefore, to a new family of compact stars, denser than the ordinary hybrid stars. In a range of parameters, one may obtain twin hybrid stars (hybrid stars with the same masses but different radii) and even triplets where three stars, with inner cores of nuclear matter, 2SC matter, and CFL matter, respectively, all have the same mass but different radii.

This review focuses on applications of the ideas of superfluidity and superconductivity in neutron stars in a broader context, ranging from the microphysics of pairing in nucleonic superfluids to macroscopic manifestations of superfluidity in pulsars. The exposition of the basics of pairing, vorticity and mutual friction can serve as an introduction to the subject. We also review some topics of recent interest, including the various types of pinning of vortices, glitches, and oscillations in neutron stars containing superfluid phases of baryonic matter.

We study the bulk viscosity of quark matter in the strong coupling regime within the two-flavor Nambu--Jona-Lasinio model. The dispersive effects that lead to non-zero bulk viscosity arise from quark-meson fluctuations above the Mott transition temperature, where meson decay into two quarks is kinematically allowed. We adopt the Kubo-Zubarev formalism and compute the equilibrium imaginary-time correlation function for pressure in the $O(1/N_c)$ power counting scheme. The bulk viscosity of matter is expressed in terms of the Lorentz components of the quark spectral function and includes multi-loop contributions which arise via re-summation of infinite geometrical series of loop diagrams. We show that the multi-loop contributions dominate the single-loop contribution close to the Mott line, whereas at high temperatures the one-loop contribution is dominant. The multi-loop bulk viscosity dominates the shear viscosity close to the Mott temperature by factors 5 to 20, but with increasing temperature the shear viscosity becomes the dominant dissipation mechanism of stresses as the one-loop contribution becomes the main source of bulk viscosity.

The Bose-Einstein condensation of $\alpha$ partciles in the multicomponent environment of dilute, warm nuclear matter is studied. We consider the cases of matter composed of light clusters with mass numbers $A\leq 4$ and matter that in addition these clusters contains $\isotope[56]{Fe}$ nuclei. We apply the quasiparticle gas model which treats clusters as bound states with infinite life-time and binding energies independent of temperature and density. We show that the $\alpha$ particles can form a condensate at low temperature $T\le 2$ MeV in such matter in the first case. When the $\isotope[56]{Fe}$ nucleus is added to the composition the cluster abundances are strongly modified at low temperatures, with an important implication that the $\alpha$ condensation at these temperatures is suppressed.

MeV-scale magnetic fields in the interiors of magnetars suppress the pairing of neutrons and protons in the $S$-wave state. In the case of a neutron condensate the suppression is the consequence of the Pauli-paramagnetism of the neutron gas, i.e., the alignment of the neutron spins along the magnetic field. The proton $S$-wave pairing is suppressed because of the Landau diamagnetic currents of protons induced by the field. The Ginzburg-Landau and BCS theories of the critical magnetic fields for unpairing are reviewed. The macrophysical implications of the suppression (unpairing) of the condensates are discussed for the rotational crust-core coupling in magnetars and the neutrino-dominated cooling era of their thermal evolution.

I review a number of recent developments in the physics of compact stars containing deconfined quark matter, including (a)~their cooling with possible phase transition from a fully gapped to a gapless phase of QCD at low temperatures and large isospin, (b)~the transport coefficients of the 2SC phase and the role played by the Aharonov-Bohm interactions between flux-tubes and unpaired fermions; (c)~rapidly rotating compact stars and spin-down and spin-up induced phase transition between hadronic and QCD matter as well as between different phases of QCD.

We study the electrical conductivity of finite-temperature crust of a warm compact star which may be formed in the aftermath of a supernova explosion or a binary neutron star merger as well as when a cold neutron star is heated by accretion of material from a companion. We focus on the temperature-density regime where plasma is in the liquid state and, therefore, the conductivity is dominated by the electron scattering off correlated nuclei. The dynamical screening of this interaction is implemented in terms of the polarization tensor computed in the hard-thermal-loop effective field theory of QED plasma. The correlations of the background ionic component are accounted for via a structure factor derived from Monte Carlo simulations of one-component plasma. With this input we solve the Boltzmann kinetic equation in relaxation time approximation taking into account the anisotropy of transport due to the magnetic field. The electrical conductivity tensor is studied numerically as a function of temperature and density for carbon and iron nuclei as well as density-dependent composition of zero-temperature dense matter in weak equilibrium with electrons. We also provide accurate fit formulas to our numerical results as well as supplemental tables which can be used in dissipative magneto-hydrodynamics simulations of warm compact stars.

The properties of $\isotope[12]{C}$, $\isotope[16]{O}$, and $\isotope[20]{Ne}$ nuclei in strong magnetic fields $B\simeq 10^{17}\,$G are studied in the context of strongly magnetized neutron stars and white dwarfs. The SKY3D code is extended to incorporate the interaction of nucleons with the magnetic field and is utilized to solve the time-independent Hartree-Fock equations with a Skyrme interaction on a Cartesian three-dimensional grid. The numerical solutions demonstrate a number of phenomena, which include a splitting of the energy levels of spin-up and -down nucleons, spontaneous rearrangment of energy levels in $\isotope[16]{O}$ at a critical field, which leads to jump-like increases of magnetization and proton current in this nucleus, and evolution of the intrinsically deformed $\isotope[20]{Ne}$ nucleus towards a more spherical shape under increasing field strength. Many of the numerical features can be understood within a simple analytical model based on the occupation by the nucleons of the lowest states of the harmonic oscillator in a magnetic field.

Electrical conductivity of finite-temperature plasma in neutron star crusts is studied for applications in magneto-hydrodynamical description of compact stars. We solve the Boltzmann kinetic equation in relaxation time approximation taking into account the anisotropy of transport due to the magnetic field, the effects of dynamical screening in the scattering matrix element and corre- lations among the nuclei. We show that conductivity has a minimum at a non-zero temperature, a low-temperature decrease and a power-law increase with increasing temperature. Selected numerical results are shown for matter composed of carbon, iron, and heavier nuclei present in the outer crusts of neutron star.

We report new simulations of cooling of compact stars containing quark cores and updated fits to the Cas A fast cooling data. Our model is built on the assumption that the transient behaviour of the star in Cas A is due to a phase transition within the dense QCD matter in the core of the star. Specifically, the fast cooling is attributed to an enhancement in the neutrino emission triggered by a transition from a fully gapped, two-flavor, red-green color-superconducting quark condensate to a superconducting crystalline or an alternative gapless, color-superconducting phase. The blue colored condensate is modeled as a Bardeen-Cooper-Schrieffer (BCS)-type color superconductor with spin-one pairing order parameter. We study the sensitivity of the fits to the phase transition temperature, the pairing gap of blue quarks and the time-scale characterizing the phase transition (the latter modelled in terms of a width parameter). Relative variations in these parameter around their best fit values larger than $10^{-3}$ spoil the fit to the data. We confirm the previous finding that the cooling curves show significant variations as a function of compact star mass, which allows one to account for dispersion in the data on the surface temperatures of thermally emitting neutron stars.

Cooling simulations of neutron stars and their comparison with the data from thermally emitting x-ray sources put constraints on the properties of axions, and by extension of any light pseudoscalar dark matter particles, whose existence has been postulated to solve the strong-CP problem of QCD. We incorporate the axion emission by pair-breaking and formation processes by $S$- and $P$-wave nucleonic condensates in a benchmark code for cooling simulations as well as provide fit formulas for the rates of these processes. Axion cooling of neutron stars has been simulated for 24 models covering the mass range 1 to 1.8 solar masses, featuring nonaccreted iron and accreted light-element envelopes, and a range of nucleon-axion couplings. The models are based on an equation state predicting conservative physics of superdense nuclear matter that does not allow for the onset of fast cooling processes induced by phase transitions to non-nucleonic forms of matter or high proton concentration. The cooling tracks in the temperature vs age plane were confronted with the (time-averaged) measured surface temperature of the central compact object in the Cas A supernova remnant as well as surface temperatures of three nearby middle-aged thermally emitting pulsars.We find that the axion coupling is limited to $f_a/10^{7}\textrm{GeV} \ge (5$--$10)$, which translates into an upper bound on axion mass $m_a \le (0.06$--$0.12)~\textrm{eV}$ for Peccei-Quinn charges of the neutron $\vert C_n \vert \sim 0.04$ and proton $\vert C_p \vert \sim 0.4$ characteristic for hadronic models of axions.

If a magnetar interior $B$-field exceeds $10^{15}$~G, it will unpair the proton superconductor in the stellar core by inducing diamagnetic currents that destroy the Cooper pair coherence. Then, the $P$-wave neutron superfluid in these non-superconducting regions will couple to the stellar plasma by scattering of protons off the quasiparticles that are confined in the cores of neutron vortices by the strong (nuclear) force. The dynamical timescales associated with this interaction span from several minutes at the crust-core interface to a few seconds in the deep core. We show that (a) the rapid crust-core coupling is incompatible with oscillation models of magnetars that completely decouple the core superfluid from the crust and (b) magnetar precession is damped by the coupling of normal fluids to the superfluid core and, if observed, needs to be forced or continuously excited by seismic activity.

We obtain the critical magnetic field required for complete destruction of $S$-wave pairing in neutron matter, thereby setting limits on the pairing and superfluidity of neutrons in the crust and outer core of magnetars. We find that for fields $B \ge 10^{17}$ G the neutron fluid is non-superfluid -- if weaker spin-1 superfluidity does not intervene -- a result with profound consequences for the thermal, rotational, and oscillatory behavior of magnetars. Because the dineutron is not bound in vacuum, cold dilute neutron matter cannot exhibit a proper BCS-BEC crossover. Nevertheless, owing to the strongly resonant behavior of the $nn$ interaction at low densities, neutron matter shows a precursor of the BEC state, as manifested in Cooper-pair correlation lengths {being} comparable to the interparticle distance. We make a systematic quantitative study of this type of BCS-BEC crossover in the presence of neutron fluid spin-polarization induced by an ultra-strong magnetic field. We evaluate the Cooper pair wave-function, quasiparticle occupation numbers, and quasiparticle spectra for densities and temperatures spanning the BCS-BEC crossover region. The phase diagram of spin-polarized neutron matter is constructed and explored at different polarizations.

We describe the microphysics, phenomenology, and astrophysical implication of a $B$-field induced unpairing effect that may occur in magnetars, if the local $B$-field in the core of a magnetar exceeds a critical value $H_{c2}$. Using the Ginzburg-Landau theory of superconductivity, we derive the $H_{c2}$ field for proton condensate taking into the correction ($\le 30\%$) which arises from its coupling to the background neutron condensate. The density dependence of pairing of proton condensate implies that $H_{c2}$ is maximal at the crust-core interface and decreases towards the center of the star. As a consequence, magnetar cores with homogenous constant fields will be partially superconducting for "medium-field" magnetars ($10^{15}\le B\le 5 \times 10^{16}$ G) whereas "strong-field" magnetars ($B>5\times 10^{16}$ G) will be void of superconductivity. The neutrino emissivity of a magnetar's core changes in a twofold manner: (i)~the $B$-field assisted direct Urca process is enhanced by orders of magnitude, because of the unpairing effect in regions where $B\ge H_{c2}$; (ii)~the Cooper-pair breaking processes on protons vanish in these regions and the overall emissivity by the pair-breaking processes is reduced by a factor of only a few.

We study the phase diagram of isospin-asymmetrical nuclear matter in the density-temperature plane, allowing for four competing phases of nuclear matter: (i) the unpaired phase, (ii) the translationally and rotationally symmetric, but isospin-asymmetrical BCS condensate, (iii) the current-carrying Larkin-Ovchinnikov-Fulde-Ferrell phase, and (iv) the heterogeneous phase-separated phase. The phase diagram of nuclear matter composed of these phases features two tri-critical points in general, as well as crossovers from the asymmetrical BCS phase to a BEC of deuterons plus a neutron gas, both for the homogeneous superfluid phase (at high temperatures) and for the heterogeneous phase (at low temperatures). The BCS-BEC type crossover in the condensate occurs as the density is reduced. We analyze in detail some intrinsic properties of these phases, including the Cooper-pair wave function, the coherence length, the occupation numbers of majority and minority nucleonic components, and the dispersion relations of quasiparticle excitations about the ground state. We show by explicit examples that the physics of the individual phases and the transition from weak to strong coupling can be well understood by tracing the behavior of these quantities.

Background: The two-flavor color superconducting (2SC) phase of quark matter is a possible constituent of the core of neutron stars. To assess its impact on the observable behavior of the star one must analyze transport properties, which in 2SC matter are controlled by the scattering of gapless fermionic modes by each other and possibly also by color-magnetic flux tubes. Purpose: We determine the electrical and thermal conductivities and the shear viscosity of 2SC matter. Methods: We use a variational formulation of transport theory, treating the strong and electromagnetic interactions via a weak coupling expansion. Results: We provide the leading order scaling of the transport coefficients with temperature and chemical potential as well as accurate fits to our numerical results. We also find that the scattering of fermions by color-magnetic flux tubes is insignificant for thermal conductivity, but may contribute to the electrical conductivity and shear viscosity in the limit of very low temperature or high magnetic field. We also estimate the transport coefficients in unpaired quark matter. Conclusions: Our calculation has set the stage for exploration of possible signatures of the presence of 2SC quark matter in neutron stars.

We construct equilibrium models of compact stars using a realistic equation of state and obtain the density range occupied by the proton superconductor in strong $B$-fields. We do so by combining the density profiles of our models with microscopic calculations of proton pairing gaps and the critical unpairing field $H_{c2}$ above which the proton type-II superconductivity is destroyed. We find that magnetars with interior homogeneous field within the range $0.1 \le B_{16}\le 2$, where $B_{16} = B/10^{16}$ G, are partially superconducting, whereas those with $B_{16} > 2$ are void of superconductivity. We briefly discuss the neutrino emissivity and superfluid dynamics of magnetars in the light of their (non)-superconductivity.

We present a simultaneous calculation of heavy single-$\Lambda$ hypernuclei and compact stars containing hypernuclear core within a relativistic density functional theory based on a Lagrangian which includes the hyperon octet and lightest isoscalar-isovector mesons which couple to baryons with density-dependent couplings. The corresponding density functional allows for SU(6) symmetry breaking and mixing in the isoscalar sector, whereby the departures in the $\sigma$-$\Lambda$ and $\sigma$-$\Sigma$ couplings away from their values implied by the SU(3) symmetric model are used to adjust the theory to the laboratory and astronomical data. We fix $\sigma$-$\Lambda$ coupling using the data on the single-$\Lambda$ hypernuclei and derive an upper bound on the $\sigma$-$\Sigma$ from the requirement that the lower bound on the maximum mass of a compact star is $2 M_{\odot}$.

Starting from equations of state of nucleonic and color-superconducting quark matter and assuming a first-order phase transition between these, we construct an equation of state of stellar matter, which contains three phases: a nucleonic phase, as well as two-flavor and three-flavor color-superconducting phases of quarks. Static sequences of the corresponding hybrid stars include massive members with masses of $\sim 2 M_{\odot}$ and radii in the range of $13< R < 16$ km. We investigate the integral parameters of rapidly rotating stars and obtain evolutionary sequences that correspond to constant rest-mass stars spinning down by electromagnetic and gravitational radiation. Physically new transitional sequences are revealed that are distinguished by a phase transition from nucleonic to color-superconducting matter for some configurations that are located between the static and Keplerian limits. The snapshots of internal structure of the star, displaying the growth or shrinkage of superconducting volume as the star's spin changes, are displayed for constant rest mass stars. We further obtain evolutionary sequences of rotating supramassive compact stars and construct pre-collapse models that can be used as initial data to simulate a collapse of color-superconducting hybrid stars to a black hole.

The recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

We discuss to which extent the modifications of the hyperon--scalar-meson coupling constants affect the equation of state (EoS) hypernuclear matter. The study is carried out within a relativistic density functional theory. The nucleonic matter is described in terms of a density-dependent parametrization of nucleon-meson couplings, whereas the hyperon--meson couplings are deduced from the octet model. We identify the parameter space of hyperon-meson couplings for which massive stellar configurations with $M\le 2.25 M_{\odot}$ exist. We also discuss the EoS at finite temperatures with and without of a trapped neutrino component and show that neutrinos stiffen the EoS and change qualitatively the composition of stellar matter.

We investigate the modification of the pion self-energy at finite temperature due to its interaction with a low-density, isospin-symmetric nuclear medium embedded in a constant magnetic background. To one loop, for fixed temperature and density, we find that the pion effective mass increases with the magnetic field. For the $\pi^{-}$, interestingly, this happens solely due to the trivial Landau quantization shift $\sim |eB|$, since the real part of the self-energy is negative in this case. In a scenario in which other charged particle species are present and undergo an analogous trivial shift, the relevant behavior of the effective mass might be determined essentially by the real part of the self-energy. In this case, we find that the pion mass decreases by $\sim 10%$ for a magnetic field $|eB|\sim m_\pi^2$, which favors pion condensation at high density and low temperatures.

We present a model of the compact star in Cassiopea A that accounts for its unusually fast cooling behavior. This feature is interpreted as an enhancement in the neutrino emission triggered by a transition from a fully gapped, two-flavor, color-superconducting phase to a crystalline phase or an alternative gapless, color-superconducting phase. By fine-tuning a single parameter -- the temperature of this transition -- a specific cooling scenario can be selected that fits the Cas A data. Such a scenario requires a massive $M\sim 2M_{\odot}$ star and is, therefore, distinctive from models invoking canonical 1.4 $M_{\odot}$ mass star with nucleonic pairing alone.

We apply thermal field theory methods to compute microscopically the nucleon self-energy arising from one-pion exchange in isospin-symmetric nuclear matter and neutron matter. A self-consistent numerical scheme is introduced and its convergence is demonstrated. The repulsive contribution from the Fock exchange diagram to the energy per nucleon in symmetric nuclear matter is obtained.

We construct an equation of state of strange quark matter in strong magnetic field within a confining model. The confinement is modeled by means of the Richardson potential for quark-quark interaction modified suitably to account for strong magnetic field. We compare our results for the equation of state and magnetization of matter to those derived within the MIT bag model. The differences between these models arise mainly due to the momentum dependence of the strong interaction between quarks in the Richardson model. Specifically, we find that the magnetization of strange quark matter in this model has much more pronounced de Haas-van Alf\'{v}en oscillations than in the MIT bag model, which is the consequence of the (static) gluon-exchange structure of the confining potential.

We study the equation of state and composition of hypernuclear matter within a relativistic density functional theory with density-dependent couplings. The parameter space of hyperon--scalar-meson couplings is explored by allowing for mixing and breaking of SU(6) symmetry, while keeping the nucleonic couplings constant fixed. The subset of equations of state, which corresponds to small values of hyperon--scalar-meson couplings allows for massive M < 2.25 M_solar compact stars; the radii of hypernuclear stars are within the range 12--14 km. We also study the equation of state of hot neutrino-rich and neutrinoless hypernuclear matter and confirm that neutrinos stiffen the equation of state and dramatically change the composition of matter by keeping the fractions of charged leptons nearly independent of the density prior to the onset of neutrino transparency. We provide piecewise polytropic fits to six representative equations of state of hypernuclear matter, which are suitable for applications in numerical astrophysics.

We study the response of a single-component pair-correlated baryonic Fermi-liquid to density, spin, and their current perturbations. A complete set of response functions is derived in the low-temperature regime both within an effective theory based on a small momentum transfer expansion and within a numerical scheme valid for arbitrary momentum transfers. A comparison of these two approaches validates the perturbative approximation within the domain of its convergence. We derive the spectral functions of collective excitations associated with the density, density-current, spin, and spin-current perturbations. The dispersion relations of density and spin fluctuations are derived and it is shown that the density fluctuations lead to exciton-like undamped bound states, whereas the spin excitations correspond to diffusive modes above the pair-breaking threshold. The contribution of the collective pair-breaking modes to the specific heat of neutron matter at subnuclear densities is computed and is shown to be comparable to that of the degenerate electron gas at not too low temperatures.

Compact stars with strong magnetic fields (magnetars) have been observationally determined to have surface magnetic fields of order of 10^{14}-10^{15} G, the implied internal field strength being several orders larger. We study the equation of state and composition of dense hypernuclear matter in strong magnetic fields in a range expected in the interiors of magnetars. Within the non-linear Boguta-Bodmer-Walecka model we find that the magnetic field has sizable influence on the properties of matter for central magnetic field B \ge 10^{17} G, in particular the matter properties become anisotropic. Moreover, for the central fields B \ge 10^{18} G, the magnetized hypernuclear matter shows instability, which is signaled by the negative sign of the derivative of the pressure parallel to the field with respect to the density, and leads to vanishing parallel pressure at the critical value B_{\rm cr} \simeq 10^{19} G. This limits the range of admissible homogeneously distributed fields in magnetars to fields below the critical value B_{\rm cr}.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

I discuss three topics in physics of massive (two solar-mass and larger) neutron stars containing deconfined quark matter: (i) the equation of state of deconfined dense quark matter and its color superconducting phases, (ii) the thermal evolution of stars with quark cores, (iii) color-magnetic flux tubes in type-II superconducting quark matter and their dynamics driven by Aharonov-Bohm interactions with unpaired fermions.

We report on a comprehensive study of the phase structure of cold, dilute nuclear matter featuring a SD condensate at non-zero isospin asymmetry, within wide ranges of temperatures and densities. We find a rich phase diagram comprising three superfluid phases, namely a Larkin-Ovchinnikov-Fulde-Ferrell phase, the ordinary BCS phase, and a heterogeneous, phase-separated BCS phase, with associated crossovers from the latter two phases to a homogeneous or phase-separated Bose-Einstein condensate of deuterons. The phase diagram contains two tricritical points (one a Lifshitz point), which may degenerate into a single tetracritical point for some degree of isospin asymmetry.

I review three problems in astrophysics of compacts stars: (i) the phase diagram of warm pair-correlated nuclear matter a sub-saturation densities at finite isospin asymmtery; (ii) the Standard Model neutrino emission from superfluid phases in neutron stars within the Landau theory of Fermi (superfluid) liquids; (iii) the beyond Standard Model physics of axionic cooling of compact stars by the Cooper pair-breaking processes.

Once formed in a supernova explosion, a neutron star cools rapidly via neutrino emission during the first 10^4-10^5 yr of its life-time. Here we compute the axion emission rate from baryonic components of a star at temperatures below their respective critical temperatures T_c for normal-superfluid phase transition. The axion production is driven by a charge neutral weak process, associated with Cooper pair breaking and recombination. The requirement that the axion cooling does not overshadow the neutrino cooling puts a lower bound on the axion decay constant f_a > 6 10^9 T_{c9}^{-1} GeV, with T_{c9} = T_c/10^9 K. This translates into a upper bound on the axion mass m_a < 10^{-3} T_{c9} eV.

Neutrino emission rate from baryonic matter in neutron stars via weak neutral vector interaction is computed up to order O(v_F^6), where v_F is the Fermi velocity in units of speed of light. The vector current polarization tensors are evaluated with full vertices which include resummed series in the particle-hole channel. The neutrino emissivity is enhanced compared to the O(v_F^4) order up of 10% for values v_F< 0.4 characteristic to baryons in compact stars.

I discuss the structure and composition of massive (two solar-mass) neutron stars containing hypernuclear and deconfined quark matter in color superconducting states. Stable configurations featuring such matter are obtained if the equation of state of hadronic matter is stiff above the saturation density, the transition to quark matter takes place at a few times the nuclear saturation density, and the repulsive vector interactions in quark matter are substantial. I also discuss our recent progress in understanding the cooling of massive compact stars with color superconducting quark cores.

The recent measurement of a $1.97\pm 0.04$ solar-mass pulsar places a stringent lower bound on the maximum mass of compact stars and therefore challenges the existence of any agents that soften the equation of state of ultra-dense matter. We investigate whether hyperons and/or quark matter can be accommodated in massive compact stars by constructing an equation of state based on a combination of phenomenological relativistic hyper-nuclear density functional and an effective model of quantum chromodynamics (the Nambu-Jona-Lasinio model). Stable configurations are obtained with $M \ge 1.97 M_{\sun}$ featuring hyper-nuclear and quark matter in color superconducting state if the equation of state of nuclear matter is stiff above the saturation density, the transition to quark matter takes place at a few times the nuclear saturation density, and the repulsive vector interactions in quark matter are substantial.

We examine the thermal evolution of a sequence of compact objects containing low-mass hadronic and high-mass quark-hadronic stars constructed from a microscopically motivated equation of state. The dependence of the cooling tracks in the temperature versus age plane is studied on the variations of the gaplessness parameter (the ratio of the pairing gap for red-green quarks to the electron chemical potential) and the magnitude of blue quark gap. The pairing in the red-green channel is modeled assuming an inhomogeneous superconducting phase to avoid tachionic instabilities and anomalies in the specific heat; the blue colored condensate is modeled as a Bardeen-Cooper-Schrieffer (BCS)-type color superconductor. We find that massive stars containing quark matter cool faster in the neutrino-cooling era if one of the colors (blue) is unpaired and/or the remaining colors (red-green) are paired in a inhomogeneous gapless superconducting state. The cooling curves show significant variations along the sequence, as the mass (or the central density) of the models is varied. This feature provides a handle for fine-tuning the models to fit the data on the surface temperatures of same-age neutron stars. In the late-time photon cooling era we observe inversion in the temperature arrangement of models, i.e., stars experiencing fast neutrino cooling are asymptotically hotter than their slowly cooling counterparts.

Magnetohydrodynamics of strongly magnetized relativistic fluids is derived in the ideal and dissipative cases, taking into account the breaking of spatial symmetries by a quantizing magnetic field. A complete set of transport coefficients, consistent with the Curie and Onsager principles, is derived for thermal conduction, as well as shear and bulk viscosities. It is shown that in the most general case the dissipative function contains five shear viscosities, two bulk viscosities, and three thermal conductivity coefficients. We use Zubarev's non-equilibrium statistical operator method to relate these transport coefficients to correlation functions of equilibrium theory. The desired relations emerge at linear order in the expansion of the non-equilibrium statistical operator with respect to the gradients of relevant statistical parameters (temperature, chemical potential, and velocity.) The transport coefficients are cast in a form that can be conveniently computed using equilibrium (imaginary-time) infrared Green's functions defined with respect to the equilibrium statistical operator.

The appearance of quark matter in the centers of compact stars has a number of astrophysical implications. In this contribution I discuss the structure and stability of compact stars, gravitational radiation from non-axisymmetric deformations, and nucleation and dynamics of vortices (flux tubes) in the color superconducting phases.

The phase diagram of charge and color neutral two-flavor color superconducting quark matter is studied including the homogeneous two-flavor superconductor (2SC) and the inhomogeneous Fulde-Ferrell (FF) phases within the Nambu-Jona-Lasinio model. The low-temperature domain T < 5 MeV of the phase diagram contains the FF phase, which borders at high temperatures to the 2SC phase. The critical temperature of phase transition from the 2SC to the unpaired state is in the range 20-30 MeV. We derive the equation of state of matter and its composition and show that matter in mature compact stars should be in the inhomogeneous FF-like superconducting state. We briefly discuss the astrophysical implications of such a phase in compact stars.

We argue that if color-superconducting quark matter exists in the core of a neutron star, it may contain a high density of flux tubes, carrying flux that is mostly color-magnetic, with a small admixture of ordinary magnetic flux. We focus on the two-flavor color-superconducting ("2SC") phase, and assume that the flux tubes are energetically stable, although this has not yet been demonstrated. The density of flux tubes depends on the nature of the transition to the color-superconducting phase, and could be within an order of magnitude of the density of magnetic flux tubes that would be found if the core were superconducting nuclear matter. We calculate the cross-section for Aharonov-Bohm scattering of gapless fermions off the flux tubes, and the associated collision time and frictional force on a moving flux tube. We discuss the other forces on the flux tube, and find that if we take in to account only the forces that arise within the 2SC core region then the timescale for expulsion of the color flux tubes from the 2SC core is of order 10^10 years.

In strong magnetic fields the transport coefficients of strange quark matter become anisotropic. We determine the general form of the complete set of transport coefficients in the presence of a strong magnetic field. By using a local linear response method, we calculate explicitly the bulk viscosities $\zperp$ and $\zpara$ transverse and parallel to the $B$-field respectively, which arise due to the non-leptonic weak processes $u+s\leftrightarrow u+d$. We find that for magnetic fields $B<10^{17}$ G, the dependence of $\zperp$ and $\zpara$ on the field is weak, and they can be approximated by the bulk viscosity for zero magnetic field. For fields $B>10^{18}$ G, the dependence of both $\zperp$ and $\zpara$ on the field is strong, and they exhibit de Haas-van Alphen-type oscillations. With increasing magnetic field, the amplitude of these oscillations increases, which eventually leads to negative $\zperp$ in some regions of parameter space. We show that the change of sign of $\zperp$ signals a hydrodynamic instability. As an application, we discuss the effects of the new bulk viscosities on the r-mode instability in rotating strange quark stars. We find that the instability region in strange quark stars is affected when the magnetic fields exceeds the value $B= 10^{17}$ G. For fields which are larger by an order of magnitude, the instability region is significantly enlarged, making magnetized strange stars more susceptible to $r$-mode instability than their unmagnetized counterparts.

Massive neutron stars may harbor deconfined quark matter in their cores. I review some recent work on the microphysics and the phenomenology of compact stars with cores made of quark matter. This includes the equilibrium and stability of non-rotating and rapidly rotating stars, gravitational radiation from deformations in their quark cores, neutrino radiation and dichotomy of fast and slow cooling, and pulsar radio-timing anomalies.

We present an equation of state and the composition of low-density supernova matter composed of light nuclei with mass number A \le 13. We work within the quasiparticle gas model, which accounts for bound states with decay time scales larger than the relevant time scale of supernova and protoneutron star evolution. The mean-field contribution is included in terms of Skyrme density functional. Deuterons, tritons, and 3H(e) nuclei appear in matter in concentrations that are substantially higher than those of heavier nuclei. We calculate the critical temperature of deuteron condensation in such matter, and demonstrate that the appearance of clusters substantially lowers the critical temperature.

We calculate the phase diagram of two-flavor quark matter in the temperature-flavor asymmetry plane in the case where there are three competing phases: the homogeneous Bardeen-Cooper-Schrieffer (BCS) phase, the unpaired phase, and a phase with broken spatial symmetry, which is here taken to be the counterpart of the Fulde-Ferrell (FF) phase in condensed matter physics. We show that the system belongs to the universality class of paramagnetic-ferromagnetic -helical systems. We vary the coupling constant of the theory, which is obtained from integrating out the gluonic degrees of freedom. In weak coupling, the FF phase is favored at arbitrary flavor asymmetries for sufficiently low temperatures; at intermediate coupling its occupancy domain is shifted towards larger asymmetries. Strong coupling features a new regime of an inhomogeneous FF state, which we identify with a current-carrying Bose-Einstein condensate of tightly bound up and down quarks. We analyze the occupation numbers and the Cooper-pair wave function and show that when the condensate momentum is orthogonal to the particle momentum the minority component contains a blocking region (breach) around the Fermi sphere in the weak-coupling limit, which engulfs more low-momentum modes as the coupling is increased, and eventually leads to a topological change in strong coupling, where the minority Fermi sphere contains either two occupied strips or an empty sphere. For non-orthogonal momenta, the blocking region is either reduced or extinct, i.e., no topological changes are observed.

The interiors of high mass compact (neutron) stars may contain deconfined quark matter in a crystalline color superconducting (CCS) state. On a basis of microscopic nuclear and quark matter equations of states we explore the internal structure of such stars in general relativity. We find that their stable sequence harbors CCS quark cores with masses M_core \le (0.78-0.82)M_{sun} and radii R_core \le 7 km. The CCS quark matter can support nonaxisymmetric deformations, because of its finite shear modulus, and can generate gravitational radiation at twice the rotation frequency of the star. Assuming that the CCS core is maximally strained we compute the maximal quadrupole moment it can sustain. The characteristic strain of gravitational wave emission $h_0$ predicted by our models are compared to the upper limits obtained by the LIGO and GEO 600 detectors. The upper limits are consistent with the breaking strain of CCS matter \sigma \le 10^{-4} and large pairing gaps \Delta \sim 50 MeV, or, alternatively, with \sigma \sim 10^{-3} and small pairing gaps \Delta \sim 15 MeV. An observationally determined value of the characteristic strain h_0 can pin down the product \sigma\Delta^2. On the theoretical side a better understanding of the breaking strain of CCS matter will be needed to predict reliably the level of the deformation of CCS quark core from first principles.

We study the long wavelength shear modes (Tkachenko waves) of triangular lattices of singly quantized vortices in neutron star interiors taking into account the mutual friction between the superfluid and the normal fluid and the shear viscosity of the normal fluid. The set of Tkachenko modes that propagate in the plane orthogonal to the spin vector are weakly damped if the coupling between the superfluid and normal fluid is small. In strong coupling, their oscillation frequencies are lower and are undamped for small and moderate shear viscosities. The periods of these modes are consistent with the observed ~100-1000 day variations in spin of PSR 1828-11.

The three-flavor crystalline color-superconducting (CCS) phase of quantum chromodynamics (QCD) is a candidate phase for the ground state of cold matter at moderate densities above the density of the deconfinement phase transition. Apart from being a superfluid, the CCS phase has properties of a solid, such as a lattice structure and a shear modulus, and hence the ability to sustain multipolar deformations in gravitational equilibrium. We construct equilibrium configurations of hybrid stars composed of nuclear matter at low, and CCS quark matter at high, densities. Phase equilibrium between these phases is possible only for rather stiff equations of state of nuclear matter and large couplings in the effective Nambu--Jona-Lasinio Lagrangian describing the CCS state. We identify a new branch of stable CCS hybrid stars within a broad range of central densities which, depending on the details of the equations of state, either bifurcate from the nuclear sequence of stars when the central density exceeds that of the deconfinement phase transition or form a new family of configurations separated from the purely nuclear sequence by an instability region. The maximum masses of our non-rotating hybrid configurations are consistent with the presently available astronomical bounds. The sequences of hybrid configurations that rotate near the mass-shedding limit are found to be more compact and thus support substantially larger spins than their same mass nuclear counterparts.

At temperatures below the critical temperature of superfluid phase transition baryonic matter emits neutrinos by breaking and recombination of Cooper pairs formed in the condensate. The strong interactions in the nuclear medium modify the weak interaction vertices and the associated neutrino loss rates. We study these modifications non-perturbatively by summing infinite series of particle-hole loops in the S-wave superfluid neutron matter. The pairing and particle-hole interactions in neutron matter are described in the framework of the BCS and Fermi-liquid theories derived from microscopic interactions. Consistent with the f-sum rule, the leading order contribution to the polarization tensor arises at O(q^2) in the small momentum transfer, q, expansion. The associated neutrino emission rate is parametrically suppressed compared to its one-loop counterpart by a factor of the order of 5\times 10^{-3}, the parameter being the baryon recoil in units of temperature.

This article provides a concise review of the problem of neutrino radiation from dense matter. The subjects addressed include quantum kinetic equations for neutrino transport, collision integrals describing neutrino radiation through charged and neutral current interactions, radiation rates from pair-correlated baryonic and color superconducting quark matter.

This chapter provides a review of microscopic theories of pairing in nuclear systems and neutron stars. Special attention is given to the mean-field BCS theory and its extensions to include effects of polarization of the medium and retardation of the interactions. Superfluidity in nuclear systems that exhibit isospin asymmetry is studied. We further address the crossover from the weak-coupling BCS description to the strong-coupling BEC limit in dilute nuclear systems. Finally, within the observational context of rotational anomalies of pulsars, we discuss models of the vortex state in superfluid neutron stars and of the mutual friction between superfluid and normal components, along with the possibility of type-I superconductivity of the proton subsystem.

This review describes the properties of hadronic phases of dense matter in compact stars. The theory is developed within the method of real-time Green's functions and is applied to study of baryonic matter at and above the saturation density. The non-relativistic and covariant theories based on continuum Green's functions and the T-matrix and related approximations to the self-energies are reviewed. The effects of symmetry energy, onset of hyperons and meson condensation on the properties of stellar configurations are demonstrated on specific examples. Neutrino interactions with baryonic matter are introduced within a kinetic theory. We concentrate on the classification, analysis and first principle derivation of neutrino radiation processes from unpaired and superfluid hadronic phases. We then demonstrate how neutrino radiation rates from various microscopic processes affect the macroscopic cooling of neutron stars and how the observed X-ray fluxes from pulsars constrain the properties of dense hadronic matter.

If deconfined quark matter exists inside compact stars, the primary cooling mechanism is neutrino radiation via the direct Urca processes d->u+e+antinu_e and u+e->d+nu_e. Below a critical temperature, T_c, quark matter forms a colour superconductor, one possible manifestation of which is a condensate of <ud> quark Cooper pairs in an electric-charge neutralising background of electrons. We compute the neutrino emission rate from such a phase, including charged pair-breaking and recombination effects, and find that on a material temperature domain below T_c the pairing-induced suppression of the neutrino emission rate is not uniformly exponential. If gapless modes are present in the condensed phase, the emissivity at low temperatures is moderately enhanced above that of completely unpaired matter. The importance of charged current pair-breaking processes for neutrino emission both in the fully gapped and partially gapped regimes is emphasised.

We study the finite temperature-density phase diagram of an attractive fermionic system that supports two-body (dimer) and three-body (trimer) bound states in free space. Using interactions characteristic for nuclear systems, we obtain the critical temperature T_c2 for the superfluid phase transition and the limiting temperature T_c3 for the extinction of trimers. The phase diagram features a Cooper-pair condensate in the high-density, low-temperature domain which, with decreasing density, crosses over to a Bose condensate of strongly bound dimers. The high-temperature, low-density domain is populated by trimers whose binding energy decreases toward the density-temperature domain occupied by the superfluid and vanishes at a critical temperature T_c3 > T_c2.

We obtain the equation of state of interacting alpha matter and the critical temperature of Bose-Einstein condensation of alpha particles within an effective scalar field theory. We start from a non-relativistic model of uniform alpha matter interacting with attractive two-body and repulsive three-body potentials and reformulate this model as a O(2) symmetric scalar \phi^6 field theory with negative quartic and positive sextic interactions. Upon restricting the Matsubara sums, near the temperature of Bose-Einstein condensation, to the zeroth order modes we further obtain an effective classical theory in three spatial dimensions. The phase diagram of the alpha matter is obtained from simulations of this effective field theory on a lattice using local Monte-Carlo algorithms.

Type-I proton superconducting cores of neutron stars break up in a magnetic field into alternating domains of superconducting and normal fluids. We examine two channels of superfluid-normal fluid friction where (i) rotational vortices are decoupled from the non-superconducting domains and the interaction is due to the strong force between protons and neutrons; (ii) the non-superconducting domains are dynamically coupled to the vortices and the vortex motion generates transverse electric fields within them, causing electronic current flow and Ohmic dissipation. The obtained dissipation coefficients are consistent with the Eulerian precession.

Dense matter in compact stars cools efficiently by neutrino emission via the direct Urca processes n -> p + e + nu and p + e -> n + nu. Below the pairing phase transition temperature Tc these processes are suppressed - at asymptotically low temperatures exponentially. We compute the emissivity of the Urca process at one loop for temperatures T < Tc, in the case where the baryons are paired in the 1S0 partial wave. The Urca process is suppressed linearly, rather than exponentially, within the temperature range 0.1 < T/Tc < 1. The charge-current Cooper pair-breaking process contributes up to 50% of the total Urca emissivity in this temperature range.

I discuss the deformed Fermi surface superconductivity (DFS) and some of its alternatives in the context of nucleonic superfluids and two flavor color superconductors that may exist in the densest regions of compact stellar objects.

We develop a model of off-mass-shell pairing correlations in nuclear systems, which is based on the meson-exchange picture of nuclear interactions. The temporal retardations in the model are generated by the Fock-exchange diagrams. The kernel of the complex gap equation for baryons is related to the in-medium spectral function of mesons, which is evaluated nonperturbatively in the random phase approximation. The model is applied to the low-density neutron matter in neutron star crusts by separating the interaction into a long-range one-pion-exchange component and a short-range component parametrized in terms of Landau Fermi liquid parameters. The resulting Eliashberg-type coupled nonlinear integral equations are solved by an iterative procedure.We find that the self-energies extend to off-shell energies of the order of several tens of MeV. At low energies the damping of the neutron pair correlations due to the coupling to the pionic modes is small, but becomes increasingly important as the energy is increased. We discuss an improved quasiclassical approximation under which the numerical solutions are obtained.

This review starts with a discussion of the hierarchy of scales, relevant to the description of superfluids in neutron stars, which motivates a subsequent elementary exposition of the Newtonian superfluid hydrodynamics. Starting from the Euler equations for a superfluid and a normal fluid we apply the tensor virial method to obtain the virial equations of the first, second, and third order and to compute their Eulerian perturbations. Special emphasis is put on the computation of perturbations of the new terms due to mutual gravitational attraction and mutual friction between the two fluids. The oscillation modes of superfluid Maclaurin spheroids are derived from the first and second order perturbed virial equations. We discuss two generic classes of oscillation modes which correspond to the co-moving and relative oscillations of two fluids. These modes decouple if the normal fluid is inviscid. We also discuss the mixing of these modes (when the normal fluid is viscous) and its effect on the dynamical and secular instabilities of the co-moving modes and their damping.

We study the perturbation modes of rotating superfluid ellipsoidal figures of equilibrium in the framework of the two-fluid superfluid hydrodynamics and Newtonian gravity. Our calculations focus on linear perturbations of background equilibria in which the two fluids move together, the total density is uniform, and the densities of the two components are proportional to one another, with ratios that are independent of position. The motions of the two fluids are coupled by mutual friction, as formulated by Khalatnikov. We show that there are two general classes of modes for small perturbations: one class in which the two fluids move together and the other in which there is relative motion between them. The former are identical to the modes found for a single fluid, except that the rate of viscous dissipation, when computed in the secular (or ``low Reynolds number'') approximation under the assumption of a constant kinematic viscosity, is diminished by a factor $f_N$, the fraction of the total mass in the normal fluid. The relative modes are completely new, and are studied in detail for a range of values for the phenomenological mutual friction coefficients, relative densities of the superfluid and normal components, and, for Roche ellipsoids, binary mass ratios. We find that there are no new secular instabilities connected with the relative motions of the two fluid components. Moreover, although the new modes are subject to viscous dissipation (a consequence of viscosity of the normal matter), they do not emit gravitational radiation at all.

We study the neutrino emissivity of strongly magnetized neutron stars due to the charged and neutral current couplings of neutrinos to baryons in strong magnetic fields. The leading order neutral current process is the one-body neutrino-pair bremsstrahlung, which does not have an analogue in the zero field limit. The leading order charged current reaction is the known generalization of the direct Urca processes to strong magnetic fields. While for superstrong magnetic fields in excess of $10^{18}$ G the direct Urca process dominates the one-body bremsstrahlung, we find that for fields on the order $10^{16}-10^{17}$ G and temperatures a few times $10^9$ K the one-body bremsstrahlung is the dominant process. Numerical computation of the resulting emissivity, based on a simple parametrization of the equation of state of the $npe$-matter in a strong magnetic field, shows that the emissivity of this reaction is of the same order of magnitude as that of the modified Urca process in the zero field limit.

We develop a neutrino transport theory, in terms of the real-time non-equilibrium Green's functions, which is applicable to physical conditions arbitrary far from thermal equilibrium. We compute the coherent neutrino radiation in cores of supernovae by evaluating the two-particle-two-hole (2p-2h) polarization function with dressed propagators. The propagator dressing is carried out in the particle-particle channel to all orders in the interaction. We show that at two loops there are two distinct sources of coherence effects in the bremsstrahlung. One is the generically off-shell intermediate state propagation, which leads to the Landau-Pomeranchuk-Migdal type suppression of radiation. We extend previous perturbative results, obtained in the leading order in quasiparticle width, by deriving the exact non-perturbative expression. A new contribution due to off-shell finial/initial baryon states is treated in the leading order in the quasiparticle width. The latter contribution corresponds to processes of higher order than second order in the virial expansion in the number of quasiparticles. At 2p-2h level, the time component of the polarization tensor for the vector transitions vanishes identically in the soft neutrino limit. Vector current thereby is conserved. The contraction of the neutral axial vector current with tensor interaction among the baryons leads to a non-vanishing contribution to the bremsstrahlung rate. These rates are evaluated numerically for finite temperature pure neutron matter at and above the nuclear saturation density.

We calculate the rate of energy radiation by bremsstrahlung of neutrino pairs by baryons in neutron stars employing a transport model in which neutrinos couple to baryons with spectral width. The coherence effects, which are included by computing the self energies with fully dressed propagators, lead to the Landau-Pomeranchuk-Migdal suppression of the neutrino emission rates in the soft neutrino and high temperature limit. A microscopic evaluation of the bremsstrahlung rate by neutrons is carried out using the Brueckner theory of nuclear matter at finite temperatures.

We study the neutron-proton pairing in nuclear matter as a function of isospin asymmetry at finite temperatures and the saturation density using realistic nuclear forces and Brueckner-renormalized single particle spectra. Our computation of the thermodynamic quantities shows that while the difference of the entropies of the superconducting and normal phases anomalously changes its sign as a function of temperature for arbitrary asymmetry, the grand canonical potential does not; the superconducting state is found to be stable in the whole temperature-asymmetry plane. The pairing gap completely disappears for density-asymmetries exceeding $\alpha_c= (n_n-n_p)/n \simeq 0.11$.

We show that the crust-core interface in neutron stars acts as a potential barrier to the peripheral neutron vortices approaching the interface in the model in which these are coupled to the proton vortex clusters. This elementary barrier because of the interaction of vortex magnetic flux with the Meissner currents set up by the crustal magnetic field at the interface. The dominant part of the force is derived from to the cluster-interface interaction. As a result of the stopping of the continuous neutron vortex current through the interface, angular momentum is stored in the superfluid layers in the vicinity of the crust-core interface during the interglitch period. Discontinuous annihilation of proton vortices on the boundary restores the neutron vortex current and spins up the observable crust on short time-scales, leading to a glitch in the spin characteristics of a pulsar.

The impact of various competing heating processes on the thermal evolution of neutron stars is investigated. We show that internal heating leads to significantly enhanced surface temperatures for pulsars of middle and old age. The heating due to thermal creep of pinned vortices and due to outward motion of proton vortices in the interior of the star leads to a better agreement with the observed data in the case of enhanced cooling. The strong pinning models are ruled out by a comparison with the cooling data on the old pulsars. For millisecond pulsars, the heating due to thermal creep of pinned vortices and chemical heating of the core have the largest impact on the surface temperatures. The angular dependence of the heating rates require two dimensional cooling simulations in general. Such a simulation is performed for a selected case in order to check the applicability of one-dimensional codes used in the past.

We consider the precession of isolated neutron stars in which superfluid is not pinned to the stellar crust perfectly. In the case of perfect pinning, Shaham (1977) showed that there are no slowly oscillatory, long-lived modes. When the assumption of perfect pinning is relaxed, new modes are found that can be long-lived, but are expected to be damped rather than oscillatory, unless the drag force on moving superfluid vortex lines has a substantial component perpendicular to the direction of relative motion. The response of a neutron star to external torques, such as the spindown torque, is also treated. We find that when computing the response of a star to perturbations, assuming perfect coupling of superfluid to normal matter from the start can miss some effects.

We derive the transport relaxation times for quasiparticle-vortex scattering processes via nuclear force, relevant for the damping of differential rotation of superfluids in the quantum liquid core of a neutron star. The proton scattering off the neutron vortices provides the dominant resistive force on the vortex lattice at all relevant temperatures in the phase where neutrons only are in the paired state. If protons are superconducting, a small fraction of hyperons and resonances in the normal state would be the dominant source of friction on neutron and proton vortex lattices at the core temperatures $T\ge 10^{7}$ K.

The phonon spectrum of Coulomb lattice in neutron star crusts above the neutron drip density is affected by the interaction with the ambient neutron Fermi-liquid. For the values of the neutron-phonon coupling constant in the range $0.1 \le \lambda \le 1$ an appreciable renormalization of the phonon spectrum occurs which can lead to a lattice instability manifested in an exponential growth of the density fluctuations. The BCS phonon exchange mechanism of superconductivity leads to neutron pairing with a gap in the neutron excitation spectrum comparable to that due to the direct nuclear interaction.

We discuss the spin evolution of pulsars in the case where a superfluid component of the star is coupled to the observable crust on long, spindown timescales. The momentum transfer from the superfluid interior results in an apparent decay of the external torque and, after a dramatic increase, to an asymptotic decrease of the generic value of the braking index, e.g. $n=3$, to values $n\sim 2.5$ if the magnetic field of the star does not decay over its lifetime. In the case where an exponential decay of the magnetic field towards a residual value occurs, the star undergoes a spin-up phase after which it could emerge in the millisecond sector of the $P$-$\dot P$ diagram.

The magnetic structure of neutron vortices in the superfluid cores of neutron stars is determined assuming that the proton condensate forms a type I superconductor. It is shown that the entrainment currents induced by the neutron vortex circulation cause the proton superconductor to break into successive domains of normal and superconducting regions. The Gibbs free-energy is found in the case in which the normal domains form cylindrical tubes coaxial with the neutron vortex. The minimum of the energy functional corresponds to a tube radius $a\sim 0.1-0.5 ~b$, where $b$ is the outer radius of the neutron vortex. The magnetic field within the tube is of the order of $ 5 \times 10^{14}$ G.

We show that the neutron star's crust-core interface acts as a potential barrier on the peripheral neutron vortices approaching this interface and thus prevents their continuous decay on it in the course of equilibrium state deceleration. The barrier arises due to the interaction of vortex magnetic flux with the Meissner currents set up by the crustal magnetic field at the interface. When the non-balanced part of the Magnus force reaches the value at which the vortices are able to annihilate at the interface, the rapid transfer of angular momentum from the superfluid spins-up the observable crust on short dynamical coupling times.

Neutrino emission from the dense hadronic component in neutron stars is subject to strong modifications due to collective effects in the nuclear medium. We implement new estimates of the neutrino emissivities of two processes operating in the nuclear medium into numerical cooling simulations of neutron stars. The first process is the modified Urca process, for which the softening of the pion exchange mode and other polarization effects as well as the neutrino emission arising from the intermediate reaction states are taken into account. The second process concerns neutrino emission through superfluid pair breaking and formation processes. It is found that the medium effects on the emissivity of the modified Urca process result in a strong density dependence, which gives a smooth crossover from the standard to the nonstandard cooling scenario for increasing star masses. For superfluid stars, the superfluid pair breaking and formation processes accelerate mildly both the standard and the nonstandard cooling scenario. This leads to a good agreement between the theoretical cooling tracks and the rather low temperatures observed for objects like PSRs 0833-45 (Vela), 0656+14, and 0630+18 (Geminga). The robustness of our findings against variations in both the underlying equation of state of baryonic matter and the used fast cooling processes is demonstrated. Hence we conclude that the two recalculated neutrino emissivities studied here enable one to reproduce theoretically most of the observed pulsar temperatures by varying the masses of neutron star models.