The role of scalar boson exchange as a mediator of the fermionic dark particle interaction and the mass of dark particle on the bulk properties of fermionic dark stars including their moment of inertia and tidal deformability are studied. We have found that the role of the attractive nature of the scalar boson exchange and the fermionic dark particle mass can control the stiffness of the fermionic dark star equation of state. By increasing the strength of scalar boson coupling, the fermionic dark star becomes more compact. As a consequence, if scalar boson exchange contribution is included the compactness of a dark star can exceed $C$=0.22. We also compare the fermionic dark stars moment of inertia and tidal deformability to those of neutron stars (with and without hyperons in neutron star core) predicted by relativistic mean field model. It is evidence that the properties of both types of stars are quite different. We also have found that the universal I-Love relation in fermionic dark stars is not affected by scalar boson exchange contribution and the fermionic dark particle mass. Possible observations of fermionic dark stars are also discussed.

We study the stability of non-relativistic polytropic stars within two modified gravity theories i.e., Beyond Horndeski gravity and Eddington inspired Born Infeld theories using the configuration entropy method. We use spatially localized bounded function of energy density as solutions from stellar effective equations to construct the corresponding configuration entropy. We use the same argument as the one used by the authors of Refs. [1,2] that the stars are stable if there is peak in configuration entropy as a function of adiabatic index curve. Specifically, the boundary between stable and unstable regions which corresponds to Chandrasekhar stability bound is indicated from the existence of the maximum peak while the most stable polytropic stars is indicated by minimum peak in the corresponding curve. We have found that the value of critical adiabatic indexes of Chandrasekhar stability bound and the most stable polytropic stars predicted by non-relativistic limit of Beyond Horndeski Gravity and Eddington inspired Born Infeld theories are different to those predicted by general relativity where the corresponding differences depend on the free parameter of both theories.

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 examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic calculations, all describing 2$M_\odot$ neutron stars. Unified EoSs for the inner-crust-core region have been built for all the phenomenological models, both relativistic and non-relativistic. Our investigation shows the existence of a strong correlation of the neutron star radii with the linear combination of the slopes of the nuclear matter incompressibility and the symmetry energy coefficients at the saturation density. Such correlations are found to be almost independent of the neutron star mass in the range $0.6\text{-}1.8M_{\odot}$. This correlation can be linked to the empirical relation existing between the star radius and the pressure at a nucleonic density between one and two times saturation density, and the dependence of the pressure on the nuclear matter incompressibility, its slope and the symmetry energy slope. The slopes of the nuclear matter incompressibility and the symmetry energy coefficients as estimated from the finite nuclei data yield the radius of a $1.4M_{\odot}$ neutron star in the range $11.09\text{-}12.86$ 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.

We study the diversities in the properties of the neutron stars arising due to the different choices for the cross-coupling between various mesons which governs the density dependence of the nuclear symmetry energy in the extended relativistic mean-field(RMF) model. For this purpose, we obtain two different families of the extended RMF model corresponding to different non-linear cross-coupling term in the isovector part of the effective Lagrangian density. The lowest order contributions for the $\delta$ mesons are also included. The different models within the same family are so obtained that they yield wide variation in the value of neutron-skin thickness in the $^{208}$Pb nucleus. These models are employed to compute the neutron star properties such as, core-crust transition density, radius and red shift at canonical mass ($1.4M_{\odot}$), tidal polarizability parameter, and threshold mass required for the enhanced cooling through direct Urca process. Most of the neutron star properties considered are significantly different(10\%-40\%) for the different families of models at a smaller neutron-skin thickness ($\sim 0.15$ fm) in the $^{208}$Pb nucleus.

By using a simple thermodynamical method we confirm the finding of Chavanis and Harko that stable Bose-Einstein condensate stars can form. However, by using a thermodynamically consistent boson equation of state, we obtain a less massive Bose-Einstein condensate star compared to the one predicted by Chavanis and Harko. We also obtain that the maximum mass of a boson star is insensitive to the change of matter temperature. However, the mass of boson star with relatively large radius depends significantly on the temperature of the boson matter.

