In this study, our main focus is to investigate the mass–radius relation and several important properties of massive neutron stars to realize the nature, behavior and evolution of these kinds of compact objects at present time. Also, we want to understand the equation-of-state of the core nuclear matter precisely with their stable equilibrium configuration. We have chosen a few massive binary pulsars and investigated on maximum attainable mass and lowest possible radius of them. We have considered Einstein–Hilbert action as [Formula: see text], where [Formula: see text] is the Ricci scalar and [Formula: see text], the trace of energy–momentum tensor with [Formula: see text] as the coupling parameter and also used modified TOV equations. The interior space-time of the spherical neutron star is matched to the exterior Schwarzschild line element at the surface of the star. The [Formula: see text]–[Formula: see text] curve can predict the maximum achievable mass is about [Formula: see text] with lowest possible radius of around [Formula: see text] km for the massive compact stars under stable equilibrium. We have obtained a clear picture of structural evolution of massive neutron stars through accelerating space-time and can put some constraints on several quantities related to them. It is depicted from our present investigation that all the derived outcomes are compatible with physically adopted regimes which reveal the viability of our current model in the context of [Formula: see text] modified gravity.