Neutron stars in generalized f(R) gravity

A gravity theory is considered with the Einstein-Hilbert Lagrangean R+aR ²+bR _μν R ^μν , R _μν being Ricci’s tensor and R the curvature scalar. The parameters a and b are taken of order 1 km². Arguments are given which suggest that the effective theory so obtained might be a fair approximation of a viable theory. A numerical integration is performed of the field equations for a free neutron gas. The result is that the star mass increases with increasing central density until about 1 solar mass and then decreases. The baryon number increases monotonically, which suggests that the theory allows stars in equilibrium with arbitrary baryon number, no matter how large.