We study the effects of dark energy (DE) anisotropic stress on features of the matter power spectrum (PS). We employ the Parametrized Post-Friedmannian (PPF) formalism to emulate an effective DE, and model its anisotropic stress properties through a two-parameter equation that governs its overall amplitude ($g_0$) and transition scale ($c_g$). For the background cosmology, we consider different equations of state to model DE including a constant $w_0$ parameter, and models that provide thawing (CPL) and freezing (nCPL) behaviors. We first constrain these parameters by using the Pantheon, BAO, $H_0$ and CMB Planck data. Then, we analyze the role played by these parameters in the linear PS. In order for the anisotropic stress not to provoke deviations larger than $10\%$ and $5\%$ with respect to the $\Lambda$CDM PS at $k \sim 0.01 \,h/\text{Mpc}$, the parameters have to be in the range $-0.30< g_0 < 0.32$, $0 \leq c_g^2 < 0.01$ and $-0.15 < g_0 < 0.16$, $0 \leq c_g^2 < 0.01$, respectively. Additionally, we compute the leading nonlinear corrections to the PS using standard perturbation theory in real and redshift space, showing that the differences with respect to the $\Lambda$CDM are enhanced, especially for the quadrupole and hexadecapole RSD multipoles.

A hundred years ago, two British expeditions measured the deflection of starlight by the sun's gravitational field, confirming the prediction made by Einstein's General theory of Relativity. One hundred years later many physicists around the world are involved in studying the consequences and use as a research tool, of the deflection of light by gravitational fields, a discipline that today receives the generic name of Gravitational Lensing. The present review aims to commemorate the centenary of Einstein's Eclipse expeditions by presenting a historical perspective of the development and milestones on gravitational light bending, covering from early XIX century speculations, to its current use as an important research tool in astronomy and cosmology.

We develop perturbation theory approaches to model the marked correlation function constructed to up-weight low density regions of the Universe, which might help distinguish modified gravity models from the standard cosmology model based on general relativity. Working within Convolution Lagrangian Perturbation Theory, we obtain that weighted correlation functions are expressible as double convolution integrals that we approximate using a combination of Eulerian and Lagrangian schemes. We find that different approaches agree within 1$\%$ on quasi non-linear scales. Compared with {\it N}-body simulations, the perturbation theory is found to provide accurate predictions for the marked correlation function of dark matter fields, dark matter halos as well as Halo Occupation Distribution galaxies down to $30$ Mpc/h. These analytic approaches help to understand the degeneracy between the mark and the galaxy bias and find a way to maximize the differences among various cosmological models.

We present a formalism to study screening mechanisms in modified theories of gravity via perturbative methods in different cosmological scenarios. We consider Einstein frame posed theories that are recast as Jordan frame theories, where a known formalism is employed, though the resulting non-linearities of the Klein-Gordon equation acquire an explicit coupling between matter and the scalar field, which is not present in Jordan frame theories. The obtained growth functions are then separated in screening and non-screened contributions to facilitate its analysis. This allows us to compare several theoretical models and to recognize patterns which can be used to differentiate models and their screening mechanisms. In particular, we find anti-screening features in the Symmetron model. In opposition, chameleon type theories, both in the Jordan and in the Einstein frame, always present a screening behaviour. Up to third order in perturbation, we find no anti-screening behaviour in theories with a Vainshtein mechanism, such as the DGP and the cubic Galileon.

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.

We present a formalism to compute Lagrangian displacement fields for a wide range of cosmologies in the context of perturbation theory up to third order. We emphasize the case of theories with scale dependent gravitational strengths, such as chameleons, but our formalism can be accommodated to other modified gravity theories. In the non-linear regime two qualitative features arise. One, as is well known, is that nonlinearities lead to a screening of the force mediated by the scalar field. The second is a consequence of the transformation of the Klein-Gordon equation from Eulerian to Lagrangian coordinates, producing frame-lagging terms that are important especially at large scales, and if not considered, the theory does not reduce to the $\Lambda$CDM model in that limit. We apply our formalism to compute the 1-loop power spectrum and the correlation function in $f(R)$ gravity by using different resummation schemes. We further discuss the IR divergences of these formalisms.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

