Hyperons in Finite and Inﬁnite Nuclear Systems

: In this work, we shortly review the role and properties of hyperons in ﬁnite and inﬁnite nuclear systems such as hypernuclei and neutron stars. Particularly, we describe different production mechanisms of hypernuclei, discuss some aspects of their γ -ray spectroscopy and their weak decay modes, and give a few strokes on their theoretical description. We reexamine also the role played by hyperons on the properties of neutron and proto-neutron stars with a special emphasis on the well-known “hyperon puzzle”, of which we discuss some of the solutions that have been proposed to tackle this problem. Finally, we review the role of hyperons on the cooling properties of newly born neutron stars and on the so-called r-mode instability.


Introduction
The presence of strange baryons, commonly known as hyperons, in hypernuclei and neutron stars permits the study of baryon interactions from an enlarged perspective, and the extension of our present knowledge of conventional nuclear physics to the SU(3)-flavor sector [1]. Hypernuclei, bound systems composed of neutrons, protons and one or more hyperons, were first observed in 1952 by Danysz and Pniewski in a balloon-flown emulsion stack where a hyperfragment was discovered [2]. Pion and proton beam production in emulsions and later in 4 He bubble chambers, where single-Λ hypernuclei were identified from the weak decay of the Λ hyperon into a proton and a π − , followed these initial cosmicray observations of hypernuclei. The advent of separated K − beams, which permitted the realization of counter experiments, lead to more systematic investigations of hypernuclei. A considerable amount of hypernuclear features, such as, e.g., the small spin-orbit strength of the hyperon-nucleon (YN) interaction or the fact that the Λ essentially retains its identity inside the nucleus, were revealed by in-flight (K − , π − ) counter experiments carried out at CERN and at Brookhaven National Laboratory (BNL). Experiments using (π + , K + ) and (K − stopped , π 0 ) reactions were conducted later at the Brookhaven AGS and KEK accelerators with higher intensities and improved energy resolution of the beams. A high-precision tool for the study of Λ-hypernuclear spectroscopy with resolutions of several hundred keV [3] is provided by the electromagnetic production of hypernuclei through the reaction (e, e K + ) carried out at the Thomas Jefferson National Laboratory (JLAB) and the Mainz Microtron Accelerator (MAMI-C). A promising new way to produce hypernuclei by using stable and unstable heavy-ion beams was proposed a few years ago by the HypHI Collaboration at FAIR/GSI [4], and it has recently allowed the observation of the Λ hyperon, and the 3 Λ H and 4 Λ H hypernuclei in a first experiment using a 6 Li beam on a 12 C target at 2 AGeV [5]. Today, thanks to the use of high-energy accelerators and modern electronic counters, more than 40 single-Λ hypernuclei, and a few double-Λ and single-Ξ ones have been identified. The existence of single-Σ hypernuclei, however, has not been experimentally confirmed yet without ambiguity, suggesting that the ΣN interaction is most likely repulsive.
In addition to hypernuclei, a big interest has being put in the study of hyperonic matter (nuclear matter with nucleonic and hyperonic degrees of freedom), especially in connection with the physics of neutron star interiors [6][7][8][9][10]. The density in the interior of neutron stars is large enough to allow for the appearance of new particles with strangeness content besides the conventional nucleons and leptons by means of weak interaction processes. Hyperons are expected to appear in neutron stars at around twice normal nuclear matter saturation density ρ 0 = 0.16 fm −3 . Neutron star properties are closely related to the underlying Equation of State (EoS) of matter at high densities. Therefore, despite the fact that hypernuclear matter is an idealized system, the theoretical determination of its EoS is an essential step towards the understanding of those neutron star properties which can be affected by the presence of hyperons. It is well known that the presence of hyperons softens the EoS and reduces the mass of neutron stars (see, e.g., [8,9]). In addition, hyperons can strongly influence also the thermal evolution and gravitational instabilities of these objects. The presence of hyperons, for instance, can modify the neutrino emissivity of dense matter and it can allow also for additional cooling mechanisms. Furthermore, hyperons dominate the bulk viscosity of matter as soon as they appear in the neutron star interior. Consequently, the emission of gravitational waves in hot and rapidly rotating neutron stars due to the so-called r-mode instability is affected also by their presence. Conversely, further constraints on the YN and hyperon-hyperon (YY) interactions can be provided by comparing the theoretical predictions for these properties with astrophysical observations.
A detailed knowledge of the EoS of hypernuclear matter over a wide range of densities is required to understand better the effect of hyperons on neutron stars. However, this is a very hard task. Two types of approaches have been traditionally used to describe baryon interactions in the nuclear medium and, to construct from them the (hyper)nuclear EoS: phenomenological and microscopic approaches. Relativistic or non-relativistic phenomenological approaches are based on effective density-dependent interactions which contain typically a certain number of parameters that are adjusted to reproduce (hyper)nuclear observables, and neutron star properties. Among the most commonly used ones, we can mention Skyrme-type interaction models and relativistic mean field (RMF) models. Several authors have used density-dependent baryon-baryon interactions based on Skyrme-type forces including hyperons to derive phenomenological EoSs of hyperonic matter [11][12][13][14][15][16][17][18][19][20][21]. Properties of nuclei and the experimental data from hypernuclei are employed within this approach to fix the in-medium nucleon-nucleon (NN), YN and YY interactions. RMF models are based on effective Lagrangian densities in which baryon-baryon interactions are described in terms of meson exchanges. A RMF description of the EoS of dense matter with hyperons is turning out today to be one of the most popular ones (see, e.g., [22][23][24][25][26][27]). The parameters of RMF models are usually determined by using, in the case of the nucleons, the properties of nuclei and nuclear bulk matter, and by employing symmetry relations and hypernuclear observables to fix the coupling constants of the hyperons with the mesons. The Quark-Meson-Coupling (QMC) model has also been employed to determine the EoS of (hyper)nuclear matter and the properties of neutron stars [28][29][30][31][32][33][34][35][36]. In the QMC model, baryons are treated as confined non-overlapping bags of three quarks where the interaction is modeled through the exchange of mesons between quarks from different bags, which are, in turn, modeled using the MIT bag model. The in-medium properties of hyperons have been also studied within the Non Linear Derivative (NLD) model, an alternative RMF approach which incorporates an explicit momentum dependence of the in-medium baryon optical potentials [37][38][39]. Microscopic approaches, on the other hand, are based on realistic two-body baryon-baryon interactions that describe the scattering data in free space. These interactions have been mainly constructed within the framework of the meson-exchange theory [40][41][42][43][44][45][46][47][48][49][50], although a new approach based on chiral perturbation theory has recently emerged as a powerful tool [51][52][53][54][55][56][57][58][59]. To obtain the EoS, one has to solve then the very complicated many-body problem. The main difficulty of this problem lies in the treatment of the repulsive core which dominates the short-range behavior of the baryon interaction. Different microscopic many-body methods have been extensively used to study nuclear matter, however, very few of them have been extended to the hypernuclear sector. To the best of our knowledge, the many-body methods extended to the strange sector include the Brueckner-Hartree-Fock (BHF) approximation [60][61][62][63][64][65][66][67] of the Brueckner-Bethe-Goldstone theory, the Hartree-Fock theory based on the soft V low k interactions [68] and the Dirac-Brueckner-Hartree-Fock (BHF) theory [69]. The Auxiliary Field Diffusion Monte Carlo method [70] was also extended to the hyperonic sector a few years ago. Very recently BHF calculations of hyperonic matter using YN interactions derived within SU(3) chiral effective field theory have also been done by the Jülich-Bonn-Munich group [71,72] and Kohno [73,74].
This work is not intended to be an exhaustive review but rather a short introduction to the several aspects of hypernuclear physics. Here, we briefly review the properties of hyperons in hypernuclei and neutron stars. Particularly, in Section 2, we shortly discuss different production mechanisms of hypernuclei, as well as some aspects of their γ-ray spectroscopy, their weak decay modes, and their theoretical description. In Section 3, we reexamine the role played by hyperons on the properties of neutron and proto-neutron stars with a special emphasis on the so-called "hyperon puzzle" and present some of the possible solutions that have been proposed to tackle it. We also reexamine in this section the effect of hyperons on the cooling properties of newly born neutron stars and on the development of the so-called r-mode instability. The manuscript is finished in Section 4 with a summary.

