The effect of pasta phases on the quark-hadron phase transition is investigated for a set of relativistic mean-field equations of state for both hadron and quark matter. The results of the full numerical solution with pasta phases are compared with those of an interpolating construction used in previous works, for which we demonstrate an adequate description of the numerical results. A one-to-one mapping of the free parameter of the construction to the physical surface tension of the quark-hadron interface is obtained for which a fit formula is given. For each pair of quark and hadron matter models the critical value of the surface tension is determined, above which the phase transition becomes close to the Maxwell construction. This result agrees well with earlier theoretical estimates. The study is extended to neutron star matter in beta equilibrium with electrons and muons and is applied to investigate the effect of pasta phases on the structure of hybrid compact stars and the robustness of a possible third family solution.

Within the relativistic mean-field framework with hadron masses and coupling constants dependent on the mean scalar field we study properties of nuclear matter at finite temperatures, baryon densities and isospin asymmetries relevant for heavy-ion collisions at laboratory energies below 2$A$ GeV. Previously constructed (KVORcut-based and MKVOR-based) models for the description of the cold hadron matter, which differ mainly by the density dependence of the nucleon effective mass and symmetry energy, are extended for finite temperatures. The baryon equation of state, which includes nucleons and $\Delta$ resonances is supplemented by the contribution of the pion gas described either by the vacuum dispersion relation or with taking into account the $s$-wave pion-baryon interaction. Distribution of the charge between components is found. Thermodynamical characteristics on $T-n$ plane are considered. The energy-density and entropy-density isotherms are constructed and a dynamical trajectory of the hadron system formed in heavy-ion collisions is described. The effects of taking into account the $\Delta$ isobars and the $s$-wave pion-nucleon interaction on pion differential cross sections, pion to proton and $\pi^-/\pi^+$ ratios are studied. The liquid-gas first-order phase transition is studied within the same models in isospin-symmetric and asymmetric systems. We demonstrate that our models yield thermodynamic characteristics of the phase transition compatible with available experimental results. In addition, we discuss the scaled variance of baryon and electric charge in the phase transition region. Effect of the non-zero surface tension on spatial redistribution of the electric charge is considered for a possible application to heavy-ion collisions at low energies.

We consider the model of composite dark matter assuming stable particles of charge $-2$ bound with primordial helium nuclei by Coulomb force in $O$He atoms. We study capture of such dark atoms in matter and propose a possibility of existence of stable $O$-enriched superheavy nuclei and $O$-nuclearites, in which heavy $O$-dark matter fermions are bound by electromagnetic forces with ordinary nuclear matter. $O$He atoms accumulation in stars and its possible effect in stellar evolution is also considered, extending the set of indirect probes for composite dark matter.

We consider a model of $O$He atomic dark matter formed by Coulomb binding of the hypothetical stable double-charged massive $O^{--}$ particles with nuclei of primordial helium. Such a dark matter can be captured by ordinary matter forming superheavy nuclei. We propose a possibility of self-bound by nuclear and electromagnetic interaction $O$-nuclearites as well as more massive gravitating ones and discuss effect of accumulation of $O$He atoms in stars and their effect in stellar evolution.

We demonstrate that the existing neutron-star cooling data can be appropriately described within "the nuclear medium cooling scenario" including hyperons under the assumption that different sources have different masses. We use a stiff equation of state of the relativistic mean-field model MKVORH$\phi$ with hadron effective couplings and masses dependent on the scalar field. It fulfills a large number of experimental constraints on the equation of state of the nuclear matter including the $2\,M_{\odot}$ lower bound for the maximum predicted neutron-star mass and the constraint for the pressure from the heavy-ion particle flow. We select appropriate $^1S_0$ proton and $\Lambda$ hyperon pairing gap profiles from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. The $^3P_2$ neutron pairing gap is assumed to be negligibly small in our scenario. The possibility of the pion, kaon and charged $\rho$-meson condensations is for simplicity suppressed. The resulting cooling curves prove to be sensitive to the value and the density dependence of the pp pairing gap and rather insensitive to the values of the $^1S_0$ neutron pairing gaps.

We investigate the robustness of third family solutions for hybrid compact stars with a quark matter core that correspond to the occurrence of high-mass twin stars against a softening of the phase transition by means of a construction that mimics the effects of pasta structures in the mixed phase. We consider a class of hybrid equations of state that exploits a relativistic mean-field model for the hadronic as well as for the quark matter phase. We present parametrizations that correspond to branches of high-mass twin star pairs with maximum masses between $2.05~M_\odot$ and $1.48~M_\odot$ having radius differences between 3.2 km and 1.5 km, respectively. When compared to a Maxwell construction with a fixed value of critical pressure $P_c$ the effect of the mixed phase construction consists in the occurrence of a region of pressures around $P_c$ belonging to the coexistence of hadronic and quark matter phases between the onset pressure at $P_H$ and the end of the transition at $P_Q$. The maximum broadening which would still allow mass twin compact stars is found to be $(P_Q-P_H)_{\rm max} \approx P_c$ for all parametrizations within the present class of models. At least the heavier of the neutron stars of the binary merger GW170817 could have been a member of the third family of hybrid stars.

