The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We present a new parameterization of quintessence potentials for dark energy based directly upon the dynamical properties of the equations of motion. Such parameterization arises naturally once the equations of motion are written as a dynamical system in terms of properly defined polar variables. We have identified two different classes of parameters, and we dubbed them as dynamical and passive parameters. The dynamical parameters appear explicitly in the equations of motion, but the passive parameters play just a secondary role in their solutions. The new approach is applied to the so-called thawing potentials and it is argued that only three dynamical parameters are sufficient to capture the evolution of the quintessence fields at late times. This work reconfirms the arbitrariness of the quintessence potentials as the recent observational data fail to constrain the dynamical parameters.

We propose a new parametrization of the background equations of motion corresponding to (canonical) single-field models of inflation, which allows a better understanding of the general properties of the solutions and of the corresponding predictions in the inflationary observables. Based on the tools of dynamical systems, the method suggests that inflation comes in two flavors: power-law and de Sitter. Power-law inflation seems to occur for a restricted type of potentials, whereas de Sitter inflation has a much broader applicability. We also show a general perturbative method, by means of series expansion, to solve the new equations of motion around the critical point of the de Sitter type, and how the method can be used for arbitrary models of de Sitter inflation. It is then argued that for the latter there are two general classes of inflationary solutions, given in terms of the behavior of the tensor-to-scalar ratio as a function of the number $N$ of $e$-folds before the end of inflation: $r \sim N^{-1}$ (Class I), or $r \sim N^{-2}$ (Class II). We give some examples of the two classes in terms of known scalar field potentials, and compare their general predictions with constraints obtained from observations.

We derive non-relativistic equations of motion for the formation of cosmological structure in a Scalar Field Dark Matter (SFDM) model corresponding to a complex scalar field endowed with a quadratic scalar potential. Starting with the full equations of motion written in the Newtonian gauge of scalar perturbations, we separate out the fields involved into relativistic and non-relativistic parts, and find the equations of motion for the latter that can be used to build up the full solution. One important assumption will also be that the SFDM field is in the regime of fast oscillations, under which its behavior is exactly that of cold dark matter. The resultant equations are quite similar to the Schr\"odinger-Poisson system of Newtonian boson stars plus relativistic leftovers. We exploit that similarity to show how to simulate, with minimum numerical effort, the formation of cosmological structure in SFDM models and others alike, and ultimately prove their viability as complete dark matter models.

In this paper we review the production of primordial black holes (PBHs) during preheating after a chaotic inflationary model. All relevant equations of motion are solved numerically in a modified version of HLattice, and we then calculate the mass variance to determine structure formation during preheating. It is found that production of PBHs can be a generic result of the model, even though the results seem to be sensitive to the values of the smoothing scale. We consider a constraint for overproduction of PBHs that could uncover some stress between inflation-preheating models and observations.

We revise the conditions for the physical viability of a cosmological model in which dark matter has bulk viscosity and also interacts with dark energy. We have also included radiation and baryonic matter components; all matter components are represented by perfect fluids, except the dark matter, that is treated as an imperfect fluid. We impose upon the model the condition of a complete cosmological dynamics that results in an either null or negative bulk viscosity, but the latter also disagrees with the Local Second Law of Thermodynamics. The model is also compared with cosmological observations at different redshifts: type Ia supernova, the shift parameter of CMB, the acoustic peak of BAO, and the Hubble parameter H(z). In general, observations consistently point out to a negative value of the bulk viscous coefficient, and in overall the fitting procedure shows no preference for the model over the standard LCDM model.

A massive, self-interacting scalar field has been considered as a possible candidate for the dark matter in the universe. We present an observational constraint to the model arising from strong lensing observations in galaxies. The result points to a discrepancy in the properties of scalar field dark matter halos for dwarf and lens galaxies, mainly because halo parameters are directly related to physical quantities in the model. This is an important indication that it becomes necessary to have a better understanding of halo evolution in scalar field dark matter models, where the presence of baryons can play an important role.

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.

Using the dynamical system approach, we describe the general dynamics of cosmological scalar fields in terms of critical points and heteroclinic lines. It is found that critical points describe the initial and final states of the scalar field dynamics, but that heteroclinic lines which give a more complete description of the evolution in between the critical points. In particular, the heteroclinic line that departs from the (saddle) critical point of perfect fluid-domination is the representative path in phase space of quintessence fields that may be viable dark energy candidates. We also discuss the attractor properties of the heteroclinic lines, and their importance for the description of thawing and freezing fields.

