In the Tolman-Oppenheimer-Volkoff model of neutron star, the spherically symmetric metrics was written with the help of auxiliary quantity "u". Specifically, "u" was established with the help of g_{rr} component of metric tensor as u = r*(1 + 1/g_{rr})/2 ("r" is the radial distance). Or, the inverse calculation is g_{rr} = -1/(1 - 2*u/r).
In the outer Schwarschild solution for the vacuum, quantity "u" is constant and also component g_{tt} of metric tensor can be given with the help of "u" as g_{tt} = 1 - 2*u/r. This is then used to gauge "u" making the deduction that the gravitational acceleration in general relativity (GR) is given as a = -(c^{2}/2)*(@g_{tt}/@r) = -c^{2}*u/r^{2} in the limit of weak field and this acceleration should be identical to its Newtonian counterpart equal to a = grad(Psi) = @(Psi)/@r = -G*M/r^{2} ("c" is speed of light, symbol "@" is used for the partial derivative, "Psi" is the Newtonian gravitational potential, "G" is the gravitational onstant, and "M" is the mass). Comparing both relations for "a", the numerical value of quantity "u" can be accounted as u = G*M/c^{2}. Or, u = r*Psi/c^{2}.
The GR was initially created as the geometric theory, including the concept of energy and potential, but not mass. However, quantity "u" has often been represented as closely related to mass "M" in the theory of neutron stars. Or, it has been directly identified to mass. (If one chooses such the system of units that G=1 and c=1, then relation u = G*M/c^{2} reduces to the apparent identity u = M.) So, the identification of "u" and "M" has been considered not only in a weak field, but it has been generalized and used in an argumentation everywhere inside the neutron star.
PROPER QUESTION: Is the representation of auxiliary metric quantity "u" as the "mass", or quantity directly related to mass (according to the relation u = G*M/c^{2}) correct? Consequently, must there always be valid u > 0 (since M > 0)? And, why we then demand, in fact, the validity of inequality |g_{rr}| > 1, which is implied by u > 0? What are the arguments, else than already mentioned, for the demand u > 0 and |g_{rr} > 1|?
I note, there was established so-called effective potential in the GR. In the radial problem (with zero angular momentum), this potential can be given as U = c^{2}*sqrt(1 - 2*u/r). We can see that U > 0 even if u < 0. Another question can be: do we need (or why do we need) two concepts of potential in GR? (The effective potential and that corresponding to relation u = r*Psi/c^{2}.)