The existence of high mass ($\sim 2M_{\odot}$) pulsars PSR J1614-2230 and PSR J0348-0432 requires the compact star matter to be described by a stiff equation of state (EoS). Presence of hyperons in neutron stars leads to a softening of the EoS that results in a decrease of the maximum neutron-star mass below the measured values of masses for PSR J1614-2230 and PSR J0348-0432 pulsars if one exploits ordinary relativistic mean-field (RMF) models (hyperon puzzle). However, within an RMF EoS with a sigma-scaled hadron effective masses and coupling constants, the maximum neutron-star mass remains above $2M_{\odot}$ even when hyperons are included. Also, other important constraints on the equation of state, e.g. the flow constraint from heavy-ion collisions are to be fulfilled. We demonstrate how a satisfactory explanation of all existing observational data for the temperature-age relation is reached within the "nuclear medium cooling" scenario with a relativistic-mean-field EoS with a sigma-scaled hadron effective masses and coupling constants including hyperons.

We extend relativistic mean-field models with hadron masses and meson-baryon coupling constants dependent on the scalar field $\sigma$, including hyperons and $\Delta(1232)$ baryons, to incorporate a possibility of the charged $\rho$ meson condensation in neutron star matter. The influence of the $\rho^-$ condensation on the equation of state proves to be strongly model dependent. In our models of one type (KVORcut-based ones) the $\rho^-$ condensation arises by a second-order phase transition above a critical density and the maximum value of the neutron star mass diminishes only slightly. The matter composition changes more significantly. In our models of other type (MKVOR*-based ones), if the system is considered at fixed density, the $\rho^-$ condensation arises by a second-order phase transition at the baryon density $n=n_{c,\rho}^{(\rm II)}$ and at a slightly higher density $n=n_{c,\rho}^{(\rm I)}$ there occurs a first-order phase transition. In a neutron star matter starting with a density $n<n_{c,\rho}^{(\rm II)}$ there appears a region of a mixed phase, or the system is described by Maxwell construction, that results in a substantial decrease of the value of the maximum neutron star mass. Nevertheless in the models under consideration the observational constraint on the maximum neutron star mass is fulfilled. Besides, in MKVOR*-based models the appearance of the $\rho^-$ condensate is accompanied by a strong rearrangement of the matter composition. Dependence of the results on a choice of the $\rho$ meson scaling functions for the effective $\rho$ meson mass and coupling constants is also investigated.

Recent evidence for high masses ($\sim 2~ M_\odot$) pulsars PSR J1614-2230 and PSR J0348-0432 requires neutron star matter to have a stiff equation of state (EoS). The thermal evolution of compact stars (CS) with stiff hadronic EoS necessitates the application of the "nuclear medium cooling" scenario with a selection of appropriate proton gap profiles together with in-medium effects (like pion softening) on cooling mechanisms in order to achieve a satisfactory explanation of all existing observational data for the temperature-age relation of CS. Here we focus on two examples from \cite{Grigorian:2016leu} for a stiff hadronic EoS without (DD2 EoS) and with (DD2vex) excluded volume correction.

A possibility of the condensation of excitations with a non-zero momentum in rectilinearly moving and rotating superfluid bosonic and fermionic (with Cooper pairing) media is considered in terms of a phenomenological order-parameter functional at zero and non-zero temperature. The results might be applicable to the description of bosonic systems like superfluid $^4$He, ultracold atomic Bose gases, charged pion and kaon condensates in rotating neutron stars, and various superconducting fermionic systems with pairing, like proton and color-superconducting components in compact stars, metallic superconductors, and neutral fermionic systems with pairing, like the neutron component in compact stars and ultracold atomic Fermi gases. Order parameters of the "mother" condensate in the superfluid and the new condensate of excitations, corresponding energy gains, critical temperatures and critical velocities are found.

Knowledge of the equation of state of the baryon matter plays a decisive role in the description of neutron stars. With an increase of the baryon density the filling of Fermi seas of hyperons and $\Delta$ isobars becomes possible. Their inclusion into standard relativistic mean-field models results in a strong softening of the equation of state and a lowering of the maximum neutron star mass below the measured values. We extend a relativistic mean-field model with scaled hadron masses and coupling constants developed in our previous works and take into account now not only hyperons but also the $\Delta$ isobars. We analyze available empirical information to put constraints on coupling constants of $\Delta$s to mesonic mean fields. We show that the resulting equation of state satisfies majority of presently known experimental constraints.

We extend the relativistic mean-field models with hadron masses and meson-baryon coupling constants dependent on the scalar $\sigma$ field, studied previously to incorporate $\Delta(1232)$ baryons. Available empirical information is analyzed to put constraints on the couplings of $\Delta$s with meson fields. Conditions for the appearance of $\Delta$s are studied. We demonstrate that with inclusion of the $\Delta$s our equations of state continue to fulfill majority of known empirical constraints including the pressure-density constraint from heavy-ion collisions, the constraint on the maximum mass of the neutron stars, the direct Urca and the gravitational-baryon mass ratio constraints.

