We compute the one-loop density power spectrum including Newtonian and relativistic contributions, as well as the primordial non-Gaussianity contributions from $f_{\rm NL}$ and $g_{\rm NL}$ in the local configuration. To this end we take solutions to the Einstein equations in the long-wavelength approximation and provide expressions for the matter density perturbation at second and third order. These solutions have shown to be complementary to the usual Newtonian cosmological perturbations. We confirm a sub-dominant effect from pure relativistic terms, manifested at scales dominated by cosmic variance, but find that a sizable effect of order one comes from $g_{\rm NL}$ values allowed by Planck-2018 constraints, manifested at scales probed by forthcoming galaxy surveys like DESI and Euclid. As a complement, we present the matter bispectrum at the tree-level including the mentioned contributions.

We present a sample of 81 Gamma-Ray Bursts (GRBs) observed by Fermi-GBM for which we compute the distance moduli and use them to constrain effective dark energy models. To overcome the circularity problem affecting the use of GRBs as distance indicators, we calibrate the Amati relation of our sample by employing a cosmology-independent technique. Specifically, the latest observational Hubble parameter data are used to approximate the cosmic expansion through a Bezier parametric curve. We subsequently obtain the distance moduli of the GRBs and include it in a suite of recent cosmological (Planck Compressed 2018 and 2012 BOSS release of BAO data) and local (Pantheon SNIa) observations of the expansion history to compute Bayesian posterior constraints for the standard cosmological model $\Lambda$CDM, $\omega$CDM, and for the CPL parametrization. We highlight the advantages that our dataset and method represent over other recent GRB data.

We study a model of inflation based on $\mathcal{N}=1$ supergravity essentially depending on one effective parameter. Under a field transformation we show that this model turns out to be equivalent to a previously studied supergravity model known to be ruled out with the original choice of the parameter. Such parameter measures the slope of the potential at observable scales. Through a Bayesian parameter estimation, it is shown how this model is compatible with recent CMB temperature measurements by {\it Planck 2018} giving rise to a simple, viable, single field model of inflation. The tensor to scalar ratio constraint is found to be $r_{0.002}<0.065$ with negative running. We discuss how observables are invariant under the field transformation which leaves unaltered the slow-roll parameters. As a consequence the use of one presentation of the model or its field-transformed version is purely a matter of convenience.

Next generation surveys will be capable of determining cosmological parameters beyond percent level. To match this precision, theoretical descriptions should look beyond the linear perturbations to approximate the observables in large scale structure. A quantity of interest is the Number density of galaxies detected by our instruments. This has been focus of interest recently, and several efforts have been made to explain relativistic effects theoretically, thereby testing the full theory. However, the results at nonlinear level from previous works are in disagreement. We present a new and independent approach to computing the relativistic galaxy number counts to second order in cosmological perturbation theory. We derive analytical expressions for the full second order relativistic observed redshift, for the angular diameter distance and for the volume spanned by a survey. Finally, we compare our results with previous works which compute the general distance-redshift relation, finding that our result is in agreement at linear order.

Dark energy models with tracker properties have gained attention due to the large range of initial conditions leading to the current value of the dark energy density parameter. A well-motivated family of these models are the so-called $\alpha$-attractors, which show the late time behavior of a cosmological constant. In the present paper we perform a model-selection analysis of a variety of $\alpha$-attractor potentials in comparison with a non-flat $\Lambda$CDM model. Specifically, we compute the Bayes Factor for the L-Model, the Oscillatory Tracker Model, the Recliner Model, and the Starobinsky Model, while considering the non-flat $\Lambda$CDM as the base model. Each model is tested through a Bayesian analysis using observations relevant to the current accelerated expansion: we employ the latest SNe Ia data, combined with cosmic clocks, the latest BOSS release of BAO data, and the Planck Compressed 2018 data. The produced Markov Chains for each model are further compared through a Bayesian evidence analysis. From the latter we conclude that the Oscillatory Tracker Model is preferred by data (even if weakly) over the non-flat $\Lambda$CDM model. Our results also suggest at the L-model is the least favoured version of the $\alpha$-attractor models considered.