This review describes the discovery of gravitational waves. We recount the journey of predicting and finding those waves, since its beginning in the early twentieth century, their prediction by Einstein in 1916, theoretical and experimental blunders, efforts towards their detection, and finally the subsequent successful discovery.

We study the spherical collapse in the Parametrized Post-Friedmannian (PPF) scheme. We use a general form of the PPF parameter related to the Poisson equation and found the equations to solve that includes a non-trivial fifth force coming from the convolution of the modified gravity term in the k-space. In order to compute a concrete model, we use the parametrization proposed by Bertschinger and Zukin. The equations of the spherical collapse are solved assuming a Gaussian density profile and we show there is no shell crossing before reaching the turn around point. We show that the fifth force does not satisfy the Birkhoff's theorem and introduces different behaviors for the density threshold $\delta_{c}$, which in this case depends on the size and shape of the initial density profile, and therefore one expects a different statistic of the collapsed objects in the universe.

We review a system of autonomous differential equations developed in our previous work [1] describing a flat cosmology filled with a barotropic fluid and a scalar field with a modified kinetic term of the form L=F(X)-V(phi). We analyze the critical points and summarize the conditions to obtain scaling solutions. We consider a set of transformations and show that they leave invariant the equations of motion for the systems in which the scaling solution is obtained, allowing to reduce the number of degrees of freedom.

Using galactic rotation curves, we test a -quantum motivated- gravity model that at large distances modifies the Newtonian potential when spherical symmetry is considered. In this model one adds a Rindler acceleration term to the rotation curves of disk galaxies. Here we consider a standard and a power-law generalization of the Rindler modified Newtonian potential that are hypothesized to play the role of dark matter in galaxies. The new, universal acceleration has to be -phenomenologically- determined. Our galactic model includes the mass of the integrated gas and stars for which we consider a free mass model. We test the model by fitting rotation curves of thirty galaxies that has been employed to test other alternative gravity models. We find that the Rindler parameters do not perform a suitable fit to the rotation curves in comparison to the Burkert dark matter profile, but the models achieve a similar fit as the NFW's profile does. However, the computed parameters of the Rindler gravity show some spread, posing the model to be unable to consistently explain the observed rotation curves.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We review the role of fluids in cosmology by first introducing them in General Relativity and then by applying them to a FRW Universe's model. We describe how relativistic and non-relativistic components evolve in the background dynamics. We also introduce scalar fields to show that they are able to yield an inflationary dynamics at very early times (inflation) and late times (quintessence). Then, we proceed to study the thermodynamical properties of the fluids and, lastly, its perturbed kinematics. We make emphasis in the constrictions of parameters by recent cosmological probes.

In order to unravel the nature of the dark matter (DM) we have proposed a particle-physics motivated model called Bound Dark Matter (BDM) that consist in DM massless particles above a threshold energy Ec that acquire mass below it due to nonperturbative methods. Therefore, the BDM model describes DM particles which are relativistic, hot dark matter (HDM) in the denser (inner) regions of galaxies and describes nonrelativistic, cold dark matter (CDM) where halo density is below rho_c = Ec^4. We test this model by fitting rotation curves from Low Surface Brightness (LSB) galaxies from The HI Nearby Galaxy Survey (THINGS). We use a particular DM cored profile that contains three parameters: a typical scale length (rs) and density (rho_0) of the halo, and a core radius (rc) stemming from the relativistic nature of the BDM model. Since the energy Ec parameterizes the phase transition due to the underlying particle physics model, it is independent on the details of galaxy and/or structure formation and therefore the DM profile parameters rs, rc, Ec are constrained, leaving only two free parameters. Through the results we agree with previous ones implying that cored profiles are preferred over the N-body motivated cuspy profiles. We also compute 2D likelihoods of the BDM parameters rc and Ec for the different galaxies and matter contents, and find an average galaxy core radius rc = 1.48kpc and a transition energy between hot and cold dark matter at Ec = 0.06 eV. The phase transition scale Ec is a new fundamental scale for our DM model well motivated theoretical origin from gauge group dynamics.

