We present an algorithm for inferring the dynamical mass of galaxy clusters directly from their respective phase-space distributions, i.e. the observed line-of-sight velocities and projected distances of galaxies from the cluster centre. Our method employs neural flows, a deep neural network capable of learning arbitrary high-dimensional probability distributions, and inherently accounts, to an adequate extent, for the presence of interloper galaxies which are not bounded to a given cluster, the primary contaminant of dynamical mass measurements. We validate and showcase the performance of our neural flow approach to robustly infer the dynamical mass of clusters from a realistic mock cluster catalogue. A key aspect of our novel algorithm is that it yields the probability density function of the mass of a particular cluster, thereby providing a principled way of quantifying uncertainties, in contrast to conventional machine learning approaches. The neural network mass predictions, when applied to a contaminated catalogue with interlopers, have a logarithmic residual scatter which goes down to 0.043 dex for the most massive clusters. This is nearly an order of magnitude improvement over the classical cluster mass scaling relation with the velocity dispersion, and outperforms recently proposed machine learning approaches. We also apply our neural flow mass estimator to a compilation of galaxy observations of some well-studied clusters with robust dynamical mass estimates, further substantiating the efficacy of our algorithm. Such sophisticated approaches would undoubtedly be relevant for robust and efficient dynamical mass inference from upcoming surveys covering unprecedented volumes of the sky.

In this paper, we study the chameleon profile in inhomogeneous density distributions and find that the fifth force in thin shell near the surface is weaker from what expected in homogeneous density distributions. Also, we check the validity of quasi-static approximation for the chameleon scalar field in the astrophysical time scales. We have investigated the rolling down behaviour of the scalar field on its effective potential inside a one solar mass red giant star by using MESA code. We have found that the scalar field is fast enough to follow the minimum of the potential. This adiabatic behaviour reduces the fifth force and extends the screened regions to lower densities where the field has smaller mass and was expected to be unscreened. As a consequence, the star evolution is similar to what expected from standard general relativity. In addition, considering the stability of star, an approximate constraint on the coupling constant $\beta$ is found.