An equation of state of cold nuclear matter with an arbitrary isotopic composition is studied within a relativistic mean-field approach with hadron masses and coupling constants depending self-consistently on the scalar mean-field. All hadron masses decrease universally with the scalar field growth, whereas meson-nucleon coupling constants can vary differently. More specifically we focus on two modifications of the KVOR model studied previously. One extension of the model (KVORcut) demonstrates that the equation of state stiffens if the increase of the scalar-field magnitude with the density is bounded from above at some value for baryon densities above the saturation nuclear density. This can be realized if the nucleon vector-meson coupling constant changes rapidly as a function of the scalar field slightly above the desired value. The other version of the model (MKVOR) utilizes a smaller value of the nucleon effective mass at the nuclear saturation density and a saturation of the scalar field in the isospin asymmetric matter induced by a strong variation of the nucleon isovector-meson coupling constant as function of the scalar field. A possibility of hyperonization of the matter in neutron star interiors is incorporated. Our equations of state fulfill majority of known empirical constraints including the pressure-density constraint from heavy-ion collisions, direct Urca constraint, gravitational-baryon mass constraint for the pulsar J0737-3039B, and the constraint on the maximum mass of the neutron stars.

Measurements demonstrate the existence of compact stars with masses in a broad range from 1.2 to 2 $M_\odot$. The most massive of these objects might be hybrid stars. To fulfill the constraint $M_{\rm max}>2~M_{\odot}$ with a reserve we exploit the stiff DD2 hadronic equation of state (EoS) without and with excluded volume (DD2vex) correction, which produce maximum neutron star masses of $M_{\rm max} = 2.43~M_{\odot}$ and $2.70~M_{\odot}$, respectively. We show that the stiffness of the EoS does not preclude an explanation of the whole set of cooling data within "nuclear medium cooling" scenario for compact stars by a variation of the star masses. We select appropriate proton gap profiles from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. However, we suppress the possibility of pion condensation. In general, the stiffer the EoS the steeper a decrease with density of the effective pion gap is required. Results are compared with previously obtained ones for the HDD EoS for which $M_{\rm max} =2.06~M_{\odot}$. The cooling of the compact star in the supernova remnant Cassiopeia A (Cas A) is explained mainly by an efficient medium modified Urca process. With the DD2vex EoS and using an effective pion gap steeper decreasing with the density and/or a proton gap shifted to smaller densities we are also able to reproduce both a strong decline compatible with ACIS-S data and HRC-S instrument data. The mass of Cas A is estimated as $1.6 - 1.9~M_{\odot}$, above the value $1.5~M_{\odot}$, which we have evaluated with the softer HDD equation of state. Different mass choices for the hottest object XMMU J173203.3-344518 are discussed. We make general remarks also on hybrid star cooling and on its dependence on the stiffness of the hadronic EoS.

We study relativistic mean-field (RMF) models including nucleons interacting with scalar, vector and iso-vector mean fields and self- and cross- mean-field interaction terms. Usually, in such a models the magnitude of the scalar field increases monotonically with the nucleon density, and the nucleon effective mass decreases. We demonstrate that the latter quantity stops to decrease and the equation of state stiffens, provided the mean-field self-interaction potential rises sharply in a narrow vicinity of the values of mean fields corresponding to nucleon densities $n> n_{*}>n_0$, where $n_0$ is the nuclear saturation density. As the result the limiting neutron star mass increases. This procedure offers a simple way to stiffen the equation of state at densities above $n_{*}$ without altering it at densities $n\le n_{0}$. The developed scheme allows an application to neutron stars of the RMF models, which are well fitted to finite nuclei but do not fulfill the experimental constraint on the limiting neutron star mass. The exemplary application of the method to the well-known FSUGold model allows us to increase the limiting neutron star mass from $1.72~M_\odot$ to $M \geq 2.01~M_\odot$.

The equation of state of cold baryonic matter is studied within a relativistic mean-field model with hadron masses and coupling constants depending on a scalar field. We demonstrate that if the effective nucleon mass stops to decrease with a density increase at densities $n>n_*>n_0$, where $n_0$ is the nuclear saturation density, the equation of state stiffens for these densities and the limiting neutron star mass increases. The stabilization of the nucleon mass can be realised if in the equation of motion for the scalar mean-field there appear a term sharply varying in a narrow vicinity of the field value corresponding to the density $n_*$. We show several possible realizations of this mechanism getting sufficiently stiff equations of state. The appearance of hyperons in dense neutron star interiors is accounted for. The obtained equations of state remain sufficiently stiff if the reduction of the $\phi$ meson mass is incorporated. Thereby, the hyperon puzzle can be resolved.

The equation of state of cold baryonic matter is studied within a relativistic mean-field model with hadron masses and coupling constants depending on the scalar field. All hadron masses undergo a universal scaling, whereas the coupling constants are scaled differently. The appearance of hyperons in dense neutron star interiors is accounted for, however the equation of state remains sufficiently stiff if a reduction of the $\phi$ meson mass is included. Our equation of state matches well the constraints known from analyses of the astrophysical data and the particle production in heavy-ion collisions.

The recent measurements of the masses of the pulsar J00737-3039B and of the companion J1756-2251 and pulsars PSR J1614-2230, PSR J0348-0432 demonstrate the existence of compact stars with masses in a broad range from 1.2 to 2 $M_\odot$. To fulfill the constraint $M_{\rm max}>2M_{\odot}$ and to demonstrate the possibility of cooling scenarios for purely hadronic and further for hybrid stars we exploit the stiff DD2 hadronic equation of state producing a maximum neutron star mass $M\simeq 2.43 M_{\odot}$. We show that the "nuclear medium cooling" scenario for neutron stars comfortably explains the whole set of cooling curves just by a variation of the star masses without the necessity for the occurrence of the direct Urca reaction. To describe the cooling data with the very stiff DD2 equation of state we select a proton gap profile from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. Fast cooling of young neutron stars like it is seen in the data for Cas A is explained with the DD2 equation of state when the following conditions are provided: the presence of an efficient medium modified Urca process, and a large proton gap at densities $n\le 2n_0$ vanishing for $n\ge (2.5 - 3) n_0$, where $n_0$ is the saturation nuclear density.