Recently a new -quantum motivated- theory of gravity has been proposed that modifies the standard Newtonian potential at large distances when spherical symmetry is considered. Accordingly, Newtonian gravity is altered by adding an extra Rindler acceleration term that has to be phenomenologically determined. Here we consider a standard and a power-law generalization of the Rindler modified Newtonian potential. The new terms in the gravitational potential are hypothesized to play the role of dark matter in galaxies. Our galactic model includes the mass of the integrated gas, and stars for which we consider three stellar mass functions (Kroupa, diet-Salpeter, and free mass model). We test this idea by fitting rotation curves of seventeen low surface brightness galaxies from The HI Nearby Galaxy Survey (THINGS). We find that the Rindler parameters do not perform a suitable fit to the rotation curves in comparison to standard dark matter profiles (Navarro-Frenk-White and Burkert) and, in addition, the computed parameters of the Rindler gravity show a high spread, posing the model as a nonacceptable alternative to dark matter.

Dark matter and dark energy are essential in the description of the late Universe, since at least the epoch of equality. On the other hand, the inflation is also necessary and demands a "dark" component, usually associated to a scalar field that dominated the dynamics and kinematics in the very early Universe. Yet, these three dark components of standard model of cosmology are independent from each other, although there are alternative models that pursue to achieve a triple unification, or at least a double. In the present work we present an update of two models that we have considered in recent years. The first is the "dark fluid" model in which dark matter and dark energy are the same thing, achieving a double unification with specific properties that exactly emulate the standard model of cosmology, given the dark degeneracy that exists in the LCDM model. The second model is given by a single F(X) scalar field Lagrangian, with which one is able to model the whole cosmological dynamics, from inflation to today, representing a triple unification model. We highlight the main properties of these models, as well as we test them against known cosmological probes.

We analyze the dynamical system defined by a universe filled with a barotropic fluid plus a scalar field with modified kinetic term of the form L = F (X) - V (phi). After a suitable choice of variables that allows us to study the phase space of the system we obtain the critical points and their stability. We find that they reduce to the ones defined for the canonical case when F (X) = X. We also study the field energy conditions to have a nonsingular bounce.

We analyze the phase space of a particular unified model of dark matter, dark energy, and inflation that we recently studied in [Phys. Rev. D 83, 063502 (2011)] whose Lagrangian is of the form L(X,phi) = F(X) - V(phi). We show that this model possesses a large set of initial conditions consistent with a successful cosmological model in which an inflationary phase is possible, followed by a matter era to end with dark energy domination. In order to understand the success of the model, we study the general features that unified dark matter (UDM) models should comply and then we analyze some particular models and find their constrictions.

