We address the 'plane of satellites problem' by studying planar configurations around two disc galaxies with no late major mergers, formed in zoom-in hydro-simulations. Due to the current lack of good quality kinematic data for M31 satellites, we use only positional information. So far, positional analyses of simulations are unable to find planes as thin and populated as the observed ones. Moreover, they miss systematicity and detail in the plane-searching techniques, as well as in the study of the properties and quality of planes, both in simulations or real data. To fill this gap, i) we extend the 4-galaxy-normal density plot method (Pawlowski et al. 2013) in a way designed to efficiently identify the best quality planes (i.e., thin and populated) without imposing extra constraints on their properties, and ii), we apply it for the first time to simulations. Using zoom-in simulations allows us to mimic MW/M31-like systems regarding the number of satellites involved as well as the galactic disc mass and morphology, in view of possible disc effects. At all timesteps analyzed in both simulations we find satellite planar configurations that are compatible, along given time intervals, with all the spatial characteristics of observed planes identified using the same methodology. However, the fraction of co-orbiting satellites within them is in general low, suggesting time-varying satellite membership. We conclude that high-quality positional planes of satellites are not infrequent in LCDM-formed disc galaxies with a quiet assembly history. Detecting kinematically-coherent, time-persistent planes demands considering the full six-dimensional phase-space information of satellites.

Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$\alpha$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[\alpha$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[\alpha$/Fe] enhanced (i.e., old) populations show an homogeneous behaviour. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the $\rm \Lambda CDM$ model, we have found that the two phases of halo mass assembly (an early, fast phase, followed by a slow one, with low mass assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing has only a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients, and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars to suffer from strong radial mixing. By linking the two-component disk concept with the two-phase halo mass assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event marking periods characterized by different physical conditions under which thick and thin disk stars were born.

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

We present a detailed analysis of the global and fine structure of four middle-mass disc galaxies obtained from simulations in a $\Lambda$CDM scenario. These objects have photometric D/T ratios in good agreement with those observed for late-type spirals, as well as kinematic properties in agreement with the observational Tully-Fisher relation. We identify the different dynamical components at z=0 on the basis of both orbital parameters and the binding energy of stars in the galaxy. In this way, we recognize a slowly rotating centrally concentrated spheroid, and two disc components supported by rotation: a thin disc with stars in nearly circular orbits, and a thick disc with orbital parameters transitional between the thin disc and the spheroid. The spheroidal component is composed mainly by old, metal-poor and {\alpha}-enhanced stars. The distribution of metals in this component shows, however, a clear bimodality with a low-metallicity peak, which could be related to a classical bulge, and a high-metallicity peak, which could be related to a pseudo-bulge. The thin disc appears in our simulations as the youngest and most metal-rich component. The radial distribution of ages and colours in this component are U-shaped: the new stars are forming in the inner regions, and then migrate through secular processes. Finally, we also find a thick disc containing about 16% of the total stellar mass and with properties that are intermediate between those of the thin disc and the spheroid. Its low-metallicity stars are {\alpha}-enhanced when compared to thin disc stars of the same metallicity. The structural parameters (e.g., the scale height) of the simulated thick discs suggest that such a component could result from the combination of different thickening mechanisms that include merger-driven processes, but also long-lived internal perturbations of the thin disc. [Abridged]

The mass assembly and star formation histories of massive galaxies identified at low redshift z in different cosmological hydrodynamical simulations, have been studied through a detailed follow-up backwards in time of their constituent mass elements (sampled by particles) of different types. Then, the configurations they depict at progressively higher zs have been analysed. The analyses show that these histories share common generic patterns, irrespective of particular circumstances. In any case, the results we have found are different depending on the particle type. The most outstanding differences follow. We have found that by z ~ 3.5 - 6, mass elements identified as stellar particles at z=0 exhibit a gaseous cosmic-web-like morphology with scales of ~ 1 physical Mpc, where the densest mass elements have already turned into stars by z ~ 6. These settings are in fact the densest pieces of the cosmic web, where no hot particles show up, and dynamically organized as a hierarchy of flow convergence regions, that is, attraction basins for mass flows. On the other hand, mass elements identified at the diffuse hot coronae surrounding massive galaxies at z = 0, do not display a clear web-like morphology at any z. Diffuse gas is heated when flow convergence regions go through contractive deformations, and most of it keeps hot and with low density along the evolution. To shed light on the physical foundations of the behaviour our analyses show up, as well as on their possible observational implications, these patterns have been confronted with some generic properties of singular flows as described by the adhesion model. We have found that these common patterns simulations show can be interpreted as a consequence of flow properties, that, moreover, could explain different generic observational results on massive galaxies or their samples. We briefly discuss some of them.[Abridged]

Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow certain scaling relations. These relations define virial planes at the halo scale, whereas at the galactic scale they define intrinsic dynamical planes that are, however, tilted relative to the virial plane. Therefore, we predict that massive galaxies must lie on fundamental planes from their formation. We briefly discuss the physical origin of the tilt in terms the physical processes underlying massive galaxy formation at high z, in the context of a two-phase galaxy formation scenario. Specifically, we have found that it lies on the different behavior of the gravitationally heated gas as compared with cold gas previously involved in caustic formation, and the mass dependence of the energy available to heat the gas.

We present SPH cosmological simulations of the formation of three disk galaxies with a detailed treatment of chemical evolution and cooling. The resulting galaxies have properties compatible with observations: relatively high disk-to-total ratios, thin stellar disks and good agreement with the Tully-Fisher and the luminosity-size relations. They present a break in the luminosity profile at 3.0 +- 0.5 disk scale lengths, while showing an exponential mass profile without any apparent breaks, in line with recent observational results. Since the stellar mass profile is exponential, only differences in the stellar populations can be the cause of the luminosity break. Although we find a cutoff for the star formation rate imposed by a density threshold in our star formation model, it does not coincide with the luminosity break and is located at 4.3 +- 0.4 disk scale lengths, with star formation going on between both radii. The color profiles and the age profiles are "U-shaped", with the minimum for both profiles located approximately at the break radius. The SFR to stellar mass ratio increases until the break, explaining the coincidence of the break with the minimum of the age profile. Beyond the break we find a steep decline in the gas density and, consequently, a decline in the SFR and redder colors. We show that most stars (64-78%) in the outer disk originate in the inner disk and afterwards migrate there. Such stellar migrations are likely the main origin of the U-shaped age profile and, therefore, of the luminosity break.

We describe an SPH model for chemical enrichment and radiative cooling in cosmological simulations of structure formation. This model includes: i) the delayed gas restitution from stars by means of a probabilistic approach designed to reduce the statistical noise and, hence, to allow for the study of the inner chemical structure of objects with moderately high numbers of particles; ii) the full dependence of metal production on the detailed chemical composition of stellar particles by using, for the first time in SPH codes, the Qij matrix formalism that relates each nucleosynthetic product to its sources; and iii) the full dependence of radiative cooling on the detailed chemical composition of gas particles, achieved through a fast algorithm using a new metallicity parameter zeta(T) that gives the weight of each element on the total cooling function. The resolution effects and the results obtained from this SPH chemical model have been tested by comparing its predictions in different problems with known theoretical solutions. We also present some preliminary results on the chemical properties of elliptical galaxies found in self-consistent cosmological simulations. Such simulations show that the above zeta-cooling method is important to prevent an overestimation of the metallicity-dependent cooling rate, whereas the Qij formalism is important to prevent a significant underestimation of the [alpha/Fe] ratio in simulated galaxy-like objects.

