We consider the thermal conductivity, shear viscosity, and momentum relaxation rates in the nucleon cores of the neutron stars. We study how the choice of the nuclear interaction and the model for three-body forces may affect these transport coefficients calculated within the Brueckner-Hartree-Fock many-body nuclear theory. We find that at relatively large densities the model dependence of the results is substantial. In addition we provide the analytical approximations which allow to incorporate our results in practical simulations.

We calculate modified Urca neutrino emission rates in the dense nuclear matter in neutron star cores. We find that these rates are strongly enhanced in the beta-stable matter in regions of the core close to the direct Urca process threshold. This enhancement can be tracked to the use of the in-medium nucleon spectrum in the virtual nucleon propagator. We describe the in-medium nucleon scattering in the non-relativistic Bruckner-Hartree-Fock framework taking into account two-body as well as the effective three-body forces, although the proposed enhancement does not rely on a particular way of the nucleon interaction treatment. Finally we suggest a simple approximate expression for the emissivity of the n-branch of the modified Urca process that can be used in the neutron stars cooling simulations with any nucleon equation of state of dense matter.

Thermal conductivity and shear viscosity of npe$\mu$ matter in non-superfluid neutron star cores are considered in the framework of Brueckner-Hartree-Fock many-body theory. We extend our previous work (Shternin, Baldo and Haensel, 2013) by analysing different nucleon-nucleon potentials and different three-body forces. We find that the use of different potentials leads up to one order of magnitude variations in the values of the nucleon contribution to transport coefficients. The nucleon contribution dominates the thermal conductivity, but for all considered models the shear viscosity is dominated by leptons.

The nuclear symmetry energy characterizes the variation of the binding energy as the neutron to proton ratio of a nuclear system is varied. This is one of the most important features of nuclear physics in general, since it is just related to the two component nature of the nuclear systems. As such it is one of the most relevant physical parameters that affect the physics of many phenomena and nuclear processes. This review paper presents a survey of the role and relevance of the nuclear symmetry energy in different fields of research and of the accuracy of its determination from the phenomenology and from the microscopic many-body theory. In recent years, a great interest was devoted not only to the Nuclear Matter symmetry energy at saturation density but also to its whole density dependence, which is an essential ingredient for our understanding of many phenomena. We analyze the nuclear symmetry energy in different realms of nuclear physics and astrophysics. In particular we consider the nuclear symmetry energy in relation to nuclear structure, astrophysics of Neutron Stars and supernovae, and heavy ion collision experiments, trying to elucidate the connections of these different fields on the basis of the symmetry energy peculiarities. The interplay between experimental and observational data and theoretical developments is stressed. The expected future developments and improvements are schematically addressed, together with most demanded experimental and theoretical advances for the next few years.

We derive a new equation of state (EoS) for neutron stars (NS) from the outer crust to the core based on modern microscopic Brueckner-Hartree-Fock (BHF) calculations using the Argonne $v_{18}$ potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use the recent Barcelona-Catania-Paris-Madrid (BCPM) nuclear energy density functional that is directly based on the same microscopic BHF calculations, and which is able to reproduce the ground-state properties of nuclei along the periodic table. The EoS of the outer crust requires the masses of neutron-rich nuclei, which are obtained through Hartree-Fock-Bogoliubov calculations with the BCPM functional when they are unknown experimentally. To compute the inner crust, Thomas-Fermi calculations in Wigner-Seitz cells are performed with the same functional. Existence of nuclear pasta is predicted in a range of average baryon densities between $\simeq$0.067 fm$^{-3}$ and $\simeq$0.0825 fm$^{-3}$, where the transition to the core takes place. The NS core is computed from the nuclear EoS of the BHF calculation assuming non-exotic constituents (core of $npe\mu$ matter). In each region of the star, we discuss the comparison of the new EoS with previous EoSes for the complete NS structure, in particular, with the Lattimer-Swesty EoS and with the Shen et al. EoS widely used in astrophysical calculations. The new microscopically derived EoS fulfills at the same time a NS maximum mass of 2~$M_\odot$ with a radius of 10 km, and a 1.5~$M_\odot$ NS with a radius of 11.7 km.

