Forces with a large radius of interaction can have a significant impact on the equation of state of matter. Low-mass neutrinos generate a long-range potential due to the exchange of neutrino pairs. We discuss a possible relationship between the neutrino masses, which determine the interaction radius of the neutrino-pair exchange potential, and the equation of state of neutron matter. Contrary to previous statements, the thermodynamic potential, when decomposed into the number of neutrino interactions, vanishes in any decomposition order, except for the interaction of two neutrons. In the one-loop approximation, long-range multiparticle neutrino interactions are stable in the infrared region for all neutrino masses and do not affect the equation of state of neutron matter or the stability of neutron stars.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Among the fermions of the Standard Model, neutrinos are the lightest particles. Their masses are at least six orders of magnitude smaller than the mass of any other charged fermion. The exchange of low-mass particles creates a long-range potential. The exchange of massless photons, e.g., leads to the Coulomb potential. Since neutrinos are fermions, long-range two-body forces can involve them through the formation of neutrino pairs [1–5]. The neutrino-pair exchange potential is similar to the van der Waals potential arising from the two-photon exchange (see, e.g., [6]). Weakly interacting light bosons provide an excellent illustration of the significant influence of weak forces with large interaction radii on the equation of state (EoS) of neutron matter and the structure of neutron stars [7, 8].
Fischbach [4] considered the effect of long-range multiparticle interactions of neutrinos on the EoS of neutron matter and concluded that the contribution of neutrino interactions to EoS diverges in the infrared region. To guarantee the finiteness of EoS and, ultimately, the stability of neutron stars, he postulated a lower limit for the neutrino masses of \(m \gtrsim 0.4\) eV. Cosmological models place an upper limit on the sum of neutrino masses of 0.13 eV [9]. According to KATRIN experiment on tritium \(\beta \) decay, the upper limit on the effective electron neutrino mass is 0.8 eV [10]. Fishbach’s estimate is close to these limits and partly intersects with them, which requires a thorough analysis of multiparticle neutrino interactions in nuclear matter. In relation to the KATRIN experiment, the mass constraint [4] is discussed in a recent paper [11].
Abada, Gavela, and Pinea [12] address the effect of multiparticle neutrino interactions in neutron matter using the standard techniques of quantum statistics (see, e.g., [13]). The authors confirm an infrared instability of EoS in each order of decomposition by the number of interactions, but conclude, nevertheless, that the total contribution of multiparticle neutrino interactions to the EOS of neutron matter is zero.
In this paper, we show that, contrary to the previous statements [4, 12], the multiparticle interactions of neutrinos are stable in the infrared region and, moreover, their contributions vanish at each term of the EoS expansion into a power series with respect to the number density. The structure of neutron stars is thereby not sensitive to the masses of neutrinos.
The effective Hamiltonian for low-energy interaction of neutrinos and neutrons is generated by the exchange of the Z boson. We consider the case of Dirac neutrinos. After averaging the neutral weak current of quarks over the neutron wave function, the effective Hamiltonian takes the form
where \({{g}_{A}}\) is the axial coupling constant of nucleons. The axial component of the weak current, which indicates the direction of the average spin of neutrons, vanishes in unpolarized matter, so the Z-boson mean field, U, is a pure vector. The massive neutrino interacts with the potential U by its left component only. In mean-field approximation, \(\langle \bar {n}(x){{\gamma }^{\mu }}(1 - {{g}_{A}}{{\gamma }_{5}})n(x)\rangle = {{g}^{{\mu 0}}}\rho \). The typical number density of neutrons is \(\rho \sim 0.4\) fm–3. The typical Fermi momentum of neutrons is a few hundred MeV, whereas the Z-boson mean-field potential equals \(U = - {{G}_{F}}\rho {\text{/}}\sqrt 2 \sim - 20\) eV. We work in the approximations of homogeneous neutron matter and flat Minkowski space. These approximations are well-founded for neutrinos with masses greater than the inverse radius of neutron stars; i.e., \(m \gtrsim 1{\text{/}}{{R}_{{\text{s}}}} \sim \) \(2 \times {{10}^{{ - 11}}}\) eV, where \({{R}_{{\text{s}}}} \sim 10\) km.
It is useful to define projection operators \(L = (1 - {{\gamma }_{5}}){\text{/}}2\), \(R = (1 + {{\gamma }_{5}}){\text{/}}2\), and \({{P}_{ \pm }} = (1 \pm \boldsymbol{\alpha} {\kern 1pt} {\mathbf{n)}}{\text{/}}2\), where \({\kern 1pt} \boldsymbol{\alpha} = {{\gamma }_{0}}\gamma \) and \({\mathbf{n}} = {\mathbf{q}}/{\text{|}}{\mathbf{q}}{\text{|}}\) is the unit vector oriented in the direction of neutrino momentum.
The effective Lagrangian of a neutrino with mass m has the form:
The Green function is defined by the quadratic form of the effective Lagrangian. In the momentum representation,
The change in the thermodynamic potential, Ω, due to an interaction is given by the well-known expression (see, e.g., [13])
In the case under consideration, \(H_{{\text{I}}}^{\lambda } = \lambda H_{{\text{I}}}^{\;}\) is the effective Hamiltonian (1) with the scaled coupling constant. In terms of the Green function,
where V is the normalization volume. This expression implies a smooth thermodynamic limit \(V \to \infty \), \(\rho = {\text{const}}\). In a realistic approach, the number of particles in a star is finite, although large. The neutrino propagator should be expanded into a power series by the neutron number density, i.e., by the parameter U, and the series should be truncated at s ~ N ≡ \({{M}_{ \odot }}{\text{/}}{{m}_{n}} = \) 1.2 × 1057, where \({{M}_{ \odot }}\) is the mass of the Sun, and \({{m}_{n}}\) is the mass of the neutron. Each term, proportional to \({{U}^{s}}\), describes the scattering of neutrino by s neutrons. If the series converges, the limit \(N \to \infty \) is well defined and the decomposition is not required.
The papers [4, 12] declare that for massless neutrinos, the individual terms of the series are proportional to \({{(U{{R}_{{\text{s}}}})}^{s}}\). If this were true, the infinite series would diverge because \(U{{R}_{{\text{s}}}} \sim {{10}^{{12}}} \gg 1\). Abada et al. further argue the transition to the limit of \(N \to \infty \) by the possibility for neutrinos to scatter several times on the same neutron. Multiple scattering involving the same neutron is possible only in higher orders of the loop expansion, whereas Abada et al. work in the one-loop approximation. In this approximation, the transition to the limit of \(N \to \infty \) is not allowed if the series does not converge.
We expand the expression (5) in a series by U. The integration by λ gives
The relations \(LR = 0\), \(L{{\gamma }_{\mu }} = {{\gamma }_{\mu }}R\) lead to the identity
In terms of the projection operators defined above, \({{\gamma }_{0}}\hat {q} = ({{q}_{0}} - \;{\text{|}}{\mathbf{q}}{\text{|}}){{P}_{ + }} + ({{q}_{0}} + \;{\text{|}}{\mathbf{q}}{\text{|}}){{P}_{ - }}\). Using the binomial formula for \({{({{\gamma }_{0}}\hat {q})}^{s}}\) and the relations \({{P}_{ + }}{{P}_{ - }} = 0\), \(P_{ \pm }^{2} = {{P}_{ \pm }}\), the right side of Eq. (7) can be simplified to give
Closing the contour of integration by \({{q}_{0}}\) in the upper half of the complex plane, we find that the integral is determined by the residues at \({{q}_{0}} = - \sqrt {{{{\mathbf{q}}}^{2}} + {{m}^{2}}} + i0\). Using the symbolic computing software package Maple,Footnote 1 it is possible to find the residues and the corresponding integrals over the momentum space for quite large s. It turns out that all the terms \(5 \leqslant s \leqslant 100\) of the series vanish identically. After regularization of the neutrino loop, the term \(s = 2\) becomes finite, whereas the terms \(s = 3,4\) vanish. For \(s = 1\), the integral over an infinitely distant contour in the upper half of the \({{q}_{0}}\)-plane cancels the contribution of the residue at \({{q}_{0}} = - \sqrt {{{{\mathbf{q}}}^{2}} + {{m}^{2}}} + i0\).
After performing the Wick rotation \({{q}_{0}} \to i\omega \) and assuming that the limit \(\tau \to - 0\) can be interchanged with the momentum integral, a more general proof can be offered. The integration by ω goes within the limits \( - \infty < \omega < \infty \). The spherical coordinate system in the Euclidean space \((\omega ,{\mathbf{q}})\) is defined by \(\omega = \eta \cos \alpha \), \({{q}_{x}} = \eta \sin \alpha \cos \beta \), \({{q}_{y}} = \eta \sin \alpha \sin \beta \cos \gamma \), and \({{q}_{z}} = \eta \sin \alpha \sin \beta \sin \gamma \). The angles are restricted by \(0 \leqslant \alpha \), \(\beta \leqslant \pi \), \(0 \leqslant \gamma \leqslant 2\pi \). The absolute value of momentum is \({\text{|}}{\mathbf{q}}{\text{|}} = \eta \sin \alpha \). The volume element is \({{d}^{4}}q = i{{\eta }^{3}}d\eta {{\sin }^{2}}\alpha d\alpha \sin \beta d\beta d\gamma \). The radial variable η takes values in the interval \((0, + \infty )\). After integration over the angles, the expression (8) becomes
All the terms in Eq. (6) vanish for \(5 \leqslant s < + \infty \), because \(\sin (\pi s) = 0\) for integer s and the integral in the η variable converges. There is \(\infty \times 0\) uncertainty for the values \(s = 1,3,4\). The ultraviolet divergence in the radial integral is eliminated by the regularization, in which case the terms \(s = 1,3,4\) similarly vanish.
The term \(s = 2\) is stable in the infrared region, as evidenced by the neutrino-pair exchange potential for zero neutrino mass [1]:
Additional contributions to the potential (10) arise from loops involving charged fermions of the Standard Model at distances closer than the electron’s Compton wavelength and from loops involving heavy bosons of the Standard Model at distances closer than the Z boson’s inverse mass.
The upper bounds on the neutrino masses [9, 10] are much lower than the effective neutrino potential \(U = - {{G}_{F}}\rho {\text{/}}\sqrt 2 \sim - 20\) eV. Under these conditions, the interaction potential (10) is modified as a result of neutrinos being trapped by the effective potential U to create a condensate. In old neutron stars with temperatures around 100 eV, only a fraction of the discrete neutrino levels are occupied in thermal equilibrium with the neutrons. The effect of neutrino condensate on the EoS of neutron matter is negligible [12, 14]. The experimental prospects of detecting the neutrino-pair exchange force in neutrino backgrounds are discussed in [15].
We examined the multiparticle contributions of the neutrino interactions to the thermodynamic potential. The infrared divergences discussed earlier in the literature are actually absent in each individual term of the decomposition (6) and, thus, in the sum. As a consequence, the limit \(N \to \infty \) is unnecessary, and keeping the number of particles large but finite is sufficient to provide the required proofs. All of the components in Eq. (6) vanish for massless and massive neutrinos, with the exception of the two-body interaction, which is negligible. Thus, under the one-loop approximation, long-range multiparticle interactions of neutrinos have no impact on the structure and stability of neutron stars.
REFERENCES
G. Feinberg and J. Sucher, Phys. Rev. 166, 1638 (1968).
S. D. H. Hsu and P. Sikivie, Phys. Rev. D 49, 4951 (1994).
J. A. Grifols, E. Masso, and R. Toldra, Phys. Lett. B 389, 563 (1996).
E. Fischbach, Ann. Phys. (N.Y.) 247, 213 (1996).
A. Segarra and J. Bernabeu, Phys. Rev. D 101, 093004 (2020).
C. Itzykson and J.-M. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980), p. 705.
M. I. Krivoruchenko, F. Simkovic, and A. Faessler, Phys. Rev. D 79, 125023 (2009).
D.-H. Wen, B.-A. Li, and L.-W. Chen, Phys. Rev. Lett. 103, 211102 (2009).
T. M. C. Abbott, M. Aguena, A. Alarcon, et al. (DES Collab.), Phys. Rev. D 105, 023520 (2022).
M. Aker, A. Beglarian, J. Behrens, et al. (KATRIN Collab.), Nat. Phys. 18, 160 (2022).
E. Fischbach, D. E. Krause, Q. le Thien, and C. Scarlett, arXiv: 2208.03790v1 [hep-ph] (2022).
As. Abada, M. B. Gavela, and O. Pinea, Phys. Lett. B 387, 315 (1996).
E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2: Theory of Condensed State, Vol. 9 of Course of Theoretical Physics (Pergamon, New York, 1980), p. 397.
A. Yu. Smirnov and F. Vissani, arXiv: 9604443v2 [hep-ph] (1996).
M. Ghosh, Yu. Grossman, W. Tangarife, X.-J. Xu, and B. Yu, arXiv: 2209.07082v2 [hep-ph] (2022).
ACKNOWLEDGMENTS
The author sincerely thanks the organizers of the ITEP Astrophysical Seminar for drawing his attention to the paper [11] and the participants of the seminar for interesting discussions. The manuscript’s preparation assistance from Nadezhda Lisitsyna is greatly appreciated.
Funding
The work was supported by the Russian Science Foundation, project no. 23-22-00307.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
The authors declare that they have no conflicts of interest.
Rights and permissions
Open Access. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Krivoruchenko, M.I. Long-range Multiparticle Interactions Induced by Neutrino Exchange in Neutron Star Matter. Jetp Lett. 117, 96–98 (2023). https://doi.org/10.1134/S0021364022603104
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1134/S0021364022603104