We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (\Sigma)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

In this letter we propose a test to detect the linearity of the dark energy equation of state, and apply it to two different Type Ia Supernova (SN Ia) data sets, Union2.1 and SNLS3. We find that: a. current SN Ia data are well described by a dark energy equation of state linear in the cosmic scale factor a, at least up to a redshift z = 1, independent of the pivot points chosen for the linear relation; b. there is no significant evidence of any deviation from linearity. This apparent linearity may reflect the limit of dark energy information extractable from current SN Ia data.

Hubble Space Telescope imaging of the galaxy cluster Abell 1689 has revealed an exceptional number of strongly lensed multiply-imaged galaxies, including high-redshift candidates. Previous studies have used this data to obtain the most detailed dark matter reconstructions of any galaxy cluster to date, resolving substructures ~25 kpc across. We examine Abell 1689 (hereafter, A1689) non-parametrically, combining strongly lensed images and weak distortions from wider field Subaru imaging, and we incorporate member galaxies to improve the lens solution. Strongly lensed galaxies are often locally affected by member galaxies, however, these perturbations cannot be recovered in grid based reconstructions because the lensing information is too sparse to resolve member galaxies. By adding luminosity-scaled member galaxy deflections to our smooth grid we can derive meaningful solutions with sufficient accuracy to permit the identification of our own strongly lensed images, so our model becomes self consistent. We identify 11 new multiply lensed system candidates and clarify previously ambiguous cases, in the deepest optical and NIR data to date from Hubble and Subaru. Our improved spatial resolution brings up new features not seen when the weak and strong lensing effects are used separately, including clumps and filamentary dark matter around the main halo. Our treatment means we can obtain an objective mass ratio between the cluster and galaxy components, for examining the extent of tidal stripping of the luminous member galaxies. We find a typical mass-to-light ratios of M/L_B = 21 inside the r<1 arcminute region that drops to M/L_B = 17 inside the r<40 arcsecond region. Our model independence means we can objectively evaluate the competitiveness of stacking cluster lenses for defining the geometric lensing-distance-redshift relation in a model independent way.

In this work we present a nonparametric approach, which works on minimal assumptions, to reconstruct the cosmic expansion of the Universe. We propose to combine a locally weighted scatterplot smoothing method and a simulation-extrapolation method. The first one (Loess) is a nonparametric approach that allows to obtain smoothed curves with no prior knowledge of the functional relationship between variables nor of the cosmological quantities. The second one (Simex) takes into account the effect of measurement errors on a variable via a simulation process. For the reconstructions we use as raw data the Union2.1 Type Ia Supernovae compilation, as well as recent Hubble parameter measurements. This work aims to illustrate the approach, which turns out to be a self-sufficient technique in the sense we do not have to choose anything by hand. We examine the details of the method, among them the amount of observational data needed to perform the locally weighted fit which will define the robustness of our reconstruction. In view of our results, we believe that our proposal offers a promising alternative for reconstructing global trends of cosmological data when there is little intuition on the relationship between the variables and we also think it even presents good prospects to generate reliable mock data points where the original sample is poor.

Cosmography provides a model-independent way to map the expansion history of the Universe. In this paper we simulate a Euclid-like survey and explore cosmographic constraints from future Baryonic Acoustic Oscillations (BAO) observations. We derive general expressions for the BAO transverse and radial modes and discuss the optimal order of the cosmographic expansion that provide reliable cosmological constraints. Through constraints on the deceleration and jerk parameters, we show that future BAO data have the potential to provide a model-independent check of the cosmic acceleration as well as a discrimination between the standard $\Lambda$CDM model and alternative mechanisms of cosmic acceleration.

Aims. We investigate the degree of improvement in dark energy constraints that can be achieved by extending Type Ia Supernova (SN Ia) samples to redshifts z > 1.5 with the Hubble Space Telescope (HST), particularly in the ongoing CANDELS and CLASH multi-cycle treasury programs. Methods. Using the popular CPL parametrization of the dark energy, w = w0 +wa(1-a), we generate mock SN Ia samples that can be projected out to higher redshifts. The synthetic datasets thus generated are fitted to the CPL model, and we evaluate the improvement that a high-z sample can add in terms of ameliorating the statistical and systematic uncertainties on cosmological parameters. Results. In an optimistic but still very achievable scenario, we find that extending the HST sample beyond CANDELS+CLASH to reach a total of 28 SN Ia at z > 1.0 could improve the uncertainty in the wa parameter by up to 21%. The corresponding improvement in the figure of merit (FoM) would be as high as 28%. Finally, we consider the use of high-redshift SN Ia samples to detect non-cosmological evolution in SN Ia luminosities with redshift, finding that such tests could be undertaken by future spacebased infrared surveys using the James Webb Space Telescope (JWST).

