Recently alternative approaches in cosmology seeks to explain the nature of dark matter as a direct result of the non-linear spacetime curvature due to different types of deformation potentials. In this context, a key test for this hypothesis is to examine the effects of deformation on the evolution of large scales structures. An important requirement for the fine analysis of this pure gravitational signature (without dark matter elements) is to characterize the position of a galaxy during its trajectory to the gravitational collapse of super clusters at low redshifts. In this context, each element in an gravitational N-body simulation behaves as a tracer of collapse governed by the process known as chaotic advection (or lagrangian turbulence). In order to develop a detailed study of this new approach we develop the COsmic LAgrangian TUrbulence Simulator (COLATUS) to perform gravitational N-body simulations based on Compute Unified Device Architecture (CUDA) for graphics processing units (GPUs). In this paper we report the first robust results obtained from COLATUS.

We study the star/galaxy classification efficiency of 13 different decision tree algorithms applied to photometric objects in the Sloan Digital Sky Survey Data Release Seven (SDSS DR7). Each algorithm is defined by a set of parameters which, when varied, produce different final classification trees. We extensively explore the parameter space of each algorithm, using the set of $884,126$ SDSS objects with spectroscopic data as the training set. The efficiency of star-galaxy separation is measured using the completeness function. We find that the Functional Tree algorithm (FT) yields the best results as measured by the mean completeness in two magnitude intervals: $14\le r\le21$ ($85.2%$) and $r\ge19$ ($82.1%$). We compare the performance of the tree generated with the optimal FT configuration to the classifications provided by the SDSS parametric classifier, 2DPHOT and Ball et al. (2006). We find that our FT classifier is comparable or better in completeness over the full magnitude range $15\le r\le21$, with much lower contamination than all but the Ball et al. classifier. At the faintest magnitudes ($r>19$), our classifier is the only one able to maintain high completeness ($>$80%) while still achieving low contamination ($\sim2.5%$). Finally, we apply our FT classifier to separate stars from galaxies in the full set of $69,545,326$ SDSS photometric objects in the magnitude range $14\le r\le21$.

Although the theoretical understanding of the nonlinear gravitational clustering has greatly advanced in the last decades, in particular by the outstanding improvement on numerical N-body simulations, the physics behind this process is not fully elucidated. The main goal of this work is the study of the possibility of a turbulent-like physical process in the formation of structures, galaxies and clusters of galaxies, by the action of gravity alone. We use simulation data from the Virgo Consortium, in ten redshift snapshots (from 0 to 10). From this we identify galaxy-sized and cluster-sized dark matter haloes, by using a FoF algorithm and applying a boundedness criteria, and study the gravitational potential energy spectra. We find that the galaxy-sized haloes energy spectrum follows closely a Kolmogorov's power law, similar to the behaviour of dynamically turbulent processes in fluids. This means that the gravitational clustering of dark matter may admit a turbulent-like representation.