Search

Context. The DESI Legacy Imaging Survey DR9, with its extensive dataset of galaxy locations and photometric redshifts, presents an opportunity to study Baryon Acoustic Oscillations (BAO) in the region covered by the ongoing DESI spectroscopic survey. Aims. We aim to investigate differences between different parts of the DR9 footprint. Furthermore, we want to measure the BAO scale for luminous red galaxies within them. Our selected redshift range of 0.6 to 0.8 corresponds to the bin in which a tension between DESI Y1 and eBOSS was found. Methods. We calculated the anisotropic two-point correlation function in a modified binning scheme to detect the BAO in DR9 data. We then use template fits based on simulations to measure the BAO scale in the imaging data. Results. Our analysis revealed the expected correlation function shape in most of the footprint areas, showing a BAO scale consistent with Planck's observations. Aside from identified mask-related data issues in the southern region of the South Galactic Cap, we also find a notable variance between the different footprints. Conclusions. We find that this variance is consistent with the difference between the DESI Y1 and eBOSS data and it supports the argument that that tension is caused by sample variance. Additionally, we also uncovered systematic biases not previously accounted for in photometric BAO studies. We emphasize the necessity of adjusting for the systematic shift in the BAO scale associated with typical photometric redshift uncertainties to ensure accurate measurements., Comment: 20 pages, 14 figures, 3 tables, accepted for publication in Astronomy and Astrophysics

We present three separate void catalogs created using a volume-limited sample of the DESI Year 1 Bright Galaxy Survey. We use the algorithms VoidFinder and V2 to construct void catalogs out to a redshift of z=0.24. We obtain 1,461 interior voids with VoidFinder, 420 with V2 using REVOLVER pruning, and 295 with V2 using VIDE pruning. Comparing our catalog with an overlapping SDSS void catalog, we find generally consistent void properties but significant differences in the void volume overlap, which we attribute to differences in the galaxy selection and survey masks. These catalogs are suitable for studying the variation in galaxy properties with cosmic environment and for cosmological studies., Comment: 17 pages, 6 figures

We measure the growth of cosmic density fluctuations on large scales and across the redshift range $0.3

We present the first robust measurement of the one-dimensional Lyman alpha (Ly$\alpha$) forest bispectrum using the complete extended Baryon Oscillation Spectroscopic Survey (eBOSS) quasar sample, corresponding to the sixteenth data release (DR16) of the Sloan Digital Sky Survey (SDSS). The measurement employs an FFT estimator over 12 redshift bins, ranging from $z=2.2$ to $z=4.4$, and extends to scales of $0.02 ~ (\mathrm{km/s})^{-1}$. The sample consists of 122,066 quasar spectra, although only the first six redshift bins contain sufficient data to extract a physical bispectrum. To validate and correct the bispectrum measurement, we use synthetic datasets generated from lognormal and 2LPT mocks. Additionally, we detect clear evidence of correlations between Si$_\mathrm{III}$ absorption lines and the Ly$\alpha$ forest within the bispectrum signal, which we describe with an extension of the model used for the analogue of 1d power spectrum signal. In this context, the pipeline developed for this study addresses the impact of instrumental and methodological systematics and is ready for application to larger spectroscopic datasets, such as those from the first year of DESI observations. Finally, A simple perturbation theory model provides a reasonable explanation of the eBOSS bispectrum, suggesting that higher-order one-dimensional statistics in the Ly$\alpha$ forest can complement cosmological model inference based on the power spectrum in future analyses., Comment: 38 pages, 19 figures

We present new Galactic reddening maps of the high Galactic latitude sky using DESI imaging and spectroscopy. We directly measure the reddening of 2.6 million stars by comparing the observed stellar colors in $g-r$ and $r-z$ from DESI imaging with the synthetic colors derived from DESI spectra from the first two years of the survey. The reddening in the two colors is on average consistent with the \cite{fitzpatrick_correcting_1999} extinction curve with $R_\mathrm{V}=3.1$. We find that our reddening maps differ significantly from the commonly used \cite{schlegel_maps_1998} (SFD) reddening map (by up to 80 mmag in $E(B-V)$), and we attribute most of this difference to systematic errors in the SFD map. To validate the reddening map, we select a galaxy sample with extinction correction based on our reddening map, and this yields significantly better uniformity than the SFD extinction correction. Finally, we discuss the potential systematic errors in the DESI reddening measurements, including the photometric calibration errors that are the limiting factor on our accuracy. The $E(g-r)$ and $E(g-r)$ maps presented in this work, and for convenience their corresponding $E(B-V)$ maps with SFD calibration, are publicly available., Comment: Submitted to the Open Journal of Astrophysics. Associated data files: https://data.desi.lbl.gov/public/papers/mws/desi_dust/y2/v1/maps/