Empirically determined values of the nuclear volume and surface symmetry energy coefficients from nuclear masses are expressed in terms of density distributions of nucleons in heavy nuclei in the local density approximation. This is then used to extract the value of the symmetry energy slope parameter $L$. The density distributions in both spherical and well deformed nuclei calculated within microscopic framework with different energy density functionals give $L = 59.0 \pm 13.0$ MeV. Application of the method also helps in a precision determination of the neutron skin thickness of nuclei that are difficult to measure accurately.

We have calculated the proton charge radius by assuming that the real proton radius is not unique and the radii are randomly distributed in a certain range. This is performed by averaging the elastic electron-proton differential cross section over the form factor cut-off. By using a dipole form factor and fitting the middle value of the cut-off to the low $Q^2$ Mainz data, we found the lowest $\chi^2/N$ for a cut-off $\Lambda=0.8203\pm 0.0003$ GeV, which corresponds to a proton charge radius $r_E=0.8333\pm 0.0004$ fm. The result is compatible with the recent precision measurement of the Lamb shift in muonic hydrogen as well as recent calculations using more sophisticated techniques. Our result indicates that the relative variation of the form factor cut-off should be around 21.5%. Based on this result we have investigated effects of the nucleon radius variation on the symmetric nuclear matter (SNM) and the neutron star matter (NSM) by considering the excluded volume effect in our calculation. The mass-radius relation of neutron star is found to be sensitive to this variation. The nucleon effective mass in the SNM as well as the equation of state of both the SNM and the NSM exhibit a similar sensitivity.

The possibility of existence of hyperons in the recently measured $2M_\odot$ pulsar PSRJ1614-2230 is explored using a diverse set of nuclear equations of state calculated within the relativistic mean-field models. Our results indicate that the nuclear equations of state compatible with heavy-ion data allow the hyperons to exist in the PSRJ1614-2230 only for significantly larger values for the meson-hyperon coupling strengths. The maximum mass configurations for these cases contain sizable hyperon fractions ($\sim 60%$) and yet masquared their counterpart composed of only nucleonic matter.

Based on the properties of the critical and the actual effective masses of sigma and omega mesons, criteria to estimate the values of the isoscalar nonlinear terms of the standard relativistic mean field model that reproduce stable equations of state in respect to particle hole excitation at high densities are derived. The relation between nuclear matter stability and the symmetric nuclear matter properties are shown. The criteria are used to analyze in a more systematic way the high-density longitudinal and transverse instabilities of some parameter sets of relativistic mean field models. The critical role of the vector and vector-scalar nonlinear terms is also discussed quantitatively.

The effect of neutrino trapping on the longitudinal dielectric function at low densities has been investigated by using different relativistic mean field models. Parameter sets G2 of Furnstahl-Serot-Tang and Z271 of Horowitz-Piekarewicz, along with the adjusted parameter sets of both models, have been used in this study. The role of the isovector adjustment and the effect of the Coulomb interaction have been also studied. The effect of the isovector adjustment is found to be more significant in the Horowitz-Piekarewicz model, not only in the neutrinoless matter, but also in the matter with neutrino trapping. Although almost independent to the variation of the leptonic fraction, the instability region of matter with neutrino trapping is found to be larger. The presence of more protons and electrons compared to the neutrinoless case is the reason behind this finding. For parameter sets with soft equation of states at low density, the appearance of a large and negative epsilon_L (q, q0 = 0) in some parts of the edge of the instability region in matter with neutrino trapping is understood as a consequence of the fact that the Coulomb interaction produced by electrons and protons interaction is larger than the repulsive isovector interaction created by the asymmetry between the proton and neutron numbers.

The instability of nuclear matter due to particle-hole excitation modes has been studied in the frame-work of several relativistic mean field (RMF) models. It is found that both the longitudinal and the transversal modes depend sensitively on the parameter sets used. The important impact of the vector and vector-scalar nonlinear terms on the stability of both modes is demonstrated. Our finding corroborates the result of previous studies, namely that certain RMF models cannot be used in high density applications. However, we show that for certain parameter sets of RMF models this shortcoming can be alleviated by adding these nonlinear terms.

