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 this work, we have investigated anisotropic neutron stars in the framework of Rastall-Rainbow gravity. All our calculations were computed using the IU-FSU realistic equation of state (EoS), in which was considered two cases: standard nucleonic composition and the one with the eight lightest baryons. From the neutron star masses and radii obtained we conclude that anisotropic pressure has significant consequences on the structure of stellar objects. In particular, when anisotropy is considered within the general relativity framework, it significantly modifies the maximum stellar mass. On the other hand, when Rastall-Rainbow gravity and anisotropy are simultaneously considered, they provide the best results for mass and radius values, including important astrophysical objects such as the LMXB NGC 6397 and the extremely massive pulsar millisecond MSP J0740 + 6620. Although the expected inclusion of hyperons in the nuclear model reproduces stellar masses smaller than those produced by standard nucleonic matter, we shown that the hyperon puzzle problem can be solved by including anisotropic effects on compact stars in the context of the Rastall-Rainbow gravity.

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.