One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $O(N^3)$ with the number of objects, can be further improved with tree methods but not enough to deal with large scale correlations of Rubin's data. However, a harmonic basis decomposition of these higher order statistics reduces the time dramatically, to scale as a two-point correlation function with the number of objects, so that the signal can be extracted in a reasonable amount of time. In this work, we aim to develop the framework to use these expansions within the Limber approximation for scalar (or spin-0) fields, such as galaxy counts, weak lensing convergence or aperture masses. We develop an estimator to extract the signal from catalogs and different phenomenological and theoretical models for its description. The latter includes halo model and standard perturbation theory, to which we add a simple effective field theory prescription based on the short range of non-locality of cosmic fields, significantly improving the agreement with simulated data. In parallel to the modeling of the signal, we develop a code that can efficiently calculate three points correlations of more than 200 million data points (a full sky simulation with Nside=4096) in $\sim$40 minutes, or even less than 10 minutes using an approximation in the searching algorithm, on a single high-performance computing node, enabling a feasible analysis for the upcoming LSST data., Comment: version 2: 31 pages, 10 figures

Measuring the growth of structure is a powerful probe for studying the dark sector, especially in light of the $\sigma_8$ tension between primary CMB anisotropy and low-redshift surveys. This paper provides a new measurement of the amplitude of the matter power spectrum, $\sigma_8$, using galaxy-galaxy and galaxy-CMB lensing power spectra of Dark Energy Spectroscopic Instrument Legacy Imaging Surveys Emission-Line Galaxies and the $\textit{Planck}$ 2018 CMB lensing map. We create an ELG catalog composed of $27$ million galaxies and with a purity of $85\%$, covering a redshift range $0 < z < 3$, with $z_{\rm mean} = 1.09$. We implement several novel systematic corrections, such as jointly modeling the contribution of imaging systematics and photometric redshift uncertainties to the covariance matrix. We also study the impacts of various dust maps on cosmological parameter inference. We measure the cross-power spectra over $f_{\rm sky} = 0.25$ with a signal-to-background ratio of up to $ 30\sigma$. We find that the choice of dust maps to account for imaging systematics in estimating the ELG overdensity field has a significant impact on the final estimated values of $\sigma_8$ and $\Omega_{\rm M}$, with far-infrared emission-based dust maps preferring $\sigma_8$ to be as low as $0.702 \pm 0.030$, and stellar-reddening-based dust maps preferring as high as $0.719 \pm 0.030$. The highest preferred value is at $\sim 3 \sigma$ tension with the $\textit{Planck}$ primary anisotropy results. These findings indicate a need for tomographic analyses at high redshifts and joint modeling of systematics., Comment: 50 pages, 24 figures (figure data can be obtained at https://zenodo.org/records/13381499)

We explore correlations between the orientations of small galaxy groups, or "multiplets", and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of 100 Mpc/h. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift z = 1. Multiplet alignment is a form of higher-order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale (< 1 Mpc/h) dynamics. Within samples of bright galaxies (BGS), luminous red galaxies (LRG) and emission-line galaxies (ELG), we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no intrinsic alignment. DESI's LRG mock galaxy catalogues created from the AbacusSummit N-body simulations produce a similar alignment signal, though with a 33% lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 Mpc/h. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples., Comment: For an accessible summary of this paper, see https://cmlamman.github.io/doc/multipletIA_summary.pdf

We infer the growth of large scale structure over the redshift range $0.4\lesssim z \lesssim 1$ from the cross-correlation of spectroscopically calibrated Luminous Red Galaxies (LRGs) selected from the Dark Energy Spectroscopic Instrument (DESI) legacy imaging survey with CMB lensing maps reconstructed from the latest Planck and ACT data. We adopt a hybrid effective field theory (HEFT) model that robustly regulates the cosmological information obtainable from smaller scales, such that our cosmological constraints are reliably derived from the (predominantly) linear regime. We perform an extensive set of bandpower- and parameter-level systematics checks to ensure the robustness of our results and to characterize the uniformity of the LRG sample. We demonstrate that our results are stable to a wide range of modeling assumptions, finding excellent agreement with a linear theory analysis performed on a restricted range of scales. From a tomographic analysis of the four LRG photometric redshift bins we find that the rate of structure growth is consistent with $\Lambda$CDM with an overall amplitude that is $\simeq5-7\%$ lower than predicted by primary CMB measurements with modest $(\sim2\sigma)$ statistical significance. From the combined analysis of all four bins and their cross-correlations with Planck we obtain $S_8 = 0.765\pm0.023$, which is less discrepant with primary CMB measurements than previous DESI LRG cross Planck CMB lensing results. From the cross-correlation with ACT we obtain $S_8 = 0.790^{+0.024}_{-0.027}$, while when jointly analyzing Planck and ACT we find $S_8 = 0.775^{+0.019}_{-0.022}$ from our data alone and $\sigma_8 = 0.772^{+0.020}_{-0.023}$ with the addition of BAO data. These constraints are consistent with the latest Planck primary CMB analyses at the $\simeq 1.6-2.2\sigma$ level, and are in excellent agreement with galaxy lensing surveys., Comment: 60 pages, 26 figures, comments welcome

We present a high-significance cross-correlation of CMB lensing maps from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) with spectroscopically calibrated luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI). We detect this cross-correlation at a significance of 38$\sigma$; combining our measurement with the Planck Public Release 4 (PR4) lensing map, we detect the cross-correlation at 50$\sigma$. Fitting this jointly with the galaxy auto-correlation power spectrum to break the galaxy bias degeneracy with $\sigma_8$, we perform a tomographic analysis in four LRG redshift bins spanning $0.4 \le z \le 1.0$ to constrain the amplitude of matter density fluctuations through the parameter combination $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4}$. Prior to unblinding, we confirm with extragalactic simulations that foreground biases are negligible and carry out a comprehensive suite of null and consistency tests. Using a hybrid effective field theory (HEFT) model that allows scales as small as $k_{\rm max}=0.6$ $h/{\rm Mpc}$, we obtain a 3.3% constraint on $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4} = 0.792^{+0.024}_{-0.028}$ from ACT data, as well as constraints on $S_8^\times(z)$ that probe structure formation over cosmic time. Our result is consistent with the early-universe extrapolation from primary CMB anisotropies measured by Planck PR4 within 1.2$\sigma$. Jointly fitting ACT and Planck lensing cross-correlations we obtain a 2.7% constraint of $S_8^\times = 0.776^{+0.019}_{-0.021}$, which is consistent with the Planck early-universe extrapolation within 2.1$\sigma$, with the lowest redshift bin showing the largest difference in mean. The latter may motivate further CMB lensing tomography analyses at $z<0.6$ to assess the impact of potential systematics or the consistency of the $\Lambda$CDM model over cosmic time., Comment: Prepared for submission to JCAP (47 pages, 13 figures)

