The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We consider a scenario where dark energy and baryons are dynamically coupled without any energy transfer. In this scenario, the background cosmology is unaffected and, at the perturbations level, the coupling only appears through the corresponding Euler equations of dark energy and baryons. We then explore some phenomenological consequences of this scenario and their signatures in several cosmological observables. In particular, we show its ability to suppress the growth of cosmic structures. We also constrain the parameters of the model with cosmological data and show that an interaction of dark energy with baryons on cosmological scales is mildly favoured.

We consider cosmological models where dark matter is charged under a dark Abelian gauge field. This new interaction is repulsive and competes with gravity on large scales and in the dynamics of galaxies and clusters. We focus on non-linear models of dark electrodynamics where the effects of the new force are screened within a K-mouflage radius that helps avoiding traditional constraints on charged dark matter models. We discuss the background cosmology of these models in a Newtonian approach and show the equivalence with relativistic Lema\^itre models where an inhomogeneous pressure due to the electrostatic interaction is present. In particular, we find that dark matter shells of different radii evolve differently as they exit their K-mouflage radii at different times, resulting in a breaking of the initial comoving evolution. In the large time regime, the background cosmology is described by comoving but inhomogeneous model with a reduced gravitational Newton constant and a negative curvature originating from the electrostatic pressure. Baryons do not directly feel the electrostatic interaction, but are influenced by the inhomogeneous matter distribution induced by the electric force. We find that shells of smaller radii evolve faster than the outer shells which feel the repulsive interaction earlier. This mimics the discrepancy between the large scale Hubble rate and the local one. Similarly, as galaxies and clusters are not screened by the new interaction, large scale global flows would result from the existence of the new dark electromagnetic interaction.

Gravitational wave (GW) oscillations occur whenever there are additional tensor modes interacting with the perturbations of the metric coupled to matter. These extra modes can arise from new spin-2 fields (as in e.g. bigravity theories) or from non-trivial realisations of the cosmological principle induced by background vector fields with internal symmetries (e.g. Yang-Mills, gaugids or multi-Proca). We develop a general cosmological framework to study such novel features due to oscillations. The evolution of the two tensor modes is described by a linear system of coupled second order differential equations exhibiting friction, velocity, chirality and mass mixing. We follow appropriate schemes to obtain approximate solutions for the evolution of both modes and show the corresponding phenomenology for different mixings. Observational signatures include modulations of the wave-form, oscillations of the GW luminosity distance, anomalous GW speed and chirality. We discuss the prospects of observing these effects with present and future GW observatories such as LIGO/VIRGO and LISA.

We perform the general relativistic stability analysis against radial oscillations of unpaired quark stars obtained using the equation of state for cold quark matter from perturbative QCD, the only free parameter being the renormalization scale. This approach consistently incorporates the effects of interactions and includes a built-in estimate of the inherent systematic uncertainties in the evaluation of the equation of state. We also take into account the constraints imposed by the recent gravitational wave event GW 170817 to the compact star masses and radii, and restrict their vibrational spectrum.

We study the effects of heavy quarks on the equation of state for cold and dense quark matter obtained from perturbative QCD, yielding observables parametrized only by the renormalization scale. We investigate the behavior of charm quark matter under the constraints of $\beta$-equilibrium and electric charge neutrality in a region of densities where perturbative QCD is in principle much more reliable. We also analyze the stability of charm stars, and find that such quark stars are unstable under radial oscillations.