Hypernuclear Physics in a Nutshell
Compared to the NN interaction, the YN and YY ones are still poorly constrained due, mainly, to the scarce amount of YN scattering data and to the complete absence of them in the YY case. The reason of this scarce amount of data should be traced back to the experimental difficulties associated with the short lifetime of hyperons and the low-beam intensity fluxes. In addition to the scattering data, information on the YN and YY interactions can be obtained, using the so-called femtoscopy technique, by measuring the correlations (in momentum space) of Yp and YY pairs produced in heavy-ion collisions [75][76][77][78][79][80][81][82][83][84][85]. The ratio of the distribution of relative momenta k * = ( p 1 − p 2 )/2 between a correlated and uncorrelated pair defines the correlation function of a given baryon pair. If the interaction of a baryon pair is attractive, then the measured correlation function will be found to be larger than one. Conversely, if the interaction of the pair is repulsive, the correlation function will take values between zero and one. The correlation function between two baryons can be theoretically expressed as [86,87] where ψ( k * , r) is the relative wave function of the baryon pair of interest and S( r) is the so-called source function, and it represents the distribution of the distance | r| at which the particles are emitted. The comparison between the theoretical correlation function and the measured one permits the testing and improvement of the existing YN and YY potentials. Lattice QCD offers also a very powerful way to derive baryon-baryon interactions, see, e.g., [88][89][90][91][92][93] for reviews. A big progress in this direction has been made in the last years by the HALQCD [94][95][96][97][98][99][100][101][102] and the NPLQCD [103][104][105] collaborations. We should note, however, that the methods employed by these two collaborations are quite different. Whereas the HALQCD collaboration follows a method to extract the different baryonbaryon potentials from the Nambu-Bethe-Salpeter wave function measured on the lattice, the NPLQCD collaboration combines calculations of correlation functions at several lightquark-mass values with low-energy effective field theory (EFT). This second approach is particularly interesting because it allows the matching of lattice QCD results with lowenergy EFT providing, in this way, the means for first predictions in the physical quark mass limit. Results for various NN, NY and YY interaction channels at a single value of the lattice spacing and of the lattice volume have been recently obtained by the HALQCD collaboration which managed to approach the region of physical masses [100][101][102]. Very recently, the NPLQCD collaboration has studied the interaction between two octet baryons for strangeness S = 0, −1, −2, −3 and −4 at low energies using larger-than-physical quark masses corresponding to a pion mass of m π ∼ 450 MeV and a kaon mass of m K ∼ 596 MeV, and have extracted the corresponding values of s-wave scattering phase shifts, low-energy scattering parameters, and binding energies of two-baryon bound systems [105]. A detailed review of the last lattice QCD developments in the strangeness sector is beyond the scope of the present paper and, therefore, the interested reader is referred to the original works of both the HALQCD and the NPLQCD collaborations for further information.
Alternative and complementary information on the YN and YY interactions can be extracted from the study of hypernuclei. The main aim of hypernuclear physics [106,107] is, in fact, to relate hypernuclear observables with the underlying YN and YY interactions. In this section, we briefly review the different production mechanisms of hypernuclei, discuss some aspects of their γ-ray spectroscopy and their weak decay modes, and we finish with a few strokes on their theoretical description.