Dark Energy Spectroscopic Instrument (DESI) uses more than 2.4 million Emission Line Galaxies (ELGs) for 3D large-scale structure (LSS) analyses in its Data Release 1 (DR1). Such large statistics enable thorough research on systematic uncertainties. In this study, we focus on spectroscopic systematics of ELGs. The redshift success rate ($f_{\rm goodz}$) is the relative fraction of secure redshifts among all measurements. It depends on observing conditions, thus introduces non-cosmological variations to the LSS. We, therefore, develop the redshift failure weight ($w_{\rm zfail}$) and a per-fibre correction ($\eta_{\rm zfail}$) to mitigate these dependences. They have minor influences on the galaxy clustering. For ELGs with a secure redshift, there are two subtypes of systematics: 1) catastrophics (large) that only occur in a few samples; 2) redshift uncertainty (small) that exists for all samples. The catastrophics represent 0.26\% of the total DR1 ELGs, composed of the confusion between O\,\textsc{ii} and sky residuals, double objects, total catastrophics and others. We simulate the realistic 0.26\% catastrophics of DR1 ELGs, the hypothetical 1\% catastrophics, and the truncation of the contaminated $1.31

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing power of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying Large-Scale Structure systematic contamination (e.g. Cosmic-Infrared Background) in the Galactic extinction maps is a likely root cause. We additionally observe that the DESI LRGs sample exhibits a complex dependency on a combination of seeing, depth, and intrinsic galaxy brightness, which is not replicated by Obiwan, suggesting discrepancies between the current simulation settings and the actual observations. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRG sample., Comment: 46 pages, 26 figures

Lyman break galaxies (LBGs) are promising probes for clustering measurements at high redshift, $z>2$, a region only covered so far by Lyman-$\alpha$ forest measurements. In this paper, we investigate the feasibility of selecting LBGs by exploiting the existence of a strong deficit of flux shortward of the Lyman limit, due to various absorption processes along the line of sight. The target selection relies on deep imaging data from the HSC and CLAUDS surveys in the $g,r,z$ and $u$ bands, respectively, with median depths reaching 27 AB in all bands. The selections were validated by several dedicated spectroscopic observation campaigns with DESI. Visual inspection of spectra has enabled us to develop an automated spectroscopic typing and redshift estimation algorithm specific to LBGs. Based on these data and tools, we assess the efficiency and purity of target selections optimised for different purposes. Selections providing a wide redshift coverage retain $57\%$ of the observed targets after spectroscopic confirmation with DESI, and provide an efficiency for LBGs of $83\pm3\%$, for a purity of the selected LBG sample of $90\pm2\%$. This would deliver a confirmed LBG density of $\sim 620$ deg$^{-2}$ in the range $2.3

We present strong constraints on the spacetime variation of the fine-structure constant $\alpha$ using the Dark Energy Spectroscopic Instrument (DESI). In this pilot work, we utilize $\sim110,000$ galaxies with strong and narrow O III $\lambda\lambda$4959,5007 emission lines to measure the relative variation $\Delta\alpha/\alpha$ in space and time. The O III doublet is arguably the best choice for this purpose owing to its wide wavelength separation between the two lines and its strong emission in many galaxies. Our galaxy sample spans a redshift range of $0

