We present a generalization of Rastall's gravity in which the conservation law of the energy-moment tensor is altered, and as a result, the trace of the energy-moment tensor is taken into account together with the Ricci scalar in the expression for the covariant derivative. Afterwards, we obtain the field equation in this theory and solve it by considering a spherically symmetric space-time. We show that the external solution has two possible classes of solutions with spherical symmetry in the vacuum in generalized Rastall's gravity. The first class of solutions is completely equivalent to the Schwarzschild solution, while the second class of solutions has the same structure as the Schwarzschild--de Sitter solution in general relativity. The generalization, in contrast to constant value $k=8\pi G$ in general relativity, has a gravitational parameter $k$ that depends on the energy density $\rho$. As an application, we perform a careful analysis of the effects of the theory on neutron stars using realistic equations of state (EoS) as inputs. Our results show that important differences on the profile of neutron stars are obtained within two representatives EoS.

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.

The possibility of modifications on general relativity is investigated. We propose an alternative theory of gravity constructed with the combination of Rastall and Rainbow theories. The hydrostatic equilibrium equations are obtained in order to test the new theory in neutron stars, whose mass-radius diagrams are obtained using modern equations of state of nuclear matter derived from relativistic mean field models and compared with the ones computed by the Tolman-Oppenheimer-Volkoff equations. We conclude that substantial modifications are obtained even for very small alterations on the two free parameters, making the reproduction of astrophysical observations an easy task.