The Consistency Relations for the Large Scale Structure provide a link between the amplitude of Baryonic Acoustic Oscillations in the squeezed bispectrum (BS) and in the power spectrum (PS). This relation depends on the large scale bias of the considered tracer, $b_\alpha$, and on the growth rate of structures, $f$. Remarkably, originating from basic symmetry principles, this relation is exact and independent on the underlying cosmological model. By analysing data from large volume simulations, both for dark matter and for haloes, we illustrate how BS and PS measurements can be used to extract $b_\alpha$ and $f$ without the need of any theoretical approximation scheme for the computation of the BS and the PS. We show that, combining measurements of the squeezed BS with the quadrupole to monopole ratios for the PS at large scales can successfully break the $b_\alpha -f$ degeneracy. We forecast that this method, applied to a Euclid-like survey, will be able to measure bias, and then the growth rate, at better than $10\%$ level, with no extra assumption.

Consistency Relations (CR) for the Large Scale Structure are exact equalities between correlation functions of different order. These relations descend from the equivalence principle and hold for primordial perturbations generated by single-field models of inflation. They are not affected by nonlinearities and hold also for biased tracers and in redshift space. We show that Baryonic Acoustic Oscillations (BAO) in the bispectrum (BS) in the squeezed limit are suppressed with respect to those in the power spectrum (PS) by a coefficient that depends on the BS configuration and on the bias parameter (and, in redshift space, also on the growth rate). We test these relations using large volume N-body simulations and show that they provide a novel way to measure large scale halo bias and, potentially, the growth rate. Since bias is obtained by comparing two directly observable quantities, the method is free from theoretical uncertainties both on the computational scheme and on the underlying cosmological model.

We present the first application to real data from the BOSS DR12 dataset of the Extractor procedure to determine the acoustic scale imprinted on Baryonic Acoustic Oscillations (BAO). We show that, being largely insensitive to the broadband shape of the Power Spectrum, this procedure requires a lower number of nuisance parameters than those used by the BOSS collaboration, For non-reconstructed data our analysis improves the accuracy on the acoustic scale by about 20 %, while for reconstructed ones we get essentially the same level of accuracy as the BOSS analysis.

We explore the deep ultraviolet (that is, short-distance) limit of the power spectrum (PS) and of the correlation function of a cold dark matter dominated Universe. While for large scales the PS can be written as a double series expansion, in powers of the linear PS and of the wavenumber $k$, we show that, in the opposite limit, it can be expressed via an expansion in powers of the form $1/k^{d+2n}$, where $d$ is the number of spatial dimensions, and $n$ is a non negative integer. The coefficients of the terms of the expansion are nonperturbative in the linear PS, and can be interpreted in terms of the probability density function for the displacement field, evaluated around specific configurations of the latter, that we identify. In the case of the Zel'dovich dynamics, these coefficients can be determined analytically, whereas for the exact dynamics they can be treated as fit, or nuisance, parameters. We confirm our findings with numerical simulations and discuss the necessary steps to match our results to those obtained for larger scales and to actual measurements.

Structure formation in 1+1 dimensions is considered, with emphasis on the effects of shell-crossing. The breakdown of the perturbative expansion beyond shell-crossing is discussed, and it is shown, in a simple example, that the perturbative series can be extended to a transseries including nonperturbative terms. The latter converges to the exact result well beyond the range of validity of perturbation theory. The crucial role of the divergences induced by shell-crossing is discussed. They provide constraints on the structure of the transseries and act as a bridge between the perturbative and the nonperturbative sectors. Then, we show that the dynamics in the deep multistreaming regime is governed by attractors. In the case of simple initial conditions, these attractors coincide with the asymptotic configurations of the adhesion model, but in general they may differ. These results are applied to a cosmological setting, and an algorithm to build the attractor solution starting from the Zel'dovich approximation is developed. Finally, this algorithm is applied to the search of `haloes' and the results are compared with those obtained from the exact dynamical equations.

We study a new procedure to measure the sound horizon scale via Baryonic Acoustic Oscillations (BAO). Instead of fitting the measured power spectrum (PS) to a theoretical model containing the cosmological informations and all the nonlinear effects, we define a procedure to project out (or to "extract") the oscillating component from a given nonlinear PS. We show that the BAO scale extracted in this way is extremely robust and, moreover, can be reproduced by simple theoretical models at any redshift. By using N-body simulations, we discuss the effect of the nonlinear evolution of the matter field, of redshift space distortions and of scale-dependent halo bias, showing that all these effects can be reproduced with sub-percent accuracy. We give a one-parameter theoretical model based on a simple (IR) modification of 1-loop perturbation theory, which reproduces the BAO scale from measurements of halo clustering in redshift space at better than $0.1\%$ level and does not need any external UV input, such as coefficients measured from N-body simulations.