We explore black hole solutions and some of its physical properties in Einstein's theory in 4D, modified by a cubic gravity term and in the presence of non-linear electrodynamics. In the context of Effective Field Theories (EFT) and under certain assumptions, these curvature and non-linear electromagnetic terms represent the first corrections to the Einstein-Maxwell theory. We obtain static and spherically symmetric generalizations to the asymptotically flat Reissner-Nordstr\"om metric using perturbative methods, showing how an asymptotic expansion solution connects with a near horizon solution for a small coupling of the curvature correction term. Finally, we discuss how these EFT corrections change the event horizon properties and also lead to measurable effects on black hole shadows and gravitational lensing around these solutions., Comment: 18 pages, 9 figures; V2, References added

We present the one-dimensional Lyman-$\alpha$ forest power spectrum measurement using the first data provided by the Dark Energy Spectroscopic Instrument (DESI). The data sample comprises $26,330$ quasar spectra, at redshift $z > 2.1$, contained in the DESI Early Data Release and the first two months of the main survey. We employ a Fast Fourier Transform (FFT) estimator and compare the resulting power spectrum to an alternative likelihood-based method in a companion paper. We investigate methodological and instrumental contaminants associated to the new DESI instrument, applying techniques similar to previous Sloan Digital Sky Survey (SDSS) measurements. We use synthetic data based on log-normal approximation to validate and correct our measurement. We compare our resulting power spectrum with previous SDSS and high-resolution measurements. With relatively small number statistics, we successfully perform the FFT measurement, which is already competitive in terms of the scale range. At the end of the DESI survey, we expect a five times larger Lyman-$\alpha$ forest sample than SDSS, providing an unprecedented precise one-dimensional power spectrum measurement., Comment: 23 pages, 23 figures, Journal version

The next generation of galaxy surveys will provide highly accurate measurements of the large-scale structure of the Universe, allowing for more stringent tests of gravity on cosmological scales. Higher order statistics are a valuable tool to study the non-Gaussianities in the matter field and to break degeneracies between modified gravity and other physical or nuisance parameters. However, understanding from first principles the behaviour of these correlations is essential to characterise deviations from General Relativity (GR), and the purpose of this work. This work uses contemporary ideas of Standard Perturbation Theory on biased tracers to characterize the three point correlation function (3PCF) at tree level for Modified Gravity models with a scale-dependent gravitational strength, and applies the theory to two specific models ($f(R)$ and DGP) that are representative for Chameleon and Vainshtein screening mechanisms. Additionally, we use a multipole decomposition, which apart from speeding up the algorithm to extract the signal from data, also helps to visualize and characterize GR deviations., Comment: 25 pages, 5 figures

We establish the existence of time-dependent solitons in a modified gravity framework, which is defined by the low energy limit of theories with a weakly broken galileon symmetry and a mass term. These are regular vacuum configurations of finite energy characterized by a single continuous parameter representing the amplitude of the scalar degree of freedom at the origin. When the central field amplitude is small the objects are indistinguishable from boson stars. In contrast, increasing the central value of the amplitude triggers the effect of higher derivative operators in the effective theory, leading to departures from the previous solutions, until the theory becomes strongly coupled and model-dependent. The higher order operators are part of the (beyond) Horndeski theory, hence the name of the compact objects. Moreover, a remnant of the galileon non-renormalization theorem guarantees that the existence and properties of these solutions are not affected by quantum corrections. Finally, we discuss the linear stability under small radial perturbations, the mass-radius relation, the compactness, the appearance of innermost stable circular orbits and photon spheres, and some astrophysical signatures (accretion disks, gravitational radiation and lensing) that may be relevant to falsify the model., Comment: 49+1 pages, 17 figures, 2 tables. To appear in JCAP

We introduce an Eulerian Perturbation Theory to study the clustering of tracers for cosmologies in the presence of massive neutrinos. Our approach is based on mapping recently-obtained Lagrangian Perturbation Theory results to the Eulerian framework. We add Effective Field Theory counterterms, IR-resummations and a biasing scheme to compute the one-loop redshift-space power spectrum. To assess our predictions, we compare the power spectrum multipoles against synthetic halo catalogues from the Quijote simulations, finding excellent agreement on scales $k\lesssim 0.25 \,h \text{Mpc}^{-1}$. One can obtain the same fitting accuracy using higher wave-numbers, but then the theory fails to give a correct estimation of the linear bias parameter. We further discuss the implications for the tree-level bispectrum. Finally, calculating loop corrections is computationally costly, hence we derive an accurate approximation wherein we retain only the main features of the kernels, as produced by changes to the growth rate. As a result, we show how FFTLog methods can be used to further accelerate the loop computations with these reduced kernels., Comment: 42 pages, 10 figures

Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for -- its applications on those scales being largely based on extrapolation and its validity sometimes questioned in the shadow of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of $N$-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel estimators beyond linear redshift distortions in two well-studied modified gravity models, chameleon $f(R)$ gravity and a braneworld model, and the potential of testing these deviations from GR using DESI. These estimators employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI., Comment: 69 pages, 36 figures, 6 tables; accepted for publication in JCAP. The modified gravity HOD galaxy catalogues used in this paper can be downloaded from http://icc.dur.ac.uk/data/, by searching for the title of this paper from the dataset table