Recently a new particle physics model called Bound Dark Matter (BDM) has been proposed in which dark matter (DM) particles are massless above a threshold energy (Ec) and acquire mass below it due to nonperturbative methods. Therefore, the BDM model describes DM particles which are relativistic, hot dark matter, in the inner regions of galaxies and describes nonrelativistic, cold dark matter, where halo density is below rho_c = Ec^4. To realize this idea in galaxies we use a particular DM cored profile that contains three parameters: a scale length (rs) and density (rho_0) of the halo, and a core radius (rc) stemming from the relativistic nature of the BDM model. We test this model by fitting rotation curves of seventeen Low Surface Brightness galaxies from The HI Nearby Galaxy Survey (THINGS). Since the energy Ec parameterizes the phase transition due to the underlying particle physics model, it is independent on the details of galaxy or structure formation and therefore the DM profile parameters rs, rc, Ec are constrained, leaving only two free parameters. The high spatial and velocity resolution of this sample allows to derive the model parameters through the numerical implementation of the chi^2-goodness-of-fit test to the mass models. We compare the fittings with those of Navarro-Frenk-White (NFW), Burkert, and Pseudo-Isothermal (ISO) profiles. Through the results we conclude that the BDM profile fits better, or equally well, than NFW, Burkert, and ISO profiles and agree with previous results implying that cored profiles are preferred over the N-body motivated cuspy profile. We also compute 2D likelihoods of the BDM parameters rc and Ec for the different galaxies and matter contents, and find an average galaxy core radius rc=300 pc and a transition energy Ec = 0.11 eV when the DM halo is the only component. In Kroupa mass model, we obtain a core rc=1.48 kpc, and energy Ec=0.06 eV.

We analyze the a set of seventeen rotation curves of Low Surface Brightness (LSB) galaxies from the The HI Nearby Galaxy Survey (THINGS) with different mass models to study the core structure and to determine a phase transition energy scale (E_c) between hot and cold dark matter, due to nonperturbative effects in the Bound Dark Matter (BDM) model. Our results agree with previous ones implying the cored profiles are preferred over the N-body motivated cuspy NFW profile. We find an average galactic core radius of r_c = 260 \times 10^{+/- 1.3} pc and a phase transition energy E_c = 0.11\times 10^{+/- 0.46} eV, that is of the same order of magnitude as the sum of the neutrino masses.

We summarize our studies on the determination of the mass of the dark matter halo, based on observations of rotation curves of test particles or of the gravitational lensing. As we show, it is not uncommon that some studies on the nature of dark matter include extra assumptions, some even on the very nature of the dark matter, what we want to determine!, and that bias the studies and the results obtained from the observation and, in some cases, imply an inconsistent system altogether.

This is a brief review of the standard model of cosmology. We first introduce the FRW models and their flat solutions for energy fluids playing an important role in the dynamics at different epochs. We then introduce different cosmological lengths and some of their applications. The later part is dedicated to the physical processes and concepts necessary to understand the early and very early Universe and observations of it.

We study some properties of the dark degeneracy, which is the fact that what we measure in gravitational experiments is the energy momentum tensor of the total dark sector, and any split into components (as in dark matter and dark energy) is arbitrary. In fact, just one dark fluid is necessary to obtain exactly the same cosmological and astrophysical phenomenology as the LCDM model. We work explicitly the first order perturbation theory and show that beyond the linear order the dark degeneracy is preserved under some general assumptions. Then, we construct the dark fluid from a collection of interacting fluids. Finally, we try to break the degeneracy with a general class of couplings to baryonic matter. Nonetheless, we show that these interactions can also be understood in the context of the LCDM model as between dark matter and baryons. For this last investigation we choose two independent parameterizations for the interactions, one inspired in electromagnetism and the other in Chameleon theories. Then, we constrain them with a joint analysis of CMB and Supernovae observational data.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We examine the problem of the collapse and fragmentation of molecular clouds with a Gaussian density distribution with high resolution, double precision numerical simulations using the GADGET-2 code. To describe the thermodynamic properties of the cloud during the collapse -to mimic the rise of temperature predicted by radiative transfer- we use a barotropic equation of state that introduces a critical density to separate the isothermal and adiabatic regimes. We discuss the effects of this critical density in the formation of multiple systems. We confirm the tendency found for Plummer and Gaussian models that if the collapse changes from isothermal to adiabatic at earlier times that occurs for the models with a lower critical density, the collapse is slowed down, and this enhances the fragments' change to survive. However, this effect happens up to a threshold density below which single systems tend to form. On the other hand, by setting a bigger initial perturbation amplitude, the collapse is faster and in some cases a final single object is formed.