We define a procedure to extract the oscillating part of a given nonlinear Power Spectrum, and derive an equation describing its evolution including the leading effects at all scales. The intermediate scales are taken into account by standard perturbation theory, the long range (IR) displacements are included by using consistency relations, and the effect of small (UV) scales is included via effective coefficients computed in simulations. We show that the UV effects are irrelevant in the evolution of the oscillating part, while they play a crucial role in reproducing the smooth component. Our "extractor" operator can be applied to simulations and real data in order to extract the Baryonic Acoustic Oscillations (BAO) without any fitting function and nuisance parameter. We conclude that the nonlinear evolution of BAO can be accurately reproduced at all scales down to $z=0$ by our fast analytical method, without any need of extra parameters fitted from simulations.

Many of the methods proposed so far to go beyond Standard Perturbation Theory break invariance under time-dependent boosts (denoted here as extended Galilean Invariance, or GI). This gives rise to spurious large scale effects which spoil the small scale predictions of these approximation schemes. By using consistency relations we derive fully non-perturbative constraints that GI imposes on correlation functions. We then introduce a method to quantify the amount of GI breaking of a given scheme, and to correct it by properly tailored counterterms. Finally, we formulate resummation schemes which are manifestly GI, discuss their general features, and implement them in the so called Time-Flow, or TRG, equations.

We study the impact of neutrino masses on the shape and height of the BAO peak of the matter correlation function, both in real and redshift space. In order to describe the nonlinear evolution of the BAO peak we run N-body simulations and compare them with simple analytic formulae. We show that the evolution with redshift of the correlation function and its dependence on the neutrino masses is well reproduced in a simplified version of the Zel'dovich approximation, in which the mode-coupling contribution to the power spectrum is neglected. While in linear theory the BAO peak decreases for increasing neutrino masses, the effect of nonlinear structure formation goes in the opposite direction, since the peak broadening by large scale flows is less effective. As a result of this combined effect, the peak decreases by $\sim 0.6 \%$ for $ \sum m_\nu = 0.15$ eV and increases by $\sim1.2 \%$ for $ \sum m_\nu = 0.3$ eV, with respect to a massless neutrino cosmology with equal value of the other cosmological parameters. We extend our analysis to redshift space and to halos, and confirm the agreement between simulations and the analytic formulae. We argue that all analytical approaches having the Zel'dovich propagator in their lowest order approximation should give comparable performances, irrespectively to their formulation in Lagrangian or in Eulerian space.

Standard cosmological perturbation theory (SPT) for the Large Scale Structure (LSS) of the Universe fails at small scales (UV) due to strong nonlinearities and to multistreaming effects. In Pietroni et al. 2011 a new framework was proposed in which the large scales (IR) are treated perturbatively while the information on the UV, mainly small scale velocity dispersion, is obtained by nonlinear methods like N-body simulations. Here we develop this approach, showing that it is possible to reproduce the fully nonlinear power spectrum (PS) by combining a simple (and fast) 1-loop computation for the IR scales and the measurement of a single, dominant, correlator from N-body simulations for the UV ones. We measure this correlator for a suite of seven different cosmologies, and we show that its inclusion in our perturbation scheme reproduces the fully non-linear PS with percent level accuracy, for wave numbers up to $k\sim 0.4\, h~{\rm Mpc^{-1}}$ down to $z=0$. We then show that, once this correlator has been measured in a given cosmology, there is no need to run a new simulation for a different cosmology in the suite. Indeed, by rescaling this correlator by a proper function computable in SPT, the reconstruction procedure works also for the other cosmologies and for all redshifts, with comparable accuracy. Finally, we clarify the relation of this approach to the Effective Field Theory methods recently proposed in the LSS context.

Time modulations at per mil level have been reported to take place in the decay constant of about 15 nuclei with period of one year (most cases) but also of about one month or one day. In this paper we give the results of the activity measurement of a 40K source and a 232Th one. The two experiments have been done at the Gran Sasso Laboratory during a period of about 500 days, above ground (40K) and underground (232Th) with a target sensitivity of a few parts over 10^5. We also give the results of the activity measurement at the time of the X-class solar flares which took place in May 2013. Briefly, our measurements do not show any evidence of unexpected time dependence in the decay rate of 40K and 232Th.

