To constrain the allowed range for the axion decay constant $f_{a}$ or, equivalently, for the axion mass $m_{a}$, we consider the cooling of a neutron star with strong proton superfluidity and normal (non-superfluid) neutrons inside its core and without strong magnetic field, by analogy with the observed supernova remnant in HESS J1731-347. For this specific case, we demonstrate that after the thermal relaxation is over, the hydrostatic structure of the neutron star can be well described with the aid of solution of Einstein field equations, applied to a sphere of fluid in hydrostatic equilibrium, derived by Tolman. The internal temperature of the neutron star is calculated assuming that the cooling occurs dominantly due to production of neutrino pairs and axions in the nn-bremsstrahlung. To impose a constraint to the axion decay constant the fact is used that the currently observed neutron star surface temperature does not deviate from the neutrino cooling scenario. For the KSVZ-axion model we find that $f_{a}>1.9\times 10^{8}$ GeV, while for the DFSZ-axion model we obtain $f_{a}>4.7\times 10^{9}$ GeV.

A brief review is given of the current state of the problem of neutrino pair emission through neutral weak currents caused by the Cooper pairs breaking and formation (PBF) in superfluid baryon matter at thermal equilibrium. The cases of singlet-state pairing with isotropic superfluid gap and spin-triplet pairing with anisotropic gap are analyzed with allowance for the anomalous weak interactions caused by superfluidity. It is shown that taking into account the anomalous weak interactions in both the vector and axial channels is very important for a correct description of neutrino energy losses through the PBF processes. The anomalous contributions lead to an almost complete suppression of the PBF neutrino emission in spin-singlet superfluids and strong reduction of the PBF neutrino losses in the spin-triplet superfluid neutron matter, which considerably slows down the cooling rate of neutron stars with superfluid cores.

Axial anomalous contributions into neutrino PBF losses due to triplet pairing of neutrons are still ignored in modeling the evolution of neutron stars. In this paper, the influence of the anomalous axial contributions onto the rate of neutron stars cooling is estimated.

The superfluid drag effect, in hydrodynamics of pulsating neutron stars, is conventionally described with the aid of the entrainment matrix relating the mass currents with the velocities of superfluid flows in the system. Equations for the entrainment matrix of a superfluid mixture of neutrons and protons are derived with allowance for the strong dependence of the energy gaps on the velocities of superfluid flows. The calculations are carried out in the frame of the Fermi-liquid theory. The equations obtained are highly nonlinear. Numerical solutions to the equations for some typical cases demonstrate that the components of the entrainment matrix possess a highly nonlinear dependence on the velocities of the two superflows simultaneously. This effect, previously ignored, can greatly influence the dynamics of neutron stars.

The observed rapid cooling of the Cassiopeia A neutron star (Cas A NS) can be interpreted as being triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). This provides a unique possibility for probing the neutron condensate in the core. Using consistent neutron star core and crust equation of state and composition, I explore the sensitivity of this interpretation to the phase state of the triplet superfluid condensate. Modeling cooling within an expected range of neutron star masses and envelope compositions, I found that the fast cooling of the Cas A NS can not be explained by the PBF processes in the conventional $^3$P$_2$ condensate with $m_{j}=0$.

In a simple model it is demonstrated that the neutron star surface temperature evolution is sensitive to the phase state of the triplet superfluid condensate. A multicomponent triplet pairing of superfluid neutrons in the core of a neutron star with participation of several magnetic quantum numbers leads to neutrino energy losses exceeding the losses from the unicomponent pairing. A phase transition of the neutron condensate into the multicomponent state triggers more rapid cooling of superfluid core in neutron stars. This makes it possible to simulate an anomalously rapid cooling of neutron stars within the minimal cooling paradigm without employing any exotic scenarios suggested earlier for rapid cooling of isolated neutron star in Cassiopeia A.

We examine the influence of the anisotropy of the superfluid energy gap and residual Fermi-liquid interactions in the triplet-correlated neutron liquid onto neutrino energy losses through neutral weak currents. The neutrino-pair emission caused by the pair breaking and formation processes and by the spin-wave decays is considered for the case of the $^{3}P_{2}$ pairing in the state with $m_{j}=0$. The simple analytical formulae are obtained. A comparison with the previous results of the average-angle approach shows that the gap anisotropy leads to quenching of the neutrino emissivity caused by the pair recombination processes on about 15% and to substantial suppression of the spin-wave decays. Residual particle-hole interactions increase the energy losses in both the channels on about 5%.