In the strong lensing regime non-parametric lens models struggle to achieve sufficient angular resolution for a meaningful derivation of the central cluster mass distribution. The problem lies mainly with cluster members which perturb lensed images and generate additional images, requiring high resolution modeling, even though we mainly wish to understand the relatively smooth cluster component. The required resolution is not achievable because the separation between lensed images is several times larger than the deflection angles by member galaxies, even for the deepest data. Here we bypass this limitation by incorporating a simple physical prior for member galaxies, using their observed positions and their luminosity scaled masses. This galaxy contribution is added to a relatively coarse Gaussian pixel grid for modeling the cluster mass distribution, extending our established WSLAP code (Diego et al. 2007). We test this new code with a simulation based on A1689, using the pixels belonging to multiply-lensed images and the observed member galaxies. Dealing with the cluster members this way leads to convergent solutions, without resorting to regularization, reproducing well the input cluster and substructures. We highlight the ability of this method to recover dark sub-components of the cluster, unrelated to member galaxies. Such anomalies can provide clues to the nature of invisible dark matter, but are hard to discover using parametrized models where substructures are defined by the visible data. With our increased resolution and stability we show, for the first time, that non-parametric models can be made sufficiently precise to locate multiply-lensed systems, thereby achieving fully self-consistent solutions without reliance on input systems from less objective means.

The recent observation that the Cosmic Microwave Background (CMB) may prefer a neutrino excess has triggered a number of works studying this possibility. The effect obtained by the non-interacting massless neutrino excess could be mimicked by some extra radiation component in the early universe, such as a cosmological gravitational wave background. Prompted by the fact that a possible candidate to source those gravitational waves would be cosmic strings, we perform a parameter fitting study with models which considers both cosmic strings and the effective number of neutrinos as free parameters, using CMB and non-CMB data. We find that there is a correlation between cosmic strings and the number of extra relativistic species, and that strings can account for all the extra radiation necessary. In fact, CMB data prefer strings at a 2sigma level, paying the price of a higher extra radiation component. CMB data also give a moderate preference for a model with ns=1. The inclusion of non-CMB data lowers both the preference for strings and for the extra relativistic species.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

We show that a large-area imaging survey using narrow-band filters could detect quasars in sufficiently high number densities, and with more than sufficient accuracy in their photometric redshifts, to turn them into suitable tracers of large-scale structure. If a narrow-band optical survey can detect objects as faint as i=23, it could reach volumetric number densities as high as 10^{-4} h^3 Mpc^{-3} (comoving) at z~1.5 . Such a catalog would lead to precision measurements of the power spectrum up to z~3-4. We also show that it is possible to employ quasars to measure baryon acoustic oscillations at high redshifts, where the uncertainties from redshift distortions and nonlinearities are much smaller than at z<1. As a concrete example we study the future impact of J-PAS, which is a narrow-band imaging survey in the optical over 1/5 of the unobscured sky with 42 filters of ~100 A full-width at half-maximum. We show that J-PAS will be able to take advantage of the broad emission lines of quasars to deliver excellent photometric redshifts, \sigma_{z}~0.002(1+z), for millions of objects.

The cosmic microwave background (CMB) is affected by the total radiation density around the time of decoupling. At that epoch, neutrinos comprised a significant fraction of the radiative energy, but there could also be a contribution from primordial gravitational waves with frequencies greater than ~ 10^-15 Hz. If this cosmological gravitational wave background (CGWB) were produced under adiabatic initial conditions, its effects on the CMB and matter power spectrum would mimic massless non-interacting neutrinos. However, with homogenous initial conditions, as one might expect from certain models of inflation, pre big-bang models, phase transitions and other scenarios, the effect on the CMB would be distinct. We present updated observational bounds for both initial conditions using the latest CMB data at small scales from the South Pole Telescope (SPT) in combination with Wilkinson Microwave Anisotropy Probe (WMAP), current measurements of the baryon acoustic oscillations, and the Hubble parameter. With the inclusion of the data from SPT the adiabatic bound on the CGWB density is improved by a factor of 1.7 to 10^6 Omega_gw < 8.7 at the 95% confidence level (C.L.), with weak evidence in favor of an additional radiation component consistent with previous analyses. The constraint can be converted into an upper limit on the tension of horizon-sized cosmic strings that could generate this gravitational wave component, with Gmu < 2 10^-7 at 95% C.L., for string tension Gmu. The homogeneous bound improves by a factor of 3.5 to 10^6 Omega_gw < 1.0 at 95% C.L., with no evidence for such a component from current data.