We examine the fully relativistic evolution of cosmic voids constituted by baryons and cold dark matter (CDM), represented by two non-comoving dust sources in a $\Lambda$CDM background. For this purpose, we consider numerical solutions of Einstein's field equations in a fluid-flow representation adapted to spherical symmetry and multiple components. We present a simple example that explores the frame-dependence of the local expansion and the Hubble flow for this mixture of two dusts, revealing that the relative velocity between the sources yields a significantly different evolution in comparison with that of the two sources in a common 4-velocity (which reduces to a Lemaitre-Tolman-Bondi model). In particular, significant modifications arise for the density contrast depth and void size, as well as in the amplitude of the surrounding over-densities. We show that an adequate model of a frame-dependent evolution that incorporates initial conditions from peculiar velocities and large-scale density contrast observations may contribute to understand the discrepancy between the local value of $H_0$ and that inferred from the CMB.

We show how to constrain inflationary models and reheating by using mixed constraints. In particular we study the physics of the reheating phase after inflation from observational constraints to the inflationary stage. We show that it is possible to determine $\omega$, the equation of state during reheating, by using the reported values of the spectral index and the {\it full} number of $e$-folds $N(n_s,\omega)= N_H(n_s)+N_{re}(n_s,w)\approx 60$, which includes the accelerated expansion and the reheating phase. We show that the reheating number of $e$-folds $N_{re}$ is quite sensitive to this equation of state. Requiring $N_{re}>0$ and a sensible value for the thermalization scale $T_{re}$, demands in general a reheating phase with $\omega\neq 0$. We exemplify the constraints with two particular examples: We show how the Starobinsky model allows only large values of $T_{re}$ if the reheating phase is dominated by dust ($w =0$), and if Primordial Black Hole production is subdominant. For the case of $N=1$ Supergravity inflation, the extra parameter of the potential provides the necessary freedom to afford lower-scale thermalization in a dust-like reheating phase and yet our method serves to determine the rest of the observable parameters.

We examine the gravitational collapse and black hole formation of multiple non--spherical configurations constructed from Szekeres dust models with positive spatial curvature that smoothly match to a Schwarzschild exterior. These configurations are made of an almost spherical central core region surrounded by a network of "pancake-like" overdensities and voids with spatial positions prescribed through standard initial conditions. We show that a full collapse into a focusing singularity, without shell crossings appearing before the formation of an apparent horizon, is not possible unless the full configuration becomes exactly or almost spherical. Seeking for black hole formation, we demand that shell crossings are covered by the apparent horizon. This requires very special fine-tuned initial conditions that impose very strong and unrealistic constraints on the total black hole mass and full collapse time. As a consequence, non-spherical non-rotating dust sources cannot furnish even minimally realistic toy models of black hole formation at astrophysical scales: demanding realistic collapse time scales yields huge unrealistic black hole masses, while simulations of typical astrophysical black hole masses collapse in unrealistically small times. We note, however, that the resulting time--mass constraint is compatible with early Universe models of primordial black hole formation, suitable in early dust-like environments. Finally, we argue that the shell crossings appearing when non-spherical dust structures collapse are an indicator that such structures do not form galactic mass black holes but virialise into stable stationary objects.

We constrain several models in Tachyonic Inflation derived from the large-$N$ formalism by considering theoretical aspects as well as the latest observational data. On the theoretical side, we assess the field range of our models by means of the excursion of the equivalent canonical field. On the observational side, we employ BK14+PLANCK+BAO data to perform a parameter estimation analysis as well as a Bayesian model selection to distinguish the most favoured models among all four classes here presented. We observe that the original potential $V \propto \textrm{sech}(T)$ is strongly disfavoured by observations with respect to a reference model with flat priors on inflationary observables. This realisation of Tachyon inflation also presents a large field range which may demand further quantum corrections. We also provide examples of potentials derived from the polynomial and the perturbative classes which are both statistically favoured and theoretically acceptable.