[Abridged] We have analysed the parameters characterising the mass, size and velocity dispersion both at the baryonic scale and at the halo scales of two samples of relaxed elliptical-like-objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. At the halo scale they have been found to satisfy virial relations; at the scale of the baryonic object the (logarithms of the) ELO stellar masses, projected stellar half-mass radii, and stellar central l.o.s. velocity dispersions define a flattened ellipsoid close to a plane (the intrinsic dynamical plane, IDP), tilted relative to the virial one, whose observational manifestation is the observed FP. The ELO samples have been found to show systematic trends with the mass scale in both, the relative content and the relative distributions of the baryonic and the dark mass ELO components, so that homology is broken in the spatial mass distribution (resulting in the IDP tilt), but ELOs are still a two-parameter family where the two parameters are correlated. The physical origin of these trends presumably lies in the systematic decrease, with increasing ELO mass, of the relative amount of dissipation experienced by the baryonic mass component along ELO stellar mass assembly. ELOs also show kinematical segregation, but it does not appreciably change with the mass scale. The non-homogeneous population of IDPs explains the role played by the virial mass to determine the correlations among intrinsic parameters. In this paper we also show that the central stellar line-of-sight velocity dispersion of ELOs, is a fair empirical estimator of the virial mass, and this explains the central role played by this quantity at determining the observational correlations.