The effects of the symmetry energy softening of the relativistic mean field (RMF) models on the properties of matter with neutrino trapping are investigated. It is found that the effects are less significant than those in the case without neutrino trapping. The weak dependence of the equation of state on the symmetry energy is shown as the main reason of this finding. Using different RMF models the dynamical instabilities of uniform matters, with and without neutrino trapping, have been also studied. The interplay between the dominant contribution of the variation of matter composition and the role of effective masses of mesons and nucleons leads to higher critical densities for matter with neutrino trapping. Furthermore, the predicted critical density is insensitive to the number of trapped neutrinos as well as to the RMF model used in the investigation. It is also found that additional nonlinear terms in the Horowitz-Piekarewicz and Furnstahl-Serot-Tang models prevent another kind of instability, which occurs at relatively high densities. The reason is that the effective sigma meson mass in their models increases as a function of the matter density.

We have studied the effects of the nucleon weak magnetism and the neutrino electromagnetic properties on the neutrino and antineutrino mean free paths for all types of neutrinos. In the calculation, we consider matters with and without neutrino trapping as the target. We have found that the difference between the two mean free paths depends not only on the neutrino energy, but also on the matter density as well as matter constituents. Compared to the nucleon weak magnetism, the neutrino electromagnetic properties are found to have negligible effects on this difference.

The sensitivity of the differential cross section of the interaction between neutrino-electron with dense matter to the possibly nonzero neutrino electromagnetic properties has been investigated. Here, the relativistic mean field model inspired by effective field theory has been used to describe non strange dense matter, both with and without the neutrino trapping. We have found that the cross section becomes more sensitive to the constituent distribution of the matter, once electromagnetic properties of the neutrino are taken into account. The effects of electromagnetic properties of neutrino on the cross section become more significant for the neutrino magnetic moment mu_nu > 10^{-10} mu_B and for the neutrino charge radius R > 10^{-5} MeV^{-1}.

The effects of vector-isovector terms predicted by relativistic mean field models based on effective field theory on the equation of state and the neutrino mean free path in neutron stars have been studied. Using a procedure similar to Ref. [1] in treating the isovector sector, a parameter set (G2*) that predicts a much smaller proton fraction than the standard one (G2) can be obtained. We found that the G2* has a softer equation of state compared to the G2 at higher densities and the disappearance of the anomalous behavior of the neutrino mean free path in neutron star does not depend on the proton fraction.

An improvement in the treatment of the isovector channel of relativistic mean field (RMF) models based on effective field theory (E-RMF) is suggested, by adding an isovector scalar (delta) meson and using a similar procedure to the one used by Horowitz and Piekarewicz to adjust the isovector-vector channel in order to achieve a softer density dependent symmetry energy of the nuclear matter at high density. Their effects on the equation of state (EOS) at high density and on the neutrino mean free path (NMFP) in neutron stars are discussed.

The mean free path of neutrino - free electron gas interaction has been calculated by taking into account the neutrino electromagnetic form factors and the possibility of neutrino oscillation. It is shown that the form factor effect becomes significant for a neutrino magnetic moment \mu_\nu > 10^{-10} mu_B and for a neutrino radius R > 10^{-6} MeV^{-1}. The mean free path is found to be sensitive to the nu_e-nu_mu and nu_e-nu_e^R transition probabilities.

The equation of state (EOS) of dense matter and neutrino mean free path (NMFP) in a neutron star have been studied by using relativistic mean field models motivated by effective field theory (ERMF). It is found that the models predict too large proton fractions, although one of the models (G2) predicts an acceptable EOS. This is caused by the isovector terms. Except G2, the other two models predict anomalous NMFP. In order to minimize the anomaly, besides an acceptable EOS, a large M* is favorable. A model with large M* retains the regularity in the NMFP even for a small neutron fraction.