Production of Hypernuclei
Several reactions can be used to produce single-Λ hypernuclei. One of them is the so-called (K − , π − ) strangeness exchange reaction where a neutron of the nucleus target hit by a K − is changed into a Λ that remains bound to the nucleus and a π − is emitted. With this reaction, it is possible to determine accurately the mass and the binding energy of the formed hypernucleus by measuring the momenta of both the incoming K − and the outgoing π − using two magnetic spectrometers with good energy resolution. Strangeness exchange reactions, initially performed at CERN, have been mainly used later at BNL in the USA, and at KEK and J-PARC in Japan. The so-called (π + , K + ) associated production reaction is another production mechanism of hypernuclei that makes use of π + beams instead of K − ones, Here, an ss pair is created from the vacuum, and a K + and a Λ are produced in the final state. The production cross section of this reaction is smaller than that of the (K − , π − ) one, but this is compensated by the fact that the intensities of the π + beams are larger than those of the K − ones. The production of hypernuclei by means of these reactions have been also performed at BNL and KEK, and later at GSI in Germany. The use of electron beams with an excellent spatial and energy resolution has allowed the electroproduction of hypernuclei by means of the (e, e K + ) reaction, providing in addition a high-precision tool for the study of hypernuclear spectroscopy, with energy resolutions of several hundred keV [3]. Electroproduction of hypernuclei is carried out at JLAB in the USA and the MAMI-C laboratory in Germany. These two laboratories are presently the only ones with the instrumental capabilities required to perform these kind of experiments. The kinematics of the elementary processes n(K − , π − )Λ, n(π + , K + )Λ and p(γ, K + )Λ underlying the three production mechanisms of single Λ-hypernuclei discussed above is shown in Figure 1, adapted from [108]. Note that the momentum transferred to the Λ is much lower in the case of the n(K − , π − )Λ reaction than in the other two. The lower the momentum transferred to the Λ is, the larger its probability of interacting with, or being bound to, the nucleus will be. In addition, the lower the momentum transferred is, the smaller the angular momentum transfer will also be and, consequently, the Λ will more easily retain the quantum numbers of the nucleon that has been eliminated in the reaction. Therefore, in the case of the n(π + , K + )Λ or p(γ, K + )Λ reactions, since the recoil momentum of the hyperon is high, the cross sections to bound states are reduced, and the produced Λ has a higher probability of escaping the nucleus.
Similar reactions can be used to produce single-Σ hypernuclei. However, as we have already said, their existence has not been yet experimentally confirmed without ambiguity. The production of double-Λ hypernuclei requires a two-step mechanism in which, first, a Ξ − is created by means of reactions such as and second, the Ξ − is captured in an atomic orbit and interacts with the nuclear core, producing two Λ's by hitting one of the protons of the nuclei This process releases an energy of about 30 MeV that, in most of the cases, is equally shared between the two Λ's, leading to the escape of one or both hyperons from the nucleus. Double-Λ hypernuclei are currently the best systems to investigate the properties of the baryon-baryon interaction in the strangeness S = −2 sector. The ΛΛ bond energy ∆B ΛΛ in double-Λ hypernuclei can be experimentally determined by measuring of the binding energies of double and single-Λ hypernuclei as A few double-Λ hypernuclei, 6 ΛΛ He, 10 ΛΛ Be and 13 ΛΛ B, have been reported in emulsion experiments. A quite large value of the ΛΛ bond energy of around 4 to 5 MeV was deduced from the subsequent analysis of these emulsion experiments. We should also note that the identification of some of these double-Λ hypernuclei was ambiguous. Therefore, careful attention should be paid when using the data from this old analysis to put any kind of constraint on the ΛΛ interaction. However, a new 6 ΛΛ He candidate with a ΛΛ bound energy ∆B ΛΛ = 1.01 ± 0.2 +0. 18 −0.11 MeV (recently corrected to ∆B ΛΛ = 0.67 ± 0.17 MeV) was observed without ambiguity in 2001 at KEK [109].
The reactions (5) and (6) can be used to produce single-Ξ hypernuclei and, in fact, a few of them have been identified. The analysis of the experimental data from reactions such as 12 Recently, however, Friedman and Gal [111] have analyzed several Ξ − p → ΛΛ two-body capture events in 12 C and 14 Ni emulsion nuclei, concluding that the Ξ-nuclear interaction is strongly attractive, with a Ξ − potential depth in nuclear matter V Ξ − ≥ 20 MeV. We should mention here the observation [112] of a deeply bound state of the Ξ − -14 N system with a binding energy of 3.87 ± 0.21 MeV [113]. This event provides the first clear evidence of a deeply bound state of this system by an attractive ΞN interaction. The latest experimental data obtained by the J-PARC E07 collaboration [114] indicate a value for the binding energy of the Ξ − in the Ξ − -14 N system of 1.27 ± 0.21 MeV. Future Ξ hypernuclei experiments are being planned at J-PARC.