We study cold dense quark matter and hybrid neutron stars with a Dyson-Schwinger quark model and various choices of the quark-gluon vertex. We obtain the equation of state of quark matter in beta equilibrium and investigate the hadron-quark phase transition in combination with a hadronic equation of state derived within the Brueckner-Hartree-Fock many-body theory. Comparing with the results for quark matter with the rainbow approximation, the Ball-Chiu ansatz and the 1BC ansatz for the quark-gluon vertex lead to a reduction of the effective interaction at finite chemical potential, qualitatively similar to the effect of our gluon propagator. We find that the phase transition and the equation of state of the quark or mixed phase and consequently the resulting hybrid star mass and radius depend mainly on a global reduction of the effective interaction due to effects of both the quark-gluon vertex and gluon propagator, but are not sensitive to the vertex ansatz.

Within a microscopic approach the structure of Neutron Stars is usually studied by modelling the homogeneous nuclear matter of the core by a suitable Equation of State, based on a many-body theory, and the crust by a functional based on a more phenomenological approach. We present the first calculation of Neutron Star overall structure by adopting for the core an Equation of State derived from the Brueckner-Hartree-Fock theory and for the crust, including the pasta phase, an Energy Density Functional based on the same Equation of State, and which is able to describe accurately the binding energy of nuclei throughout the mass table. Comparison with other approaches is discussed. The relevance of the crust Equation of state for the Neutron Star radius is particularly emphasised.

We calculate the effective masses of neutrons and protons in dense nuclear matter within the microscopic Brueckner-Hartree-Fock many-body theory and study the impact on the neutrino emissivity processes of neutron stars. We compare results based on different nucleon-nucleon potentials and nuclear three-body forces. Useful parametrizations of the numerical results are given. We find substantial in-medium suppression of the emissivities, strongly dependent on the interactions.

We calculate thermal conductivity and shear viscosity of nucleons in dense nuclear matter of neutron star cores in the non-relativistic Brueckner-Hartree-Fock framework. Nucleon-nucleon interaction is described by the Argonne v18 potential with addition of the Urbana IX effective three-body forces. We find that this three body force leads to decrease of the kinetic coefficients with respect to the two-body case. The results of calculations are compared with electron and muon transport coefficients as well as with the results of other authors.

The thermomagnetic evolution of the young neutron star in Cassiopea A is studied by considering fast neutrino emission processes. In particular, we consider neutron star models obtained from the equation of state computed in the framework of the Brueckner-Bethe-Goldstone many-body theory and variational methods, and models obtained with the Akmal-Pandharipande-Ravenhall equation of state. It is shown that it is possible to explain a fast cooling regime as the one observed in the neutron star in Cassiopea A if the Joule heating produced by dissipation of the small-scale magnetic field in the crust is taken into account. We thus argue that it is difficult to put severe constraints on the superfluid gap if the Joule heating is considered.

We compare a set of equations of state derived within microscopic many-body approaches, and study their predictions as far as phenomenological data on nuclei from heavy ion collisions, and astrophysical observations on neutron stars are concerned. We find that all the data, taken together, put strong constraints not easy to be fulfilled accurately. However, no major discrepancies are found among the selected equations of state and with respect to the data. The results provide an estimate of the uncertainty on the theoretical prediction at a microscopic level of the nuclear equation of state.

We analyze the nuclear matter correlation properties in terms of the pair correlation function. To this aim we systematically compare the results for the variational method in the Lowest Order Constrained Variational (LOCV) approximation and for the Bruekner-Hartree-Fock (BHF) scheme. A formal link between the Jastrow correlation factor of LOCV and the Defect Function (DF) of BHF is established and it is shown under which conditions and approximations the two approaches are equivalent. From the numerical comparison it turns out that the two correlation functions are quite close, which indicates in particular that the DF is approximately local and momentum independent. The Equations of State (EOS) of Nuclear Matter in the two approaches are also compared. It is found that once the three-body forces (TBF) are introduced the two EOS are fairly close, while the agreement between the correlation functions holds with or without TBF.

In this paper we study the dynamics of flavor transformation for neutrinos propagating in the very dense environment of astrophysical compact objects as Type II supernova in post collapse phase and proto-neutron stars. The analysis is based on the generalized Boltzmann equation, which incorporates the neutrinoneutrino interaction. We focus the analysis on the possible collective flavor dynamics, which can displays bipolar and synchronized flavor oscillations, and the associated transition from single-angle to multi-angle regime. To solve the kinetic equations we use an expansion in Legendre polynomials and a two-flavor scheme. On the basis of the numerical simulations we argue that neutrino flavor dynamics is suppressed by the presence of the electron component up to a certain distance from the neutrino sphere, where a transition to multi-angle regime occurs. This distance will move towards the neutrino sphere as the neutrino emission epoch proceeds.