In this work we introduce two new polynomial parametrizations of dark energy and explore their correlation properties. The parameters to fit are the equation of state values at z=0 and z=0.5, which have naturally low correlation and have already been shown to improve the popular Chevallier-Polarski-Linder (CPL) parametrization. We test our models with low redshift astronomical probes: type Ia supernovae and baryon acoustic oscillations (BAO), in the form of both current and synthetic data. Specifically, we present simulations of measurements of the radial and transversal BAO scales similar to those expected in a BAO high precision spectroscopic redshift survey similar to EUCLID. According to the Bayesian deviance information criterion (DIC), which penalizes large errors and correlations, we show that our models perform better than the CPL re-parametrization proposed by Wang (in terms of z=0 and z=0.5). This is due to the combination of a lower correlation and smaller relative errors. The same holds for a frequentist perspective: our Figure-of-Merit is larger for our parametrizations.

Model independent reconstructions of dark energy have received some attention. The approach that addresses the reconstruction of the dimensionless coordinate distance and its two first derivatives using a polynomial fit in different redshift windows is well developed \cite{DalyDjorgovski1,DalyDjorgovski2,DalyDjorgovski3}. In this work we offer new insights into the problem by focusing on two types of observational probes: SNeIa and GRBs. Our results allow to highlight some of the intrinsic weaknesses of the method. One of the directions we follow is to consider updated observational samples. Our results indicate than conclusions on the main dark energy features as drawn from this method are intimately related to the features of the samples themselves (which are not quite ideal). This is particularly true of GRBs, which manifest themselves as poor performers in this context. In contrast to original works, we conclude they cannot be used for cosmological purposes, and the state of the art does not allow to regard them on the same quality basis as SNeIa. The next direction we contribute to is the question of how the adjusting of some parameters (window width, overlap, selection criteria) affect the results. We find again there is a considerable sensitivity to these features. Then, we try to establish what is the current redshift range for which one can make solid predictions on dark energy evolution. Finally, we strengthen the former view that this model is modest in the sense it provides only a picture of the global trend. But, on the other hand, we believe it offers an interesting complement to other approaches given that it works on minimal assumptions.

Using real and synthetic Type Ia SNe (SNeIa) and baryon acoustic oscillations (BAO) data representing current observations forecasts, this paper investigates the tension between those probes in the dark energy equation of state (EoS) reconstruction considering the well known CPL model and Wang's low correlation reformulation. In particular, here we present simulations of BAO data from both the the radial and transverse directions. We also explore the influence of priors on Omega_m and Omega_b on the tension issue, by considering 1-sigma deviations in either one or both of them. Our results indicate that for some priors there is no tension between a single dataset (either SNeIa or BAO) and their combination (SNeIa+BAO). Our criterion to discern the existence of tension (sigma-distance) is also useful to establish which is the dataset with most constraining power; in this respect SNeIa and BAO data switch roles when current and future data are considered, as forecasts predict and spectacular quality improvement on BAO data. We also find that the results on the tension are blind to the way the CPL model is addressed: there is a perfect match between the original formulation and that by the low correlation optimized, but the errors on the parameters are much narrower in all cases of our exhaustive exploration, thus serving the purpose of stressing the convenience of this reparametrization.

We study the possibility of detecting oscillating patterns in the equation of state (EoS) of the dark energy using different cosmological datasets. We follow a phenomenological approach and study three different oscillating models for the EoS, one of them periodic and the other two damped (proposed here for the first time). All the models are characterised by the amplitude value, the centre and the frequency of oscillations. In contrast to previous works in the literature, we do not fix the value of the frequency to a fiducial value related to the time extension of chosen datasets, but consider a discrete set of values, so to avoid arbitrariness and try and detect any possible time period in the EoS. We test the models using a recent collection of SNeIa, direct Hubble data and Gamma Ray Bursts data. Main results are: I. even if constraints on the amplitude are not too strong, we detect a trend of it versus the frequency, i.e. decreasing (and even negatives) amplitudes for higher frequencies; II. the centre of oscillation (which corresponds to the present value of the EoS parameter) is very well constrained, phantom behaviour is excluded at $1\sigma$ level and trend which is in agreement with the one for the amplitude appears; III. the frequency is hard to constrain, showing similar statistical validity for all the values of the discrete set chosen, but the best fit of all the scenarios considered is associated with a period which is in the redshift range depicted by our cosmological data. The "best" oscillating models are compared with $\Lambda$CDM using dimensionally consistent a Bayesian approach based information criterion and the conclusion reached is the non existence of significant evidence against dark energy oscillations.

We propose a model based on a DBI action for the unification of dark matter and dark energy. This is supported by the results of the study of its background behavior at early and late times, and reinforced by the analysis of the evolution of perturbations. We also perform a Bayesian analysis to set observational constraints on the parameters of the model using type Ia SN, CMB shift and BAO data. Finally, to complete the study we investigate its kinematics aspects, such as the effective equation of state parameter, acceleration parameter and transition redshift. Particularizing those parameters for the best fit one appreciates that an effective phantom is preferred.