Within the "nuclear medium cooling" scenario of neutron stars all reliably known temperature - age data, including those of the central compact objects in the supernova remnants of Cassiopeia A and XMMU-J1732, can be comfortably explained by a set of cooling curves obtained by variation of the star mass within the range of typical observed masses. The recent measurements of the high masses of the pulsars PSR J1614-2230 and PSR J0348-0432 on the one hand, and of the low masses for PSR J0737-3039B and the companion of PSR J1756-2251 on the other, provide independent proof for the existence of neutron stars with masses in a broad range from $\sim 1.2$ to 2 $M_\odot$. The values $M>2 M_{\odot}$ call for sufficiently stiff equations of state for neutron star matter. We investigate the response of the set of neutron star cooling curves to a stiffening of the nuclear equation of state so that maximum masses of about $2.4 M_\odot$ would be accessible and to a deconfinement phase transition from such stiff nuclear matter in the outer core to color superconducting quark matter in the inner core. Without readjustment of cooling inputs the mass range required to cover all cooling data for the stiff DD2 equation of state should include masses of $2.426 M_\odot$ for describing the fast cooling of CasA while the existence of a quark matter core accelerates the cooling so that CasA cooling data are described with a hybrid star of mass $1.674 M_\odot$.

A possibility of the condensation of excitations with non-zero momentum in moving superfluid media is considered in terms of the Ginzburg-Landau model. The results might be applicable to the superfluid $^4$He, ultracold atomic Bose gases, various superconducting and neutral systems with pairing, like ultracold atomic Fermi gases and the neutron component in compact stars. The order parameters, the energy gain, and critical velocities are found.

We demonstrate that stability of r-modes in young rapidly rotating pulsars might be explained if one takes into account strong medium modifications of the nucleon-nucleon interaction because of the softening of pionic degrees of freedom in dense nucleon matter. Presence of the efficient direct Urca processes is not required. Within our model the most rapidly rotating observed young pulsar PSR J0537-6910 should have the mass $\geq 1.8M_{\odot}$.

Recently, Elshamounty et al. performed a reanalysis of the surface temperature of the neutron star in the supernova remnant Cassiopeia A on the basis of Chandra data measured during last decade, and added a new data point. We show that all reliably known temperature data of neutron stars including those belonging to Cassiopea A can be comfortably explained in our "nuclear medium cooling" scenario of neutron stars. The cooling rates account for medium-modified one-pion exchange in dense matter, polarization effects in the pair-breaking-formation processes operating on superfluid neutrons and protons paired in the 1S_0 state, and other relevant processes. The emissivity of the pair-breaking-formation process in the 3P_2 state is a tiny quantity within our scenario. Crucial for a successful description of the Cassiopeia A cooling proves to be the thermal conductivity from both, the electrons and nucleons, being reduced by medium effects. Moreover, we exploit an EoS which stiffens at high densities due to an excluded volume effect and is capable of describing a maximum mass of 2.1 M_sun, thus including the recent measurements of PSR J1614-2230 and PSR J0348+0432.

A role of a particle-particle p-wave spin interaction in Fermi liquids with s-wave pairing is studied. Depending on the sign of the interaction there arises either the new exciton collective mode below the pair-breaking threshold or the diffusive excitation mode above the threshold. The Landau parameters which control the interaction strength are evaluated for various systems: the dilute fermion gases, degenerate electron liquid, metals, atomic nuclei and neutron matter. The interaction removes also the square-root singularity in the phase space of pair breaking processes. It is shown how these effects influence the neutrino emissivity in the neutron Cooper-pair recombinations in neutron stars.

We demonstrate that the high-quality cooling data observed for the young neutron star in the supernova remnant Cassiopeia A over the past 10 years--as well as all other reliably known temperature data of neutron stars--can be comfortably explained within the "nuclear medium cooling" scenario. The cooling rates of this scenario account for medium-modified one-pion exchange in dense matter and polarization effects in the pair-breaking formations of superfluid neutrons and protons. Crucial for the successful description of the observed data is a substantial reduction of the thermal conductivity, resulting from a suppression of both the electron and nucleon contributions to it by medium effects. We also find that possibly in as little as about ten years of continued observation, the data may tell whether or not fast cooling processes are active in this neutron star.

The Larkin-Migdal approach to a cold superfluid Fermi liquid is generalized for a non-equilibrium system. The Schwinger-Keldysh diagram technique is applied. The developed formalism is applicable to the pairing in the states with arbitrary angular momenta. We consider the white body radiation problem by calculating probabilities of different direct reactions from a piece of a fermion superfluid. The closed diagram technique is formulated in terms of the full Green's functions for systems with the pairing correlation. The cutting rules are used to classify the diagrams representing one-nucleon, two-nucleon, etc. processes in the matter. The important role of multi-piece diagrams for the vector-current conservation is demonstrated. In the case of equilibrated systems, dealing with dressed Green's functions, we demonstrate correspondence between calculations in the Schwinger-Kadanoff-Baym-Keldysh formalism and the ordinary Matsubara technique. As an example we consider neutrino radiation from the neutron pair breaking and formation processes in case of a singlet pairing. Necessary correlation effects are included. The in-medium renormalization of normal and anomalous vertices is performed.