Recently we have studied in great detail a model of Hybrid Natural Inflation (HNI) by constructing two simple effective field theories. These two versions of the model allow inflationary energy scales as small as the electroweak scale in one of them or as large as the Grand Unification scale in the other, therefore covering the whole range of possible energy scales. In any case the inflationary sector of the model is of the form $V(\phi)=V_0 \left(1+a \cos(\phi/f)\right)$ where $0\leq a<1$ and the end of inflation is triggered by an independent waterfall field. One interesting characteristic of this model is that the slow-roll parameter $\epsilon(\phi)$ is a non-monotonic function of $\phi$ presenting a {\it maximum} close to the inflection point of the potential. Because the scalar spectrum $\mathcal{P}_s(k)$ of density fluctuations when written in terms of the potential is inversely proportional to $\epsilon(\phi)$ we find that $\mathcal{P}_s(k)$ presents a {\it minimum} at $\phi_{min}$. The origin of the HNI potential can be traced to a symmetry breaking phenomenon occurring at some energy scale $f$ which gives rise to a (massless) Goldstone boson. Non-perturbative physics at some temperature $T<f$ might occur which provides a potential (and a small mass) to the originally massless boson to become the inflaton (a pseudo-Nambu-Goldstone boson). Thus the inflaton energy scale $\Delta$ is bounded by the symmetry breaking scale, $\Delta\equiv V_H^{1/4} <f.$ To have such a well defined origin and hierarchy of scales in inflationary models is not common. We use this property of HNI to determine bounds for the inflationary energy scale $\Delta$ and for the tensor-to-scalar ratio $r$.

Oscillating scalar fields are useful to model a variety of matter components in the universe. One or more scalar fields participate in the reheating process after inflation, while at much lower energies scalar fields are robust dark matter candidates. Pertaining structure formation in these models, it is well known that inhomogeneities of the Klein-Gordon field are unstable above the characteristic De Broglie wavelength. In this paper we show that such instability implies the existence of a threshold amplitude for the collapse of primordial fluctuations. We use this threshold to correctly predict the cut--off scale of the matter power spectrum in the scalar field dark matter model. Furthermore, for a Klein-Gordon field during reheating we show that this same threshold allows for abundant production of structure (oscillons but not necessarily black holes). Looking at the production of Primordial Black Holes (PBHs) in this scenario we note that the sphericity condition yields a much lower probability of PBH formation at the end of inflation. Remarkably, even after meeting such stringent condition, we find that PBHs may be overproduced during reheating. We finally constrain the epochs at which an oscillating Klein-Gordon field could dominate the early universe.