[Abridged] We have analysed the mass and velocity distributions of two samples of relaxed elliptical-like-objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. Our analysis shows that they are embedded in extended, massive dark matter haloes, and they also have an extended corona of hot diffuse gas. Dark matter haloes have experienced adiabatic contraction along their assembly process. The relative ELO dark- to bright-mass content and space distributions show broken homology, and they are consistent with observational results on the dark matter fraction at the central regions, as well as on the gradients of the mass-to-light ratio profiles for boxy ellipticals, as a function of their stellar masses. These results indicate that massive ellipticals miss stars at their central regions, as compared with less massive ones. Our simulations indicate that these missing baryons could be found beyond the virial radii as a hot, diffuse plasma. The projected stellar mass profiles of our ELOs can be well fit by the S\'ersic (1968) formula. Parameters of these fits reproduce observational correlations of the shape parameter with size, luminosity and velocity dispersion. The total mass density profiles show a power-law behaviour over a large r/r_{vir} interval, consistent with data on massive lens ellipticals at shorter radii. The velocity dispersion profiles show kinematical segregation, with no systematic mass dependence and a positive anisotropy, roughly independent of the radial distance outside the central regions. The consistency with observations strongly suggests that they could also describe important intrinsic characteristics of real ellipticals, as well as some of their properties recently inferred from observational data.

Combined with other CMB experiments, the WMAP survey provides an accurate estimate of the baryon density of the Universe. In the framework of the standard Big Bang Nucleosynthesis (BBN), such a baryon density leads to predictions for the primordial abundances of $^{4}$He and D in good agreement with observations. However, it also leads to a significant discrepancy between the predicted and observed primordial abundance of $^{7}$Li. Such a discrepancy is often termed as 'the lithium problem'. In this paper, we analyze this problem in the framework of scalar-tensor theories of gravity. It is shown that an expansion of the Universe slightly slower than in General Relativity before BBN, but faster during BBN, solves the lithium problem and leads to $^4$He and D primordial abundances consistent with the observational constraints. This kind of behavior is obtained in numerous scalar-tensor models, both with and without a self-interaction potential for the scalar field. In models with a self-interacting scalar field, the convergence towards General Relativity is ensured without any condition, thanks to an attraction mechanism which starts to work during the radiation-dominated epoch.

We present results of a study on the evolution of the parameters characterizing the structure and dynamics of the relaxed elliptical-like objects (ELOs) identified at z=0, z=1 and z=1.5 in a set of hydrodynamical, self-consistent simulations operating in the context of a concordance cosmological model. The values of the stellar mass, the stellar half-mass radius and the stellar mean-square velocity have been measured in each ELO and found to populate, at any z, a flattened ellipsoid close to a plane (the dynamical plane, DP). Our simulations indicate that, at the intermediate zs considered, individual ELOs evolve, increasing the values of these parameters as a consequence of on-going mass assembly, but, nevertheless, their DP is roughly preserved within its scatter, in agreement with observations of the Fundamental Plane of ellipticals at different zs. We briefly discuss how this lack of significant dynamical and structural evolution in ELO samples arises, in terms of the two different phases operating in the mass aggregation history of their dark matter halos. According with our simulations, most dissipation involved in ELO formation takes place at the early violent phase, causing the stellar mass, the stellar half-mass radius and the stellar mean-square velocity parameters to settle down to the DP, and, moreover, the transformation of most of the available gas into stars. In the subsequent slow phase, ELO stellar mass growth preferentially occurs through non-dissipative processes, so that the DP is preserved and the ELO star formation rate considerably decreases. These results hint, for the first time, to a possible way of explaining, in the context of cosmological simulations, different apparently paradoxical observational results on ellipticals.