The influence of the admixture of the $^{3}F_{2}$ state onto collective spin oscillations and neutrino emission processes in the triplet superfluid neutron liquid is studied in the BCS approximation. The eigen mode of spin oscillations with $\omega\simeq\sqrt{58/35}\Delta$ is predicted to exist in the triplet superfluid neutron condensate besides the already known mode $\omega\simeq\Delta/\sqrt{5}$. Excitation of the high-frequency spin oscillations in the condensate occurs through the tensor interactions between quasiparticles. Neutrino energy losses through neutral weak currents are found to consist of three separate contributions caused by a recombination of broken Cooper pairs and by weak decays of the collective modes of spin oscillations. Neutrino decays of the low-frequency spin waves can play an important role in the cooling scenario of neutron stars. Weak decays of the high-frequency oscillations that occur only if the tensor forces are taken into account in the pairing interactions does not modify substantially the total energy losses. Simple expressions are suggested for the total neutrino emissivity.

The eigen mode of spin oscillations with $\omega\simeq \sqrt{58/35}\Delta$ is predicted to exist besides already known spin waves with $\omega \simeq\Delta /\sqrt{5}$ in the triplet superfluid neutron condensate in the inner core of neutron stars. The new mode is kinematically able to decay into neutrino pairs through neutral weak currents. The problem is considered in BCS approximation for the case of $^{3}P_{2}-^{3}F_{2}$ pairing with a projection of the total angular momentum $m_{j}=0$ which is conventionally considered as preferable one at supernuclear densities.

The linear response of a superfluid neutron liquid onto external vector field is studied for the case of $^{3}P_{2}-\,^{3}F_{2}$ pairing. The consideration is limited to the case when the wave-length of the perturbation is large as compared to the coherence length in the superfluid matter and the transferred energy is small in comparison with the gap amplitude. The obtained results are used to analyse the collisionless phonon-like excitations of the condensate of superfluid neutrons. In particular, we analyze the case of neutron condensation into the state with $m_{j}=0$ which is conventionally considered as the preferable one in the bulk matter of neutron stars. Zero sound (if it exists) is found to be anisotropic and undergoes strong decrement below some temperature threshold depending substantially on the intensity of Fermi-liquid interactions.

Neutrino energy losses through neutral weak currents in the triplet-spin superfluid neutron liquid are studied for the case of condensate involving several magnetic quantum numbers. Low-energy excitations of the multicomponent condensate in the timelike domain of the energy and momentum are analyzed. Along with the well-known excitations in the form of broken Cooper pairs, the theoretical analysis predicts the existence of collective waves of spin density at very low energy. Because of a rather small excitation energy of spin waves, their decay leads to a substantial neutrino emission at the lowest temperatures, when all other mechanisms of neutrino energy loss are killed by a superfluidity. Neutrino energy losses caused by the pair recombination and spin-wave decays are examined in all of the multicomponent phases that might represent the ground state of the condensate, according to modern theories, and for the case when a phase transition occurs in the condensate at some temperature. Our estimate predicts a sharp increase in the neutrino energy losses followed by a decrease, along with a decrease in the temperature, that takes place more rapidly than it would without the phase transition. We demonstrate the important role of the neutrino radiation caused by the decay of spin waves in the cooling of neutron stars.

The linear response of a neutron spin-triplet superfluid onto external weak axial-vector field is studied for the case of $^{3}P_{2}$ pairing with a projection of the total angular momentum $m_{j}=0$. The problem is considered in the BCS approximation discarding Fermi-liquid effects. The anomalous axial-vector vertices of neutron quasiparticles possess singularities at some frequencies which specify existence of undamped spin-density waves in the Cooper condensate. The spin waves are of a low excitation energy and are kinematically able to decay into neutrino pairs through neutral weak currents. We evaluate the neutrino emissivity from the spin wave decays in the bulk neutron superfluid in old neutron stars. This calculation predicts significant energy losses from within a neutron star at lowest temperatures when all other mechanisms of neutrino emission are killed by the neutron and proton superfluidity.