In this paper we study in detail the dynamics of flavor transformation for neutrinos propagating in the very dense environment of astrophysical compact objects as Type II supernova in post collapse phase and proto-neutron stars. The analysis is based on the formalism by Strack and Burrows, who introduced the generalized Boltzmann equation for Wigner phase space density. In appropriate limits the formalism reduces to the usual evolution equations for the wave functions or for density matrix elements of Liouville equation. We incorporate the most important aspects of neutrino propagation physics: the phenomenology of standard oscillations with MSW resonance induced by ordinary matter, collective behavior due to self-interaction, which can produce bipolar and synchronized flavor oscillations, whose relevance is recognized by recent literature, the combination of these effects with collisions including scattering, emission and absorption of neutrinos.

We review our knowledge on the properties of the nuclear medium that have been studied, along many years, on the basis of many-body theory, laboratory experiments and astrophysical observations. First we consider the realm of phenomenological laboratory data and astrophysical observations, and the hints they can give on the characteristics that the nuclear medium should possess. The analysis is based on phenomenological models, that however have a strong basis on physical intuition and an impressive success. More microscopic models are also considered, and it is shown that they are able to give invaluable information on the nuclear medium, in particular on its Equation of State. The interplay between laboratory experiments and astrophysical observations are particularly stressed, and it is shown how their complementarity enriches enormously our insights into the structure of the nuclear medium. We then introduce the nucleon-nucleon interaction and the microscopic many-body theory of nuclear matter, with a critical discussion about the different approaches and their results. The Landau Fermi Liquid theory is introduced and briefly discussed. As illustrative example, we discuss neutron matter at very low density, and it is shown how it can be treated within the many-body theory. A section is dedicated to the pairing problem. The connection with nuclear structure is then discussed, on the basis of the Energy Density Functional method. The possibility to link the physics of exotic nuclei and the astrophysics of neutron stars is particularly stressed. Finally we discuss the thermal properties of the nuclear medium, in particular the liquid-gas phase transition and its connection with the phenomenology on heavy ion reactions and the cooling evolution of neutron stars. The presentation has been taken for non-specialists and possibly for non-nuclear physicists.

We study the hadron-quark phase transition in the interior of neutron stars. For the hadronic sector, we use a microscopic equation of state involving nucleons and hyperons derived within the Brueckner-Hartree-Fock many-body theory with realistic two-body and three-body forces. For the description of quark matter, we employ the Dyson-Schwinger approach and compare with the MIT bag model. We calculate the structure of neutron star interiors comprising both phases and find that with the Dyson-Schwinger model, the hadron-quark phase transition takes place only when hyperons are excluded, and that a two-solar-mass hybrid star is possible only if the nucleonic equation of state is stiff enough.

We study the possible collective plasma modes which can affect neutron-star thermodynamics and different elementary processes in the baryonic density range between nuclear saturation ($\rho_0$) and $3\rho_0$. In this region, the expected constituents of neutron-star matter are mainly neutrons, protons, electrons and muons ($npe\mu$ matter), under the constraint of beta equilibrium. The elementary plasma excitations of the $pe\mu$ three-fluid medium are studied in the RPA framework. We emphasize the relevance of the Coulomb interaction among the three species, in particular the interplay of the electron and muon screening in suppressing the possible proton plasma mode, which is converted into a sound-like mode. The Coulomb interaction alone is able to produce a variety of excitation branches and the full spectral function shows a rich structure at different energy. The genuine plasmon mode is pushed at high energy and it contains mainly an electron component with a substantial muon component, which increases with density. The plasmon is undamped for not too large momentum and is expected to be hardly affected by the nuclear interaction. All the other branches, which fall below the plasmon, are damped or over-damped.

We study the collective density modes which can affect neutron-star thermodynamics in the baryonic density range between nuclear saturation ($\rho_0$) and $3\rho_0$. In this region, the expected constituents of neutron-star matter are mainly neutrons, protons and electrons ($npe$ matter), under the constraint of beta equilibrium. The elementary excitations of this $npe$ medium are studied in the RPA framework. We emphasize the effect of Coulomb interaction, in particular the electron screening of the proton plasmon mode. For the treatment of the nuclear interaction, we compare two modern Skyrme forces and a microscopic approach. The importance of the nucleon effective mass is observed.