Neutrino emission in processes of breaking and formation of nucleon Cooper pairs is calculated in the framework of the Larkin-Migdal and the Leggett approaches to the description of superfluid Fermi liquids at finite temperatures. We explain peculiarities of both approaches and explicitly demonstrate that they lead to the same expression for the emissivity in pair breaking and formation processes.

Solutions of hydrodynamical equations are presented for the equation of state of the Var der Waals type allowing for the first order phase transition. Attention is focused on description of the hadron-quark phase transition in heavy ion collisions. It is shown that fluctuations dissolve and grow as if the fluid is effectively very viscous. Even in spinodal region germs are growing slowly due to viscosity and critical slowing down. This prevents enhancement of fluctuations in the near-critical region, which is frequently considered as a signal of the critical point in heavy ion collisions.

Neutrino emission in processes of breaking and formation of neutron and proton Cooper pairs is calculated within the Larkin-Migdal-Leggett approach for a superfluid Fermi liquid. We demonstrate explicitly that the Fermi-liquid renormalization respects the Ward identity and assures the weak vector current conservation. The systematic expansion of the emissivities for small temperatures and nucleon Fermi velocity, v_{F,i}, i=n,p, is performed. Both neutron and proton processes are mainly controlled by the axial-vector current contributions, which are not strongly changed in the superfluid matter. Thus, compared to earlier calculations the total emissivity of processes on neutrons paired in the 1S_0 state is suppressed by a factor ~(0.9-1.2) v_{F,n}^2. A similar suppression factor (~v_{F,p}^2) arises for processes on protons.

A new scheme for testing nuclear matter equations of state (EsoS) at high densities using constraints from neutron star phenomenology and a flow data analysis of heavy-ion collisions is suggested. An acceptable EoS shall not allow the direct Urca process to occur in neutron stars with masses below $1.5~M_{\odot}$, and also shall not contradict flow and kaon production data of heavy-ion collisions. Compact star constraints include the mass measurements of 2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807, of 2.0 +/- 0.1 M_sun from the innermost stable circular orbit for 4U 1636-536, the baryon mass - gravitational mass relationships from Pulsar B in J0737-3039 and the mass-radius relationships from quasiperiodic brightness oscillations in 4U 0614+09 and from the thermal emission of RX J1856-3754. This scheme is applied to a set of relativistic EsoS constrained otherwise from nuclear matter saturation properties with the result that no EoS can satisfy all constraints simultaneously, but those with density-dependent masses and coupling constants appear most promising.

We consider the early cooling evolution of strongly magnetized strange stars in a CFL phase with high gap Delta > 100 MeV. We demonstrate how this model may explain main features of the gamma-ray burst phenomena and also yield a strong star kick. The mechanism is based on beaming of neutrino emission along the magnetic vortex lines. We show that for sufficiently high initial temperatures T_0 ~ 30 - 60 MeV and surface magnetic fields B_s ~ 10^{15} - 10^{17} G, the energy release within the narrow beam is up to 10^{52} erg with a magnetic field dependent time scale between 10^{-2} s (for a smaller magnetic field) to 10 s. The above mechanism together with the parity violation of the neutrino-producing weak interaction processes in a magnetic field allow for the strange star kick. The higher the magnetic field the larger is the stars kick velocity. These velocities may cover the same range as observed pulsar kick velocities.

We study the dependence of the cooling of isolated neutron stars on the magnitude of the $3P_2$ neutron gap. It is demonstrated that our ``nuclear medium cooling'' scenario is in favor of a suppressed value of the $3P_2$ neutron gap.

A QCD-motivated dynamical-quasiparticle model with parameters adjusted to reproduce the lattice-QCD equation of state is extrapolated from region of high temperatures and moderate baryonic densities to the domain of high baryonic densities and zero temperature. The resulting equation of state matched with realistic hadronic equations of state predicts a phase transition into the quark phase at higher densities than those reachable in neutron star interiors. This excludes the possibility of the existence of hybrid (hadron-quark) stars. Pure quark stars are possible and have low masses, small radii and very high central densities. Similar results are obtained for a simple bag model with massive quarks, fitted to reproduce the same lattice results. Self-bound quark matter is also excluded within these models. Uncertainties in the present extrapolation re discussed. Comparison with standard bag models is made.

We study the dependence of the cooling of isolated neutron stars on the magnitude of the $3P_2$ neutron gap. It is demonstrated that our ``Nuclear medium cooling scenario'' is in favor of a suppressed value of the $3P_2$ neutron gap.