Neutrino emission due to the pair breaking and formation processes in the bulk triplet superfluid in neutron stars is investigated with taking into account of anomalous weak interactions. We consider the problem in the BCS approximation discarding Fermi-liquid effects. In this approach we derive self-consistent equations for anomalous vector and axial-vector vertices of weak interactions taking into account the $^{3}P_{2}- ^{3}F_{2}$ mixing. Further we simplify the problem and consider the pure $^{3}P_{2}$ pairing with $m_{j}=0$, as is adopted in the minimal cooling paradigm. As was expected because of current conservation we have obtained a large suppression of the neutrino emissivity in the vector channel. More exactly, the neutrino emission through the vector channel vanishes in the nonrelativistic limit $V_F=0$. The axial channel is also found to be moderately suppressed. The total neutrino emissivity is suppressed by a factor of $1.9\times10^{-1}$ relative to original estimates using bare weak vertices.

The linear response of a nonrelativistic superfluid baryon system on an external weak field is investigated while taking into account of the Fermi-liquid interactions. We generalize the theory developed by Leggett for a superfluid Fermi-liquid at finite temperature to the case of timelike momentum transfer typical of the problem of neutrino emission from neutron stars. A space-like kinematics is also analysed for completeness and compared with known results. We use the obtained response functions to derive the neutrino energy losses caused by recombination of broken pairs in the electrically neutral superfluid baryon matter. We find that the dominant neutrino radiation occurs through the axial-vector neutral currents. The emissivity is found to be of the same order as in the BCS approximation, but the details of its temperature dependence are modified by the Fermi-liquid interactions. The role of electromagnetic correlations in the pairing case of protons interacting with the electron background is discussed in the conclusion.

We examine the effective interaction of nonrelativistic fermions with an external vector field in superfluid systems. In contrast to the complicated vertex equation, usually used in this case, we apply the approach which does not employ an explicit form of the pairing interaction. This allows to obtain a simple analytic expression for the vertex function only in terms of the order parameter and other macroscopic parameters of the system. We use this effective vertex to analyze the linear response function of the superfluid medium at finite temperatures. At the time-like momentum transfer, the imaginary part of the response function is found to be proportional to the fourth power of small Fermi velocity, i.e. the energy losses through vector currents are strongly suppressed. As an application, we calculate the neutrino energy losses through neutral weak currents caused by the pair recombination in the superfluid neutron matter at temperatures lower than the critical one for S-wave pairing. This approach confirms a strong suppression of the neutrino energy losses as predicted in Ref.[4].

Neutrino emission caused by Cooper pairing of baryons in neutron stars is recalculated by accurately taking into account for conservation of the vector weak current. The vector current contribution to the neutrino emissivity is found to be several orders of magnitude smaller than that obtained before by different authors. Therefore, the neutrino energy losses due to singlet-state pairing of baryons can in practice be neglected in simulations of neutron star cooling. This makes negligible the neutrino radiation from pairing of protons or hyperons. The neutrino radiation from triplet pairing takes place through axial weak currents. For these states, when the total momentum projection is $m_{j}=0$, the vanishing of the vector weak current contribution results in the suppression of the neutrino energy losses by about 25%. The neutrino emissivity due to triplet pairing with $| m_{j}| =2$ is suppressed by about a factor of 3, caused by the collective contribution of spin-density fluctuations in the condensate.

Neutrino emission caused by singlet Cooper pairing of baryons in neutron stars is recalculated by accurately taking into account for conservation of the vector weak currents. The neutrino emissivity via the vector weak currents is found to be several orders of magnitude smaller than that obtained before by different authors. This makes unimportant the neutrino radiation from singlet pairing of protons or hyperons.

The process of $\nu \bar{\nu}$ radiation from a relativistic plasma of electrons and positrons is studied within the Random Phase Approximation.

At temperatures less than a few MeV, the efficiency of the neutrino-pair bremsstrahlung from nucleons interracting with the condensed pion field is much less than previously estimated by other authors. The physical reason for this is a periodic structure of the developed pion condensate. The repeated interactions of a nucleon with the periodic pion field modify the nucleon spectrum, which is split into two bands. The energy gap between the bands is about a few MeV or larger, even if the amplitude of the pion field is small compared with the pion mass.

In the mean field approximation, we study the effects of weak magnetism and pseudoscalar interaction in the neutrino energy losses caused by the direct Urca processes on relativistic nucleons in the degenerate baryon matter. Our formula for the neutrino energy losses incorporates the effects of nucleon recoil, parity violation, weak magnetism, and pseudoscalar interaction. For numerical testing of our formula, we use a self-consistent relativistic model of the multicomponent baryon matter. We found that, due to weak magnetism effects, relativistic emissivities increase by approximately 40-50%, while the pseudoscalar interaction only slightly suppresses the energy losses, approximately by 5%.