We study the structure of hybrid stars with the Field Correlator Method, extended to the zero temperature limit, for the quark phase. For the hadronic phase, we use the microscopic Brueckner-Hartree- Fock many-body theory. The comparison with the neutron star mass phenomenology puts serious constraints on the currently adopted values of the gluon condensate $G_2 \simeq 0.006-0.007 \rm {GeV^4}$, and the large distance static $Q \bar Q$ potential.

We explore the relevance of confinement in quark matter models for the possible quark core of neutron stars. For the quark phase, we adopt the equation of state (EoS) derived with the Field Correlator Method, extended to the zero temperature limit. For the hadronic phase, we use the microscopic Brueckner-Hartree-Fock many-body theory. We find that the currently adopted value of the gluon condensate $G_2 \simeq 0.006-0.007 \rm {GeV^4}$, which gives a critical temperature $T_c \simeq 170 \rm MeV$, produces maximum masses which are only marginally consistent with the observational limit, while larger masses are possible if the gluon condensate is increased.

Supernovae of Type II is a phenomenon that occurs at the end of evolution of massive stars when the iron core of the star exceeds a mass limit. After collapse of the core under gravity the shock wave alone does not succeed in expelling the mass of the star and in this sense the role of neutrinos is the most important mechanism to do so. During the emission of neutrinos flavor conversion is possible, related the phenomenon of oscillations, which however depends directly on the particular density profile of the medium. In this paper we present results of numerical simulations of neutrino flavor conversion in protoneutron stars and after collapse. The probabilities of survival for a given flavor in a complete three-flavors framework is estimated through an algorithm which conserves unitarity to a high degree of accuracy. The sensitivity of the results to the different adopted models for the protoneutron star structure is examined in detail demonstrating how the neutrino signal could be used to check the validity of models.

Recently, it has been shown that the standard Nambu-Jona-Lasinio (NJL) model is not able to reproduce the correct QCD behavior of the gap equation at large density, and therefore a different cutoff procedure at large momenta has ben proposed. We found that, even with this density dependent cutoff procedure, the pure quark phase in neutron stars (NS) interiors is unstable, and we argue that this could be related to the lack of confinement in the original NJL model.

The quark matter equation of state (EOS) derived from the standard Nambu - Jona-Lasinio (NJL) model is soft enough to render neutron stars (NS) unstable at the onset of the deconfined phase, and no pure quark matter can be actually present in its interior. Since this is a peculiarity of the NJL model, we have studied a modified NJL model with a momentum cut-off which depends on the density. This procedure, which improves the agreement between QCD and NJL model at large density, modifies the standard NJL equation of state, and then it is potentially relevant for the stability analysis of neutron stars. We show that also within this approach, the NS instability still persists, and that the vacuum pressure, as a signal of quark confinement, has a fundamental role for the NS stability. In this respect, our conclusions point to a relationship between confinement and NS stability.

We study the structure of protoneutron stars within the finite temperature Brueckner-Bethe-Goldstone many-body theory. If nucleons, hyperons, and leptons are present in the stellar core, we find that neutrino trapping stiffens considerably the equation of state, because hyperon onsets are shifted to larger baryon density. However, the value of the critical mass turns out to be smaller than the ``canonical'' value 1.44 $M_\odot$. We find that the inclusion of a hadron-quark phase transition increases the critical mass and stabilizes it at about 1.5--1.6 $M_\odot$.

We study the hadron-quark phase transition in the interior of protoneutron stars. For the hadronic sector, we use a microscopic equation of state involving nucleons and hyperons derived within the finite-temperature Brueckner-Bethe-Goldstone many-body theory, with realistic two-body and three-body forces. For the description of quark matter, we employ the MIT bag model both with a constant and a density-dependent bag parameter. We calculate the structure of protostars with the equation of state comprising both phases and find maximum masses below 1.6 solar masses. Metastable heavy hybrid protostars are not found.

We study the structure of newly born neutron stars (protoneutron stars) within the finite temperature Brueckner-Bethe-Goldstone theoretical approach including also hyperons. We find that for purely nucleonic stars both finite temperature and neutrino trapping reduce the value of the maximum mass. For hyperonic stars the effect is reversed, because neutrino trapping shifts the appearance of hyperons to larger baryon density and stiffens considerably the equation of state.

