Strange quark matter in beta equilibrium at high densities is studied in a quark confinement model. Two equations of state are dynamically generated for the {\it same} set of model parameters used to describe the nucleon: one corresponds to a chiral restored phase with almost massless quarks and the other to a chiral broken phase. The chiral symmetric phase saturates at around five times the nuclear matter density. Using the equation of state for this phase, compact bare quark stars are obtained with radii and masses in the ranges $R\sim 5 - 8$ km and $M\sim M_\odot$. The energy per baryon number decreases very slowly from the center of the star to the periphery, remaining above the corresponding values for the iron or the nuclear matter, even at the edge. Our results point out that strange quark matter at very high densities may not be absolutely stable and the existence of an energy barrier between the two phases may prevent the compact quarks stars to decay to hybrid stars.

Equations of state for strange quark matter in beta equilibrium at high densities are used to investigate the structure (mass and radius) of compact objects. The chromodielectric model is used as a general framework for the quark interactions, which are mediated by chiral mesons, $\sigma$ and $\vec \pi$, and by a confining chiral singlet dynamical field, $\chi$. Using a quartic potential for $\chi$, two equations of state for the same set of model parameters are obtained, one with a minimum at around the nuclear matter density $\rho_0$ and the other at $\rho \sim 5 \rho_0$. Using the latter equation of state in the Tolman-Oppenheimer-Volkoff equations we found solutions corresponding to compact objects with $R\sim 5 - 8$ km and $M\sim M_\odot$. The phenomenology of recently discovered X-ray sources is compatible with the type of quark stars that we have obtained.

We introduce in this work a class of relativistic models for nuclear matter and neutron stars which exhibits a parameterization, through mathematical constants, of the non-linear meson-baryon couplings. For appropriate choices of the parameters, it recovers current QHD models found in the literature: Walecka and the Zimanyi and Moszkowski models (ZM and ZM3). For other choices of the parameters, the new models give very interesting and new physical results.We have obtained numerical values for the maximum mass and radius of a neutron star sequence, red-shift, hyperon populations, radial distribution of particles, among other relevant static properties of these stellar objects. The phenomenology of neutron stars in ZM models is presented and compared to the phenomenology obtained in other versions of the Walecka model. We have found that the ZM3 model is too soft and predicts a very small maximum neutron star mass, $\sim0.72 M_\odot$. A strong similarity between the results of ZM-like models and those with exponential couplings is noted. Sensibility of the results to the specific choice of the values for the binding energy and saturation density is pointed out. Finally, we discuss the very intense scalar condensates found in the interior of neutron stars which may lead to negative effective masses.