We study the differences and equivalences between the non-perturbative description of the evolution of cosmic structure furnished by the Szekeres dust models (a non-spherical exact solution of Einstein's equations) and the dynamics of Cosmological Perturbation Theory (CPT) for dust sources in a $\Lambda$CDM background. We show how the dynamics of Szekeres models can be described by evolution equations given in terms of "exact fluctuations" that identically reduce (at all orders) to evolution equations of CPT in the comoving isochronous gauge. We explicitly show how Szekeres linearised exact fluctuations are specific (deterministic) realisations of standard linear perturbations of CPT given as random fields but, as opposed to the latter perturbations, they can be evolved exactly into the full non-linear regime. We prove two important results: (i) the conservation of the curvature perturbation (at all scales) also holds for the appropriate approximation of the exact Szekeres fluctuations in a $\Lambda$CDM background, and (ii) the different collapse morphologies of Szekeres models yields, at nonlinear order, different functional forms for the growth factor that follows from the study of redshift space distortions. The metric based potentials used in linear CPT are computed in terms of the parameters of the linearised Szekeres models, thus allowing us to relate our results to linear CPT results in other gauges. We believe that these results provide a solid starting stage to examine the role of non-perturbative General Relativity in current cosmological research.

We take a pragmatic, model independent approach to single field slow-roll canonical inflation by imposing conditions, not on the potential, but on the slow-roll parameter $\epsilon(\phi)$ and its derivatives $\epsilon^{\prime }(\phi)$ and $\epsilon^{\prime\prime }(\phi)$, thereby extracting general conditions on the tensor-to-scalar ratio $r$ and the running $n_{sk}$ at $\phi_{H}$ where the perturbations are produced, some $50$ $-$ $60$ $e$-folds before the end of inflation. We find quite generally that for models where $\epsilon(\phi)$ develops a maximum, a relatively large $r$ is most likely accompanied by a positive running while a negligible tensor-to-scalar ratio implies negative running. The definitive answer, however, is given in terms of the slow-roll parameter $\xi_2(\phi)$. To accommodate a large tensor-to-scalar ratio that meets the limiting values allowed by the Planck data, we study a non-monotonic $\epsilon(\phi)$ decreasing during most part of inflation. Since at $\phi_{H}$ the slow-roll parameter $\epsilon(\phi)$ is increasing, we thus require that $\epsilon(\phi)$ develops a maximum for $\phi > \phi_{H}$ after which $\epsilon(\phi)$ decrease to small values where most $e$-folds are produced. The end of inflation might occur trough a hybrid mechanism and a small field excursion $\Delta\phi_e\equiv |\phi_H-\phi_e |$ is obtained with a sufficiently thin profile for $\epsilon(\phi)$ which, however, should not conflict with the second slow-roll parameter $\eta(\phi)$. As a consequence of this analysis we find bounds for $\Delta \phi_e$, $r_H$ and for the scalar spectral index $n_{sH}$. Finally we provide examples where these considerations are explicitly realised.

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3--dimensional networks of cold dark matter structures (over--densities and/or density voids) undergoing "pancake" collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities of structures that evolved, from linear initial data at the last scattering surface, to fully non--linear 10--20 Mpc. scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained -- but fully relativistic non--linear and non--perturbative -- description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.

We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.

We study tachyon inflation within the large-$N$ formalism, which takes a prescription for the small Hubble flow slow--roll parameter $\epsilon_1$ as a function of the large number of $e$-folds $N$. This leads to a classification of models through their behaviour at large $N$. In addition to the perturbative $N$ class, we introduce the polynomial and exponential classes for the $\epsilon_1$ parameter. With this formalism we reconstruct a large number of potentials used previously in the literature for Tachyon Inflation. We also obtain new families of potentials form the polynomial class. We characterize the realizations of Tachyon Inflation by computing the usual cosmological observables up to second order in the Hubble flow slow--roll parameters. This allows us to look at observable differences between tachyon and canonical single field inflation. The analysis of observables in light of the Planck 2015 data shows the viability of some of these models, mostly for certain realization of the polynomial and exponential classes.

We examine the relation between the dynamics of Lema\^{\i}tre-Tolman-Bondi (LTB) dust models (with and without $\Lambda$) and the dynamics of dust perturbations in two of the more familiar formalisms used in cosmology: the metric based Cosmological Perturbation Theory (CPT) and the Covariant Gauge Invariant (GIC) perturbations. For this purpose we recast the evolution of LTB models in terms of a covariant and gauge invariant formalism of local and non-local "exact fluctuations " on a Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) background defined by suitable averages of covariant scalars. We examine the properties of these fluctuations, which can be defined for a confined comoving domain or for an asymptotic domain extending to whole time slices. In particular, the non-local density fluctuation provides a covariant and precise definition for the notion of the "density contrast ". We show that in their linear regime these LTB exact fluctuations (local and non-local) are fully equivalent to the conventional cosmological perturbations in the synchronous-comoving gauge of CPT and to GIC perturbations. As an immediate consequence, we show the time-invariance of the spatial curvature perturbation in a simple form. The present work may provide important theoretical connections between the exact and perturbative (linear or no-linear) approach to the dynamics of dust sources in General Relativity.