We study the hadron-quark phase transition in the interior of neutron stars (NS). For the hadronic sector, we use a microscopic equation of state (EOS) involving nucleons and hyperons derived within the Brueckner-Bethe-Goldstone many-body theory, with realistic two-body and three-body forces. For the description of quark matter, we employ both the MIT bag model with a density dependent bag constant, and the color dielectric model. We calculate the structure of NS interiors with the EOS comprising both phases, and we find that the NS maximum masses are never larger than 1.7 solar masses, no matter the model chosen for describing the pure quark phase.

We discuss the high-density nuclear equation of state within the Brueckner-Hartree-Fock approach. Particular attention is paid to the effects of nucleonic three-body forces, the presence of hyperons, and the joining with an eventual quark matter phase. The resulting properties of neutron stars, in particular the mass-radius relation, are determined. It turns out that stars heavier than 1.3 solar masses contain necessarily quark matter.

We explore the relevance of color superconductivity inside a possible quark matter core for the bulk properties of neutron stars. For the quark phase we use an Nambu--Jona-Lasinio (NJL) type model, extended to include diquark condensates. For the hadronic phase, a microscopic many-body model is adopted, with and without strangeness content. In our calculations, a sharp boundary is assumed between the hadronic and the quark phases. For NJL model parameters fitted to vacuum properties we find that no star with a pure quark core does exist. Nevertheless the presence of color superconducting phases can lower the neutron star maximum mass substantially. In some cases, the transition to quark matter occurs only if color superconductivity is present. Once the quark phase is introduced, the value of the maximum mass stays in any case below the value of two solar masses.