We use the field theoretical model to perform relativistic calculations of neutrino energy losses caused by the direct Urca processes on nucleons in the degenerate baryon matter. By our analysis, the direct neutron decay in the superdense nuclear matter under beta equilibrium is open only due to the isovector meson fields, which create a large energy gap between protons and neutrons in the medium. Our expression for the neutrino energy losses, obtained in the mean field approximation, incorporates the effects of nucleon recoil, parity violation, weak magnetism, and pseudoscalar interaction. For numerical testing of our formula, we use a self-consistent relativistic model of the multicomponent baryon matter. The relativistic emissivity of the direct Urca reactions is found substantially larger than predicted in the non-relativistic approach. We found that, due to weak magnetism effects, relativistic emissivities increase by approximately 40-50%, while the pseudoscalar interaction only slightly suppresses the energy losses, approximately by 5%.

We derive a relativistic expression for neutrino energy losses caused by the direct Urca processes in degenerate baryon matter of neutron stars. We use two different ways to calculate the emissivity caused by the reactions to our interest. First we perform a standard calculation by Fermi's ''golden'' rule. The second calculation, resulting in the same expression, is performed with the aid of polarization functions of the medium. Our result for neutrino energy losses strongly differs from previous non-relativistic results. We also discuss nonconservation of the baryon vector current in reactions through weak charged currents in the medium, when the asymmetry between protons and neutrons is considered. The above effects, not discussed in the literature before, substantially modify the polarization functions responsible for the induced weak charged currents in baryon matter.

We have investigated the bound states of an electron and positron in superstrong magnetic fields typical for neutron stars. The complete relativistic problem of positronium in a strong magnetic field has not been succesfully solved up to now. In particular, we have studied the positronium when it moves relativistically across the magnetic field. A number of problems which deal with the pulsar magnetosphere, as well as the evolution of protoneutron stars, could be considered as a field for application.

We present a method to calculate $\nu \bar{\nu}$ energy losses from neutron star crusts, which automatically takes into account for collective effects, and allows to calculate the total emissivity without a separate consideration of particular processes. We show that the formula we obtain describe the known results for the emissivity due to plasmon decay and Bremsstrahlung from degenerate electrons, when one of this processes dominates. In the case of low temperatures, our formula gives a suppression of the electron vector weak-current contribution to $\nu \bar{\nu}$ Bremsstrahlung, due to the collective effects discussed in this paper.

The neutrino emission due to formation and breaking of Cooper pairs of protons in superconducting cores of neutron stars is considered with taking into account the electromagnetic coupling of protons to ambient electrons. It is shown that collective response of electrons to the proton quantum transition contributes coherently to the complete interaction with a neutrino field and enhances the neutrino-pair production. Our calculation shows that the contribution of the vector weak current to the $\nu \bar{\nu}$ emissivity of protons is much larger than that calculated by different authors without taking into account the plasma effects. Partial contribution of the pairing protons to the total neutrino radiation from the neutron star core is very sensitive to the critical temperatures for the proton and neutron pairing. We show domains of these parameters where the neutrino radiation, caused by a singlet-state pairing of protons is dominating.

We have considered collective effects in $\nu \bar{\nu}$ synchrotron radiation from an ultrarelativistic degenerate electron gas in neutron stars with strong magnetic fields. For this problem we apply a calculation method which explicitly makes use of the fact that the radiating electron moves semi-classically, but takes into account the interaction among particles in a quantum way. First we apply this method to calculate $\nu \bar{\nu}$ synchrotron radiation by an ultrarelativistic electron in vacuum and we compare this result with that obtained previously by other techniques. When a degenerate plasma is considered, we show that collective effects lead to an essential enhancement (about three times) of the vector weak-current contribution to neutrino pair emissivity.

We calculate the emissivity for the direct URCA process in strongly magnetized, degenerate matter in neutron stars, under $\beta $-equilibrium. We show that, if the magnetic field is large enough for protons and electrons to be confined to the ground Landau levels, the field-free threshold condition on proton concentration no longer holds, and direct URCA reactions are open for an arbitrary proton concentration. Direct URCA process leads to an early phase of fast neutron star cooling. This circumstance allows us to constrain the initial magnetic field inside observed pulsars.