In the forthcoming large volume galaxy surveys higher order statistics will provide complementary information to the usual two point statistics. Low variance estimators of the Three Point Correlation Function (3CPF) of discrete data count triangle configurations with vertices mixing data and random catalogues. Large density random catalogues are used to reduce the shot noise, which leads to a computational cost of one or two orders of magnitude more than the pure data histogram. In this paper, we explore time reductions of the isotropic 3PCF random sampling terms in periodic boxes without using random catalogues. In the first approach, based on Hamilton's construction of his famous two point estimator, we use an ad-hoc two point correlation term, while for the second procedure we construct the operators from a geometrical viewpoint, using two sides and their opening angle to describe the 3PCF triangle configurations. We map the last result to the three triangle side basis either numerically or analytically, and show that the latter approach performs best when applied to synthetic data. Moreover, we elaborate on going beyond periodic boxes, discuss other low variance n-point estimators and present useful 3PCF visualization schemes., Comment: minor changes

It is commonly believed that a successful period of inflation driven by a single or several scalar fields requires a specific hierarchy of masses given by $M_{inf}\ll H \ll M_{heavy}$, where $M_{inf}$ can correspond to several or a single light field and $M_{heavy}\,\,\,$ corresponds to any heavy field that might be integrated out if it satisfies suitable conditions. This is at the heart of the so called $\eta$-problem in inflation, since large contributions to the masses of the inflatons might spoil the slow-roll conditions required for inflation. We show that, while this is an unavoidable conclusion in single field inflation, in multifield inflation, heavy fields as defined above, may be fully responsible for a successful period of what we call fat slow-roll inflation. Moreover we show that in this scenario, the turning rate of the inflationary trajectory, $\Omega/H$, is larger than one. Thus, the $\eta$-problem is evaded with large turns in fat inflation. Depending on the perturbations' mass spectra, cosmological predictions will differ either slightly or largely with respect to those of the single field case. We illustrate this scenario in a concrete example in Type IIB string flux compactifications, where a probe D5-brane moving along the angular and radial directions in a warped throat drives fat D5-brane natural inflation. An instantaneous superplanckian decay constant can be defined, consistent with our low energy approximations, thanks to the strong warping of the geometry. We compute the cosmological observables, which differ from the single field case. We also discuss fat inflation in the context of recently proposed swampland de Sitter conjectures., Comment: 41 pages, 13 figures. Minor typos fixed, references added, appendix A expanded. To appear in JCAP

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

Starting from the Gleyzes-Langlois-Piazza-Vernizzi action, we derive the most general effective theory that is invariant under internal shifts and a $\mathbb{Z}_2$ mirror symmetry in the scalar sector. Contrary to what one may think, this model presents a dark matter tracker previous to the dark energy domination. We show that, in an empty universe and to linear order in perturbations, the scalar mode clusters in exactly the same way as standard nonrelativistic cold dark matter. This also holds for the subsector of the theory where the speed of propagation of gravitational waves equals that of light, in agreement with the recent multimessenger observation. However, the inclusion of standard model particles introduces nontrivial couplings of the gravitational scalar mode to baryons, modifying their clustering properties. We argue that no arrangement of the parameters of the model can reduce the extra scalar to precisely behave as cold dark matter., Comment: 11 pages, no figures. To appear in Phys. Rev. D

The direct detection of gravitational waves opens new perspectives for measuring properties of gravitationally bound compact objects. It is then important to investigate black holes and neutron stars in alternative theories of gravity, since they can have features that make them observationally distinguishable from their General Relativity (GR) counterparts. In this work, we examine a special case of vector Galileons, a vector-tensor theory of gravity with interesting cosmological properties, which consists of a one parameter modification of the Einstein-Maxwell action. Within this theory, we study configurations describing asymptotically flat, spherically symmetric black holes and neutron stars. The set of black hole solutions in this theory is surprisingly rich, generalising results found in GR or in related scalar-tensor theories. We investigate the properties and conserved charges of black holes, using both analytical and numerical techniques, highlighting configurations that are more compact than in GR. We then study properties of neutron stars, showing how the vector profile can influence the star internal structure. Depending on properties of matter and fields inside the star, neutron stars can be more massive than in GR, and they can be even more compact than Schwarzschild black holes, making these objects observationally interesting. We also comment on possible extensions of our configurations to magnetically charged or rotating configurations., Comment: 16 pages plus appendices, 9 figures

