We present an analysis of the Counts-in-Cells (CiC) statistics of subhalos in the publically available Illustris-TNG cosmological simulations (TNG100-1, TNG100-3 and TNG300-3), in redshifts ranging from $z = 0$ to $z=5$. We fit the CiC statistics to two models: the gravitational quasi-equilibrium distribution (GQED) and the negative binomial distribution (NBD). We analyzed differences between ``full'' runs and their associated ``dark-only'' counterparts, in terms of GQED and NBD, and also by the use of gaussian kernel density estimates. We discuss the behaviour of the best-fit parameters $\bar{N}$ (average number of subhalos per cell), and the clustering parameters $b$ (GQED) and $g$ (NBD), in terms of redshifts and comoving cell sizes. We found that ``full'' and ``dark-only'' TNG runs follow similar CiC distributions. However, there are quantifiable distinctions between these sets of runs in the evolution and scaling of the best-fit parameters. We compare our results with previous numerical, observational and ``mock'' catalogs, and found similar scaling and evolutionary trends in all samples, up to factors in the amplitude of parameter values, possibly regulated by different magnitude cutoffs in the observational samples. Our results suggest the overall validity of the GQED and NBD models for describing the large-scale statistics of extended, multi-component gravitational systems and subject to complex baryonic physics. Our analysis indicate that gravithermodynamics is generally suitable for describing gravitational clustering in the presence of baryonic physics and, with proper refinements, it may provide additional theoretical guidance for the modelling of bias.

We revisit the "two-component virial theorem" (2VT) in the light of recent theoretical and observational results related to the "dark matter"(DM) problem. This modification of the virial theorem offers a physically meaningful framework to investigate possible dynamical couplings between the baryonic and DM components of extragalactic systems. In particular, we examine the predictions of the 2VT with respect to the "acceleration-discrepancy relation" (ADR). Considering the combined data (composed of systems supported by rotation and by velocity dispersion), we find that: (i) the overall behavior of the 2VT is consistent with the ADR; and (ii) the 2VT predicts a nearly constant behavior in the lower acceleration regime, as suggested in recent data on dwarf spheroidals. We also briefly comment on possible differentiations between the 2VT and some modified gravity theories.

We report on constraints to the cosmological jerk parameter ($j$) and to possible variability of the fine-structure constant ($\Delta \alpha/\alpha$) based on stochastic quintessence models of dark energy, discussed by Chongchitnan and Efstathiou (2007). We confirm the results by these authors in the sense that many viable solutions can be obtained, obeying current observational constraints in low redshifts. We add the observables $j$ and $\Delta \alpha/\alpha$ to this conclusion. However, we find peculiarities that may produce, in the nearby Universe, potential observational imprints in future cosmological data. We conclude, for redshifts $z \lesssim 3$, that: {\it (i) } $j(z)$ fluctuates due to the stochasticity of the models, reaching an amplitude of up to $5\%$ relatively to the $\Lambda$CDM model value ($j_{\Lambda {\rm CDM}}=1$); and {\it (ii)} by contrasting two distinct ("extreme") types of solutions, variabilities in $\alpha(z)$, linked to a linear coupling ($\zeta$) between the dark energy and electromagnetic sectors, are weakly dependent on redshift, for couplings of the order $|\zeta| \sim 10 ^{-4}$, even for large variations in the equation of state parameter at relatively low redshifts. Nonlinear couplings produce an earlier and steeper onset of the evolution in $\Delta \alpha/\alpha (z)$, but can still accommodate the data for weak enough couplings.

Stellar and galactic systems are objects in dynamical equilibrium that are composed of ordinary baryonic matter hypothetically embedded in extended dominant dark matter halos. Our aim is to investigate the scaling relations and dissipational features of these objects over a wide range of their properties, taking the dynamical influence of the dark matter component into account. We study the physical properties of these self-gravitating systems using the two-component virial theorem in conjunction with data that embrace a wide range of astrophysical systems. We find that the scaling relations defined by the properties of these objects admit a dark-to-luminous density ratio parameter as a natural requirement in this framework. We also probe dissipational effects on the fundamental surface defined by the two-component virial theorem and discuss their relations with respect to the region devoid of objects in the data distribution. Our results indicate complementary contributions of dissipation and dark matter to the orign of scaling relations in astrophysical systems.

In this paper we present a study of the cosmic flow from the point of view of how clusterings at different dynamical regimes in an expanding universe evolve according to a `coarse-grained' partitioning of their ranked energy distribution. By analysing a Lambda-CDM cosmological simulation from the Virgo Project, we find that cosmic flows evolve in an orderly sense, when tracked from their coarse-grained energy cells, even when nonlinearities are already developed. We show that it is possible to characterize scaling laws for the Pairwise Velocity Distribution in terms of the energy cells, generally valid at the linear and nonlinear clustering regimes.

We present an analysis of the behaviour of the `coarse-grained' (`mesoscopic') rank partitioning of the mean energy of collections of particles composing virialized dark matter halos in a Lambda-CDM cosmological simulation. We find evidence that rank preservation depends on halo mass, in the sense that more massive halos show more rank preservation than less massive ones. We find that the most massive halos obey Arnold's theorem (on the ordering of the characteristic frequencies of the system) more frequently than less massive halos. This method may be useful to evaluate the coarse-graining level (minimum number of particles per energy cell) necessary to reasonably measure signatures of `mesoscopic' rank orderings in a gravitational system.

