In this work we extend the perturbation theory for modified gravity (MG) in two main aspects. First, the construction of matter overdensities from Lagrangian displacement fields is shown to hold in a general framework, allowing us to find Standard Perturbation Theory (SPT) kernels from known Lagrangian Perturbation Theory (LPT) kernels. We then develop a theory of biased tracers for generalized cosmologies, extending already existing formalisms for $\Lambda$CDM. We present the correlation function in Convolution-LPT and the power spectrum in SPT for $\Lambda$CDM, $f(R)$ Hu-Sawicky, and DGP braneworld models. Our formalism can be applied to many generalized cosmologies and to facilitate it, we are making public a code to compute these statistics. We further study the relaxation of bias with the use of a simple model and of excursion set theory, showing that in general the bias parameters have smaller values in MG than in General Relativity.

Scalar fields have been used as candidates for dark matter in the universe, from axions with masses $\sim10^{-5}$eV until ultra-light scalar fields with masses $\sim10^{-22}$eV. Axions behave as cold dark matter while the ultra-light scalar fields galaxies are Bose-Einstein condensate drops. The ultra-light scalar fields are also called scalar field dark matter model. In this work we study rotation curves for low surface brightness spiral galaxies using two scalar field models: the Gross-Pitaevskii Bose-Einstein condensate in the Thomas-Fermi approximation and a scalar field solution of the Klein-Gordon equation. We also used the zero disk approximation galaxy model where photometric data is not considered, only the scalar field dark matter model contribution to rotation curve is taken into account. From the best-fitting analysis of the galaxy catalog we use, we found the range of values of the fitting parameters: the length scale and the central density. The worst fitting results (values of $\chi^2_{red}$ much greater than 1, on the average) were for the Thomas-Fermi models, i.e., the scalar field dark matter is better than the Thomas-Fermi approximation model to fit the rotation curves of the analysed galaxies. To complete our analysis we compute from the fitting parameters the mass of the scalar field models and two astrophysical quantities of interest, the dynamical dark matter mass within 300 pc and the characteristic central surface density of the dark matter models. We found that the value of the central mass within 300 pc is in agreement with previous reported results, that this mass is $\approx 10^{7}$ $M_\odot/$pc$^2$, independent of the dark matter model. And, on the contrary, the value of the characteristic central surface density do depend on the dark matter model.

In this work we study rotation curves of spiral galaxies using a model of dark matter based on a scalar-tensor theory of gravity. We show how to estimate the scalar field dark matter parameters using a sample of observed rotation curves.

In this work we present a model of dark matter based on scalar-tensor theory of gravity. With this scalar field dark matter model we study the non-linear evolution of the large scale structures in the universe. The equations that govern the evolution of the scale factor of the universe are derived together with the appropriate Newtonian equations to follow the non-linear evolution of the structures. Results are given in terms of the power spectrum that gives quantitative information on the large-scale structure formation. The initial conditions we have used are consistent with the so called concordance $\Lambda$CDM model.

We study the large-scale structure formation in the Universe in the frame of scalar-tensor theories as an alternative to general relativity. We review briefly the Newtonian limit of non-minimally coupled scalar-tensor theories and the evolution equations of the $N$-body system that is appropriate to study large-scale structure formation in the Universe. We compute the power-spectrum of the universe at present epoch and show how the large-scale structure depends on the scalar field contribution.

We present numerical $N$-body simulation studies of large-scale structure formation. The main purpose of these studies is to analyze the several models of dark matter and the role they played in the process of large-scale structure formation. We analyze the standard and more successful case, i.e., the cold dark matter with cosmological constant ($\Lambda$CDM). We compare the results of this model with the corresponding results of other alternative models, in particular, the models that can be built from the Newtonian limit of alternative theories of gravity like scalar-tensor theories. An specific model is the one that considers that the scalar field is non-minimally coupled to the Ricci scalar in the Einstein-Hilbert Lagrangian that gives, in the Newtonian limit an effective gravitational force that is given by two contributions: the standard Newtonian potential plus a Yukawa potential that comes from a massive scalar field. Comparisons of the models are done by analyzing the snapshots of the $N$-body system at z=0 for several values of the SF parameters.

We study bar formation in galactic disks as a consequence of the collision of two spiral galaxies under the influence of a potential which is obtained from the Newtonian limit of a scalar--tensor theory of gravity. We found that dynamical effects depend on parameters ($\alpha$, $\lambda$) of the theory. In particular, we observe that the bar is shorter for weaker tidal perturbations, which in turn corresponds to smaller values of $\lambda$ used in our numerical experiments.

We present simulations within the framework of scalar-tensor theories, in the Newtonian limit, to investigate the influence of massive scalar fields on the dynamics of the collision of two equal spherical clouds. We employ a SPH code modified to include the scalar field to simulate two initially non-rotating protogalaxies that approach each other, and as a result of the tidal interaction, intrinsic angular momentum is generated. We have obtained sufficient large values of J/M to suggest that intrinsic angular momentum can be the result of tidal interactions.

We build a spherical halo model for galaxies using a general scalar-tensor theory of gravity in its Newtonian limit. The scalar field is described by a time-independent Klein-Gordon equation with a source that is coupled to the standard Poisson equation of Newtonian gravity. Our model, by construction, fits both the observed rotation velocities of stars in spirals and a typical luminosity profile. As a result, the form of the new Newtonian potential, the scalar field, and dark matter distribution in a galaxy are determined. Taking into account the constraints for the fundamental parameters of the theory (lambda,alpha), we analyze the influence of the scalar field in the dark matter distribution, resulting in shallow density profiles in galactic centers.