We report on a study of the parameters characterizing the mass and velocity distributions of two samples of relaxed elliptical-like objects (ELOs) identified, at z=0, in a set of self-consistent hydrodynamical simulations operating in the context of a concordance cosmological model. Star formation (SF) has been phenomenologically implemented in the simulations in the framework of the turbulent sequential scenario through a threshold gas density and an efficiency parameter. Each ELO sample is characterized by the values these parameters take. We have found that the (logarithms of the) ELO stellar masses, projected stellar half-mass radii, and stellar central line-of-sight (LOS) velocity dispersions define dynamical fundamental planes (FPs). Zero points depend on the particular values that the SF parameters take, while slopes do not change. The ELO samples have been found to show systematic trends with the mass scale in both the relative content and the relative distributions of the baryonic and the dark mass ELO components. The physical origin of these trends lies in the systematic decrease, with increasing ELO mass, of the relative dissipation experienced by the baryonic mass component along ELO mass assembly, resulting in a tilt of the dynamical FP relative to the virial plane. The dynamical FPs shown by the two ELO samples are consistent with that shown by the SDSS elliptical sample in the same variables, with no further need for any relevant contribution from stellar population effects to explain the observed tilt. These effects could, however, have contributed to the scatter of the observed FP, as the dynamical FPs have been found to be thinner than the observed one. The results we report on hint, for the first time, at a possible way to understand the tilt of the observed FP in a cosmological context.

We present results of a study of the stellar age distributions in the sample of elliptical-like objects (ELOs) identified at z=0 in four simulations operating in the context of a concordance cosmological model. The simulations show that the formation of most stars in each ELO of the sample is a consequence of violent dynamical events, either fast multiclump collapse at high z, or mergers at lower z. This second way can explain the age spread as well as the dynamical peculiarities observed in some ellipticals, but its relative weight is never dominant and decreases as the ELO mass at the halo scale, $M_{vir}$, increases, to such an extent that some recent mergers contributing an important fraction to the total ELO mass can possibly contribute only a small fraction of new born stars. More massive objects have older means and narrower spreads in their stellar age distributions than less massive ones. The ELO sample shows also a tight correlation between $M_{vir}$ and the central stellar l.o.s. velocity dispersion, $\sigma_{los}$. This gives a trend of the means and spreads of ELO stellar populations with $\sigma_{los}$ that is consistent, even quantitatively, with the age effects observationally detected in the stellar populations of elliptical galaxies. Therefore, these effects can be explained as the observational manifestation of the intrinsic correlations found in the ELO sample between $M_{vir}$ and the properties of the stellar age distribution, on the one hand, and $M_{vir}$ and $\sigma_{los}$, on the other hand. These correlations hint, for the first time, at a possible way to reconcile age effects in ellipticals, and, particularly, the increase of $\alpha / <Fe>$ ratios with $\sigma_{los}$, with the hierarchical clustering paradigm.

We present results of an analysis of the structural and kinematical properties of a sample of elliptical-like objects (ELOs) identified in four hydrodynamical, self-consistent simulations run with the DEVA code (Serna et al. 2003). Star formation has been implemented in the code through a simple phenomenological parameterization, that takes into account stellar physics processes only implicitly through the values of a threshold gas density, $\rho_{\rm g,thres}$, and an efficiency parameter, $c_*$. The four simulations operate in the context of a $\Lambda$CDM cosmological model consistent with observations, resolve ELO mass assembly at scales up to $\simeq$ 2 kpc, and differ in the values of their star formation parameters. Stellar masses, projected half-mass radii and central l.o.s. velocity dispersions, $\sigma_{\rm los, 0}$, have been measured on the ELO sample and their values compared with data from the Sloan digital sky survey. For the first time in self-consistent simulations, a good degree of agreement has been shown, including the Faber-Jackson and the $D_n - \sigma_{\rm los, 0}$ relations, among others, but only when particular values of the $\rho_{\rm g,thres}$ and $c_*$ parameters are used. This demostrates the effect that the star formation parameterization has on the ELO mass distribution. Additionally, our results suggest that it is not strictly necessary, at the scales resolved in this work, to appeal to energy sources other than gravitational (as for example supernovae feedback effects) to account for the structure and kinematics of large ellipticals.