In this work a supersymmetric cosmological model is analyzed in which we consider a general superfield action of a homogeneous scalar field supermultiplet interacting with the scale factor in a supersymmetric FRW model. There appear fermionic superpartners associated with both the scale factor and the scalar field, and classical equations of motion are obtained from the super-Wheeler-DeWitt equation through the usual WKB method. The resulting supersymmetric Einstein-Klein-Gordon equations contain extra radiation and stiff matter terms, and we study their solutions in flat space for different scalar field potentials. The solutions are compared to the standard case, in particular those corresponding to the exponential potential, and their implications for the dynamics of the early Universe are discussed in turn.

We study the consequences of further modification of $f(R,R_{\mu\nu} R^{\mu\nu}, R_{\mu\nu\sigma\rho} R^{\mu\nu\sigma\rho}) / f(R)$-theories by means of the Dirac-Born-Infeld deformation procedure, which amounts to the replacement of $f$ by $\lambda(\sqrt{1+2f/\lambda}-1)$ (the free parameter $\lambda$ fixes an additional energy scale). We pay special attention to the definition of masses of the linearized propagating degrees of freedom since these are important to judge about the stability of the linearization around vacuum background spaces. In this context we discuss the subtleties associated with expanding $f(R,R_{\mu\nu} R^{\mu\nu},R_{\mu\nu\sigma\rho} R^{\mu\nu\sigma\rho})$-Lagrangians around maximally symmetric spaces of constant curvature, as well as with equivalence of the linearized Lagrangian to a scalar-tensor theory. Investigation of the consequences of applying the Dirac-Born-Infeld strategy to further modify quadratic theories, on the stability of de Sitter vacuum, as well as its impact on the cosmological dynamics, is the main concern of this paper.

The dynamics of a cosmological model fueled by scalar field dark matter with a cosh-like potential plus a cosmological constant is investigated in detail. It is revealed that the late-time attractor is always the de Sitter solution, and that, depending on the values of the free parameters, the oscillating solution of the scalar field -- modeling cold dark matter -- mediates between some early stage (say, the radiation-dominated solution) and the accelerating de Sitter attractor.

We review the attractor properties of the simplest chaotic model of inflation, in which a minimally coupled scalar field is endowed with a quadratic scalar potential. The equations of motion in a flat Friedmann-Robertson-Walker universe are written as an autonomous system of equations, and the solutions of physical interest appear as critical points. This new formalism is then applied to the study of inflation dynamics, in which we can go beyond the known slow-roll formalism of inflation.

Dark energy and dark matter are only indirectly measured via their gravitational effects. It is possible that there is an exchange of energy within the dark sector, and this offers an interesting alternative approach to the coincidence problem. We consider two broad classes of interacting models where the energy exchange is a linear combination of the dark sector densities. The first class has been previously investigated, but we define new variables and find a new exact solution, which allows for a more direct, transparent and comprehensive analysis. The second class has not been investigated in general form before. We give general conditions on the parameters in both classes to avoid unphysical behavior (such as negative energy densities).

We propose a hybrid inflation model with a complex waterfall field which contains an interaction term that breaks the U(1) global symmetry associated to the waterfall field charge. We show that the asymmetric evolution of the real and imaginary parts of the complex field during the phase transition at the end of inflation translates into a charge asymmetry. The latter strongly depends on the vev of the waterfall field, which is well constrained by diverse cosmological observations.

Equilibrium configurations for a self-gravitating scalar field with self-interaction are constructed. The corresponding Schr\"odinger-Poisson (SP) system is solved using finite differences assuming spherical symmetry. It is shown that equilibrium configurations of the SP system are late-time attractor solutions for initially quite arbitrary density profiles, which relax and virialize through the emission of scalar field bursts; a process dubbed gravitational cooling. Among other potential applications, these results indicate that scalar field dark matter models (in its different flavors) tolerate the introduction of a self-interaction term in the SP equations. This study can be useful in exploring models in which dark matter in galaxies is not point-like.

