A new family of nonrelativistic, Newtonian, non-quantum equilibrium configurations describing galactic halos is introduced, by considering strange quark matter conglomerates with masses larger than about 8 GeV as new possible components of the dark matter. Originally introduced to explain the state of matter in neutron stars, such conglomerates may also form in the high-density and temperature conditions of the primordial Universe and then decouple from ordinary baryonic matter, providing the fundamental components of dark matter for the formation of pristine gravitational potential wells and the subsequent evolution of cosmic structures. The obtained results for halo mass and radius are consistent with the rotational velocity curve observed in the Galaxy. Additionally, the average density of such dark matter halos is similar to that derived for halos of dwarf spheroidal galaxies, which can therefore be interpreted as downscaled versions of larger dark matter distributions around Milky Way-sized galaxies and hint for a common origin of the two families of cosmic structures.

Recent observations of the Virgo cluster and the Local Group suggested that some galaxies are flowing out from their parent cluster. This may be the signature that dark energy (DE) acts significantly also on small cosmological scales. By means of direct N-body simulations we performed several simulations, in which the effect of DE and gravity are taken into account, aiming to determine whether DE can produce an outflow of galaxies compatible with observations. Comparing the different simulations, our results suggest that the observed outflow of galaxies is likely due to the local effect of DE.

We discuss how the centres of galaxy clusters evolve in time, showing the results of a series of direct N-body simulations. In particular, we followed the evolution of a galaxy cluster with a mass $M_{clus} \simeq 10^{14} $M$_{\odot}$ in different configurations. The dynamical evolution of the system leads in all the cases to the formation of dense and massive sub-structures in the cluster centre, that form in consequence of a series of collisions and merging among galaxies travelling in the cluster core. We investigate how the structural properties of the main merging product depends on the characteristics of those galaxies that contributed to its formation.

We performed a series of direct N-body simulations with the aim to follow the dynamical evolution of a galaxy cluster (GC) ($M_{clus}\simeq 10^{14} M_{\odot}$) in different environment. The results show the formation of heavy sub-structures in the cluster centre in consequence of multiple merging among the innermost galaxies. Moreover we investigate the dynamics of super-massive black holes (SMBHs) residing in the centre of galaxies that form the most massive sub-structure.

Recently, observational data and high precision mapping of the local velocity field of Local Group and Virgo cluster have revealed a linear velocity-distance relation of the outermost galaxies, properly referred to as Local Hubble Flow. By means of direct N-body method, we performed several simulations in which a galaxy cluster undergoes the action of the Dark Energy force and of the gravitational one induced by the gas. We reproduced the so-called Hubble diagrams, to highlight the outflow of the galaxies lying in the external region of the cluster. Our preliminary results suggest that the observed outflow of galaxies is likely due to the local effect of Dark Energy. Furthermore, the accuracy of the N-body method used, allows us to follow the merging process among some galaxies with the aim to reproduce the formation of a single compact object in the centre of the cluster.

In this paper we study the equilibrium configurations of anisotropic self-gravitating fermions, by extending to general relativity the solutions obtained in a previous paper. This treatment also generalizes to anisotropic systems the relativistic self-gravitating Fermi gas model, by considering different degrees of anisotropy. We discuss some important characteristics of the models and the obtained density profiles, and generalize the relation between the anisotropy and the mass of particles in the relativistic regime. These relativistic models may also be applied to the study of superdense neutron stars with anisotropic pressure or super-Chandrasekhar white dwarfs generated by the presence of a magnetic field.

Systems of self-gravitating fermions constitute a topic of great interest in astrophysics, due to the wide field of applications. In this paper, we consider the gravitational equilibrium of spherically symmetric Newtonian models of collisionless semidegenerate fermions. We construct numerical solutions by taking into account the effects of the anisotropy in the distribution function and considering the prevalence of tangential velocity. In this way, our models generalize the solutions obtained for isotropic Fermi-Dirac statistics. We also extend the analysis to equilibrium configurations in the classical regime and in the fully degenerate limit, recovering, for different levels of anisotropy, hollow equilibrium configurations obtained in Maxwellian regime. Moreover, in the limit of full degeneracy, we find a direct expression relating the anisotropy with the mass of the particles composing the system.

We discuss the possibility to analyze the problem of gravothermal catastrophe in a new way, by obtaining thermodynamical equations to apply to a selfgravitating system. By using the King distribution function in the framework of statistical mechanics we treat the globular clusters evolution as a sequence of quasi-equilibrium thermodynamical states.

We summarize some observational comparison concerning the features of globular clusters (GCs) population in connection to the evolution of King models. We also make a comparison with some extragalactic GCs systems, in order to underline the effects of the main body on the dynamical evolution.