We revisit the integral formulation (or Green's function approach) of Einstein's equations in the context of braneworlds. The integral formulation has been proposed independently by several authors in the past, based on the assumption that it is possible to give a reinterpretation of the local metric field in curved spacetimes as an integral expression involving sources and boundary conditions. This allows one to separate source-generated and source-free contributions to the metric field. As a consequence, an exact meaning to Mach's Principle can be achieved in the sense that only source-generated (matter fields) contributions to the metric are allowed for; universes which do not obey this condition would be non-Machian. In this paper, we revisit this idea concentrating on a Randall-Sundrum-type model with a non-trivial cosmology on the brane. We argue that the role of the surface term (the source-free contribution) in the braneworld scenario may be quite subtler than in the 4D formulation. This may pose, for instance, an interesting issue to the cosmological constant problem.

We present a new set of dissipationless N-body simulations aiming to better understand the pure dynamical aspects of the ``Fundamental Plane'' (FP) of elliptical galaxies. We have extended our previous hierarchical merger scheme by considering the Hernquist profile for the initial galaxy model. Two-component Hernquist galaxy models were also used to study the effect of massive dark halos on the end-products characteristics. We have also performed new collapse simulations including initial spin. We found that the one-component Hernquist mergers give results similar to those found for the one-component King models, namely both were able to build-up small scatter FP-like correlations with slopes consistent with what is found for the near-infrared FP of nearby galaxies. The two-component models also reproduce a FP-like correlation, but with a significantly steeper slope. This is in agreement with what has been found for elliptical galaxies at higher redshift (0.1 $<$ z $<$ 0.6). We discuss some structural properties of the simulated galaxies and their ability to build-up FP-like correlations. We confirm that collapses generally do not follow a FP-like correlation regardless of the initial spin. We suggest that the evolution of gradients in the gravitational field of the merging galaxies may be the main ingredient dictating the final non-homology property of the end products.

Following on from the numerical work of Capelato, de Carvalho & Carlberg (1995, 1997), where dissipationless merger simulations were shown to reproduce the "Fundamental Plane" (FP) of elliptical galaxies, we investigate whether the end products of pure, spherically symmetric, one-component dissipationless {\it collapses} could also reproduce the FP. Past numerical work on collisionless collapses have addressed important issues on the dynamical/structural characteristics of collapsed equilibrium systems. However, the study of collisionless collapse in the context of the nature of the FP has not been satisfactorily addressed yet. Our aim in this paper is to focus our attention on the resulting collapse of simple one-component spherical models with a range of different initial virial coefficients. We find that the characteristic correlations of the models are compatible with virialized, centrally homologous systems. Our results strengthen the idea that merging may be a fundamental ingredient in forming non-homologous objects.

Motivated by present indirect evidences that galaxies are surrounded by dark matter halos, we investigate whether their physical properties can be described by a formulation of the virial theorem which explicitly takes into account the gravitational potential term representing the interaction of the dark halo with the barionic or luminous component. Our analysis shows that the application of such a ``two-component virial theorem'' not only accounts for the scaling relations displayed, in particular, by elliptical galaxies, but also for the observed properties of all virialized stellar systems, ranging from globular clusters to galaxy clusters.

We investigate the dynamical properties of the cluster of galaxies A3558 (Shapley 8). Studying a region of one square degree ($\sim$ 3 Mpc$^2$) centered on the cluster cD galaxy, we have obtained a statistically complete photometric catalog with positions and magnitudes of 1421 galaxies (down to a limiting magnitude of $B \sim 21$). This catalog has been matched to the recent velocity data obtained by Mazure et al. (1997) and from the literature, yielding a radial velocity catalog containing 322 galaxies. Our analysis shows that the position/velocity space distribution of galaxies shows significant substructure. A central bimodal core detected previously in preliminary studies is confirmed by using the Adaptive Kernel Technique and Wavelet Analysis. We show that this central bimodal subtructure is nevertheless composed of a projected feature, kinematically unrelated to the cluster, plus a group of galaxies probably in its initial merging phase into a relaxed core. The cD velocity offset with respect to the average cluster redshift, reported earlier by several authors, is completely eliminated as a result of our dynamical analysis. The untangling of the relaxed core component also allows a better, more reliable determination of the central velocity dispersion, which in turn eliminates the ``$\beta$-problem'' for A3558. The cluster also shows a ``preferential'' distribution of subclumps coinciding with the direction of the major axis position angle of the cD galaxy and of the central X-ray emission ellipsoidal distribution, in agreement with an anisotropic merger scenario.