The universal character of the gravitational interaction provided by the equivalence principle motivates a geometrical description of gravity. The standard formulation of General Relativity \`a la Einstein attributes gravity to the spacetime curvature, to which we have grown accustomed. However, this perception has masked the fact that two alternative, though equivalent, formulations of General Relativity in flat spacetimes exist, where gravity can be fully ascribed either to torsion or to non-metricity. The latter allows a simpler geometrical formulation of General Relativity that is oblivious to the affine spacetime structure. Generalisations along this line permit to generate teleparallel and symmetric teleparallel theories of gravity with exceptional properties. In this work we explore modified gravity theories based on non-linear extensions of the non-metricity scalar. After presenting some general properties and briefly studying some interesting background cosmologies (including accelerating solutions with relevance for inflation and dark energy), we analyse the behaviour of the cosmological perturbations. Tensor perturbations feature a re-scaling of the corresponding Newton's constant, while vector perturbations do not contribute in the absence of vector sources. In the scalar sector we find two additional propagating modes, hinting that $f(Q)$ theories introduce, at least, two additional degrees of freedom. These scalar modes disappear around maximally symmetric backgrounds because of the appearance of an accidental residual gauge symmetry corresponding to a restricted diffeomorphism. We finally discuss the potential strong coupling problems of these maximally symmetric backgrounds caused by the discontinuity in the number of propagating modes.

We reconsider the construction of general derivative self-interactions for a massive Proca field. The constructed Lagrangian is such that the vector field propagates at most three degrees of freedom, thus avoiding the ghostly nature of a fourth polarisation. The construction makes use of the well-known condition for constrained systems of having a degenerate Hessian. We briefly discuss the casuistry according to the nature of the existing constraints algebra. We also explore various classes of interesting new interactions that have been recently raised in the literature. For the sixth order Lagrangian that satisfies the constraints by itself we prove its topological character, making such a term irrelevant. There is however a window of opportunity for exploring other classes of fully-nonlinear interactions that satisfy the constraint algebra by mixing terms of various order.

We investigate protoneutron star matter using the state-of-the-art perturbative equation of state for cold and dense QCD in the presence of a fixed lepton fraction in which both electrons and neutrinos are included. Besides computing the modifications in the equation of state due to the presence of trapped neutrinos, we show that stable strange quark matter has a more restricted parameter space. We also study the possibility of nucleation of unpaired quark matter in the core of protoneutron stars by matching the lepton-rich QCD pressure onto a hadronic equation of state, namely TM1 with trapped neutrinos. Using the inherent dependence of perturbative QCD on the renormalization scale parameter, we provide a measure of the uncertainty in the observables we compute.

To understand the self-sustenance of subcritical turbulence in spectrally stable shear flows, we performed direct numerical simulations of homogeneous shear turbulence for different aspect ratios of the flow domain and analyzed the dynamical processes in Fourier space. There are no exponentially growing modes in such flows and the turbulence is energetically supported only by the linear growth of perturbation harmonics due to the shear flow non-normality. This non-normality-induced, or nonmodal growth is anisotropic in spectral space, which, in turn, leads to anisotropy of nonlinear processes in this space. As a result, a transverse (angular) redistribution of harmonics in Fourier space appears to be the main nonlinear process in these flows, rather than direct or inverse cascades. We refer to this type of nonlinear redistribution as the nonlinear transverse cascade. It is demonstrated that the turbulence is sustained by a subtle interplay between the linear nonmodal growth and the nonlinear transverse cascade that exemplifies a well-known bypass scenario of subcritical turbulence. These two basic processes mainly operate at large length scales, comparable to the domain size. Therefore, this central, small wave number area of Fourier space is crucial in the self-sustenance; we defined its size and labeled it as the vital area of turbulence. Outside the vital area, the nonmodal growth and the transverse cascade are of secondary importance. Although the cascades and the self-sustaining process of turbulence are qualitatively the same at different aspect ratios, the number of harmonics actively participating in this process varies, but always remains quite large. This implies that the self-sustenance of subcritical turbulence cannot be described by low-order models.

A large class of modified theories of gravity used as models for dark energy predict a propagation speed for gravitational waves which can differ from the speed of light. This difference of propagations speeds for photons and gravitons has an impact in the emission of gravitational waves by binary systems. Thus, we revisit the usual quadrupolar emission of binary system for an arbitrary propagation speed of gravitational waves and obtain the corresponding period decay formula. We then use timing data from the Hulse-Taylor binary pulsar and obtain that the speed of gravitational waves can only differ from the speed of light at the percentage level. This bound places tight constraints on dark energy models featuring an anomalous propagations speed for the gravitational waves.

In this work we revisit the construction of theories for a massive vector field with derivative self-interactions such that only the 3 desired polarizations corresponding to a Proca field propagate. We start from the decoupling limit by constructing healthy interactions containing second derivatives of the Stueckelberg field with itself and also with the transverse modes. The resulting interactions can then be straightforwardly generalized beyond the decoupling limit. We then proceed to a systematic construction of the interactions by using the Levi-Civita tensors. Both approaches lead to a finite family of allowed derivative self-interactions for the Proca field. This construction allows us to show that some higher order terms recently introduced as new interactions trivialize in 4 dimensions by virtue of the Cayley-Hamilton theorem. Moreover, we discuss how the resulting derivative interactions can be written in a compact determinantal form, which can also be regarded as a generalization of the Born-Infeld lagrangian for electromagnetism. Finally, we generalize our results for a curved background and give the necessary non-minimal couplings guaranteeing that no additional polarizations propagate even in the presence of gravity.

We extend previous results on healthy derivative self-interactions for a Proca field to the case of a set of massive vector fields. We obtain non-gauge invariant derivative self-interactions for the vector fields that maintain the appropriate number of propagating degrees of freedom. In view of the potential cosmological applications, we restrict to interactions with an internal rotational symmetry. We provide a systematical construction order by order in derivatives of the fields and making use of the antisymmetric Levi-Civita tensor. We then compare with the one single vector field case and show that the interactions can be broadly divided into two groups, namely the ones obtained from a direct extension of the generalized Proca terms and genuine multi-Proca interactions with no correspondence in the single Proca case. We also discuss the curved spacetime version of the interactions to include the necessary non-minimal couplings to gravity. Finally, we explore the cosmological applications and show that there are three different vector fields configurations giving rise to isotropic solutions. Two of them have already been considered in the literature and the third one, representing a combination of the first two, is new and offers unexplored cosmological scenarios.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

Spacetime with general linear vector distortion is introduced. Thus, the torsion and the nonmetricity of the affine connection are assumed to be proportional to a vector field (and not its derivatives). The resulting two-parameter family of non-Riemannian geometries generalises the conformal Weyl geometry and some other interesting special cases. Taking into account the leading order quadratic curvature correction to the Einstein-Hilbert action results uniquely in the one-parameter extension of the Starobinsky inflation known as the alpha-attractor. The most general quadratic curvature action introduces, in addition to the canonical vector kinetic term, novel ghost-free vector-tensor interactions.

We study a general class of gravitational theories formulated in the Palatini approach and derive the equations governing the evolution of tensor perturbations. In the absence of torsion, the connection can be solved as the Christoffel symbols of an auxiliary metric which is non-trivially related to the space-time metric. We then consider background solutions corresponding to a perfect fluid and show that the tensor perturbations equations (including anisotropic stresses) for the auxiliary metric around such a background take an Einstein-like form. This facilitates the study in a homogeneous and isotropic cosmological scenario where we explicitly establish the relation between the auxiliary metric and the space-time metric tensor perturbations. As a general result, we show that both tensor perturbations coincide in the absence of anisotropic stresses.

We generalize the ultraviolet sector of gravitation via a Born-Infeld action using lessons from massive gravity. The theory contains all of the elementary symmetric polynomials and is treated in the Palatini formalism. We show how the connection can be solved algebraically to be the Levi-Civita connection of an effective metric. The non-linearity of the algebraic equations yields several branches, one of which always reduces to General Relativity at low curvatures. We explore in detail a {\it minimal} version of the theory, for which we study solutions in the presence of a perfect fluid with special attention to the cosmological evolution. In vacuum we recover Ricci-flat solutions, but also an additional physical solution corresponding to an Einstein space. The existence of two physical branches remains for non-vacuum solutions and, in addition, the branch that connects to the Einstein space in vacuum is not very sensitive to the specific value of the energy density. For the branch that connects to the General Relativity limit we generically find three behaviours for the Hubble function depending on the equation of state of the fluid, namely: either there is a maximum value for the energy density that connects continuously with vacuum, or the energy density can be arbitrarily large but the Hubble function saturates and remains constant at high energy densities, or the energy density is unbounded and the Hubble function grows faster than in General Relativity. The second case is particularly interesting because it could offer an interesting inflationary epoch even in the presence of a dust component. Finally, we discuss the possibility of avoiding certain types of singularities within the minimal model.

We study the Horndeski vector-tensor theory that leads to second order equations of motion and contains a non-minimally coupled abelian gauge vector field. This theory is remarkably simple and consists of only 2 terms for the vector field, namely: the standard Maxwell kinetic term and a coupling to the dual Riemann tensor. Furthermore, the vector sector respects the U(1) gauge symmetry and the theory contains only one free parameter, M^2, that controls the strength of the non-minimal coupling. We explore the theory in a de Sitter spacetime and study the presence of instabilities and show that it corresponds to an attractor solution in the presence of the vector field. We also investigate the cosmological evolution and stability of perturbations in a general FLRW spacetime. We find that a sufficient condition for the absence of ghosts is M^2>0. Moreover, we study further constraints coming from imposing the absence of Laplacian instabilities. Finally, we study the stability of the theory in static and spherically symmetric backgrounds (in particular, Schwarzschild and Reissner-Nordstr\"om-de Sitter). We find that the theory, quite generally, do have ghosts or Laplacian instabilities in regions of spacetime where the non-minimal interaction dominates over the Maxwell term. We also calculate the propagation speed in these spacetimes and show that superluminality is a quite generic phenomenon in this theory.

A screening mechanism for conformal vector-tensor modifications of general relativity is proposed. The conformal factor depends on the norm of the vector field and makes the field to vanish in high dense regions, whereas drives it to a non-null value in low density environments. Such process occurs due to a spontaneous symmetry breaking mechanism and gives rise to both the screening of fifth forces as well as Lorentz violations. The cosmology and local constraints are also computed.

We study the cosmological evolution for a universe in the presence of a continuous tower of massive scalar fields which can drive the current phase of accelerated expansion of the universe and, in addition, can contribute as a dark matter component. The tower consists of a continuous set of massive scalar fields with a gaussian mass distribution. We show that, in a certain region of the parameter space, the {\it heavy} modes of the tower (those with masses much larger than the Hubble expansion rate) dominate at early times and make the tower behave like the usual single scalar field whose coherent oscillations around the minimum of the potential give a matter-like contribution. On the other hand, at late times, the {\it light} modes (those with masses much smaller than the Hubble expansion rate) overcome the energy density of the tower and they behave like a perfect fluid with equation of state ranging from 0 to -1, depending on the spectral index of the initial spectrum. This is a distinctive feature of the tower with respect to the case of quintessence fields, since a massive scalar field can only give acceleration with equation of state close to -1. Such unique property is the result of a synergy effect between the different mass modes. Interestingly, we find that, for some choices of the spectral index, the tower tracks the matter component at high redshifts (or it can even play the role of the dark matter) and eventually becomes the dominant component of the universe and give rise to an accelerated expansion.

We show that a longitudinal gauge degree of freedom for a vector field is equivalent to a Pais-Uhlenbeck scalar field. With the help of this equivalence, we can determine natural interactions of this field with scalars and fermions. Since the theory has a global U(1) symmetry, we have the usual conserved current of the charged fields, thanks to which the dynamics of the scalar field is not modified by the interactions. We use this fact to consistently quantize the theory even in the presence of interactions. We argue that such a degree of freedom can only be excited by gravitational effects like the inflationary era of the early universe and may play the role of dark energy in the form of an effective cosmological constant whose value is linked to the inflation scale.

n this paper we study the effect of biasing on the power spectrum at large scales. We show that even though non-linear biasing does introduce a white noise contribution on large scales, the $P(k)\propto k^n$ behavior of the matter power spectrum on large scales may still be visible and above the white noise for about one decade. We show, that the Kaiser biasing scheme which leads to linear bias of the correlation function on {\em large} scales, also generates a linear bias of the {\rm power spectrum} on rather small scales. This is a consequence of the divergence on small scales of the pure Harrison-Zeldovich spectrum. However, biasing becomes k-dependent when we damp the underlying power spectrum on small scales. We also discuss the effect of biasing on the baryon acoustic oscillations.

We discuss an extended version of electromagnetism in which the usual gauge fixing term is promoted into a physical contribution that introduces a new scalar state in the theory. This new state can be generated from vacuum quantum fluctuations during an inflationary era and, on super-Hubble scales, gives rise to an effective cosmological constant. The value of such a cosmological constant coincides with the one inferred from observations as long as inflation took place at the electroweak scale. On the other hand, the new state also generates an effective electric charge density on sub-Hubble scales that produces both vorticity and magnetic fields with coherent lengths as large as the present Hubble horizon.