In a recent paper, Erik Verlinde has developed the interesting possibility that spacetime and gravity may emerge from the entangled structure of an underlying microscopic theory. In this picture, dark matter arises as a response to the standard model of particle physics from the delocalized degrees of freedom that build up the dark energy component of the Universe. Dark matter physics is then regulated by a characteristic acceleration scale $a_0$, identified with the radius of the (quasi)-de Sitter universe we inhabit. For a point particle matter source, or outside an extended spherically symmetric object, MOND's empirical fitting formula is recovered. However, Verlinde's theory critically departs from MOND when considering the inner structure of galaxies, differing by a factor of 2 at the centre of a regular massive body. For illustration, we use the eight classical dwarf spheroidal satellites of the Milky Way. These objects are perfect testbeds for the model given their approximate spherical symmetry, measured kinematics, and identified missing mass. We show that, without the assumption of a maximal deformation, Verlinde's theory can fit the velocity dispersion profile in dwarf spheroidals with no further need of an extra dark particle component. If a maximal deformation is considered, the theory leads to mass-to-light ratios that are marginally larger than expected from stellar population and formation history studies. We also compare our results with the recent phenomenological interpolating MOND function of McGaugh {\it et al}, and find a departure that, for these galaxies, is consistent with the scatter in current observations., Comment: 12 pages, 5 figures, 1 table. To appear in MNRAS

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Black hole configurations offer insights on the non-linear aspects of gravitational theories, and can suggest testable predictions for modifications of General Relativity. In this work, we examine exact black hole configurations in vector-tensor theories, originally proposed to explain dark energy by breaking the Abelian symmetry with a non-minimal coupling of the vector to gravity. We are able to evade the no-go theorems by Bekenstein on the existence of regular black holes in vector-tensor theories with Proca mass terms, and exhibit regular black hole solutions with a profile for the longitudinal vector polarization, characterised by an additional charge. We analytically find the most general static, spherically symmetric black hole solutions with and without a cosmological constant, and study in some detail their features, such as how the geometry depends on the vector charges. We also include angular momentum, and find solutions describing slowly-rotating black holes. Finally, we extend some of these solutions to higher dimensions., Comment: 13 pages, v2: references added, minor typos corrected

In the context of a cubic Galileon model in which the Vainshtein mechanism suppresses the scalar field interactions with matter, we study low-density stars with slow rotation and static relativistic stars. We develop an expansion scheme to find approximated solutions inside the Vainshtein radius, and show that deviations from General Relativity (GR), while considering rotation, are also suppressed by the Vainshtein mechanism. In a quadratic coupling model, in which the scalarisation effect can significantly enhance deviations from GR in normal scalar tensor gravity, the Galileon term successfully suppresses the large deviations away from GR. Moreover, using a realistic equation of state, we construct solutions for a relativistic star, and show that deviations from GR are more suppressed for higher density objects. However, we found that the scalar field solution ceases to exist above a critical density, which roughly corresponds to the maximum mass of a neutron star. This indicates that, for a compact object described by a polytropic equation of state, the configuration that would collapse into a black hole cannot support a non-trivial scalar field., Comment: 32 pages, 11 figures, published in JCAP

We study linear perturbations around time dependent spherically symmetric solutions in the Lambda_3 massive gravity theory, which self-accelerate in the vacuum. We find that the dynamics of the scalar perturbations depend on the coordinate choice for the background solutions. For particular choices of coordinates there is a symmetry enhancement, leaving no propagating scalar degrees of freedom at linear order in perturbations. In contrast, any other coordinate choice propagates a single scalar mode. We find that the Hamiltonian of this scalar mode is unbounded from below for all self-accelerating solutions, signalling an instability., Comment: references added; accepted for publication in JCAP

Starting from the general Horndeski action, we derive the most general effective theory for scalar perturbations around flat space that allows us to screen fifth forces via the Vainshtein mechanism. The effective theory is described by a generalization of the Galileon Lagrangian, which we use to study the stability of spherically symmetric configurations exhibiting the Vainshtein effect. In particular, we discuss the phenomenological consequences of a scalar-tensor coupling that is absent in the standard Galileon Lagrangian. This coupling controls the superluminality and stability of fluctuations inside the Vainshtein radius in a way that depends on the density profile of a matter source. Particularly we find that the vacuum solution is unstable due to this coupling., Comment: 6 pages, matches version published by PRD, references added

Massive gravity is a good theoretical laboratory to study modifications of General Relativity. The theory offers a concrete set-up to study models of dark energy, since it admits cosmological self-accelerating solutions in the vacuum, in which the size of the acceleration depends on the graviton mass. Moreover, non-linear gravitational self-interactions, in the proximity of a matter source, manage to mimic the predictions of linearised General Relativity, hence agreeing with solar-system precision measurements. In this article, we review our work in the subject, classifying, on one hand, static solutions, and on the other hand, self-accelerating backgrounds. For what respects static solutions we exhibit black hole configurations, together with other solutions that recover General Relativity near a source via the Vainshtein mechanism. For the self-accelerating solutions we describe a wide class of cosmological backgrounds, including an analysis of their stability., Comment: 17 pages, 4 figures. Invited review to CQG focus issue on massive gravity