We consider cosmologies in which a dark-energy scalar field interacts with cold dark matter. The growth of perturbations is followed beyond the linear level by means of the time-renormalization-group method, which is extended to describe a multi-component matter sector. Even in the absence of the extra interaction, a scale-dependent bias is generated as a consequence of the different initial conditions for baryons and dark matter after decoupling. The effect is enhanced significantly by the extra coupling and can be at the 2-3 percent level in the range of scales of baryonic acoustic oscillations. We compare our results with N-body simulations, finding very good agreement.

We present fully nonlinear consistency relations for the squeezed bispectrum of Large Scale Structure. These relations hold when the matter component of the Universe is composed of one or more species, and generalize those obtained in Peloso, Pietroni '13 and Kehagias, Riotto '13 in the single species case. The multi-species relations apply to the standard dark matter + baryons scenario, as well as to the case in which some of the fields are auxiliary quantities describing a particular population, such as dark matter halos or a specific galaxy class. If a large scale velocity bias exists between the different populations new terms appear in the consistency relations with respect to the single species case. As an illustration, we discuss two physical cases in which such a velocity bias can exist: (1) a new long range scalar force in the dark matter sector (resulting in a violation of the equivalence principle in the dark matter-baryon system), and (2) the distribution of dark matter halos relative to that of the underlying dark matter field.

We study the gravitational clustering of spherically symmetric overdensities and the statistics of the resulting dark matter halos in the "symmetron model", in which a new long range force is mediated by a $Z_2$ symmetric scalar field. Depending on the initial radius of the overdensity, we identify two distinct regimes: for small initial radii the symmetron mediated force affects the spherical collapse at all redshifts; for initial radii larger than some critical size this force vanishes before collapse because of the symmetron screening mechanism. In both cases overdensities collapse earlier than in the $\Lambda$CDM and statistically tend to form more massive dark matter halos. Regarding the halo-mass function of these objects, we observe order one departures from standard $\Lambda$CDM predictions. The formalism developed here can be easily applied to other models where fifth-forces participate to the dynamics of the gravitational collapse.

We discuss the constraints imposed on the nonlinear evolution of the Large Scale Structure (LSS) of the universe by galilean invariance, the symmetry relevant on subhorizon scales. Using Ward identities associated to the invariance, we derive fully nonlinear consistency relations between statistical correlators of the density and velocity perturbations, such as the power spectrum and the bispectrum. These relations are valid up to O (f_{NL}^2) corrections. We then show that most of the semi-analytic methods proposed so far to resum the perturbative expansion of the LSS dynamics fail to fulfill the constraints imposed by galilean invariance, and are therefore susceptible to non-physical infrared effects. Finally, we identify and discuss a nonperturbative semi-analytical scheme which is manifestly galilean invariant at any order of its expansion.

A new computational scheme for the nonlinear cosmological matter power spectrum (PS) is presented. Our method is based on evolution equations in time, which can be cast in a form extremely convenient for fast numerical evaluations. A nonlinear PS is obtained in a time comparable to that needed for a simple 1-loop computation, and the numerical implementation is very simple. Our results agree with N-body simulations at the percent level in the BAO range of scales, and at the few-percent level up to $k ~ 1$ h/Mpc at $z >= 0.5$, thereby opening the possibility of applying this tool to scales interesting for weak lensing. We clarify the approximations inherent to this approach as well as its relations to previous ones, such as the Time Renormalization Group, and the multi-point propagator expansion. We discuss possible lines of improvements of the method and its intrinsic limitations by multi streaming at small scales and low redshifts.

We consider a dark energy fluid with arbitrary sound speed and equation of state and discuss the effect of its clustering on the cold dark matter distribution at the non-linear level. We write the continuity, Euler and Poisson equations for the system in the Newtonian approximation. Then, using the time renormalization group method to resum perturbative corrections at all orders, we compute the total clustering power spectrum and matter power spectrum. At the linear level, a sound speed of dark energy different from that of light modifies the power spectrum on observationally interesting scales, such as those relevant for baryonic acoustic oscillations. We show that the effect of varying the sound speed of dark energy on the non-linear corrections to the matter power spectrum is below the per cent level, and therefore these corrections can be well modelled by their counterpart in cosmological scenarios with smooth dark energy. We also show that the non-linear effects on the matter growth index can be as large as 10-15 per cent for small scales.

