We propose a density-dependent function for the attractive interaction in the original van der Waals model to correctly describe the flow constraint at the high-density regime of the symmetric nuclear matter. After a generalization to asymmetric nuclear matter, it was also possible to study the stellar matter regime from this new model. The mass-radius relation for neutron stars under $\beta$-equilibrium is found to agree with recent X-ray observations. The neutron star masses supported against gravity, obtained from some parametrizations of the model, are in the range of $(1.97-2.07)M_{\odot}$, compatible with observational data from the PSR J0348+0432 pulsar. Furthermore, we verify the reliability of the model in predicting tidal deformabilities of the binary system related to the GW170817 neutron star merger event and find a full agreement with the new bounds obtained by the LIGO/Virgo collaboration.

In the present work we used five different versions of the quark-meson coupling (QMC) model to compute astrophysical quantities related to the GW170817 event and to neutron star cooling process. Two of the models are based on the original bag potential structure and three versions consider a harmonic oscillator potential to confine the quarks. The bag-like models also incorporate the pasta phase used to describe the inner crust of neutron stars. We show that the pasta phase always play a minor or negligible role in all studies. Moreover, while no clear correlation between the models that satisfy the GW170817 constraints and the slope of the symmetry energy is found, a clear correlation is observed between the slope and the fact that the cooling is fast or slow, i.e., fast (slow) cooling is related to higher (lower) values of the slope. We did not find one unique model that can describe, at the same time, GW170817 constraints and give a perfect description of the possible cooling processes.

Using the nuclear equation of states for a large variety of relativistic and non-relativistic force parameters, we calculate the static and rotating masses and radii of neutron stars. From these equation of states, we also evaluate the properties of rotating neutron stars, such as rotational and gravitational frequencies, moment of inertia, quadrupole deformation parameter, rotational ellipcity and gravitational wave strain amplitude. The estimated gravitational wave strain amplitude of the star is found to be $\sim 10^{-23}$.