We study the Hybrid Natural Inflation (HNI) model and some of its realisations in the light of recent CMB observations, mainly Planck temperature and WMAP-9 polarization, and compare with the recent release of BICEP2 dataset. The inflationary sector of HNI is essentially given by the potential $V(\phi) = V_0(1+a\cos (\frac{\phi}{f} ) )$, where $a$ is a positive constant smaller or equal to one and $f$ is the scale of (pseudo Nambu-Goldstone) symmetry breaking. We show that to describe the HNI model realisations we only need two observables; the spectral index $n_s$, the tensor-to-scalar ratio, and a free parameter in the amplitude of the cosine function $a$. We find that in order to make the HNI model compatible with the BICEP2 observations, we require a large positive running of the spectra. We find that this could over-produce primordial black holes in the most consistent case of the model. This situation could be aleviated if, as recently argued, the BICEP2 data do not correspond to primordial gravitational waves.

We show how the non-linearity of general relativity generates a characteristic non-Gaussian signal in cosmological large-scale structure that we calculate at all perturbative orders in a large scale limit. Newtonian gravity and general relativity provide complementary theoretical frameworks for modelling large-scale structure in $\Lambda$CDM cosmology; a relativistic approach is essential to determine initial conditions which can then be used in Newtonian simulations studying the non-linear evolution of the matter density. Most inflationary models in the very early universe predict an almost Gaussian distribution for the primordial metric perturbation, $\zeta$. However, we argue that it is the Ricci curvature of comoving-orthogonal spatial hypersurfaces, $R$, that drives structure formation at large scales. We show how the non-linear relation between the spatial curvature, $R$, and the metric perturbation, $\zeta$, translates into a specific non-Gaussian contribution to the initial comoving matter density that we calculate for the simple case of an initially Gaussian $\zeta$. Our analysis shows the non-linear signature of Einstein's gravity in large-scale structure.

We analyse Hybrid Natural Inflation in view of the recent results for the tensor index reported by BICEP2. We find that it predicts a large running of the scalar spectrum which is potentially detectable by large scale structure and $21\, \mathrm{cm} $ observations. The running of the running is also relatively large becoming close to $10^{-2}$. Along the way, we find general consistency relations at which observables are subject if the slow-roll approximation is imposed. Failure to satisfy these equations by the values obtained for the observables in surveys would be a failure of the slow-roll approximation itself.

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation, we compute the probability distribution of curvature fluctuations, paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter space.

The goal of the present paper is to set initial conditions for structure formation at non-linear order, consistent with general relativity, while also allowing for primordial non-Gaussianity. We use the non-linear continuity and Raychauduri equations, which together with the non-linear energy constraint determine the evolution of the matter density fluctuation in general relativity. We solve this equations at first and second order in a perturbative expansion, recovering and extending previous results derived in the matter-dominated limit and in the Newtonian regime. We present a second-order solution for the comoving density contrast in a $\Lambda$CDM universe, identifying non-linear contributions coming from the Newtonian growing mode, primordial non-Gaussianity and intrinsic non-Gaussianity, due to the essential non-linearity of the relativistic constraint equations. We discuss the application of these results to initial conditions in N-body simulations, showing that relativistic corrections mimic a non-zero non-linear parameter $f_{\rm NL}$.

We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic "second order Poisson equation" is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field.

In this article we perform a second order perturbation analysis of the gravitational metric theory of gravity $ f(\chi) = \chi^{3/2} $ developed by Bernal et al. (2011). We show that the theory accounts in detail for two observational facts: (1) the phenomenology of flattened rotation curves associated to the Tully-Fisher relation observed in spiral galaxies, and (2) the details of observations of gravitational lensing in galaxies and groups of galaxies, without the need of any dark matter. We show how all dynamical observations on flat rotation curves and gravitational lensing can be synthesised in terms of the empirically required metric coefficients of any metric theory of gravity. We construct the corresponding metric components for the theory presented at second order in perturbation, which are shown to be perfectly compatible with the empirically derived ones. It is also shown that under the theory being presented, in order to obtain a complete full agreement with the observational results, a specific signature of Riemann's tensor has to be chosen. This signature corresponds to the one most widely used nowadays in relativity theory. Also, a computational program, the MEXICAS (Metric EXtended-gravity Incorporated through a Computer Algebraic System) code, developed for its usage in the Computer Algebraic System (CAS) Maxima for working out perturbations on a metric theory of gravity, is presented and made publicly available.

