The Fisher gAlaxy suRvey cOde ($\texttt{FARO}$) is a new public Python code that computes the Fisher matrix for galaxy surveys observables. The observables considered are the linear multitracer 3D galaxy power spectrum, the linear convergence power spectrum for weak lensing, and the linear multitracer power spectrum for the correlation between galaxy distribution and convergence. The code allows for tomographic and model-independent analysis in which, for scale-independent growth, the following functions of redshift $A_a (z) \equiv \sigma_{8}(z) \, b_a (z)$, $R(z) \equiv \sigma_{8}(z) \, f(z)$, $L(z) \equiv \Omega_{m} \, \sigma_{8}(z) \, \Sigma (z)$ and $E(z) \equiv H(z)/H_0$, together with the function of scale $\hat{P}(k)$, are taken as free parameters in each redshift and scale bins respectively. In addition, a module for change of variables is provided to project the Fisher matrix on any particular set of parameters required. The code is built to be as fast as possible and user-friendly. As an application example, we forecast the sensitivity of future galaxy surveys like DESI, Euclid, J-PAS and LSST and compare their performance on different redshift and scale ranges.

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)

The next generation of galaxy surveys will allow us to test one of the most fundamental assumptions of the standard cosmology, i.e., that gravity is governed by the general theory of relativity (GR). In this paper we investigate the ability of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) to constrain GR and its extensions. Based on the J-PAS information on clustering and gravitational lensing, we perform a Fisher matrix forecast on the effective Newton constant, $\mu$, and the gravitational slip parameter, $\eta$, whose deviations from unity would indicate a breakdown of GR. Similar analysis is also performed for the DESI and Euclid surveys and compared to J-PAS with two configurations providing different areas, namely an initial expectation with 4000 $\mathrm{deg}^2$ and the future best case scenario with 8500 $\mathrm{deg}^2$. We show that J-PAS will be able to measure the parameters $\mu$ and $\eta$ at a sensitivity of $2\% - 7\%$, and will provide the best constraints in the interval $z = 0.3 - 0.6$, thanks to the large number of ELGs detectable in that redshift range. We also discuss the constraining power of J-PAS for dark energy models with a time-dependent equation-of-state parameter of the type $w(a)=w_0+w_a(1-a)$, obtaining $\Delta w_0=0.058$ and $\Delta w_a=0.24$ for the absolute errors of the dark energy parameters.

We analyze the sensitivity of galaxy and weak-lensing surveys to detect preferred directions in the gravitational interaction. We consider general theories of gravity involving additional vector degrees of freedom with non-vanishing spatial components in the background. We use a model-independent parametrization of the perturbations equations in terms of four effective parameters, namely, the standard effective Newton constant $G_{eff}$ and slip parameter $\gamma$ for scalar modes and two new parameters $\mu_Q$ and $\mu_h$ for vector and tensor modes respectively, which are required when preferred directions are present. We obtain the expressions for the multipole galaxy power spectrum in redshift space and for the weak-lensing shear, convergence and rotation spectra in the presence of preferred directions. By performing a Fisher matrix forecast analysis, we estimate the sensitivity of a future Euclid-like survey to detect this kind of modification of gravity. We finally compare with the effects induced by violations of statistical isotropy in the primordial power spectrum and identify the observables which could discriminate between them.

We consider the problem of parametrizing modified gravity theories that include an additional vector field in the sub-Hubble regime within the quasi-static approximation. We start from the most general set of second order equations for metric and vector field perturbations and allow for both temporal and spatial components of the background vector field. We find that in the case in which dark matter obeys standard conservation equations, eight parameters are needed to fully characterize the theory. If dark matter vorticity can be neglected, the number of independent parameters is reduced to four. In addition to the usual scale and redshift dependence, the effective parameters have an additional angular dependence induced by the preferred direction set by the background vector. In the considered sub-Hubble regime, we show that this angular dependence appears only through even multipoles and generates anisotropies in the growth function which translate into anisotropies in the galaxy and lensing convergence power spectra. The angular dependence generated by the preferred direction is different from that induced by redshift space distortions and could be disentangled in the data collected by future galaxy surveys.

It is well-known that an extremely accurate parametrization of the growth function of matter density perturbations in $\Lambda$CDM cosmology, with errors below $0.25 \%$, is given by $f(a)=\Omega_{m}^{\gamma} \,(a)$ with $\gamma \simeq 0.55$. In this work, we show that a simple modification of this expression also provides a good description of growth in modified gravity theories. We consider the model-independent approach to modified gravity in terms of an effective Newton constant written as $\mu(a,k)=G_{eff}/G$ and show that $f(a)=\beta(a)\Omega_{m}^{\gamma} \,(a)$ provides fits to the numerical solutions with similar accuracy to that of $\Lambda$CDM. In the time-independent case with $\mu=\mu(k)$, simple analytic expressions for $\beta(\mu)$ and $\gamma(\mu)$ are presented. In the time-dependent (but scale-independent) case $\mu=\mu(a)$, we show that $\beta(a)$ has the same time dependence as $\mu(a)$. As an example, explicit formalae are provided in the DGP model. In the general case, for theories with $\mu(a,k)$, we obtain a perturbative expansion for $\beta(\mu)$ around the General Relativity case $\mu=1$ which, for $f(R)$ theories, reaches an accuracy below $1 \%$. Finally, as an example we apply the obtained fitting functions in order to forecast the precision with which future galaxy surveys will be able to measure the $\mu$ parameter.

In the context of f(R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star's edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic $R^2$ and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable f(R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the f(R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.