Context. The external regions of galaxy clusters may be under strong influence of the dark energy, discovered by observations of the SN Ia at redshift z < 1. Aims. The presence of the dark energy in the gravitational equilibrium equation, with the Einstein $\Lambda$ term, contrasts the gravity, making the equilibrium configuration more extended in radius. Methods. In this paper we derive and solve the kinetic equation for an equilibrium configuration in presence of the dark energy, by considering the Newtonian regime, being the observed velocities of the galaxies inside a cluster largely smaller than the light velocity. Results. The presence of the dark energy in the gravitational equilibrium equation leads to wide regions in the $W_0-\rho_{\Lambda}$ diagram where the equilibrium solutions are not permitted, due to the prevalence of the effects of the dark energy on the gravity.

We obtain a solution for the hydrodynamic outflow of the polytropic gas from the gravitating center, in presence of the uniform Dark Energy (DE). The antigravity of DE is enlightening the outflow and make the outflow possible at smaller initial temperature, at the same density. The main property of the wind in presence of DE is its unlimited acceleration after passing the critical point. In application of this solution to the winds from galaxy clusters we suggest that collision of the strongly accelerated wind with another galaxy cluster, or with another galactic cluster wind could lead to the formation of a highest energy cosmic rays.

{We consider the Coma cluster of galaxies as a gravitationally bound physical system embedded in the perfectly uniform static dark energy background as implied by the $\Lambda$CDM cosmology.} {We ask if the density of dark energy is high enough to affect the structure of a large rich cluster of galaxies?} {We use recent observational data on the cluster together with our theory of local dynamical effects of dark energy, including the zero-gravity radius $R_{\rm ZG}$ of the local force field as the key parameter.} {1) Three masses are defined which characterize the structure of a regular cluster: the matter mass $M_{\rm M}$, the dark-energy effective mass $M_{\rm DE}$ ($<0$) and the gravitating mass $M_{\rm G}$ ($= M_{\rm M} + M_{\rm DE}$). 2) A new matter density profile is suggested which reproduces well the observational data for the Coma cluster in the radius range from 1.4 Mpc to 14 Mpc and takes into account the dark energy background. 3) Using this profile, we calculate upper limits for the total size of the Coma cluster, $R \le R_{\rm ZG} \approx 20$ Mpc, and its total matter mass, $M_{\rm M} \la M_{\rm M}(R_{\rm ZG}) = 6.2 \times 10^{15} M_{\odot}$.} {The dark energy antigravity affects strongly the structure of the Coma cluster at large radii $R \ga 14$ Mpc and should be taken into account when its total mass is derived.}

We solve the equation of the equilibrium of the gravitating body, with a polytropic equation of state of the matter $P=K\rho^{\gamma}$, with $\gamma=1+1/n$, in the frame of the Newtonian gravity, with non-zero cosmological constant $\Lambda$. We consider the cases with $n=1,\,\,1.5,\,\,3$ and construct series of solutions with a fixed value of $\Lambda$. For each value of $n$, the non-dimensional equation of the static equilibrium has a family of solutions, instead of the unique solution of the Lane-Emden equation at $\Lambda=0$. The equilibrium state exists only for central densities $\rho_0$ larger than the critical value $\rho_c$. There are no static solutions at $\rho_0 < \rho_c$. We find the values of $\rho_c$ for each value of $n$ and show that the presence of dark energy decrease the dynamic stability of the configuration. We apply our results for analyzing the possibility of existence of equilibrium states for cluster of galaxies in the present universe with non-zero $\Lambda$.

We derive a variational principle for the dynamical stability of a cluster as a gas sphere in a box. Newtonian clusters are always dynamically stable and, for relativistic clusters, the relation between dynamical and thermodynamical instabilities is analyzed. The boundaries between dynamically and thermodynamically stable and unstable models are found numerically for relativistic stellar systems with different cut off parameters. A criterion based on binding energy curve is used for determination of the boundary of dynamical stability.

We introduce a one-dimensional toy model of globular clusters. The model is a version of the well-known gravitational sheets system, where we take additionally into account mass and energy loss by evaporation of stars at the boundaries. Numerical integration by the "exact" event-driven dynamics is performed, for initial uniform density and Gaussian random velocities. Two distinct quasi-stationary asymptotic regimes are attained, depending on the initial energy of the system. We guess the forms of the density and velocity profiles which fit numerical data extremely well and allow to perform an independent calculation of the self-consistent gravitational potential. Some power-laws for the asymptotic number of stars and for the collision times are suggested.