Semi-analytical methods, based on Eulerian perturbation theory, are a promising tool to follow the time evolution of cosmological perturbations at small redshifts and at mildly nonlinear scales. All these schemes are based on two approximations: the existence of a smoothing scale and the single-stream approximation, where velocity dispersion of the dark matter fluid, as well as higher moments of the particle distributions, are neglected. Despite being widely recognized, these two assumptions are, in principle, incompatible, since any finite smoothing scale gives rise to velocity dispersion and higher moments at larger scales. We describe a new approach to perturbation theory, where the Vlasov and fluid equations are derived in presence of a finite coarse-graining scale: this allows a clear separation between long and short distance modes and leads to a hybrid approach where the former are treated perturbatively and the effect of the latter is encoded in external source terms for velocity, velocity dispersion, and all the higher order moments, which can be computed from N-body simulations. We apply the coarse-grained perturbation theory to the computation of the power spectrum and the cross-spectrum between density and velocity dispersion, and compare the results with N-body simulations, finding good agreement.

One of the nicest results in cosmological perturbation theory is the analytical resummaton of the leading corrections at large momentum, which was obtained by Crocce and Scoccimarro for the propagator. Using an exact evolution equation, we generalize this result, by showing that a class of next-to-leading corrections can also be resummed at all orders in perturbation theory. The new corrections modify the propagator by a few percent in the Baryonic Acoustic Oscillation range of scales, and therefore cannot be neglected in resummation schemes aiming at an accuracy compatible with future generation galaxy surveys. Similar tools can be employed to derive improved approximations for the Power Spectrum.

We address the issue of the cosmological bias between matter and galaxy distributions, looking at dark-matter haloes as a first step to characterize galaxy clustering. Starting from the linear density field at high redshift, we follow the centre of mass trajectory of the material that will form each halo at late times (proto-halo). We adopt a fluid-like description for the evolution of perturbations in the proto-halo distribution, which is coupled to the matter density field via gravity. We present analytical solutions for the density and velocity fields, in the context of renormalized perturbation theory. We start from the linear solution, then compute one-loop corrections for the propagator and the power spectrum. Finally we analytically resum the propagator and we use a suitable extension of the time-renormalization-group method (Pietroni 2008) to resum the power spectrum. For halo masses M<10^{14} Msol/h our results at z=0 are in good agreement with N-body simulations. Our model is able to predict the halo-matter cross spectrum with an accuracy of 5 per cent up to k = 0.1 h/Mpc approaching the requirements of future galaxy redshift surveys.

We apply the time-renormalization group approach to study the effect of primordial non-Gaussianities in the non-linear evolution of cosmological dark matter density perturbations. This method improves the standard perturbation approach by solving renormalization group-like equations governing the dynamics of gravitational instability. The primordial bispectra constructed from the dark matter density contrast and the velocity fields represent initial conditions for the renormalization group flow. We consider local, equilateral and folded shapes for the initial non-Gaussianity and analyze as well the case in which the non-linear parameter f_{NL} parametrizing the strength of the non-Gaussianity depends on the momenta in Fourier space through a power-law relation, the so-called running non-Gaussianity. For the local model of non-Gaussianity we compare our findings for the power-spectrum with those of recent N-body simulations and find that they accurately fit the N-body data up to wave-numbers k \sim 0.25 h/Mpc at z=0. We also present predictions for the (reduced) matter bispectra for the various shapes of non-Gaussianity.

Future large scale structure observations are expected to be sensitive to small neutrino masses, of the order of 0.05 eV or more. However, forecasts are based on the assumption that by the time at which these datasets will be available, the non-linear spectrum in presence of neutrino mass will be predicted with an accuracy at least equal to the neutrino mass effect itself, i.e. about 3%. Motivated by these considerations, we present the computation of the non-linear power spectrum of LambdaCDM models in the presence of massive neutrinos using the Renormalization Group time-flow approach, which amounts to a resummation of perturbative corrections to the matter power spectrum to all orders. We compare our results with those obtained with other methods, i.e. linear theory, one-loop perturbation theory and N-body simulations and show that the time-RG method improves the one-loop method in fitting the N-body data, especially in determining the suppression of the matter power spectrum when neutrino are massive with respect to the linear power spectrum.