We study the hadron-quark phase transition in the interior of neutron stars (NS's). We calculate the equation of state (EOS) of hadronic matter using the Brueckner-Bethe-Goldstone formalism with realistic two-body and three-body forces, as well as a relativistic mean field model. For quark matter we employ the MIT bag model constraining the bag constant by using the indications coming from the recent experimental results obtained at the CERN SPS on the formation of a quark-gluon plasma. We find necessary to introduce a density dependent bag parameter, and the corresponding consistent thermodynamical formalism. We calculate the structure of NS interiors with the EOS comprising both phases, and we find that the NS maximum masses fall in a relatively narrow interval, $1.4 M_\odot \leq M_{\rm max} \leq 1.7 M_\odot$. The precise value of the maximum mass turns out to be only weakly correlated with the value of the energy density at the assumed transition point in nearly symmetric nuclear matter.

Massive neutron stars (NS) are expected to possess a quark core. While the hadronic side of the NS equation of state (EOS) can be considered well established, the quark side is quite uncertain. While calculating the EOS of hadronic matter we have used the Brueckner-Bethe-Goldstone formalism with realistic two-body and three-body forces, as well as a relativistic mean field model. For quark matter we employ the MIT bag model constraining the bag constant by exploiting the recent experimental results obtained at CERN on the formation of a quark-gluon plasma. We calculate the structure of NS interiors with the EOS comprising both phases, and we find that the NS maximum masses fall in a relatively narrow interval, $1.45 M_\odot \leq M_{\rm max} \leq 1.65 M_\odot$, near the lower limit of the observational range.

We make a first calculation of eigenfrequencies of radial pulsations of neutron stars with quark cores in a general relativistic formalism given by Chandrasekhar. The equations of state (EOS) used to estimate such eigenfrequencies have been derived by taking proper care of the hadron-quark phase transition. The hadronic EOS's have been obtained in the framework of the Brueckner-Hartree-Fock and relativistic mean field theories, whereas the quark EOS has been derived within the MIT bag model. We find that the periods of oscillations of neutron stars with a quark core show a kink, which is associated with the presence of a mixed phase region. Also, oscillation periods show significant differences between ordinary neutron stars and neutron stars with dynamical quark phases.

The Bethe-Brueckner-Goldstone many-body theory of the Nuclear Equation of State is reviewed in some details. In the theory, one performs an expansion in terms of the Brueckner two-body scattering matrix and an ordering of the corresponding many-body diagrams according to the number of their hole-lines. Recent results are reported, both for symmetric and for pure neutron matter, based on realistic two-nucleon interactions. It is shown that there is strong evidence of convergence in the expansion. Once three-body forces are introduced, the phenomenological saturation point is reproduced and the theory is applied to the study of neutron star properties. One finds that in the interior of neutron stars the onset of hyperons strongly softens the Nuclear Equation of State. As a consequence, the maximum mass of neutron stars turns out to be at the lower limit of the present phenomenological observation.

In the framework of the non-relativistic Brueckner-Bethe-Goldstone theory, we derive a microscopic equation of state for asymmetric and $\beta$-stable matter containing $\Sigma^-$ and $\Lambda$ hyperons. We mainly study the effects of three-body forces (TBFs) among nucleons on the hyperon formation and the equation of state (EoS). We find that, when TBFs are included, the stellar core is almost equally populated by nucleons and hyperons. The resulting EoS, which turns out to be extremely soft, has been used in order to calculate the static structure of neutron stars. We obtain a value of the maximum mass of 1.26 solar masses (1 solar mass $M_o \simeq 1.99 \cdot 10^{33} g$). Stellar rotations increase this value by about 12%.

In the framework of the Brueckner-Bethe-Goldstone theory, we determine a fully microscopic equation of state for asymmetric and $\beta$-stable nuclear matter containing $\sim$ and $\la$ hyperons. We use the Paris and the new Argonne $Av_{18}$ two-body nucleon interaction, whereas the nucleon-hyperon interaction is described by the Njimegen soft-core model. We stress the role played by the three-body nucleon interaction, which produces a strong repulsion at high densities. This enhances enormously the hyperon population, and produces a strong softening of the equation of state, which turns out almost independent on the nucleon-nucleon interaction. We use the new equation of state in order to calculate the structure of static neutron stars. We obtain a maximum mass configuration with $M_{\rm max}$ = 1.26 (1.22) when the Paris ($Av_{18}$) nucleon potential is adopted. Central densities are about 10 times normal nuclear matter density. Stellar rotations, treated within a perturbative approach, increase the value of the limiting mass by about 12%.

We study the properties of static and rotating neutron stars adopting non-relativistic equations of state (EOS) for asymmetric nuclear matter based on the Brueckner-Hartree-Fock (BHF) scheme. The BHF calculation, with the continuous choice for the single particle potential, appears to be very close to the full EOS, which includes the three-hole line contribution calculated by solving the Bethe-Fadeev equations within the gap choice for the single particle potential. Three-body forces are included in order to reproduce the correct saturation point for nuclear matter. A comparison with fully relativistic many-body calculations of nuclear matter EOS is made.

We calculate static properties of non-rotating neutron stars (NS's) using a microscopic equation of state (EOS) for asymmetric nuclear matter, derived from the Brueckner-Bethe-Goldstone many-body theory with explicit three-body forces. We use the Argonne AV14 and the Paris two-body nuclear force, implemented by the Urbana model for the three-body force. We obtain a maximum mass configuration with $ M_{max} = 1.8 M_{\sun}$ ($M_{max} = 1.94 M_{\sun}$) when the AV14 (Paris) interaction is used. They are both consistent with the observed range of NS masses. The onset of direct Urca processes occurs at densities $n \geq 0.65~fm^{-3}$ for the AV14 potential and $n \geq 0.54~fm^{-3}$ for the Paris potential. Therefore, NS's with masses above $M^{Urca} = 1.4 M_{\sun}$ for the AV14 and $M^{Urca} = 1.24 M_{\sun}$ for the Paris potential can undergo very rapid cooling, depending on the strength of superfluidity in the interior of the NS. The comparison with other microscopic models for the EOS shows noticeable differences.

We calculate static properties of non-rotating neutron stars (NS's) using a microscopic equation of state (EOS) for asymmetric nuclear matter. The EOS is computed in the framework of the Brueckner--Bethe--Goldstone many--body theory. We introduce three-body forces in order to reproduce the correct saturation point of nuclear matter. A microscopic well behaved EOS is derived. We obtain a maximum mass configuration with $M_{max} = 1.8 M_\odot$, a radius $R = 9.7$ km and a central density $n_c = 1.34~fm^{-3}$. We find the proton fraction exceeds the critical value $x^{Urca}$, for the onset of direct Urca processes, at densities $n \geq 0.45~fm^{-3}$. Therefore, in our model, NS's with masses above $M^{Urca} = 0.96 M_\odot$ can undergo very rapid cooling depending on whether or not nucleon superfluidity in the interior of the NS takes place. A comparison with other microscopic models for the EOS is done, and neutron star structure is calculated for these models too.