γ-ray Spectroscopy of Hypernuclei
Excited states of hypernuclei can be produced when a nucleon in a p or a higher shell is replaced by a hyperon. The energy of the hypernuclear excited states can be released either by emitting nucleons or, sometimes, γ-rays. The analysis of hypernuclear excited states with very good energy resolution has been possible thanks to the detection of γray transitions in single-Λ hypernuclei. The construction of large-acceptance germanium detectors, dedicated to hypernuclear γ-ray spectroscopy, has overcome some of the initial technical difficulties found in the application of γ-ray spectroscopy to hypernuclei. These difficulties were mostly associated with the detection efficiency of γ-ray measurements and with the necessity of covering a large solid angle with γ-ray detectors. Several weak points in hypernuclear γ-ray spectroscopy, however, still persist such as, for instance, the fact that the observation of γ-rays is mostly limited to the low excitation region, maybe up to the Λ p-shell. The reason is that a number of single-particle Λ states are bound in heavy Λ hypernuclei with a potential depth of ∼28 MeV but the energy levels of many single-particle states are above the neutron and proton emission thresholds. Another weak point is clearly the fact that the γ-ray transition only measures the energy difference between two states and this single energy information is not enough to fully identify the two levels. This problem can be solved, of course, by measuring two γ-rays in coincidence. We show in Figure 2 the energy of a Λ hyperon in the single-particle states s, p, d, f and g of several hypernuclei, deduced from emulsion, (K − , π − ) and (π + , K + ) reactions, as a function of the mass number to the power −2/3. The value of ∼28 MeV extrapolated at A −2/3 = 0 is usually interpreted as the binding energy of a single Λ hyperon in infinite symmetric nuclear matter at saturation density, and it is used to fix the parameters of the majority of the models of the hyperonic EoS. Systematic spectroscopic studies of single Λ hypernuclei indicate that the ΛN interaction is attractive [115]. Energy of a Λ hyperon in the single-particle states s, p, d, f and g of several hypernuclei as a function of A −2/3 deduced from emulsion, (K − , π − ) and (π + , K + ) reactions. The lines are drawn just to help the reader.
To finish this section, in Figure 3, we show, as an example, the level scheme and γ-ray transitions of 16 Λ O identified and determined by γ-ray spectroscopy using the (K − , π − ) reaction and the germanium detector array Hyperball at BNL [116]. The twin peaks observed confirm the hypernuclear fine structure for the (1 − → 1 − ) and (1 − → 0 − ) transitions in 16 Λ O. We note that the small spacing between the twin peaks is due to the spin dependence of the ΛN interaction.

Weak Decay of Hypernuclei
The so-called mesonic weak decay is the main decay mode of the Λ hyperon in free space, where ∼60% of the times the Λ decays into a proton and a π − , and ∼40% into a neutron and a π 0 . When the Λ is bound in the nucleus, however, this mode is strongly suppressed by the Pauli exclusion principle since the momentum of the outgoing nucleon (∼100 MeV/c) is smaller than the typical Fermi momentum of a nucleon in the nucleus (∼270 MeV/c). Consequently, in hypernuclei (specially in medium and heavy ones), the dominant decay mode becomes the non-mesonic one where the Λ interacts with one (or more) of the surrounding nucleons. The weak decay of hypernuclei has been mainly studied within the framework of meson-exchange models [117,118] and, more recently, using the effective field theory [119,120]. In [121,122], the interested reader can find two comprehensive reviews on the theoretical aspects of hypernuclear weak decay.
In Figure 4, we present the weak decay rate Γ (expressed in units of the decay rate of the Λ in free space) as function of the total number of particles A + 1. The figure has been adapted from the original one in [123]. The dot, dashed and solid lines show, respectively, the theoretical predictions of the mesonic Γ M , non-mesonic Γ N M and total Γ T decay rates. The curves labeled Γ 1 and Γ 2 correspond to the contributions of one-nucleon and twonucleon induced decay modes to the non-mesonic decay rate (see Equations (10) and (11)). Experimental values of the total and non-mesonic decay rates are given, respectively, by the squares and circle marks. As it can be seen in the figure, the mesonic decay mode gets blocked as A increases, while the non-mesonic decay increases up to a saturation value of the order of the free decay, reflecting the short-range nature of the weak ∆S = 1 baryon-baryon interaction. 10 100 Total number of particles A+1  Figure 4. Weak decay rate Γ as a function of the total number of particles in units of the weak decay rate of the Λ in free space Γ free Λ . Dot, dashed and solid lines show, respectively, the theoretical predictions of the mesonic Γ M , non-mesonic Γ N M and total Γ T decay rates. Dot-dashed lines labeled Γ 1 and Γ 2 display the contributions of one-nucleon and two-nucleon induced decay modes to the non-mesonic decay rate (see Equations (10) and (11)). Experimental values of the total and nonmesonic decay rates are given by the squares and circle marks, respectively. Figure adapted from the original one in [123].