Nonlinear effects are crucial in order to compute the cosmological matter power spectrum to the accuracy required by future generation surveys. Here, a new approach is presented, in which the power spectrum, the bispectrum and higher order correlations, are obtained -- at any redshift and for any momentum scale -- by integrating a system of differential equations. The method is similar to the familiar BBGKY hierarchy. Truncating at the level of the trispectrum, the solution of the equations corresponds to the summation of an infinite class of perturbative corrections. Compared to other resummation frameworks, the scheme discussed here is particularly suited to cosmologies other than LambdaCDM, such as those based on modifications of gravity and those containing massive neutrinos. As a first application, we compute the Baryonic Acoustic Oscillation feature of the power spectrum, and compare the results with perturbation theory, the halo model, and N-body simulations. The density-velocity and velocity-velocity power spectra are also computed, showing that they are much less contaminated by nonlinearities than the density-density one. The approach can be seen as a particular formulation of the renormalization group, in which time is the flow parameter.

Inflationary observables, like the power spectrum, computed at one- and higher-order loop level seem to be plagued by large infra-red corrections. In this short note, we point out that these large infra-red corrections appear only in quantities which are not directly observable. This is in agreement with general expectations concerning infra-red effects.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

Considering the general Lagrangian of k-essence models, we study and classify them through variables connected to the fluid equation of state parameter w_\kappa. This allows to find solutions around which the scalar field describes a mixture of dark matter and cosmological constant-like dark energy, an example being the purely kinetic model proposed by Scherrer. Making the stronger assumption that the scalar field Lagrangian is exactly constant along solutions of the equation of motion, we find a general class of k-essence models whose classical trajectories directly describe a unified dark matter/dark energy (cosmological constant) fluid. While the simplest case of a scalar field with canonical kinetic term unavoidably leads to an effective sound speed c_s=1, thereby inhibiting the growth of matter inhomogeneities, more general non-canonical k-essence models allow for the possibility that c_s << 1 whenever matter dominates.

Scalar-Tensor theories of gravity can be formulated in different frames, most notably, the Einstein and the Jordan one. While some debate still persists in the literature on the physical status of the different frames, a frame transformation in Scalar-Tensor theories amounts to a local redefinition of the metric, and then should not affect physical results. We analyze the issue in a cosmological context. In particular, we define all the relevant observables (redshift, distances, cross-sections, ...) in terms of frame-independent quantities. Then, we give a frame-independent formulation of the Boltzmann equation, and outline its use in relevant examples such as particle freeze-out and the evolution of the CMB photon distribution function. Finally, we derive the gravitational equations for the frame-independent quantities at first order in perturbation theory. From a practical point of view, the present approach allows the simultaneous implementation of the good aspects of the two frames in a clear and straightforward way.

Renormalization Group (RG) techniques have been successfully employed in quantum field theory and statistical physics. Here we apply RG methods to study the non-linear stages of structure formation in the Universe. Exact equations for the power spectrum, the bispectrum, and all higher order correlation functions can be derived for any underlying cosmological model. A remarkable feature of the RG flow is the emergence of an intrinsic UV cutoff, due to dark matter velocity dispersion, which improves the convergence of the equations at small scales. As a consequence, the method is able to follow the non-linear evolution of the power-spectrum down to zero redshift and to length-scales where perturbation theory fails. Our predictions accurately fit the results of $N$-body simulations in reproducing the ``Baryon Acoustic Oscillations'' features of the power-spectrum, which will be accurately measured in future galaxy surveys and will provide a probe to distinguish among different dark energy models.

A semi-analytic approach to the computation of the non-linear power-spectrum of dark matter density fluctuations is proposed. The method is based on the Renormalization Group technique and can be applied to any underlying cosmological model. Our prediction on the baryonic acoustic oscillations in a LambdaCDM model accurately fits the results of N-body simulations down to zero redshift, where perturbation theory fails.