We investigate vector contributions to the Lagrangian of $\Lambda_3-$massive gravity in the decoupling limit, the less explored sector of this theory. The main purpose is to understand the stability of maximally symmetric %self-accelerating vacuum solutions. Around self-accelerating configurations, vector degrees of freedom become strongly coupled since their kinetic terms vanish, so their dynamics is controlled by higher order interactions. Even in the decoupling limit, the vector Lagrangian contains an infinite number of terms. We develop a systematic method to covariantly determine the vector Lagrangian at each order in perturbations, fully manifesting the symmetries of the system. We show that, around self-accelerating solutions, the structure of higher order $p$-form Galileons arise, avoiding the emergence of a sixth BD ghost mode. However, a careful analysis shows that there are directions along which the Hamiltonian is unbounded from below. This instability can be interpreted as one of the available fifth physical modes behaving as a ghost. Therefore, we conclude that self-accelerating configurations, in the decoupling limit of $\Lambda_3$-massive gravity, are generically unstable., Comment: 16 pages, 2 figures

We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background metric is flat while the Stuckelberg fields have a non-trivial profile. Around this background, at linear order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these findings., Comment: 9 pages, no figures

We study static, spherically symmetric solutions in a recently proposed ghost-free model of non-linear massive gravity. We focus on a branch of solutions where the helicity-0 mode can be strongly coupled within certain radial regions, giving rise to the Vainshtein effect. We truncate the analysis to scales below the gravitational Compton wavelength, and consider the weak field limit for the gravitational potentials, while keeping all non-linearities of the helicity-0 mode. We determine analytically the number and properties of local solutions which exist asymptotically on large scales, and of local (inner) solutions which exist on small scales. We find two kinds of asymptotic solutions, one of which is asymptotically flat, while the other one is not, and also two types of inner solutions, one of which displays the Vainshtein mechanism, while the other exhibits a self-shielding behaviour of the gravitational field. We analyse in detail in which cases the solutions match in an intermediate region. The asymptotically flat solutions connect only to inner configurations displaying the Vainshtein mechanism, while the non asymptotically flat solutions can connect with both kinds of inner solutions. We show furthermore that there are some regions in the parameter space where global solutions do not exist, and characterise precisely in which regions of the phase space the Vainshtein mechanism takes place., Comment: 21 pages, 7 figures, published version

We explore the possibility of realising self-accelerated expansion of the Universe taking into account the vector components of a massive graviton. The effective action in the decoupling limit contains an infinite number of terms, once the vector degrees of freedom are included. These can be re-summed in physically interesting situations, which result in non-polynomial couplings between the scalar and vector modes. We show there are self-accelerating background solutions for this effective action, with the possibility of having a non-trivial profile for the vector fields. We then study fluctuations around these solutions and show that there is always a ghost, if a background vector field is present. When the background vector field is switched off, the ghost can be avoided, at the price of entering into a strong coupling regime, in which the vector fluctuations have vanishing kinetic terms. Finally we show that the inclusion of a bare cosmological constant does not change the previous conclusions and it does not lead to a ghost mode in the absence of a background vector field., Comment: 23 pages, 2 figures

We investigate strong coupling effects in a covariant massive gravity model, which is a candidate for a ghost free non-linear completion of Fierz-Pauli. We analyse the conditions to recover general relativity via Vainshtein mechanism in the weak field limit, and find three main cases depending on the choice of parameters. In the first case, the potential is such that all non-linearities disappear and the vDVZ discontinuity cannot be avoided. In the second case, the Vainshtein mechanism allows to recover general relativity within a macroscopic radius from a source. In the last case, the strong coupling of the scalar graviton completely shields the massless graviton, and weakens gravity when approaching the source. In the second part of the paper, we explore new exact vacuum solutions, that asymptote de Sitter or anti de Sitter space depending on the parameter choice. The curvature of the space is proportional to the mass of the graviton, thus providing a cosmological background which may explain the present day acceleration in terms of the graviton mass. Moreover, by expressing the potential for non-linear massive gravity in a convenient form, we also suggest possible connections with a higher dimensional framework., Comment: 20 pages, no figures. Version accepted by PRD

We study spherically symmetric solutions in a covariant massive gravity model, which is a candidate for a ghost-free non-linear completion of the Fierz-Pauli theory. There is a branch of solutions that exhibits the Vainshtein mechanism, recovering General Relativity below a Vainshtein radius given by $(r_g m^2)^{1/3}$, where $m$ is the graviton mass and $r_g$ is the Schwarzschild radius of a matter source. Another branch of exact solutions exists, corresponding to Schwarzschild-de Sitter spacetimes where the curvature scale of de Sitter space is proportional to the mass squared of the graviton., Comment: 5 pages, 1 figure, version accepted by PRL