Theoretical Description of Hypernuclei
A simple theoretical description of a hypernucleus consists of an ordinary nucleus with a hyperon sitting in a single particle state of an effective hyperon-nucleus mean field potential. Based on this simple description, several approaches have been followed to derive the properties of hyperons in finite nuclei. Traditionally, Woods-Saxon potentials have been used, for instance, to describe in a shell model picture the single-particle properties of the Λ from medium to heavy hypernuclei [124][125][126][127]. To improve the overall fit of the Λ single-particle energies, non-localities and density-dependent effects have been included in non-relativistic Hartree-Fock calculations with Skyrme type YN interactions [13][14][15][16][17][18][19][20][21]. Relativistic mean field theory [24,[128][129][130][131][132][133][134][135][136][137] and Dirac phenomenology [138,139] have been also employed to perform hypernuclear structure calculations. Several hypernuclear structure studies based on ab initio approaches do also exist in the literature [140][141][142][143][144][145][146][147]. The single-particle properties of the Λ in the hypernucleus are derived in these studies from effective YN G-matrices built from bare YN interactions which describe the scarce scattering data in free space. A Quantum Monte Carlo calculation of single-and double-Λ hypernuclei has also been recently done using two-and three-body forces between the Λ and the nucleons [70,148]. We would like to note that the NPLQCD collaboration has been able to obtain the binding energies of the light hypernuclei including 3 Λ He, 4 Λ H and 4 ΛΛ He [149]. The quality of the description of hypernuclei in most of these approaches relies on the validity of the mean field picture. Correlations induced by the YN interaction can, however, change substantially this picture and, therefore, should not be ignored. While many authors have extensively studied the correlations of nucleons in nuclear matter and finite nuclei, those of hyperons have not received so much attention so far. The effect of the Λ correlations in nuclear matter, beyond the mean field description, was studied for the fist time by Robertson and Dickhoff [150] using the Green's function formalism. These authors calculated the spectral function and quasi-particle parameters of the Λ finding results qualitatively similar to those of the nucleons. They showed that the Λ is, in general, less correlated than the nucleons. A few years ago, the author of the present review studied the spectral function of the Λ hyperon in finite nuclei [151], showing, in agreement with the work of Robertson and Dickhoff, that the Λ is less correlated than the nucleons, and confirming the idea that it maintains its identity inside the nucleus. The results of this study showed also that in hypernuclear production reactions, the Λ hyperon is formed mostly in a quasi-free state.
As an example of a theoretical calculation of hypernuclei, we briefly describe here a microscopic method that allows one to determine the single-particle bound states of a Λhyperon in finite nuclei. This method starts with the construction of all the YN G-matrices which describe the interaction between a hyperon and a nucleon in infinite nuclear matter. To this end, the coupled-channel Bethe-Goldstone equation is solved. These G-matrices are then used to obtain the YN G-matrices in finite nuclei through the following integral equation: which expresses the finite nuclei G-matrices, G FN , in terms of the nuclear matter ones, G, and the difference between the finite-nucleus and the nuclear-matter propagators, written schematically as (Q/E) FN − (Q/E). This difference, which accounts for the relevant intermediate particle-particle states has been shown to be quite small and, thus, in all practical calculations, G FN can be well approximated by truncating the expansion of Equation (12) up to second order in the nuclear matter G-matrices. Therefore, we have Using then G FN as an effective YN interaction, one can obtain the Λ self-energy in the BHF approximation (see diagram (a) of Figure 5). This approximation can be split into the sum of two contributions: the one shown by diagram (b), which originates from the first-order term on the right-hand side of Equation (13), and that of diagram (c), which stands for the so-called two-particle-one-hole (2p1h) correction, where the intermediate particle-particle propagator has to be viewed as the difference of propagators appearing in Equation (13). Solving finally the Schrödinger equation with the real part of the Λ selfenergy, it is then possible to determine, as mentioned before, the different Λ single-particle bound states. Further details of this method can be found, e.g., in [146,151]. As an example of the application of this method, in Table 1, we show the energies of the Λ single-particle bound states in several hypernuclei. The results have been obtained with the NLO13 [58] and the NLO19 [59] chiral YN interaction of the Jülich-Bonn-Munich group for different values of the cutoff of the interaction including both contributions (first order term and 2p1h correction) to the Λ self-energy. We note that, due to technicalities of this method, it can be only applied to hypernuclei consisting of a closed-shell nuclear core plus a Λ sitting in a single-particle state. The values reported are to be compared with the experimental separation energies for the corresponding hypernuclei. However, since experimental data for the particular hypernuclei we consider do not always exist, the comparison is done with the closest representative hypernuclei for which experimental information is available. As it can be seen, in general, there is an underbinding of light hypernuclei such as 5 Λ He and 13 Λ C while the description of medium and heavy hypernuclei improves. Note, however, that neither the NLO13 interaction nor the NLO19 one yield a quantitative description of all medium and heavy hypernuclei. Whereas the NLO19 interaction describes reasonably well 17 Λ O, 41 Λ Ca and 91 Λ Zr, it seems to overbind 209 Λ Pb. For the latter, the predictions of the NLO13 interaction are more in line with the experiment. We note that the results of the calculation agree with the experimental fact that the spin-orbit splitting of the p-, d-, f -and g-states is very small. The interested reader is referred to [147] for further detail of these results.

Hyperons and Neutron Stars
The presence of hyperons in the interior of neutron stars has been considered by many authors for more than 60 years since the seminal work of Ambartsumyan and Saakyan [154]. Different phenomenological or microscopic models have been used to describe the neutron star matter EoS. All the works show that hyperons may appear in the neutron star interior at densities around 2 − 3ρ 0 . The reason for their appearance is simply that, at such densities, the nucleon chemical potential is large enough to make the conversion of nucleons into hyperons energetically favorable. As a result of this conversion, the Fermi pressure exerted by nucleons is relieved and, therefore, the EoS becomes softer. Consequently, the mass of the star, and particularly, its maximum value M max is reduced. How much the EoS is softened, and how much M max is reduced depends on the attractive or repulsive character of the YN and YY interactions. In general, attractive (repulsive) interactions lead to an earlier (later) onset and larger (smaller) concentration of hyperons and, thus, to a stronger (more moderate) softening of the EoS and a larger (smaller) reduction of M max . It is well known (see, e.g., [66,67]), however, that hyperons equalize the effect of different nucleonic interactions through several compensation mechanisms: a stiffer nucleonic EoS will lead to an earlier onset of hyperons enhancing in this way the softening due to their presence. Conversely, a retarded onset of a certain hyperon species will favor the appearance of other species leading also to a softer EoS. As a result, M max is surprisingly quite insensitive to the pure nucleonic EoS, and even to the details of the YN and YY interactions (see, e.g., Figure 2 in [66] or Figure 3 in [67]).

The Hyperon Puzzle and Some Possible Solutions
Although the presence of hyperons in neutron stars seems to be energetically unavoidable, the strong softening of the EoS associated with the onset of hyperons (notably in microscopic models) leads to values of M max not compatible with observations. This controversy is known in the literature as the "hyperon puzzle", and is currently a subject of intensive research. The discrepancy between the theory and observations became more dramatic after the measurements in the last decade of unusually high masses of the millisecond pulsars PSR J1903+0327 (1.667 ± 0.021M ) [155], PSR J1614-2230 (1.928 ± 0.017M ) [156], PSR J0348+0432 (2.01 ± 0.04M ) [157], and the most massive one observed up to now PSR J0740+6620 (M > 2.14 +0. 10 −0.09 M ) [158], which rule out almost all currently proposed EoS with hyperons (both microscopic and phenomenological).
To solve this puzzle, a mechanism (or mechanisms) is necessary that could provide the additional repulsion needed to make the EoS stiffer and, therefore, M max compatible with the current observational limits. Three possible mechanisms that could provide such additional repulsion are: (i) more repulsive hyperon-hyperon interactions driven by either repulsive vector mesons exchanges [22,23,[25][26][27] or less-attractive scalar σ meson exchange [24] (ii) repulsive hyperonic three-body forces [159][160][161][162][163][164][165][166][167] and (iii) a phase transition to deconfined quark matter at densities below the hyperon threshold [168][169][170][171][172][173][174][175][176][177][178][179]. The possible appearance of other hadronic degrees of freedom, such as the ∆ isobar or meson condensates, that can push the onset of hyperons to higher densities, has been also considered. An interesting possibility to circumvent the problem is the so-called two-families scenario proposed in [180,181], in which stars made of hadrons are stable only up to (1.5 − 1.6)M while most massive compact stars are entirely made of strange quark matter. In the following, we briefly review some of these possible solutions.

Hyperon-Hyperon Repulsion
This possible solution has been explored mostly in the context of RMF models [22][23][24][25][26][27], and it is based on the well-known fact that, in meson-exchange models of nuclear forces, while the scalar meson σ is responsible for the intermediate-range attraction, vector mesons generate repulsion at short distances. It has been argued that if the interaction between two hyperons, driven by vector mesons, is repulsive enough or the attraction mediated through the exchange of the σ meson is weak enough, then the EoS could be sufficiently stiff to reconcile the current high pulsar mass observations with the existence of hyperons in neutron stars. However, the strength of meson exchanges should be modified consistently with hypernuclear data which requires, at least, the ΛN interaction to be attractive and suitably tuned to the hypernuclear data [115]. Such tuning is not required if the repulsive vector meson interactions act only among the hyperons through the exchange of the strange φ vector meson (which couples only to hyperons). In this way, the onset of hyperons is shifted to higher densities and it is possible to obtain neutron stars with maximum masses larger than 2M and a significant amount of hyperons in their interior.

Hyperonic Three-Body Forces
It is well known that three-nucleon forces are fundamental to reproduce accurately the properties of few-nucleon systems as well as the empirical saturation point of symmetric nuclear matter in non-relativistic many-body approaches. It seems natural, therefore, to suggest that three-body forces of the type NNY, NYY and YYY, involving one or more hyperons, could provide, as in the case of three-nucleon forces, the additional repulsion needed at high densities to make the EoS stiff enough, solving in this way the hyperon puzzle. This idea was suggested even before the observation of neutron stars with ∼2M (see e.g., [159]), and it has been explored by a number of authors in the last years [160][161][162][163][164][165][166][167]. However, no general consensus has been reached yet regarding the role played by the hyperonic three-body forces in solving the hyperon puzzle. A multi-Pomeron exchange potential (MPP) model to introduce a universal three-body repulsion among three baryons in the hyperonic matter EoS was proposed in [162,163,165]. This universal three-body repulsive potential was based on the extended soft core (ESC) baryon-baryon interaction of the Nijmegen group [48][49][50]. The strength of the MPP was determined by analyzing the nucleus-nucleus scattering with the use of a G-matrix folding potential derived from the ESC interaction complemented with the MPP and a three-nucleon attractive part, added phenomenologically in order to reproduce the nuclear saturation properties. The results of those works [162,163,165] showed that when the MPP contribution was taken into account, universally for all baryons, neutron star radii R at a typical mass 1.5M were predicted to be around 12.3-13.1 km [166], and a maximum mass of ∼2.2M was obtained. This result for the maximum mass is in contradiction with that reported in [161] where the case of a universal three-body repulsion was also analyzed. The authors of [161] used a model based on the BHF approach of hypernuclear matter using the Argonne V18 NN potential [182] and the Nijmegen YN soft core NSC89 one [44] supplemented with additional simple phenomenological density-dependent contact terms to establish numerically lower and upper limits to the effect of hyperonic three-body forces on the maximum mass of neutron stars. Assuming that the strength of these forces was either smaller than or as large as the pure nucleonic ones, the results reported in [161] showed that, although the employed hyperonic three-body forces stiffened the EoS, they were, however, unable to provide the repulsion needed to make the predicted maximum masses compatible with the recent observations. In [164], a Monte Carlo calculation of pure neutron matter with a nonvanishing Λ-hyperon concentration was carried out including NN, NNN, ΛN and NNΛ two-and three-body interactions. In particular, the NNΛ force used in this work was tuned in order to provide a reasonable description of the measured Λ separation energy of several hypernuclei [70]. The authors of this work concluded that, with the model they considered, the presence of hyperons in the core of neutron stars could not be satisfactorily established and, consequently, according to these authors, there is no clear incompatibility with astrophysical observations when the Λ is included. However, one should note, that the presence of protons, necessary to establish the correct β-equilibrium inside neutron stars and, thus, a proper treatment of nuclear matter was neglected in their calculation. Very recently, the authors [167] have studied the effects of chiral hyperonic three-body forces on neutron stars, showing that the inclusion of a moderate repulsive NNΛ force leads already to an EoS stiff enough such that the resulting neutron star maximum mass is compatible with the largest currently measured massed. These authors have also examined the effect of this NNΛ force on the separation energy of a Λ in some hypernuclei obtaining a good description of the experimental data. This is in agreement with the results of [183], where also a moderate hyperonic three-body force was found to be enough to reproduce the binding energies of single-Λ hypernuclei, but contrary to the results reported in [70] where a very repulsive NNΛ was required to describe the data. In conclusion, it seems that although hyperonic three-body forces offer an interesting microscopic solution to the hyperon puzzle, the uncertainties associated to these forces are still too large to allow for a definite conclusion.

Quark Matter Phase Transition below the Hyperon Threshold
An early phase transition from hadronic matter to deconfined quark matter at densities below the hyperon threshold could provide the solution to the hyperon puzzle, as it has been suggested by several authors. In this case, massive stars could actually be hybrid stars with a stiff quark matter core. This solution, however, leads to a new question: can quarks provide the sufficient repulsion required to produce a 2M neutron star? In [168,169], the maximum mass was found to be in a relatively narrow interval, ∼1.4M < M max <∼ 1.7M which is incompatible with the current observational limit of 2M . To obtain values of M max larger than 2M , the interaction among quarks should fulfill two important and necessary conditions. First of all, it should be significantly repulsive, for example, in vector channels in order to guarantee that the EoS is stiff enough. Second, it should be strongly attractive in certain channels, leading to color superconductivity, so that the deconfined quark-matter phase is energetically favorable with respect to the hadronic one. Current theoretical descriptions of quark matter at high density rely on phenomenological models, which are constrained using the few available experimental information on high density baryonic matter from heavy-ion collisions. Several models of hybrid stars have been proposed in the recent years with the necessary properties to generate a 2M (see, e.g., [171][172][173][174][175][176][177][178][179]). Very recently, Shahrbaf et al. [178,179] have applied, for the first time, the finite-range polynomial interpolation method of Masuda et al. [184,185] for constructing a transition between hadronic and quark matter phases. The predicted maximum mass of the hybrid star is about 2.2M in agreement with current observations. The observation of 2M neutron stars, on the other hand, may also help constrain better the models of hybrid and strange stars (compact stars completely made of deconfined u, d and s quark matter), and improve our present understanding of the hadron-quark phase transition.

∆ Isobar and Kaon Condensation in Neutron Stars
An alternative possible solution to the hyperon puzzle that has been also considered is the appearance of other hadronic degrees of freedom, such as for instance the ∆ isobar or meson condensates, that can push the onset of hyperons to higher densities. The presence of the ∆ isobar in neutron stars was usually neglected because in former calculations its threshold density was found to be larger than the typical densities of neutron star cores. A few years ago, it has been shown [186], however, that the onset density of the ∆ isobar depends crucially on the slope parameter L of the nuclear symmetry energy. Using recent experimental constraints and a state-of-the-art EoS, the authors of [186] have shown that the ∆ isobar could actually appear at densities below the hyperon threshold. However, they found that, as soon as the ∆ appears, the EoS becomes also considerably soft and, as a consequence, M max is reduced to values below the current observational limits, giving rise to what they have named the "∆ puzzle". Very recently, using an RMF approach, the authors of [187] found, however, that the presence of the ∆ isobar in the neutron star interior is compatible with the observation of 2 M millisecond pulsars, provided that the couplings of the ∆ to the σ and ω meson fields are at least 10% stronger than the corresponding ones of the nucleons. The possibility that pions or kaons form Bose-Einstein condensates in the interior of neutron stars has also been extensively considered in the literature [188][189][190][191][192][193]. However, as in the case of the hyperons or the ∆ isobar, the appearance of a meson condensate induces also a strong softening of the EoS reducing the value of M max below the current observational limits too.

Effect of Hyperons on Proto-Neutron Stars
Thermal effects and neutrino trapping affect the properties of newly born neutron stars during the first tens of seconds after their formation. In particular, the composition and the overall stiffness of the EoS of the star are strongly influenced by these two effects. Matter becomes more proton rich, the number of muons is significantly reduced and the onset of hyperons is shifted to higher densities [194][195][196][197][198]. Furthermore, the number of strange particles is on average smaller, and the EoS is stiffer in comparison with the cold and neutrino-free case.
A very important consequence of the neutrino trapping in dense matter is the possibility of having metastable neutron stars and a delayed formation of a "low mass" (M = 1 − 2M ) black hole. To illustrate this, we show in Figure 6 the gravitational mass M G of the star as a function of its baryonic mass M B (proportional to the total number of baryons in the star) obtained by the authors of [197]. When hyperons are present in the star (Figure 6a), the deleptonization lowers the range of gravitational masses that can be supported by the EoS from about 1.59M to about 1.28M (see dotted horizontal lines in the figure). The neutron star baryonic mass can be considered constant during the evolution from the initial proto-neutron star configuration to the final neutrino-free one because most of the matter accretion on the forming neutron star happens in a very early stage after its birth (t < 1 s). For this particular calculation, proto-neutron stars born with gravitational masses between 1.28M and 1.59M (a baryonic mass between 1.40M and 1.72M ) will be stabilized by neutrino trapping effects long enough to carry out nucleosynthesis accompanying a Type-II supernova explosion. After neutrinos leave the star, the EoS is softened and it cannot support anymore the star against its own gravity and the newborn star collapses then to a black hole [194][195][196]. Conversely, when only nucleons are considered to be the relevant baryonic degrees of freedom (Figure 6b), no metastability occurs and a black hole is unlikely to be formed during the deleptonization because the gravitational mass increases during this stage which happens at (almost) constant baryonic mass. If a black hole were to form from a star with only nucleons, it is much more likely to form during the post-bounce accretion stage.

Hyperons and Neutron Star Cooling
The cooling of a newly born hot neutron star is driven in a first stage by the neutrino emission from the interior, and later by the emission of photons at the surface. Depending on whether the number of involved baryons is one or two, neutrino emission reactions can be classified as fast or slow processes, respectively. The simplest possible neutrino emission reaction is the well-known direct Urca process n → p + l +ν l , p + l → n + ν l . (14) This is a fast mechanism which, due to momentum conservation, is only possible when the proton fraction exceeds a critical value x DURCA ∼ 11% to 15% [199]. Other neutrino reactions that lead to medium or slow cooling scenarios, but that are operative at any density and proton fraction, are the modified Urca processes: the bremsstrahlung: or the Cooper pair formation: this last operating only when the temperature of the star drops below the critical temperature for neutron superfluidity or proton superconductivity. The presence of hyperons can affect the cooling of neutron stars because they can modify neutrino emissivities allowing for additional cooling mechanisms such as, for example, the direct and modified hyperonic Urca processes. These additional neutrino emission reactions, however, can lead to temperatures at the surface of the star that are much lower than those observed, unless they are suppressed by hyperon pairing gaps. The study of hyperon superfluidity becomes then of particular interest because it could play a key role in the thermal history of neutron stars. Nevertheless, whereas the presence of superfluid neutrons in the inner crust of neutron stars, and that of superfluid neutrons together with superconducting protons in their quantum fluid interior are well established and have been the subject of many studies, a quantitative estimation of the hyperon pairing has not received so much attention, and just few calculations exist in the literature [200][201][202][203][204][205][206].

Hyperons and R-Modes
It is well known that if a neutron star rotates with a spin frequency above the so-called Kepler frequency Ω K , i.e., the absolute maximum rotational frequency of a neutron star, then matter is ejected from the star's equator [207,208]. Different types of perturbations, however, can lead to instabilities that prevent the star from reaching rotational frequencies as high as Ω K and set more stringent limits on their rotation [209]. One of these possible instabilities which is particularly interesting is the so-called r-mode instability [210,211]. It is a toroidal mode of oscillation whose restoring force is the Coriolis force which leads to the emission of gravitational waves in hot and rapidly rotating neutron stars through the Chandrasekhar-Friedman-Schutz mechanism [212][213][214]. The r-mode grows with the emission of gravitational waves whereas it is damped by dissipation mechanisms. The r-mode becomes unstable when the driving time of the gravitational radiation is shorter than the damping time associated to viscous processes. A rapidly rotating neutron star could transfer in this case a significant fraction of its angular momentum and rotational energy to the emitted gravitational waves. Bulk (ξ) and shear (η) viscosities are the main dissipation mechanisms of r-and other pulsation modes. At high temperatures (T > 10 9 K), the bulk viscosity is the dominant mechanism and, therefore, it is the most important one for hot young neutron stars. It results from the variations in the pressure and the density induced by the oscillation mode which drives the star away from β-equilibrium. The weak interaction tries then to reestablish the equilibrium with the consequent dissipation of energy. If no hyperons or other exotic components are present in the neutron star interior, the bulk viscosity is mainly determined by the direct and modified Urca processes. However, as soon as hyperons appear, new mechanisms such as strong interaction reactions weak non-leptonic hyperon reactions or direct and modified hyperonic Urca (see Equations (18) and (19)) contribute to the bulk viscosity and dominate it for densities above 2-3 times the nuclear saturation density. Hyperon bulk viscosity has been considered by several authors, see, e.g., [215][216][217][218][219][220][221][222][223][224][225][226][227][228][229].
The time dependence of an r-mode oscillation is given by e iωt−t/τ(Ω,T) where ω is the frequency of the mode and 1 τ(Ω, T) is an overall timescale which describes both its exponential growth due to the gravitational wave emission and its decay due to the viscous damping, being τ GW , τ ξ and τ η the time scales associated, respectively, to the gravitational wave emission and the bulk and shear viscosity dampings. If τ GW is shorter than both τ ξ and τ η , the mode will exponentially grow, whereas in the opposite case, it will be quickly damped away. For each star at a given temperature T, one can define a critical angular velocity Ω c as the smallest root of the equation 1 τ(Ω c , T) = 0 .
This equation defines the boundary of what is usually called the r-mode instability region. If the angular velocity of a neutron star is smaller than its corresponding Ω c , then the star is stable against the r-mode instability. Conversely, a star with an angular velocity larger than Ω c will develop an instability that will cause a rapid loss of angular momentum through gravitational radiation until its angular velocity falls below the critical value. As an example, on panel (a) of Figure 7, we show the r-mode instability region for a pure nucleonic (black solid line) and a hyperonic (red dashed line) star with 1.27M both described using the BHF theory with realistic baryon interactions. On panel (b) of the figure, we show the contributions from direct and modified nucleonic Urca processes as well as from the weak non-leptonic process n + n ↔ p + Σ − included in the calculation of the bulk viscosity. As it is clearly seen, the r-mode instability is smaller for the hyperonic star. The reason is simply the increase of the bulk viscosity due to the presence of hyperons which makes the damping of the r-mode more efficient. Contributions from the direct and modified nucleonic Urca processes as well as from the weak non-leptonic process n + n ↔ p + Σ − are included.

Summary
In this work, we have briefly reviewed several topics related with the physics of hypernuclei and hyperons in neutron stars. In particular, we have revised the different production mechanisms of hypernclei as well as some aspects of their γ-ray spectroscopy and their weak decay modes. We have also given a few strokes on their theoretical description. We have discussed also the main effects of hyperons on the properties of neutron and proto-neutron stars with an emphasis on the well known "hyperon puzzle", the problem of the strong softening of the EoS of dense matter due to the appearance of hyperons which leads to maximum masses of compact stars that are not compatible with the recent observations of approximately 2M millisecond pulsars. We have shortly reexamined some of the different solutions that have been proposed to tackle this problem. The first of the mechanisms we have revised consists of the inclusion of a repulsive hyperonhyperon interaction through the exchange of vector mesons and it has been mainly explored in the context of RMF models. However, since presently there are not enough experimental data to constrain the YN and YY interaction accurately, it is not clear whether the large repulsion invoked in such models is realistic. The second one requires the inclusion of repulsive hyperonic three-body forces. However, at present, it is still an open issue whether these forces can, by themselves, solve completely the hyperon puzzle or not, although, it seems that even if they cannot provide the full answer they can contribute to it in an important way. The third possible solution to the problem we have reviewed is to consider the possibility of a phase transition to deconfined quark matter at densities below the hyperon threshold. However, the description of the quark phase via phenomenological models also suffers from uncertainties. Most models of hybrid stars unanimously agree that to construct massive neutron stars, it is required both a sufficiently stiff hadronic EoS as well as a color superconducting quark phase with strong interaction among quarks to provide sufficient repulsion. We have also briefly discussed a possible solution to the hyperon puzzle by invoking the appearance of other hadronic degrees of freedom such as the ∆ isobar or meson condensates that push the onset of hyperons to higher densities. Finally, we have also revised how the presence of hyperons can affect the cooling of neutron stars and the r-mode instability window through modifications of the microscopic input of the weak interaction rates and transport coefficients, such as the bulk viscosity, of dense matter.
Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Conflicts of Interest:
The author declares no conflict of interest.