A host of dark energy models and non-standard cosmologies predict an enhanced Hubble rate in the early Universe: perfectly viable models, which satisfy Big Bang Nucleosynthesis (BBN), cosmic microwave background and general relativity tests, may nevertheless lead to enhancements of the Hubble rate up to many orders of magnitude. In this paper we show that strong bounds on the pre-BBN evolution of the Universe may be derived, under the assumption that dark matter is a thermal relic, by combining the dark matter relic density bound with constraints coming from the production of cosmic-ray antiprotons by dark matter annihilation in the Galaxy. The limits we derive can be sizable and apply to the Hubble rate around the temperature of dark matter decoupling. For dark matter masses lighter than 100 GeV, the bound on the Hubble-rate enhancement ranges from a factor of a few to a factor of 30, depending on the actual cosmological model, while for a mass of 500 GeV the bound falls in the range 50-500. Uncertainties in the derivation of the bounds and situations where the bounds become looser are discussed. We finally discuss how these limits apply to some specific realizations of non-standard cosmologies: a scalar-tensor gravity model, kination models and a Randall-Sundrum D-brane model.

The two most popular candidates for dark energy, i.e. a cosmological constant and quintessence, are very difficult to distinguish observationally, mostly because the quintessence field does not have sizable fluctuations. We study a scalar field model for dark energy in which the scalar field is invariant under reflection symmetry, phi -> -phi. Under general assumptions, there is a phase transition at late times (z < 0.5). Before the phase transition, the field behaves as a cosmological constant. After the phase transition, a time-dependent scalar condensate forms, the field couples with dark matter and develops sizable perturbations tracking those of dark matter. The background cosmological evolution is in agreement with existing observations, but might be distinguished from that of a cosmological constant by future Supernovae surveys. The growth of cosmological perturbations carries the imprint of the phase transition, however a non-linear approach has to be developed in order to study it quantitatively.

A scalar field interacting differently with dark matter and baryons may explain why their ratio is of order unity today. We provide three working examples, checking them against the observations of CMB, Large Scale Structure, supernovae Ia, and post-newtonian tests of gravity. Such a scenario could make life much easier for supersymmetric dark matter candidates.

Scalar-tensor theories of gravity provide a consistent framework to accommodate an ultra-light quintessence scalar field. While the equivalence principle is respected by construction, deviations from General Relativity and standard cosmology may show up at nucleosynthesis, CMB, and solar system tests of gravity. After imposing all the bounds coming from these observations, we consider the expansion rate of the universe at WIMP decoupling, showing that it can lead to an enhancement of the dark matter relic density up to few orders of magnitude with respect to the standard case. This effect can have an impact on supersymmetric candidates for dark matter.

In Scalar-Tensor theories of gravity, the expansion rate of the universe in the past may differ considerably from the standard one. After imposing the constraints coming from nucleosynthesis, CMB, type Ia supernovae, and solar system tests of General Relativity, we consider the expansion rate of the universe at WIMP decoupling, showing that it can lead to an enhancement of the dark matter relic density of some orders of magnitude with respect to the standard case. This effect may have a deep impact on most popular candidates for dark matter, as the supersymmetric neutralino.

Can local fluctuations of a ``Quintessence'' scalar field play a dynamical role in the gravitational clustering and cosmic structure formation process? We address this question in the general framework of scalar-tensor theories of gravity. Non-linear energy density perturbations, both in the scalar field and matter component, and linear metric perturbations are accounted for in the perturbed Einstein's equations. We derive the Newtonian limit of the relevant equations for clustering in scalar-tensor cosmologies. We then specialize to non-linear perturbations of the ``Extended Quintessence'' model of Dark Energy; in such a model, a non-minimally coupled scalar field is thought to be responsible for driving the present accelerated phase of the Universe expansion. The interplay between Dark Energy and Dark Matter is displayed in the equations governing the growth of structure in the Universe.

Non-linear gravitational clustering in a universe dominated by dark energy, modelled by a `quintessence' scalar field, and cold dark matter with space-time varying mass is studied. Models of this type, where the variable mass is induced by dependence on the scalar field, as suggested by string theory or extra-dimensions, have been proposed as a viable solution of the coincidence problem. A general framework for the study of the non-linear phases of structure formation in scalar field cosmologies is provided, starting from a general relativistic treatment of the combined dark matter-dark energy system. As a first application, the mildly non-linear evolution of dark matter perturbations is obtained by a straightforward extension of the Zel'dovich approximation. We argue that the dark energy fluctuations may play an active role in cosmological structure formation if the scalar field effective potential develops a temporary spinodal instability during the evolution.

Brane world models with `large' extra dimensions with radii in the r_l ~ 0.01- 0.1 mm range and smaller ones at r_s < (1 TeV)^(-1) have the potential to solve the cosmic coincidence problem, i.e. the apparently fortuitous equality between dark matter and dark energy components today. The main ingredient is the assumption of a stabilization mechanism fixing the total volume of the compact submanifold, but allowing for shape deformations. The latter are associated with phenomenologically safe ultra-light scalar fields. Bulk fields Casimir energy naturally plays the role of dark energy, which decreases in time because of expanding r_l. Stable Kaluza Klein states may play the role of dark matter with increasing, O(1/r_s), mass. The cosmological equations exhibit attractor solutions in which the global equation of state is negative, the ratio between dark energy and dark matter is constant and the observed value of the ratio is obtained for two large extra dimensions. Experimental searches of large extra dimensions should take into account that, due to the strong coupling between dark matter and radii dynamics, the size of the large extra dimensions inside the galactic halo may be smaller than the average value.