We study a unification model for dark energy, dark matter, and inflation with a single scalar field with non canonical kinetic term. In this model the kinetic term of the Lagrangian accounts for the dark matter and dark energy, and at early epochs a quadratic potential accounts for slow roll inflation. The present work is an extension to the work by Bose and Majumdar [1] with a more general kinetic term that was proposed by Chimento in [2]. We demonstrate that the model is viable at the background and linear perturbation levels.

We present a scenario in which a scalar field dark energy is coupled to the trace of the energy momentum tensor of the baryonic matter fields. In the slow-roll regime, this interaction could give rise to the cosmological features of dark matter. We work out the cosmological background solutions and fit the parameters of the model using the Union 2 supernovae data set. Then, we develop the cosmological perturbations up to linear order, and we find that the perturbed variables have an acceptable behavior, in particular the density contrast of baryonic matter grows similar to that in the $\Lambda$CDM model for a suitable choice of the strength parameter of the coupling.

Attractor solutions that give dynamical reasons for dark energy to act like the cosmological constant, or behavior close to it, are interesting possibilities to explain cosmic acceleration. Coupling the scalar field to matter or to gravity enlarges the dynamical behavior; we consider both couplings together, which can ameliorate some problems for each individually. Such theories have also been proposed in a Higgs-like fashion to induce gravity and unify dark energy and dark matter origins. We explore restrictions on such theories due to their dynamical behavior compared to observations of the cosmic expansion. Quartic potentials in particular have viable stability properties and asymptotically approach general relativity.

There are two observations of galaxies that can offer some insight into the nature of the dark matter (DM), namely the rotation curves and the gravitational lensing. While the first one can be studied using the Newtonian limit, the second one is completely relativistic. Each one separately can not determine the nature of DM, but both together give us key information about this open problem. In this work we use a static and spherically symmetric metric to model the DM halo in a galaxy or in a galaxy cluster. The metric contains two free functions, one associated with the distribution of mass and the other one with the gravitational potential. We use galactic, typical rotation curves to univocally determine the kinematics of the halos. We compute separately the mass functions for a perfect fluid and a scalar field, and demonstrate that both models can be fitted to the observations, though with different masses. We then employ lensing to discriminate between these models. This procedure represents a test of models using two measurements: rotation curves and lensing.

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 $\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.

We are considering the cosmological consequences of an induced gravity theory coupled to the minimal standard model of particle physics. The non-minimal coupling parameter between gravity and the Higgs field must then be very large, yielding some new cosmological consequences for the early Universe and new constraints on the Higgs mass. As an outcome, new inflation is only possible for very special initial conditions producing first a short contraction era after which an inflationary expansion automatically follows; a chaotic inflationary scenario is successfully achieved. The contrast of density perturbations required to explain the seed of astronomic structures are obtained for very large values of the Higgs mass ($M_{H} >> G_{F}^{-1/2}$), otherwise the perturbations have a small amplitude; in any case, the spectral index of scalar perturbations agrees with the observed one.

We investigate the cosmological consequences of a theory of induced gravity in which the scalar field is identified with the Higgs field of the first symmetry breaking of a mi\-ni\-mal SU(5) GUT. The mass of the X-boson determines a great value for the coupling constant of gravity-particle physics. Because of this fact, a "slow" rollover dynamics for the Higgs field is not possible in a "new" inflation scenario and, moreover, a contraction era for the scale factor in the early universe exists, \-after\- which inflation follows automatically; "chaotic" inflation is performed without problems. Inflation is successfully achieved due to the relationship among the masses of particle physics at that scale: The Higgs-, X-boson- and Planck-masses. As a result the particle physics pa\-ra\-me\-ter $\lambda$ is not fine-tuned as usual in order to predict acceptable values of re-heating temperature and density and gravitational wave perturbations. Moreover, if the coherent Higgs oscillations didn't decay they could explain the missing mass problem of cosmology.