We propose the disappearance of "the hyperon puzzle" in neutron star (NS) by invoking two new-physics prescriptions: modified gravity theory and braneworld scenario. By assuming that NS lives on a $3$-brane within a $5d$ empty AdS bulk, gravitationally governed by Eddington-inspired Born-Infeld (EiBI) theory, the field equations can be effectively cast into the usual Einstein's with "apparent" anisotropic energy-momentum tensor. Solving the corresponding brane-TOV equations numerically, we study its mass-radius relation. It is known that the appearance of finite brane tension $\lambda$ reduces the compactness of the star. The compatibility of the braneworld results with observational constraints of NS mass and radius can be restored in our model by varying the EiBI's coupling constant, $\kappa$. We found that within the astrophysically-accepted range of parameters ($0<\kappa<6\times10^6\text{m}^2$ and $\lambda\gg1~ \text{MeV}^4$) the NS can have mass $\sim2.1~ \text{M}_\odot$ and radius $\sim10$ km.

We investigate the effect of the inclusion of relativistic Coulomb terms in a confined-isospin-density-dependent-mass (CIDDM) model of strange quark matter (SQM). We found that if we include Coulomb term in scalar density form, SQM equation of state (EOS) at high densities is stiffer but if we include Coulomb term in vector density form is softer than that of standard CIDDM model. We also investigate systematically the role of each term of the extended CIDDM model. Compared with what was reported in Ref.~\cite {ref:isospin}, we found the stiffness of SQM EOS is controlled by the interplay among the the oscillator harmonic, isospin asymmetry and Coulomb contributions depending on the parameter's range of these terms. We have found that the absolute stable condition of SQM and the mass of 2 $M_\odot$ pulsars can constrain the parameter of oscillator harmonic $\kappa_1$ $\approx 0.53$ in the case Coulomb term excluded. If the Coulomb term is included, for the models with their parameters are consistent with SQM absolute stability condition, the $2.0~M_{\odot}$ constraint more prefer the maximum mass prediction of model with scalar Coulomb term than that of model with vector Coulomb term. On contrary, the high densities EOS predicted by model with vector Coulomb is more compatible with recent pQCD result \cite{ref:pressure} than that predicted by model with scalar Coulomb. Furthermore, we also observed the quark composition in a very high density region depends quite sensitively on the kind of Coulomb term used.

We investigate the mass-radius relation of neutron star (NS) with hyperons inside its core by using the Eddington-inspired Born-Infeld (EiBI) theory of gravity. The equation of state of the star is calculated by using the relativistic mean field model under which the standard SU(6) prescription and hyperons potential depths are used to determine the hyperon coupling constants. We found that, for $4\times 10^{6}~\rm{m^2}~\lesssim~\kappa \lesssim~6\times 10^{6}~\rm{m^2}$, the corresponding NS mass and radius predicted by the EiBI theory of gravity is compatible with observational constraints of maximum NS mass and radius. The corresponding $\kappa$ value is also compatible with the $\kappa$ range predicted by the astrophysical-cosmological constraints. We also found that the parameter $\kappa$ could control the size and the compactness of a neutron star.