We show that within a recently developed nonlocal chiral quark model the critical density for a phase transition to color superconducting quark matter under neutron star conditions can be low enough for these phases to occur in compact star configurations with masses below 1.3 M_solar. We study the cooling of these objects in isolation for different values of the gravitational mass. Our equation of state (EoS) allows for 2SC quark matter with a large quark gap \~100 MeV for u and d quarks of two colors that coexists with normal quark matter within a mixed phase in the hybrid star interior. We argue that, if the phases with unpaired quarks were allowed, the corresponding hybrid stars would cool too fast. If they occured for M < 1.3 M_solar, as it follows from our EoS, one could not appropriately describe the neutron star cooling data existing by today. We discuss a "2SC+X" phase, as a possibility to have all quarks paired in two-flavor quark matter under neutron star constraints, where the X-gap is of the order of 10 keV - 1 MeV. Density independent gaps do not allow to fit the cooling data. Only the presence of an X-gap that decreases with increase of the density could allow to appropriately fit the data in a similar compact star mass interval to that following from a purely hadronic model. This scenario is suggested as an alternative explanation of the cooling data in the framework of a hybrid star model.

We study relativistic mean-field models with hadron masses and coupling constants depending self-consistently on a scalar meson field. We demonstrate that by the field redefinition some models can be equivalently transformed into each other. Thereby the large variety of scaling functions for masses and couplings can be reduced to a restricted set of functions with a constrained dependence on a scalar field. We show how by choosing properly the latter scaling functions one may stiffen or soften the equation of state at high densities and simultaneously increase the threshold density for the direct Urca process without any change of the description of nuclear matter close to the saturation density. The stiffening of the equation of state might be motivated by recent neutron star mass measurements, whereas the increase of the threshold density for the direct Urca process is required by the analysis of neutron star cooling data. We demonstrate that if a rho meson is included in a mean-field model as a non-Abelian gauge boson, then there is a possibility for a charged rho-meson condensation in dense nuclear matter. We show that such a novel phase can be realized in neutron star interiors already for sufficiently low densities, typically $\sim 3\div 4 n_0$, where $n_0$ is the nuclear saturation density. In the framework of the relativistic mean field model the new phase arises in a second-order phase transition. The appearance of a $\rho^-$ condensate significantly alters the proton fraction in a neutron star but changes moderately the equation of state. The neutrino emissivity of the processes involving a $\rho^-$ meson condensate is estimated.

We show that within a recently developed nonlocal, chiral quark model the critical densities for a phase transition to color superconducting quark matter under neutron star conditions can be low enough that these phases occur in compact star configurations with masses below 1.4 M_sun. We study the cooling of these objects in isolation for different values of the gravitational mass and thus different composition and structure of the interior. Our equation of state allows for a 2SC phase with a large quark gap Delta ~ 100 MeV for u and d quarks of two colors, a normal quark matter phase and their coexistence in a mixed phase within the hybrid star interior. We argue that, if the phases with unpaired quarks were allowed, the corresponding hybrid stars would cool too fast to describe the neutron star cooling data existing by today. We incorporate other attractive channels permitting a weak pairing of the residual quarks which remained unpaired in the 2SC phase and demonstrate that the model does not contradict the cooling data if the weak pairing gaps are of the order of 0.1 MeV.

We study the cooling of isolated neutron stars. The main cooling regulators are introduced: EoS, thermal transport, heat capacity, neutrino and photon emissivity, superfluid nucleon gaps. Neutrino emissivity includes main processes. A strong impact of medium effects on the cooling rates is demonstrated. With taking into account of medium effects in reaction rates and in nucleon superfluid gaps modern experimental data can be well explained.

Thermal fluctuations of the color superconducting order parameter in dense quark matter are investigated in terms of the phenomenological Ginzburg - Landau approach. Our estimates show that fluctuations of the di-quark gap may strongly affect some of thermodynamic quantities even far below and above the critical temperature. If the critical temperature $T_c$ of the di-quark phase transition were rather high ($\gsim (50\div 70)$ MeV) one could expect a manifestation of fluctuations of the di-quark gap in the course of heavy ion collisions (above $T_c$). For $T_c \sim 50 $MeV color superconducting fluctuations may also affect an initial stage of the hybrid star evolution.

Fluctuations of the color superconducting order parameter in dense quark matter at finite temperatures are investigated in terms of the phenomenological Ginzburg - Landau approach. Our estimates show that fluctuations of the di-quark gap may strongly affect some of thermodynamic quantities even far below and above the critical temperature. If the critical temperature of the di-quark phase transition were rather high one could expect a manifestation of fluctuations of the di-quark gap in the course of heavy ion collisions.

Nuclear medium effects in the neutrino cooling of neutron stars through the exotic reaction channel \gamma \gamma --> \pi^0--> \nu_R \bar{\nu_L} (\nu_L \bar{\nu_R}) are incorporated. Throughout the paper we discuss different possibilities of right-handed neutrinos, massive left-handed neutrinos and standard massless left-handed neutrinos (reaction is then allowed only with medium modified vertices). It is demonstrated that multi-particle effects suppress the rate of this reaction channel by 6-7 orders of magnitude that does not allow to decrease existing experimental upper limit on the corresponding \pi^0\nu\bar{\nu} coupling. Other possibilities of the manifestation of the given reaction channel in differente physical situations, e.g. in the quark color superconducting cores of some neutron stars, are also discussed. We demonstrate that in the color-flavor-locked superconducting phase for temperatures T < (0.1-10) MeV (depending on the effective pion mass and the decay width) the process is feasibly the most efficient neutrino cooling process, although the absolute value of the reaction is rather small.