Phantom energy can be visualized as a scalar field with a (non-canonical) negative kinetic energy term. We use the dynamical system formalism to study the attractor behavior of a cosmological model containing a phantom scalar field $\phi$ endowed with an exponential potential of the form $V(\phi)=V_0 \exp(-\lambda \kappa \phi)$, and a perfect fluid with constant equation of state $\gamma$; the latter can be of the phantom type too. As in the canonical case, three characteristic solutions can be identified. The scaling solution exists but is either unstable or of no physical interest. Thus, there are only two stable critical points which are of physical interest, corresponding to the perfect fluid and scalar field dominated solutions, respectively. The most interesting case arises for $0> \gamma > -\lambda^2/3$, which allows the coexistence of the three solutions. The main features of each solution are discussed in turn.

Using numerical techniques, we study the collapse of a scalar field configuration in the Newtonian limit of the spherically symmetric Einstein--Klein--Gordon (EKG) system, which results in the so called Schr\"odinger--Newton (SN) set of equations. We present the numerical code developed to evolve the SN system and topics related, like equilibrium configurations and boundary conditions. Also, we analyze the evolution of different initial configurations and the physical quantities associated to them. In particular, we readdress the issue of the gravitational cooling mechanism for Newtonian systems and find that all systems settle down onto a 0--node equilibrium configuration.

Dark matter in the universe seems to be one of the most important puzzles science has to face in this moment. In this essay we point out that dark matter could be a spin-0 fundamental interaction of Nature rather than a simple particle. From this hypothesis follows that dark matter behaves just as standard cold dark matter at cosmological level while still in good agreement with observations at galactic scales. This new interaction could be one of the scalar fields predicted by higher-dimensional theories.

The curvaton was introduced recently as a distinct inflationary mechanism for generating adiabatic density perturbations. Implicit in that scenario is that the curvaton offers a new mechanism for reheating after inflation, as it is a form of energy density not diluted by the inflationary expansion. We consider curvaton reheating in the context of a braneworld inflation model, {\em steep inflation}, which features a novel use of the braneworld to give a new mechanism for ending inflation. The original steep inflation model featured reheating by gravitational particle production, but the inefficiency of that process brings observational difficulties. We demonstrate here that the phenomenology of steep inflation is much improved by curvaton reheating.

In this letter, we develop a Newtonian approach to the collapse of galaxy fluctuations of scalar field dark matter under initial conditions inferred from simple assumptions. The full relativistic system, the so called Einstein-Klein-Gordon, is reduced to the Schr\"odinger-Newton one in the weak field limit. The scaling symmetries of the SN equations are exploited to track the non-linear collapse of single scalar matter fluctuations. The results can be applied to both real and complex scalar fields.

Ultra-light scalar fields provide an interesting alternative to WIMPS as halo dark matter. In this paper we consider the effect of embedding a supermassive black hole within such a halo, and estimate the absorption probability and the accretion rate of dark matter onto the black hole. We show that the accretion rate would be small over the lifetime of a typical halo, and hence that supermassive central black holes can coexist with scalar field halos.

In this paper, we study some interesting properties of a spherically symmetric oscillating soliton star made of a real time-dependent scalar field which is called an oscillaton. The known final configuration of an oscillaton consists of a stationary stage in which the scalar field and the metric coefficients oscillate in time if the scalar potential is quadratic. The differential equations that arise in the simplest approximation, that of coherent scalar oscillations, are presented for a quadratic scalar potential. This allows us to take a closer look at the interesting properties of these oscillating objects. The leading terms of the solutions considering a quartic and a cosh scalar potentials are worked in the so called stationary limit procedure. This procedure reveals the form in which oscillatons and boson stars may be related and useful information about oscillatons is obtained from the known results of boson stars. Oscillatons could compete with boson stars as interesting astrophysical objects, since they would be predicted by scalar field dark matter models.

A unified model is developed within the context of the braneworld paradigm, where a single scalar field can act as both the inflaton field in the very early universe and also as strong, self-interacting dark matter in the post-inflationary universe. Reheating proceeds due to the overproduction and subsequent evaporation of primordial black holes. Observational constraints, most notably from gravitational waves, are satisfied if the probability of PBH formation is sufficiently high.

We present a general description of the scalar field dark matter (SFDM) hypothesis in the cosmological context. The scenario of structure formation under such a hypothesis is based on Jeans instabilities of fluctuations of the scalar field. It is shown that it is possible to form stable long lived objects consisting of a wide range of typical galactic masses around $10^{12}M_{\odot}$ once the parameters of the effective theory are fixed with the cosmological constraints. The energy density at the origin of such an object is smooth as it should.

We present a scenario for galaxy formation based on the hypothesis of scalar field dark matter. We interpret galaxy formation through the collapse of a scalar field fluctuation. We find that a cosh potential for the self-interaction of the scalar field provides a reasonable scenario for galactic formation, which is in agreement with cosmological observations and phenomenological studies in galaxies.

