“…With the same idea some new monotonicity formulas (even with non-harmonic functions) were then found to study static and sub-static metrics in general relativity [4,7,[11][12][13]. Recently, this level set approach with harmonic functions applied to certain nonparabolic Riemannian 3-manifolds with nonnegative scalar curvature have produced several results, we mention the sharp comparisons about the rate of decay of an energy-like quantity and the area of the level sets of the minimal positive Green function, obtained by Munteanu and Wang in [41] which, as an application, allowed Chodosh and Li [20] to prove a conjecture of Schoen on the stable minimal hypersurfaces in the Euclidean space R 4 and some new proofs of the well-known positive mass theorem [6,16]. Then, considering linearly growing harmonic functions (similarly to [16], which was influenced by a pioneering work of Stern [43]), some asymptotically flat versions of the spacetime positive mass theorem have been proven in [15,31], while the monotonicity formula in [6] has been extended in [3] to the nonlinear potential theoretic setting, replacing the harmonic functions by p-harmonic functions, to obtain a simpler proof of the Riemannian Penrose inequality with a single black hole.…”