We present numerical N-body simulation studies of large-scale structure formation. The main purpose of these studies is to analyze the several models of dark matter and the role they played in the process of large-scale structure formation. We analyze in this work a flat cold dark matter dominated model known as the Santa Barbara cluster. We compare the results for this model using the standard Newtonian limit of general relativity with the corresponding results of using the Newtonian limit of scalar-tensor theories. An specific model is the one that considers that the scalar field is non-minimally coupled to the Ricci scalar in the Einstein-Hilbert Lagrangian. Comparisons of the models are done showing results of rotation curves, density profiles, and velocity dispersions for halos formed at z=0. We analyze, in particular, the Santa Barbara cluster and its possible equation of state.

We present $N$-body cosmlogical simulations in the framework of the Newtonian limit of scalar-tensor theories of gravity. The scalar field is described by a modified Helmholtz equation with a source that is coupled to the standard Poisson equation of Newtonian gravity. The effective gravitational force is given by two contributions: the standard Newtonian potential plus a Yukawa potential stemming from massive scalar fields. In particular, we consider simulations of $\Lambda$CDM models and compute the density and velocity profiles of the most massive groups found at z=0.

We present $\Lambda$CDM $N$-body cosmological simulations in the framework of a static general scalar-tensor theory of gravity. Due to the influence of the non-minimally coupled scalar field, the gravitational potential is modified by a Yukawa type term, yielding a new structure formation dynamics. We present some preliminary results and, in particular, we compute the density and velocity profiles of the most massive group.

We computed flat rotation curves from scalar-tensor theories in their weak field limit. Our model, by construction, fits a flat rotation profile for velocities of stars. As a result, the form of the scalar field potential and DM distribution in a galaxy are determined. By taking into account the constraints for the fundamental parameters of the theory $(\lambda, \alpha)$, it is possible to obtain analytical results for the density profiles. For positive and negative values of $\alpha$, the DM matter profile is as cuspy as NFW's.

We use the Newtonian limit of a general scalar-tensor theory around a background field to study astrophysical effects. The gravitational theory modifies the standard Newtonian potential by adding a Yukawa term to it, which is quantified by two theoretical parameters: $\lambda$, the lenghtscale of the gravitational interaction and its strength, $\alpha$. Within this formalism we firstly present a numerical study on the formation of bars in isolated galaxies. We have found for positive $\alpha$ that the modified gravity destabilizes the galactic discs and leads to rapid bar formation in isolated galaxies. Values of $\lambda$ in the range $\approx 8$ -- 14 kpc produce strongest bars in isolated models. Then, we extent this work by considering tidal effects due to interacting galaxies. We send two spirals to collide and study the bar properties of the remnant. We characterize the bar kinematical properties in terms of our parameters ($\lambda, \alpha$).

The joint influence of numerical parameters such as the number of particles N, the gravitational softening length $\epsilon$ and the time-step $\Delta t$ is investigated in the context of galaxy simulations. For isolated galaxy models we have performed a convergence study and estimated the numerical parameters ranges for which the relaxed models do not deviate significantly from its initial configuration. By fixing N, we calculate the range of the mean interparticle separation $\lambda(r)$ along the disc radius. We have found that in the simulations with N=1310720 particles $\lambda$ varies by a factor of 6, and the corresponding final Toomre's parameters Q change by only about 5 per cent. By decreasing N, the $\lambda$ and Q ranges broaden. Large $\epsilon$ and small N cause an earlier bar formation. For a given set of parameters the disc heating is smaller with the Plummer softening than with the spline softening. For galaxy collision models numerical simulations indicate that the properties of the formed bars strongly depend upon the selection of N and $\epsilon$. Large values of the gravitational softening parameter and a small number of particles results in the rapid formation of a well defined, slowly rotating bar. On the other hand, small values of $\epsilon$ produce a small, rapidly rotating disc with tightly wound spiral arms, and subsequently a weak bar emerges. We have found that by increasing N, the bar properties converge and the effect of the softening parameter diminishes. Finally, in some cases short spiral arms are observed at the ends of the bar that change periodically from trailing to leading and vice-versa - the wiggle.

We present a formulation for potential-density pairs to describe axisymmetric galaxies in the Newtonian limit of scalar-tensor theories of gravity. The scalar field is described by a modified Helmholtz equation with a source that is coupled to the standard Poisson equation of Newtonian gravity. The net gravitational force is given by two contributions: the standard Newtonian potential plus a term stemming from massive scalar fields. General solutions have been found for axisymmetric systems and the multipole expansion of the Yukawa potential is given. In particular, we have computed potential-density pairs of galactic disks for an exponential profile and their rotation curves.

A family of potential--density pairs has been found for spherical halos and bulges of galaxies in the Newtonian limit of scalar--tensor theories of gravity. The scalar field is described by a Klein--Gordon equation with a source that is coupled to the standard Poisson equation of Newtonian gravity. The net gravitational force is given by two contributions: the standard Newtonian potential plus a term stemming from massive scalar fields. General solutions have been found for spherical systems. In particular, we compute potential--density pairs of spherical galactic systems, and some other astrophysical quantities that are relevant to generating initial conditions for spherical galaxy simulations.