In recent papers we have proposed that the dark matter of the Universe could be from scalar field origin. In this letter, we find that if the scale of renormalization of the model is of order of the Planck Mass, then a scalar field $\Phi $ endowed with the scalar potential $V=V_{o}[\cosh {(\lambda \sqrt{\kappa_{0}}\Phi)}-1] can be a reliable model for dark matter in galaxies. The predicted scattering cross section fits the value required for self-interacting dark matter and the additional degree of freedom of the theory is of order of hundreds of TeV.

This work is a review of the last results of research on the Scalar Field Dark Matter model of the Universe at cosmological and at galactic level. We present the complete solution to the scalar field cosmological scenario in which the dark matter is modeled by a scalar field $\Phi$ with the scalar potential $V(\Phi)=V_{0}(cosh {(\lambda \sqrt{\kappa_{0}}\Phi)}-1)$ and the dark energy is modeled by a scalar field $\Psi$, endowed with the scalar potential $\tilde{V}(\Psi)= \tilde{V_{0}}(\sinh{(\alpha \sqrt{\kappa_{0}}\Psi)})^{\beta}$, which together compose the 95% of the total matter energy in the Universe. The model presents successfully deals with the up to date cosmological observations, and is a good candidate to treat the dark matter problem at the galactic level.

We present the complete solution to a 95% scalar field cosmological model in which the dark matter is modeled by a scalar field $\Phi$ with the scalar potential $V(\Phi)=V_{o}[ \cosh {(\lambda \sqrt{\kappa_{o}}\Phi)}-1]$ and the dark energy is modeled by a scalar field $\Psi$, endowed with the scalar potential $\tilde{V}(\Psi)=\tilde{V_{o}}[ \sinh {(\alpha \sqrt{\kappa_{o}}\Psi)}] ^{\beta}$. This model has only two free parameters, $\lambda$ and the equation of state $\omega_{\Psi}$. With these potentials, the fine tuning and the cosmic coincidence problems are ameliorated for both dark matter and dark energy and the models agrees with astronomical observations. For the scalar dark matter, we clarify the meaning of a scalar Jeans lenght and then the model predicts a suppression of the Mass Power Spectrum for small scales having a wave number $k > k_{min,\Phi}$, where $k_{min,\Phi} \simeq 4.5 h {\rm Mpc}^{-1}$ for $\lambda \simeq 20.28$. This last fact could help to explain the dearth of dwarf galaxies and the smoothness of galaxy core halos. From this, all parameters of the scalar dark matter potential are completely determined. The dark matter consists of an ultra-light particle, whose mass is $m_{\Phi} \simeq 1.1\times 10^{-23} {\rm eV}$ and all the success of the standard cold dark matter model is recovered. This implies that a scalar field could also be a good candidate as the dark matter of the Universe.

Continuing with previous works, we present a cosmological model in which dark matter and dark energy are modeled by scalar fields $\Phi $ and $\Psi$, respectively, endowed with the scalar potentials $V(\Phi)=V_{o}[ \cosh {(\lambda \sqrt{\kappa_{o}}\Phi)}-1] $ and $\tilde{V}(\Psi)=\tilde{V_{o}}[ \sinh {(\alpha \sqrt{\kappa_{o}}\Psi)}] ^{\beta}$. This model contains 95% of scalar field. We obtain that the scalar dark matter mass is $m_{\Phi}\sim 10^{-26}eV.$ The solution obtained allows us to recover the success of the standard CDM. The implications on the formation of structure are reviewed. We obtain that the minimal cutoff radio for this model is $r_{c}\sim 1.2 kpc.$

In this paper we propose a quintessence model with the potential $V(\Phi )=V_{o}[ \sinh {(\alpha \sqrt{\kappa_{o}}\Delta \Phi})] ^{\beta}$, which asymptotic behavior corresponds to an inverse power-law potential at early times and to an exponential one at late times. We demonstrate that this is a tracker solution and that it could have driven the Universe into its current inflationary stage. The exact solutions and the description for a complete evolution of the Universe are also given. We compare such model with the current cosmological observations.

We investigate the hypothesis that the scalar field is the dark matter and the dark energy in the Cosmos, wich comprises about 95% of the matter of the Universe. We show that this hypothesis explains quite well the recent observations on type Ia supernovae.