It is a puzzle why the densities of dark matter and dark energy are nearly equal today when they scale so differently during the expansion of the universe. This conundrum may be solved if there is a coupling between the two dark sectors. In this paper we assume that dark matter is made of cold relics with masses depending exponentially on the scalar field associated to dark energy. Since the dynamics of the system is dominated by an attractor solution, the dark matter particle mass is forced to change with time as to ensure that the ratio between the energy densities of dark matter and dark energy become a constant at late times and one readily realizes that the present-day dark matter abundance is not very sensitive to its value when dark matter particles decouple from the thermal bath. We show that the dependence of the present abundance of cold dark matter on the parameters of the model differs drastically from the familiar results where no connection between dark energy and dark matter is present. In particular, we analyze the case in which the cold dark matter particle is the lightest supersymmetric particle.

The so called `cosmic coincidence' problem seems to suggest a coupling between dark energy and dark matter. A possibility could be given by dark matter particles with masses depending exponentially on the scalar field associated to dark energy. If the latter also has an exponential potential, attractor solutions exist with constant Omega_{DM} and Omega_{DE} and negative effective equation of state. Due to the very low mass of the dark energy field, its coupling to baryons should be strongly suppressed. Here we present a natural realization of this scenario in a model with large extra dimensions.

Matsumoto and Yoshimura have recently argued that the number density of heavy particles in a thermal bath is not necessarily Boltzmann-suppressed for T << M, as power law corrections may emerge at higher orders in perturbation theory. This fact might have important implications on the determination of WIMP relic densities. On the other hand, the definition of number densities in a interacting theory is not a straightforward procedure. It usually requires renormalization of composite operators and operator mixing, which obscure the physical interpretation of the computed thermal average. We propose a new definition for the thermal average of a composite operator, which does not require any new renormalization counterterm and is thus free from such ambiguities. Applying this definition to the model of Matsumoto and Yoshimura we find that it gives number densities which are Boltzmann-suppressed at any order in perturbation theory. We discuss also heavy particles which are unstable already at T=0, showing that power law corrections do in general emerge in this case.

Scalar fields with inverse power-law effective potentials may provide a negative pressure component to the energy density of the universe today, as required by cosmological observations. In order to be cosmologically relevant today, the scalar field should have a mass $m_\phi = O(10^{-33} {\mathrm eV})$, thus potentially inducing sizable violations of the equivalence principle and space-time variations of the coupling constants. Scalar-tensor theories of gravity provide a framework for accommodating phenomenologically acceptable ultra-light scalar fields. We discuss non-minimally coupled scalar-tensor theories in which the scalar-matter coupling is a dynamical quantity. Two attractor mechanisms are operative at the same time: one towards the tracker solution, which accounts for the accelerated expansion of the Universe, and one towards general relativity, which makes the ultra-light scalar field phenomenologically safe today. As in usual tracker-field models, the late-time behavior is largely independent on the initial conditions. Strong distortions in the cosmic microwave background anisotropy spectra as well as in the matter power spectrum are expected.

Recent data point in the direction of a cosmological constant dominated universe. We investigate the role of supersymmetric QCD with N_f < N_c as a possible candidate for dynamical cosmological constant (``quintessence''). When N_f > 1, the multiscalar dynamics is fully taken into account, showing that a certain degree of flavor symmetry in the initial conditions is necessary in order to attain the required late-time ``tracking'' attractor behavior. A possible solution to the problem of long-range interactions of gravitational strength is proposed. Finally we show how, as a consequence of the interaction with the Standard Model fields, the early cosmological evolution of the scalar fields may be dramatically altered, and the approach to the tracking behavior may become much more efficient.