We describe DEVA, a multistep AP3M-like-SPH code particularly designed to study galaxy formation and evolution in connection with the global cosmological model. This code uses a formulation of SPH equations which ensures both energy and entropy conservation by including the so-called $\bn h$ terms. Particular attention has also been paid to angular momentum conservation and to the accuracy of our code. We find that, in order to avoid unphysical solutions, our code requires that cooling processes must be implemented in a non-multistep way. We detail various cosmological simulations which have been performed to test our code and also to study the influence of the $\bn h$ terms. Our results indicate that such correction terms have a non-negligible effect on some cosmological simulations, especially on high density regions associated either to shock fronts or central cores of collapsed objects. Moreover, they suggest that codes paying a particular attention to the implementation of conservation laws of physics at the scales of interest, can attain good accuracy levels in conservation laws with limited computational resources.

In this paper, we analyze the conditions for convergence toward General Relativity of scalar-tensor gravity theories defined by an arbitrary coupling function $\alpha$ (in the Einstein frame). We show that, in general, the evolution of the scalar field $(\phi)$ is governed by two opposite mechanisms: an attraction mechanism which tends to drive scalar-tensor models toward Einstein's theory, and a repulsion mechanism which has the contrary effect. The attraction mechanism dominates the recent epochs of the universe evolution if, and only if, the scalar field and its derivative satisfy certain boundary conditions. Since these conditions for convergence toward general relativity depend on the particular scalar-tensor theory used to describe the universe evolution, the nucleosynthesis bounds on the present value of the coupling function, $\alpha_0$, strongly differ from some theories to others. For example, in theories defined by $\alpha \propto \mid\phi\mid$ analytical estimates lead to very stringent nucleosynthesis bounds on $\alpha_0$ ($\lesssim 10^{-19}$). By contrast, in scalar-tensor theories defined by $\alpha \propto \phi$ much larger limits on $\alpha_0$ ($\lesssim 10^{-7}$) are found.

We present a method which enables exact solutions to be found for at homogeneous and isotropic scalar-tensor cosmologies with an arbitrary $\omega(\Phi)$ function and satisfying the general perfect fluid state equation $P=(\gamma-1)\rho c^2$. This method has been used to analyze a wide range of asymptotic analytical solutions at early and late times for different epochs in the cosmic history: false vacuum inflationary models, vacuum and radiation-dominated models, and matter-dominated models. We also describe the qualitative behavior of models at intermediary times and give exact solutions at any time for some particular scalar-tensor theories.

We study the stability of rotating collisionless self-gravitating spherical systems by using high resolution N-body experiments on a Connection Machine CM-5. We added rotation to Ossipkov-Merritt (hereafter OM) anisotropic spherical systems by using two methods. The first method conserves the anisotropy of the distribution function defined in the OM algorithm. The second method distorts the systems in velocity space. We then show that the stability of systems depends both on their anisotropy and on the value of the ratio between the total kinetic energy and the rotational kinetic energy. We also test the relevance of the stability parameters introduced by Perez et al. (1996) for the case of rotating systems.

In a previous paper (Serna & Alimi 1996), we have pointed out the existence of some particular scalar-tensor gravity theories able to relax the nucleosynthesis constraint on the cosmic baryonic density. In this paper, we present an exhaustive study of primordial nucleosynthesis in the framework of such theories taking into account the currently adopted observational constraints. We show that a wide class of them allows for a baryonic density very close to that needed for the universe closure. This class of theories converges soon enough towards General Relativity and, hence, is compatible with all solar-system and binary pulsar gravitational tests. In other words, we show that primordial nucleosynthesis does not always impose a very stringent bound on the baryon contribution to the density parameter.

