“…About selecting a typical integrand among various options as mentioned in the Appendix, we must point out that as our starting point is motivated from the application point of view, we adopt the following properties to propose the functional form of the integrand in Equation (18): directional dependency which is encoded in the first derivative of the underlying field and also inspired by the definition of crossing statistics (Ryden 1988); intuitively, our suggestion belongs to the moment and cumulant definition of density field which is more reasonable compared to other complicated functions; taking other typical functions into account, namely, f (|∇δ ( r, s) | n , ∇ 2 δ ( r, s) , ∇ i ∇ j δ ( r, s) ∇ i ∇ j δ ( r, s) ,...) is in principle allowed but it turns out that the higher derivative the higher computational time consuming and even opens new room for the higher value of numerical uncertainty. Generally, the shear tensor ij and characteristic radius of local extrema (R * ∼ σ 1 /σ 2 ), are relevant when we are dealing with the local extrema (Bardeen et al 1986;Vafaei Sadr & Movahed 2021). As long as our purpose is focusing on directional dependency, we do not need to examine the extrema condition expressed by the second derivatives; consequently, a reasonable choice is adopting the first derivative of the density field.…”