A collision of orbifold planes in eleven dimensions has been proposed as an explanation of the hot big bang. When the two planes are close to each other, the winding membranes become the lightest modes of the theory, and can be effectively described in terms of fundamental strings in a ten dimensional background. Near the brane collision, the eleven-dimensional metric is an Euclidean space times a 1+1-dimensional Milne universe. However, one may expect small perturbations to lead into a more general Kasner background. In this paper we extend the previous classical analysis of winding membranes to Kasner backgrounds, and using the Hamiltonian equations, solve for the wave function of loops with circular symmetry. The evolution across the singularity is regular, and explained in terms of the excitement of higher oscillation modes. We also show there is finite particle production and unitarity is preserved., Comment: 28 pages, 10 figures

By studying perturbations about the vacuum, we show that Horava gravity suffers from two different strong coupling problems, extending all the way into the deep infra-red. The first of these is associated with the principle of detailed balance and explains why solutions to General Relativity are typically not recovered in models that preserve this structure. The second of these occurs even without detailed balance and is associated with the breaking of diffeomorphism invariance, required for anisotropic scaling in the UV. Since there is a reduced symmetry group there are additional degrees of freedom, which need not decouple in the infra-red. Indeed, we use the Stuckelberg trick to show that one of these extra modes become strongly coupled as the parameters approach their desired infra-red fixed point. Whilst we can evade the first strong coupling problem by breaking detailed balance, we cannot avoid the second, whatever the form of the potential. Therefore the original Horava model, and its "phenomenologically viable" extensions do not have a perturbative General Relativity limit at any scale. Experiments which confirm the perturbative gravitational wave prediction of General Relativity, such as the cumulative shift of the periastron time of binary pulsars, will presumably rule out the theory., Comment: 11 pages

We study the evolution of non-periodic cosmic string loops containing Y-junctions, such as may form during the evolution of a network of (p,q) cosmic superstrings. We set up and solve the Nambu-Goto equations of motion for a loop with junctions, focusing attention on a specific static and planar initial loop configuration. After a given time, the junctions collide and the Nambu-Goto description breaks down. We also study the same loop configuration in a U(1)xU(1) field theory model that allows composite vortices with corresponding Y-junctions. We show that the field theory and Nambu-Goto evolution are remarkably similar until the collision time. However, in the field theory evolution a new phenomenon occurs: the composite vortices can unzip, producing in the process new Y-junctions, whose separation may grow significantly, destabilizing the configuration. In particular, an initial loop with two Y-junctions may evolve to a configuration with six Y-junctions (all distant from each other). Setting up this new configuration as an initial condition for Nambu Goto strings, we solve for its evolution and establish conditions under which it is stable to the decay mode seen in the field theory case. Remarkably, the condition closely matches that seen in the field theory simulations, and is expressed in terms of simple parameters of the Nambu-Goto system. This implies that there is an easy way to understand the instability in terms of which region of parameter space leads to stable or unstable unzippings., Comment: 16 pages, 11 figures, typos corrected

We consider a magnetic Bianchi I braneworld, embedded in between two Schwarzschild-AdS spacetimes, boosted equal amounts in opposite directions and compare them to the analagous solution in four-dimensional General Relativity. The efficient dissipation of anisotropy on the brane is explicitly demonstrated, a process we dub braneworld isotropization. From the bulk point of view, we attribute this to anisotropic energy being carried into the bulk by hot gravitons leaving the brane. From the brane point of view this can be interpreted in terms of the production of particles in the dual CFT. We explain how this result enables us to gain a better understanding of the behaviour of anisotropic branes already studied in the literature. We also show how there is evidence of particles being over-produced, and comment on how this may ultimately provide a possible observational signature of braneworlds., Comment: 21 pages, 3 figures

We consider one of the simplest time-dependent backgrounds in M-theory, describing the shrinking away of the M-theory dimension with the other spatial dimensions static. As the M-theory dimension becomes small, the situation becomes well-described by string theory in a singular cosmological background where the string coupling tends to zero but the $\alpha'$-corrections become large, near the cosmic singularity. We compute these $\alpha'$-corrections, both for the background and for linearized perturbations, in heterotic string theory, and show they may be reproduced by a map from eleven-dimensional M-theory., Comment: 22 pages, 4 figures

One of the simplest time-dependent solutions of M theory consists of nine-dimensional Euclidean space times 1+1-dimensional compactified Milne space-time. With a further modding out by Z_2, the space-time represents two orbifold planes which collide and re-emerge, a process proposed as an explanation of the hot big bang. When the two planes are near, the light states of the theory consist of winding M2-branes, describing fundamental strings in a particular ten-dimensional background. They suffer no blue-shift as the M theory dimension collapses, and their equations of motion are regular across the transition from big crunch to big bang. In this paper, we study the classical evolution of fundamental strings across the singularity in some detail. We also develop a simple semi-classical approximation to the quantum evolution which allows one to compute the quantum production of excitations on the string and implement it in a simplified example., Comment: 38 pages, 19 figures