We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed Particle Hydrodynamics) method while, unlike most previous implementations of SPH, gravitational forces are computed by a PP (Particle-Particle) approach. another important feature of this code is that hydrodynamics equations optionally include the correction terms appearing when $h(t,\br)$ is not constant. Our code has been implemented by using the data parallel programming model on CM5(CM). PPASPH has been applied to study the importance of adaptative smoothing correction terms on the entropy conservation. We confirm Hernquist's (1993) interpretation of the entropy violation observed in previous SPH simulations as a result of having neglected these terms. An improvement on the entropy conservation is not found by merely considering larger numbers of particles or different $N_S$ choices. The correct continuum description is only obtained if the $\bn h$ correction terms are included. Otherwise, the entropy conservation is always rather poor as compared to that found for the total energy.

We present a detailed calculation of the light element production in the framework of Scalar-Tensor (ST) theories of gravitation. The coupling function \omega has been described by an appropriate form which reproduces all the possible asymptotic behaviors at early times of viable ST cosmological models with a monotonic \omega(\Phi). This form gives an exact representation for most of the particular theories proposed in the literature, but also a first-order approximation to many other theories. In most of ST theories, the nucleosynthesis bounds lead to cosmological models which do not significantly differ from the standard FRW ones. We have found however a particular class of ST theories where the expansion rate of the universe during nucleosynthesis can be very different from that found in GR, while the present value of \omega is high enough to ensure compatibility with solar-system experiments. In the framework of this class of theories, right primordial abundances can be obtained for a baryon density range much wider (2.8 \lesssim \eta_{10} \lesssim 58.7) than in GR. Consequently, the usual constraint on the baryon contribution to the density parameter of the universe can be drastically relaxed (0.01 \lesssim \Omega_{b0} \lesssim 1.38) by considering these gravity theories. This is the first time that a ST theory is found to be compatible both with primordial nucleosynthesis and solar-system experiments while implying cosmological models significantly different from the FRW ones.

We analyze the qualitative behaviors of scalar-tensor cosmologies with an arbitrary monotonic $\omega(\Phi)$ function. In particular, we are interested on scalar-tensor theories distinguishable at early epochs from General Relativity (GR) but leading to predictions compatible with solar-system experiments. After extending the method developed by Lorentz-Petzold and Barrow, we establish the conditions required for convergence towards GR at $t\rightarrow\infty$. Then, we obtain all the asymptotic analytical solutions at early times which are possible in the framework of these theories. The subsequent qualitative evolution, from these asymptotic solutions until their later convergence towards GR, has been then analyzed by means of numerical computations. From this analysis, we have been able to establish a classification of the different qualitative behaviors of scalar-tensor cosmological models with an arbitrary monotonic $\omega(\Phi)$ function.

We propose a new hierarchical method which uses dynamical arguments to find and describe substructures in galaxy clusters. This method (hereafter h--method or h--analysis) uses a hierarchical clustering analysis to determine the relationship between galaxies according to their relative binding energies. We have tested from N-body simulations, the two following features of the proposed method: 1) It extracts subgroups which are much more stable during the cluster evolution than those given by other techniques. 2) There exists a reasonable similarity between the structures found when only the coordinates ($x, y, v_z$) provided by "observations" are considered, and those found by using the six phase--space coordinates We have applied this method to two Abell clusters: ABCG151 and ABCG2670. Our results imply that ABCG151 is separated into two clusters, one of them is again divided into two subclusters. ABCG2670 has no subclustering. Our method allows however to extract the most bound galaxies in its dynamical core.

We have performed a series of N-body experiments on Connection Machine-5 in order to simulate the formation of galaxy clusters gravitationally dominated by a massive dark background. In accordance with previous authors we find an extremely inhomogeneous evolution where subcondensations are continually formed and merged. The final distribution of galaxies is more centrally condensed than that of dark matter particles. We have analyzed the origin of this galaxy-dark matter segregation and also the origin of subclustering leading to this effect. We have then analysed, analytically and numerically, the dependence of final segregation on the physical parameters characterizing the model of protocluster. We also find that such a segregation persists even when inelastic encounters of dark halos around galaxies are taken into account. We conclude that this effect cannot be in general avoided in any hierarchical clustering scenario.