Possibility of structured mixed phases at first order phase transitions in neutron stars is re-examined by taking into account charge screening and surface effects. Transition from the hadron $npe$ phase to the quark phase is studied. Two possibilities, the mixed phase and two separate phases given by the double-tangent (Maxwell) construction are considered. Inhomogeneous profiles of the electric potential and their contribution to the energy are analytically calculated. The electric field configurations determine the droplet size and the geometry of structures.

Kaon polarization operator in dense baryonic matter of arbitrary isotopic composition is calculated including s- and p-wave kaon-baryon interactions. The regular part of the polarization operator is extracted from the realistic kaon-nucleon interaction based on the chiral and 1/N_c expansion. Contributions of the Lambda(1116), Sigma(1195), Sigma*(1385) resonances are taken explicitly into account in the pole and regular terms with inclusion of mean-field potentials. The baryon-baryon correlations are incorporated and fluctuation contributions are estimated. Results are applied for K- in neutron star matter. Within our model a second-order phase transition to the s-wave K- condensate state occurs at $rho_c \gsim 4 \rho_0$ once the baryon-baryon correlations are included. We show that the second-order phase transition to the p-wave $K^-$ condensate state may occur at densities $\rho_c \sim 3\div 5 \rho_0$ in dependence on the parameter choice. We demonstrate that a first-order phase transition to a proton-enriched (approximately isospin-symmetric) nucleon matter with a p-wave K- condensate can occur at smaller densities, $\rho\lsim 2 \rho_0$. The transition is accompanied by the suppression of hyperon concentrations.

Possibility of structured mixed phases at first order phase transitions is examined with taking into account of charge screening and surface effects. Hadron-quark phase transition in dense neutron star interior is considered, as concrete example.

This review demonstrates that neutrino emission from dense hadronic component in neutron stars is subject of strong modifications due to collective effects in the nuclear matter. With the most important in-medium processes incorporated in the cooling code an overall agreement with available soft X ray data can be easily achieved. With these findings so called "standard" and "non-standard" cooling scenarios are replaced by one general "nuclear medium cooling scenario" which relates slow and rapid neutron star coolings to the star masses (interior densities). In-medium effects take important part also at early hot stage of neutron star evolution decreasing the neutrino opacity for less massive and increasing for more massive neutron stars. A formalism for calculation of neutrino radiation from nuclear matter is presented that treats on equal footing one-nucleon and multiple-nucleon processes as well as reactions with resonance bosons and condensates. Cooling history of neutron stars with quark cores is also discussed.

The behaviour of the magnetic field inside the superconducting quark matter core of a neutron star is investigated in the framework of the Ginzburg-Landau theory. We take into account the simultaneous coupling of the diquark condensate field to the usual magnetic and to the gluomagnetic gauge fields. We solve the problem for three different physical situations: a semi-infinite region with a planar boundary, a spherical region, and a cylindrical region. We show that Meissner currents near the quark core boundary effectively screen the external static magnetic field.

We study the consequences of superconducting quark cores (with color-flavor-locked phase as representative example) for evolution of temperature profiles and the cooling curves in quark-hadron hybrid stars and in hypothetical self-bounded objects having no a hadron shell (quark core neutron stars). The quark gaps are varied from 0 to $\Delta_q =50$ MeV. For hybrid stars we find time scales of $1\div5$, $5\div10$ and $50\div100$ years for the formation of a quasistationary temperature distribution in the cases $\Delta_q =0$, 0.1 MeV and $\gsim$ 1 MeV, respectively. These time scales are governed by the heat transport within quark cores for large diquark gaps ($\Delta \gsim$ 1 MeV) and within the hadron shell for small diquark gaps ($\Delta \lsim 0.1$ MeV). For quark core neutron stars we find a time scale $\simeq 300$ years for the formation of a quasistationary temperature distribution in the case $\Delta \gsim$ 10 MeV and a very short one for $\Delta \lsim$ 1 MeV. If hot young compact objects will be observed they can be interpreted as manifestation of large gap color superconductivity. Depending on the size of the pairing gaps, the compact star takes different paths in the ${lg}T_s $ vs. ${lg} t$ diagram where $T_s$ is the surface temperature. Compared to the corresponding hadronic model which well fits existing data the model for the hybrid neutron star (with a large diquark gap) shows too fast cooling. The same conclusion can be drawn for the corresponding self-bound objects.

This review demonstrates that neutrino emission from dense hadronic component in neutron stars is subject of strong modifications due to collective effects in the nuclear matter. With the most important in-medium processes incorporated in the cooling code an overall agreement with available soft $X$ ray data can be easily achieved. With these findings so called "standard" and "non-standard" cooling scenarios are replaced by one general "nuclear medium cooling scenario" which relates slow and rapid neutron star coolings to the star masses (interior densities). In-medium effects take important part also at early hot stage of neutron star evolution decreasing the neutrino opacity for less massive and increasing for more massive neutron stars. A formalism for calculation of neutrino radiation from nuclear matter is presented that treats on equal footing one-nucleon and multiple-nucleon processes as well as reactions with resonance bosons and condensates.