Currently, most of the numerical simulations of structure formation use Newtonian gravity. When modelling pressureless dark matter, or `dust', this approach gives the correct results for scales much smaller than the cosmological horizon, but for scenarios in which the fluid has pressure this is no longer the case. In this article, we present the correspondence of perturbations in Newtonian and cosmological perturbation theory, showing exact mathematical equivalence for pressureless matter, and giving the relativistic corrections for matter with pressure. As an example, we study the case of scalar field dark matter which features non-zero pressure perturbations. We discuss some problems which may arise when evolving the perturbations in this model with Newtonian numerical simulations and with CMB Boltzmann codes.

Using dimensional analysis techniques we present an extension of Newton's gravitational theory built under the assumption that Milgrom's acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton's gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

We study the possibility of reheating the universe through the evaporation of primordial black holes created at the end of inflation. This is shown to allow for the unification of inflation and dark matter under the dynamics of a single scalar field. We determine the necessary conditions to recover the standard Big Bang by the time of nucleosynthesis after reheating through black holes.

We show that the Modified Newtonian Dynamics (MOND) regime can be fully recovered as the weak-field limit of a particular theory of gravity formulated in the metric approach. This is possible when Milgrom's acceleration constant is taken as a fundamental quantity which couples to the theory in a very consistent manner. As a consequence, the scale invariance of the gravitational interaction is naturally broken. In this sense, Newtonian gravity is the weak-field limit of general relativity and MOND is the weak-field limit of that particular extended theory of gravity.

We show that the Modified Newtonian Dynamics (MOND) regime can be fully recovered as the weak-field limit of a particular theory of gravity formulated in the metric approach. This is possible when Milgrom's acceleration constant is taken as a fundamental quantity which couples to the theory in a very consistent manner. As a consequence, the scale invariance of the gravitational interaction is naturally broken. In this sense, Newtonian gravity is the weak-field limit of general relativity and MOND is the weak-field limit of that particular extended theory of gravity. We also prove that a Noether's symmetry approach to the problem yields a conserved quantity coherent with this relativistic MONDian extension.

In the pursuit of a general formulation for a modified gravitational theory at the non-relativistic level and as an alternative to the dark matter hypothesis, we construct a model valid over a wide variety of astrophysical scales. Through the inclusion of Milgrom's acceleration constant into a gravitational theory, we show that very general formulas can be constructed for the acceleration felt by a particle. Dimensional analysis shows that this inclusion naturally leads to the appearance of a mass-length scale in gravity, breaking its scale invariance. A particular form of the modified gravitational force is constructed and tested for consistency with observations over a wide range of astrophysical environments, from solar system to extragalactic scales. We show that over any limited range of physical parameters, which define a specific class of astrophysical objects, the dispersion velocity of a system must be a power law of its mass and size. These powers appear linked together through a natural constraint relation of the theory. This yields a generalised gravitational equilibrium relation valid for all astrophysical systems. A general scheme for treating spherical symmetrical density distributions is presented, which in particular shows that the fundamental plane of elliptical galaxies, the Newtonian virial equilibrium, the Tully-Fisher and the Faber-Jackson relations, as well as the scalings observed in local dwarf spheroidal galaxies, are nothing but particular cases of that relation when applied to the appropriate mass-length scales. We discuss the implications of this approach for a modified theory of gravity and emphasise the advantages of working with the force, instead of altering Newton's second law of motion, in the formulation of a gravitational theory.

This thesis begins with a study of the origin of cosmological fluctuations with special attention to those cases in which the non-Gaussian correlation functions are large. The analysis shows that perturbations from an almost massless auxiliary field generically produce large values of the non-linear parameter f_NL. The effects of including non-Gaussian correlation functions in the statistics of cosmological structure are explored by constructing a non-Gaussian probability distribution function (PDF). Such PDF is derived for the comoving curvature perturbation from first principles in the context of quantum field theory, with n-point correlation functions as the only input. The non-Gaussian PDF is then used to explore two important problems in the physics of primordial black holes (PBHs): First, to compute non-Gaussian corrections to the number of PBHs generated from the primordial curvature fluctuations. The second application concerns new cosmological observables. The formation of PBHs is known to depend on two main physical characteristics: the strength of the gravitational field produced by the initial curvature inhomogeneity and the pressure gradient at the edge of the curvature configuration. We account for the probability of finding these configurations by using two parameters: The amplitude of the inhomogeneity and its second radial derivative, evaluated at the centre of the configuration. The implications of the derived probability for the fraction of mass in the universe in the form of PBHs are discussed.

In this paper we derive the probability of the radial profiles of spherically symmetric inhomogeneities in order to provide an improved estimation of the number density of primordial black holes (PBHs). We demonstrate that the probability of PBH formation depends sensitively on the radial profile of the initial configuration. We do this by characterising this profile with two parameters chosen heuristically: the amplitude of the inhomogeneity and the second radial derivative, both evaluated at the centre of the configuration. We calculate the joint probability of initial cosmological inhomogeneities as a function of these two parameters and then find a correspondence between these parameters and those used in numerical computations of PBH formation. Finally, we extend our heuristic study to evaluate the probability of PBH formation taking into account for the first time the radial profile of curvature inhomogeneities.

We present a ballistic description of the formation and propagation of the working surface of a relativistic jet. Using simple laws of conservation of mass and linear momentum at the working surface, we obtain a full description of the working surface flow parametrised by the initial velocity and mass injection rate. This allows us to compute analytically the energy release at any time in the working surface. We compare this model with the results obtained numerically through a new hydrodynamical code applied to the propagation of a relativistic fluid in one dimension in order to test the limits of our study. Finally, we compare our analytical results with observed light curves of five long gamma ray bursts and show that our model is in very good agreement with observations using simple periodic variations of the injected velocity profiles. This simple method allows us to recover initial mass discharge and energy output ejected during the burst.

This paper explores the consequences of non-Gaussian cosmological perturbations for the formation of primordial black holes (PBHs). A non-Gaussian probability distribution function (PDF) of curvature perturbations is presented with an explicit contribution from the three-point correlation function to linear order. The consequences of this non-Gaussian PDF for the large perturbations that form PBHs are then studied. Using the observational limits for the non-Gaussian parameter $f_{NL}$, new bounds to the mean amplitude of curvature perturbations are derived in the range of scales relevant for PBH formation.

We show how to obtain the probability density function for the amplitude of the curvature perturbation, R, produced during an epoch of slow-roll, single-field inflation, working directly from n-point correlation functions of R. These n-point functions are the usual output of quantum field theory calculations, and as a result we bypass approximate statistical arguments based on the central limit theorem. Our method can be extended to deal with arbitrary forms of non-Gaussianity, appearing at any order in the n-point hierarchy. We compute the probability density for the total smoothed perturbation within a Hubble volume, \epsilon, and for the spectrum of \epsilon. When only the two-point function is retained, exact Gaussian statistics are recovered. When the three-point function is taken into account, we compute explicitly the leading slow-roll correction to the Gaussian result.

Self-similar solutions to the problem of a strong imploding relativistic shock wave are calculated. These solutions represent the relativistic generalisation of the Newtonian Gouderley-Landau-Stanyukovich problem of a strong imploding spherical shock wave converging to a centre. The solutions are found assuming that the pre-shocked flow has a uniform density, and are accurate for sufficiently large times after the formation of the shock wave.

We present a ballistic description of the propagation of the working surface of a relativistic jet. Using simple laws of conservation of mass and linear momentum at the working surface, we obtain a full description of the jet flow parametrised by the initial velocity and mass injection. This analysis will soon be applied to particular cases of time-dependent injection of mass and velocity into the jet.

The radio/optical alignment effect for small powerful radio galaxies has been shown to be produced by shock waves formed by the interaction of the head of the jet and/or cocoon with clouds embedded in the interstellar/intergalactic medium. We present here preliminary results of analytical and numerical solutions that have been made to account for the production of implosive shock waves induced by embedding cold clouds in the radio lobe of expanding powerful radio sources.