Effects of the propagation of particles, which have a finite life-time and an according width in their mass spectrum, are discussed in the context of transport description. First, the importance of coherence effects (Landau-Pomeranchuk-Migdal effect) on production and absorption of field quanta in non-equilibrium dense matter is considered. It is shown that classical diffusion and Langevin results correspond to re-summation of certain field-theory diagrams formulated in terms of full non-equilibrium Green's functions. Then the general properties of broad resonances in dense and hot systems are discussed in the framework of a self-consistent and conserving Phi-derivable method of Baym at the examples of the rho-meson in hadronic matter and the pion in dilute nuclear matter. Further we address the problem of a transport description that properly accounts for the damping width of the particles. The Phi-derivable method generalized to the real-time contour provides a self-consistent and conserving kinetic scheme. We derive a generalized expression for the non-equilibrium kinetic entropy flow, which includes corrections from fluctuations and mass-width effects. In special cases an H-theorem is proved. Memory effects in collision terms give contributions to the kinetic entropy flow that in the Fermi-liquid case recover the famous bosonic type T^3 ln T correction to the specific heat of liquid Helium-3. At the example of the pion-condensate phase transition in dense nuclear matter we demonstrate important part played by the width effects within the quantum transport.

It is proposed that a rearrangement of single-particle degrees of freedom may occur in a portion of the quantum fluid interior of a neutron star. Such a rearrangement is associated with the pronounced softening of the spin-isospin collective mode which, under increasing density, leads to pion condensation. Arguments and estimates based on fundamental relations of many-body theory show that one realization of this phenomenon could produce very rapid cooling of the star via a direct nucelon Urca process displaying a $T^5$ dependence on temperature.

The effects of the propagation of particles which have a finite life-time and an according width in their mass spectrum are discussed in the context of transport descriptions. In the first part the coupling of soft photon modes to a source of charged particles is studied in a classical model which can be solved completely in analytical terms. The solution corresponds to a re-summation of certain field theory diagrams. The general properties of broad resonances in dense finite temperature systems are discussed at the example of the $\rho$-meson in hadronic matter. The second part addresses the problem of transport descriptions which also account for the damping width of the particles. The Kadanoff--Baym equation after gradient approximation together with the $\Phi$-derivable method of Baym provides a self-consistent and conserving scheme. Memory effects appearing in collision term diagrams of higher order are discussed. We derive a generalized expression for the nonequilibrium kinetic entropy flow, which includes corrections from fluctuations and mass-width effects. In special cases an $H$-theorem is proved. Memory effects in collision terms provide contributions to the kinetic entropy flow that in the Fermi-liquid case recover the famous bosonic type $T^3 \ln T$ correction to the specific heat of liquid Helium-3.

The effect of color superconductivity on the cooling of quark stars and neutron stars with large quark cores is investigated. Various known and new quark-neutrino processes are studied. As a result, stars being in the color flavor locked (CFL) color superconducting phase cool down extremely fast. Quark stars with no crust cool down too rapidly in disagreement with X-ray data. The cooling of stars being in the N_f =2 color superconducting (2SC) phase with a crust is compatible with existing X-ray data. Also the cooling history of stars with hypothetic pion condensate nuclei and a crust does not contradict the data.

Within the real-time formulation of nonequilibrium field theory, generalized transport equations are derived avoiding the standard quasiparticle approximation. They permit to include unstable particles into the transport scheme. In order to achieve a self-consistent, conserving and thermodynamically consistent description, we generalize the Baym's $\Phi$-functional method to genuine nonequilibrium processes. The developed transport description naturally includes all those quantum features already inherent in the corresponding equilibrium limit. Memory effects appearing in collision term diagrams of higher order are discussed. The variational properties of $\Phi$-functional permit to derive a generalized expression for the non-equilibrium kinetic entropy flow, which includes fluctuations and mass width effects. In special cases an $H$-theorem is demonstrated implying that the entropy can only increase with time. Memory effects in the kinetic terms provide corrections to the kinetic entropy flow that in equilibrium limit recover the famous bosonic type $T^3 \ln T$ correction to the specific heat of Fermi liquids like Helium-3.

Production of negative strangeness in reactions of inelastic anti-neutrino scattering on a nucleus provides information on the modification of strange degrees of freedom in nuclear matter. We calculate cross-sections of the reaction channels $\bar{\nu}_{e(\mu)}\to e^+(\mu^+) + K^{-}$ and $\bar{\nu}_{e(\mu)} + p \to \Lb + e^+(\mu^+)$ and investigate their sensitivity to the medium effects. In particular, we consider effects induced by the presence of a low-energy excitation mode in the $K^-$ spectrum, associated with the correlated $\Lambda$-particle and proton-hole states, and the renormalization of the weak interaction in medium. In order to avoid double counting, various contributions to anti-neutrino scattering are classified with the help of the optical theorem, formulated within the non-equilibrium Green's function technique.

We consider massive photon decay reactions via intermediate states of electron-electron-holes and proton-proton-holes into neutrino-antineutrino pairs in the course of neutron star cooling. These reactions may become operative in hot neutron stars in the region of proton pairing where the photon due to the Higgs-Meissner effect acquires an effective mass $m_{\gamma}$ that is small compared to the corresponding plasma frequency. The contribution of these reactions to neutrino emissivity is calculated; it varies with the temperature and the photon mass as $T^{3/2}m_{\gamma}^{7/2} e^{-m_{\gamma}/T}$ for $T < m_{\gamma}$. Estimates show that these processes appear as extra efficient cooling channels of neutron stars at temperatures $T \simeq (10^9-10^{10})$ K.

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.