We determine the Hubble constant $H_0$ precisely ($2.3\%$ uncertainty) in a manner independent of cosmological model through Gaussian process regression, using strong lensing and supernova data. Strong gravitational lensing of a variable source can provide a time-delay distance $D_{\Delta t}$ and angular diameter distance to the lens $D_{\rm{d}}$. These absolute distances can anchor Type Ia supernovae, which give an excellent constraint on the shape of the distance-redshift relation. Updating our previous results to use the H0LiCOW program's milestone dataset consisting of six lenses, four of which have both $D_{\Delta t}$ and $D_{\rm{d}}$ measurements, we obtain $H_0=72.8_{-1.7}^{+1.6}\rm{\ km/s/Mpc}$ for a flat universe and $H_0=77.3_{-3.0}^{+2.2}\rm{\ km/s/Mpc}$ for a non-flat universe. We carry out several consistency checks on the data and find no statistically significant tensions, though a noticeable redshift dependence persists in a particular systematic manner that we investigate. Speculating on the possibility that this trend of derived Hubble constant with lens distance is physical, we show how this can arise through modified gravity light propagation, which would also impact the weak lensing $\sigma_8$ tension.

Fast Radio Bursts (FRBs) are bright radio transients with millisecond duration at cosmological distances. Since compact dark matter/objects (COs) could act as lenses and cause split of this kind of very short duration signals, Mu$\rm{\tilde{n}}$oz et al. (2016) has proposed a novel method to probe COs with lensing of FRBs. In this Letter, we for the first time apply this method to real data and give constraints of the nature of COs with currently available FRB observations. We emphasize the information from dynamic spectra of FRBs is quite necessary for identifying any lensed signals and find no echoes in the existing data. The null search gives a constraint comparable to that from galactic wide binaries, though the methods of redshift inference from dispersion measure would impact a little. Furthermore, we make an improved forecast basing on the distributions of real data for the ongoing and upcoming telescopes. Finally, we discuss the situation where one or more lensed signals will be detected. In such a case, the parameter space of COs can be pinned down very well since the lens mass can be directly determined through the observed flux ratio and time delay between split images.

Unlike the electromagnetic radiation from astrophysical objects, gravitational waves (GWs) from binary star mergers have much longer wavelengths and are coherent. For ground-based GW detectors, when the lens object between the source and the earth has mass $\sim 1-10^5M_\odot$, the diffraction effect should be considered since the chirping wavelengths are comparable to the scale of the barrier (its Schwarzschild radius). The waveform will thus be distorted as the fringes. In this work, we show that signals from the third-generation GW detectors like the Einstein Telescope (ET) would be a smoking gun for probing the nature of compact dark matter (CDM) or primordial black holes. Detection of the lensing effects becomes harder when the lens mass is smaller. ET is more sensitive than LIGO, the constraint is available for CDM mass $>5M_\odot$ while LIGO can only detect the mass $>100M_\odot$. For a null search of the fringes, one-year observation of ET can constrain the CDM density fraction to $\sim10^{-2}-10^{-5}$ in the mass range $M_{\rm{CDM}}=10M_\odot-100M_\odot$.

General relativity reproduces main current cosmological observations, assuming the validity of cosmic distance duality relation (CDDR) at all scales and epochs. However, CDDR is poorly tested in the redshift interval between the farthest observed Type Ia supernovae (SN Ia) and that of the Cosmic Microwave background (CMB). We present a new idea of testing the validity of CDDR, through the multiple measurements of high-redshift quasars. Luminosity distances are derived from the relation between the UV and X-ray luminosities of quasars, while angular diameter distances are obtained from the compact structure in radio quasars. This will create a valuable opportunity where two different cosmological distances from the same kind of objects at high redshifts are compared. Our constraints are more stringent than other currently available results based on different observational data and show no evidence for the deviation from CDDR at $z\sim 3$. Such accurate model-independent test of fundamental cosmological principles can become a milestone in precision cosmology.

The knowledge about the black hole mass function (BHMF) and its evolution would help to understand the origin of the BHs and how BH binaries formed at different stages of the history of the Universe. We demonstrate the ability of future third generation gravitational wave (GW) detector -- the Einstein Telescope (ET) to infer the slope of the BHMF and its evolution with redshift. We perform the Monte Carlo simulation of the measurements of chirp signals from binary BH systems (BBH) that could be detected by ET, including the BH masses and their luminosity distances ($d_L$). We use the mass of a primary black hole in each binary system to infer the BHMF as a power-law function with slope parameter as $\alpha$. Taking into account the bias that could be introduced by the uncertainty of measurements and by the selection effect, we carried out the numerical tests and find that only one thousand of GW events registered by ET ($\sim1\%$ amount of its yearly detection rate) could accurately infer the $\alpha$ with a precision of $\alpha\sim0.1$. Furthermore, we investigate the validity of our method to recover a scenario where $\alpha$ evolves with redshift as $\alpha(z) = \alpha_0 + \alpha_1\frac{z}{1+z}$. Taking a thousand of GW events and using $d_L$ as the redshift estimator, our tests show that one could infer the value of evolving parameter $\alpha_1$ accurately with the uncertainty level of $\sim0.5$. Our numerical tests verify the reliability of our method. The uncertainty levels of the inferred parameters can be trusted directly for the several sets of the parameter we assumed, yet shouldn't be treated as a universal level for the general case.

We propose a new model-independent strategy to calibrate the distance relation in Type Ia supernova (SN) observations and to probe the intrinsic properties of SNe Ia, especially the absolute magnitude $M_B$, basing on strong lensing observations in the upcoming Large Synoptic Survey Telescope (LSST) era. The strongly lensed quasars can provide the Time Delay Distances (TDD) and the Angular Diameter Distances (ADD) to the lens galaxies. These absolute distance measurements can model-independently anchor the SNe Ia at cosmological distances. We simulated 55 high-quality lensing systems with $5\%$ uncertainties for both TDD and ADD measurements basing on future observation conditions. For the time delay distances and the angular diameter distances as the calibration standards, the calibrated $1\sigma$ uncertainties of $M_{B}$ are approximately 0.24 mag and 0.03 mag, respectively. Besides, we also consider an evolving distance relation, for example, caused by the cosmic opacity. In this case, the $1\sigma$ uncertainties of $M_B$ calibrated with TDD and ADD are approximately 0.31 mag and 0.06 mag, respectively. The results show that the ADD method will be a promising tool for calibrating supernovae.

Strongly lensed quasar systems with time delay measurements provide "time delay distances", which are a combination of three angular diameter distances and serve as powerful tools to determine the Hubble constant $H_0$. However, current results often rely on the assumption of the $\Lambda$CDM model. Here we use a model-independent method based on Gaussian process to directly constrain the value of $H_0$. By using Gaussian process regression, we can generate posterior samples of unanchored supernova distances independent of any cosmological model and anchor them with strong lens systems. The combination of a supernova sample with large statistics but no sensitivity to $H_0$ with a strong lens sample with small statistics but $H_0$ sensitivity gives a precise $H_0$ measurement without the assumption of any cosmological model. We use four well-analyzed lensing systems from the state-of-art lensing program H0LiCOW and the Pantheon supernova compilation in our analysis. Assuming the Universe is flat, we derive the constraint $H_0=72.2 \pm 2.1\,$km/s/Mpc, a precision of $2.9\%$. Allowing for cosmic curvature with a prior of $\Omega_{k}=[-0.2,0.2]$, the constraint becomes $H_0=73.0_{-3.0}^{+2.8}\,$km/s/Mpc.

The cosmic distance duality relation (CDDR) is a fundamental assumption in cosmological studies. Given the redshift $z$, it relates luminosity distance $D^L$ with angular diameter distance $D^A$ through $(1+z)^2D^A/D^L\equiv1$. Many efforts have been devoted to test the CDDR with various observational approaches. However, to the best of our knowledge, those methods are always affected by cosmic opacity which could make the CDDR violated due to the non-conservation of photon number. Such mechanism is more related with astropartical physics. In this work, to directly study the nature of space-time, i.e., to disentangle it from astropartical physics, we propose a new strategy to test the CDDR with strong lensing providing $D^A$ and gravitational wave providing $D^L$. It is known that the propagation of gravitational wave is unaffected by cosmic opacity. We demonstrate distances from optical lensing observation are also opacity-free. These two kinds of distance measurements make it possible to test any metric theories of gravity where photons travel along null geodesics. Our results show the constraints on deviations of the CDDR will be very competitive to current techniques.

Strongly lensed quasars with time-delay measurements are well known to provide the "time-delay distances" $D_{\Delta t}=(1+z_L)D_LD_S/D_{LS}$ and the angular diameter distances to lens galaxies $D_L$. These two kinds of distances give stringent constraints on cosmological parameters. In this work, we explore a different use of time-delay observables: Under the assumption of a flat Universe, strong lensing observations can accurately measure the angular diameter distances to sources $D_S$. The corresponding redshifts of quasars may be up to $z_S\sim4$ according to the forecast. The high-redshift distances would sample the Hubble diagram between SNe Ia and CMB, cosmological-model-independently providing direct information on the evolution of the nature of our Universe, for example, the dark energy Equation-of-State parameter $w(z)$. We apply our method to the existing lensing system SDSS 1206+4332 and get $D_S=2388_{-978}^{+2632}Mpc$ at $z_S=1.789$. We also make a forecast for the era of LSST. The uncertainty of $D_S$ depends on the redshifts of lens and source, the uncertainties of $D_{\Delta t}$ and $D_L$, and the correlation between $D_{\Delta t}$ and $D_L$ as well. Larger correlation would result in tighter $D_S$ determination.

Assuming the $\Lambda$CDM model, the CMB and BAO observations indicate a very flat Universe. Model-independent measurements are therefore worth studying. Time delays measured in lensed quasars provide the time delay distances. When compared with the luminosity distances from Supernova Ia observation, the measurements can provide the curvature information under the Distance Sum Rule of FLRW metric. This method is limited by the low redshifts of SNe Ia. In this work, we propose gravitational waves from the Einstein Telescope as standard sirens which reach higher redshifts covering the redshift range of lensed quasars from Large Synoptic Survey Telescope, could provide much more stringent constraints on the curvature. We first consider a conservative case where only 100 gravitational waves with electromagnetic counterparts are available, the $1\sigma$ uncertainty for the curvature parameter $\Omega_k$ is 0.057. In an optimistic case with 1000 signals available, then $\Omega_k$ uncertainty is 0.027. Combining with SNe Ia from Dark Energy Survey, $\Omega_k$ can be further constrained to 0.027 and 0.018, respectively.

The next generation of ground-based gravitational wave (GW) detectors, e.g. the Einstein Telescope, is expected to observe a significant number of strongly lensed GW events as predicted in many previous papers. However, all these works ignored the impact of the Earth's rotation on this prediction. Multiple lensed images arrive at the Earth at different time, thus the ground-based detector has different responses to the lensed images due to different orientations of the detector relative to the GW source direction. Therefore the amplitudes of the GW signal from different images are modulated appropriately, in addition to the lensing magnification. In order to assess this effect, we performed Monte Carlo simulations to calculate the event rate of lensed GW signals. Our conclusion is that the Earth's rotation has a non-negligible impact on the event rate of lensed GW image. The updated event rates decrease by factors of $\sim40\%, \sim20\%, \sim10\%$, for NS-NS, BH-NS, BH-BH systems respectively.

Gravitational wave predicted by General Relativity is the transverse wave of spatial strain. Several gravitational waveform signals from binary black holes and from a binary neutron star system accompanied by electromagnetic counterparts have been recorded by advanced LIGO and advanced Virgo. In analogy to light, the spatial fringes of diffraction and interference should also exist as the important features of gravitational waves. We propose that observational detection of such fringes could be achieved through gravitational lensing of continuous gravitational waves. The lenses would play the role of the diffraction barriers. Considering peculiar motions of the observer, the lens and the source, the spatial amplitude variation of diffraction or interference fringes should be detectable as an amplitude modulation of monochromatic gravitational signal.

Strong lens time delays have been widely used in cosmological studies, especially to infer $H_0$. The upcoming LSST will provide several hundred well measured time delays from the light curves of lensed quasars. However, due to the inclination of the finite AGN accretion disc and the differential magnification of the coherent temperature fluctuations, the microlensing by the stars can lead to changes in the actual time delay on the light-crossing time scale of the emission region $\sim days$. We first study how this would change the uncertainty of $H_0$ in the LSST era, assuming the microlensing time delays can be well estimated. We adopt 1/3, 1 and 3 days respectively as the typical microlensing time delay uncertainties. The relative uncertainty of $H_0$ will be enlarged to $0.47\%$, $0.51\%$, $0.76\%$, respectively from the one without microlensing impact $0.45\%$. Then, due to the lack of understandings on the quasar models and microlensing patterns, we also test the reliability of the results if one neglects this effect in the analysis. The biases of $H_0$ will be $0.12\%$, $0.22\%$ and $0.70\%$, respectively, suggesting that 1 day is the cut-off for a robust $H_0$ estimate.

The purpose of this work is to investigate the prospects of using the future standard siren data without redshift measurements to constrain cosmological parameters. With successful detections of gravitational wave (GW) signals an era of GW astronomy has begun. Unlike the electromagnetic domain, GW signals allow direct measurements of luminosity distances to the sources, while their redshifts remain to be measured by identifying electromagnetic counterparts. This leads to significant technical problems for almost all possible BH-BH systems. It is the major obstacle to cosmological applications of GW standard sirens. In this paper, we introduce the general framework of using luminosity distances alone for cosmological inference. The idea is to use the prior knowledge of the redshift probability distribution for coalescing sources from the intrinsic merger rates assessed with population synthesis codes. Then the posterior probability distributions for cosmological parameters can be calculated. We demonstrate the performance of our method on the simulated mock data and show that the luminosity distance measurement would enable an accurate determination of cosmological parameters up to $20%$ uncertainty level. We also find that in order to infer $H_0$ to 1% level with flat $\Lambda$CDM model, we need about $10^5$ events.

The cold dark matter scenario of hierarchical large-scale structure formation predicts the existence of abundant subhalos around large galaxies. However, the number of observed dwarf galaxies is far from this theoretical prediction, suggesting that most of the subhalos could be dark or quite faint. Gravitational lensing is a powerful tool to probe the mass distribution directly irrespective of whether it is visible or dark. Time delay anomalies in strongly lensed quasar systems are complementary to flux ratio anomalies in probing dark matter substructure in galaxies. Here we propose that lensed gravitational waves detected by the third-generation ground detectors with quite accurate time delay measurements could be a much better tool for this study than conventional techniques. Combined with good quality images of lensed host galaxies identified by the electromagnetic counterpart measurements, lensed GW signals could make the systematic errors caused by dark matter substructures detectable at several percent levels, depending on their mass functions, internal distribution of subhalos and lensing system configuration.

Model-independent constraints on the spatial curvature are not only closely related to important problems such as the evolution of the Universe and properties of dark energy, but also provide a test of the validity of the fundamental Copernican principle. In this paper, with the distance sum rule in the Friedmann-Lema\^{i}tre-Robertson-Walker metric, we achieve model-independent measurements of the spatial curvature from the latest type Ia supernovae and strong gravitational lensing (SGL) observations. We find that a spatially closed Universe is preferred. Moreover, by considering different kinds of velocity dispersion and subsample, we study possible factors which might affect model-independent estimations for the spatial curvature from SGL observations. It is suggested that the combination of observational data from different surveys might cause a systematic bias and the tension between the spatially flat Universe and SGL observations is alleviated when the subsample only from the Sloan Lens ACS Survey is used or a more complex treatment for the density profile of lenses is considered.

The standard siren approach of gravitational wave cosmology appeals to the direct luminosity distance estimation through the waveform signals from inspiralling double compact binaries, especially those with electromagnetic counterparts providing redshifts. It is limited by the calibration uncertainties in strain amplitude and relies on the fine details of the waveform. The Einstein Telescope is expected to produce $10^4-10^5$ gravitational wave detections per year, $50-100$ of which will be lensed. Here we report a waveform-independent strategy to achieve precise cosmography by combining the accurately measured time delays from strongly lensed gravitational wave signals with the images and redshifts observed in the electromagnetic domain. We demonstrate that just 10 such systems can provide a Hubble constant uncertainty of $0.68\%$ for a flat Lambda Cold Dark Matter universe in the era of third generation ground-based detectors.

We present a new model-independent strategy for testing the Friedmann-Lema\^{\i}tre-Robertson-Walker metric and constraining cosmic curvature, based on future time delay measurements of strongly lensed quasar-elliptical galaxy systems from the Large Synoptic Survey Telescope and supernova observations from the Dark Energy Survey. The test only relies on geometric optics. It is independent of the energy contents of the universe and the validity of the Einstein equation on cosmological scales. The study comprises two levels: testing the FLRW metric through the Distance Sum Rule and determining/constraining cosmic curvature. We propose an effective and efficient (redshift) evolution model for performing the former test, which allows us to concretely specify the violation criterion for the FLRW Distance Sum Rule. If the FLRW metric is consistent with the observations, then, on the second level, the cosmic curvature parameter will be constrained to $\sim0.057$ or $\sim0.041$ ($1\sigma$), depending on the availability of high-redshift supernovae, much more stringent than current model-independent techniques. We also show that the bias in the time delay method might be well controlled, leading to robust results. The proposed method is a new independent tool for both testing the fundamental assumptions of homogeneity and isotropy in cosmology and for determining cosmic curvature. It is complementary to cosmic microwave background plus baryon acoustic oscillation analyses, which normally assume a cosmological model with dark energy domination in the late-time universe.

We propose a new model-independent measurement strategy for the propagation speed of gravitational waves (GWs) based on strongly lensed GWs and their electromagnetic (EM) counterparts. This can be done in two ways: by comparing arrival times of GWs and their EM counterparts and by comparing the time delays between images seen in GWs and their EM counterparts. The lensed GW-EM event is perhaps the best way to identify an EM counterpart. Conceptually, this method does not rely on any specific theory of massive gravitons or modified gravity. Its differential setting (i.e., measuring the difference between time delays in GW and EM domains) makes it robust against lens modeling details (photons and GWs travel in the same lensing potential) and against internal time delays between GW and EM emission acts. It requires, however, that the theory of gravity is metric and predicts gravitational lensing similar to general relativity. We expect that such a test will become possible in the era of third-generation gravitational-wave detectors, when about 10 lensed GW events would be observed each year. The power of this method is mainly limited by the timing accuracy of the EM counterpart, which for kilonovae is around $10^4$ s.Thisuncertaintycanbesuppressedbya factor of $\sim 10^{10}$, if strongly lensed transients of much shorter duration associated with the GW event can be identified. Candidates for such short transients include short gamma-ray burst and fast radio bursts.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW ($H_0$ Lenses in COSMOGRAIL's Wellspring), a program that aims to measure $H_0$ with $<3.5\%$ uncertainty from five lens systems (B1608+656, RXJ1131-1231, HE0435-1223, WFI2033-4723 and HE1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.5$ tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

Model-independent estimations for the spatial curvature not only provide a test for the fundamental Copernican principle assumption, but also can effectively break the degeneracy between curvature and dark energy properties. In this paper, we propose to achieve model-independent constraints on the spatial curvature from observations of standard candles and standard clocks, without assuming any fiducial cosmology and other priors. We find that, for the popular Union2.1 type Ia supernovae (SNe Ia ) observations, the spatial curvature is constrained to be $\Omega_K=-0.045_{-0.172}^{+0.176}$. For the latest joint light-curve analysis (JLA) of SNe Ia observations, we obtain $\Omega_K=-0.140_{-0.158}^{+0.161}$. It is suggested that these results are in excellent agreement with the spatially flat Universe. Moreover, compared to other approaches aiming for model-independent estimations of spatial curvature, this method also achieves constraints with competitive precision.

The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.

Under very general assumptions of metric theory of spacetime, photons traveling along null geodesics and photon number conservation, two observable concepts of cosmic distance, i.e. the angular diameter and the luminosity distances are related to each other by the so-called distance duality relation (DDR) $D^L=D^A(1+z)^2$. Observational validation of this relation is quite important because any evidence of its violation could be a signal of new physics. In this paper we introduce a new method to test DDR based on strong gravitational lensing systems and type Ia supernovae under a flat universe. The method itself is worth attention, because unlike previously proposed techniques, it does not depend on all other prior assumptions concerning the details of cosmological model. We tested it using a new compilation of strong lensing systems and JLA compilation of type Ia supernovae and found no evidence of DDR violation. For completeness, we also combined it with previous cluster data and showed its power on constraining DDR. It could become a promising new probe in the future in light of forthcoming massive strong lensing surveys and because of expected advances in galaxy cluster modlelling.

We present our study on cosmic opacity, which relates to changes in photon number as photons travel from the source to the observer. Cosmic opacity may be caused by absorption/scattering due to matter in the universe, or by extragalactic magnetic fields that can turn photons into unobserved particles (e.g. light axions, chameleons, gravitons, Kaluza-Klein modes), and it is crucial to correctly interpret astronomical photometric measurements like type Ia supernovae observations. On the other hand, the expansion rate at different epochs, i.e. the observational Hubble parameter data $H(z)$, are obtained from differential ageing of passively evolving galaxies or from baryon acoustic oscillations and thus are not affected by cosmic opacity. In this work, we first construct opacity-free luminosity distances from $H(z)$ determinations, taking correlations between different redshifts into consideration for our error analysis. Moreover, we let the light-curve fitting parameters, accounting for distance estimation in type Ia supernovae observations, free to ensure that our analysis is authentically cosmological-model-independent and gives a robust result. Any non-zero residuals between these two kinds of luminosity distances can be deemed as an indication of the existence of cosmic opacity. While a transparent universe is currently consistent with the data, our results show that strong constraints on opacity (and consequently on physical mechanisms that could cause it) can be obtained in a cosmological-model-independent fashion.

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $\gamma'$ of the total mass density profile $\rho_{tot}\propto r^{-\gamma'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $\gamma'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

The time delays between point-like images in gravitational lens systems can be used to measure cosmological parameters. The number of lenses with measured time delays is growing rapidly; the upcoming \emph{Large Synoptic Survey Telescope} (LSST) will monitor $\sim10^3$ strongly lensed quasars. In an effort to assess the present capabilities of the community to accurately measure the time delays, and to provide input to dedicated monitoring campaigns and future LSST cosmology feasibility studies, we have invited the community to take part in a "Time Delay Challenge" (TDC). The challenge is organized as a set of "ladders," each containing a group of simulated datasets to be analyzed blindly by participating teams. Each rung on a ladder consists of a set of realistic mock observed lensed quasar light curves, with the rungs' datasets increasing in complexity and realism. The initial challenge described here has two ladders, TDC0 and TDC1. TDC0 has a small number of datasets, and is designed to be used as a practice set by the participating teams. The (non-mandatory) deadline for completion of TDC0 was the TDC1 launch date, December 1, 2013. The TDC1 deadline was July 1 2014. Here we give an overview of the challenge, we introduce a set of metrics that will be used to quantify the goodness-of-fit, efficiency, precision, and accuracy of the algorithms, and we present the results of TDC0. Thirteen teams participated in TDC0 using 47 different methods. Seven of those teams qualified for TDC1, which is described in the companion paper II.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

We propose a new consistent method to test of the distance-duality (DD) relation which related angular diameter distances (DA) to the luminosity distances (DL) in a cosmology-independent way. In order to avoid any bias brought by redshift incoincidence between galaxy clusters and Type Ia Supernave (SNe Ia), as well as to ensure the integrity of the galaxy clusters samples, we obtain the luminosity distance of a certain SN Ia point at the same redshift of the corresponding galaxy cluster by interpolating from the nearby SNe Ia. With the observational data at the same redshifts of the angular diameter distances from the complete 38 galaxy cluster sample for the spherical model and the corrected luminosity distances interpolated from the Union2 set, we find that $\eta \equiv {D_L}{(1+z)}^{-2}/{D_A}=1$ is satisfied within $2\sigma$ confidence level for various parameterizations of $\eta(z)$, which are more stringent than previous testing results without considering redshift bias.

In this paper, we quantify the ability of a future measurement of the Sandage-Loeb test signal from the Cosmic-Dynamic-Experiment-like spectrograph to constrain some popular modified gravity theories including DGP brane-world scenario, $f(R)$ modified gravity and $f(T)$ gravity theory. We find that the Sandage-Loeb test measurements are able to break degeneracies between model parameters markedly and thus greatly improve cosmological constraints for all concerned modified gravity theories when combined with the latest observations of CMB--shift parameter. However, they yield almost the same degeneracy directions between model parameters as that from the distance ratio data derived from the latest observations of the cosmic microwave background and baryonic acoustic oscillations (CMB/BAO). Moreover, for the $f(R)$ modified gravity, the Sandage-Loeb test could provide completely different bounded regions in model parameters space as compared to CMB/BAO and thus supplement strong complementary constraints.

We use the newly published 28 observational Hubble parameter data ($H(z)$) and current largest SNe Ia samples (Union2.1) to test whether the universe is transparent. Three cosmological-model-independent methods (nearby SNe Ia method, interpolation method and smoothing method) are proposed through comparing opacity-free distance modulus from Hubble parameter data and opacity-dependent distance modulus from SNe Ia . Two parameterizations, $\tau(z)=2\epsilon z$ and $\tau(z)=(1+z)^{2\epsilon}-1$ are adopted for the optical depth associated to the cosmic absorption. We find that the results are not sensitive to the methods and parameterizations. Our results support a transparent universe.

In this Letter, we propose a new generalized Ricci dark energy (NGR) model to unify Ricci dark energy (RDE) and XCDM. Our model can distinguish between RDE and XCDM by introducing a parameter $\beta$ called weight factor. When $\beta=1$, NGR model becomes the usual RDE model. The XCDM model is corresponding to $\beta=0$. Moreover, NGR model permits the situation where neither $\beta=1$ nor $\beta=0$. We then perform a statefinder analysis on NGR model to see how $\beta$ effects the trajectory on the $r-s$ plane. In order to know the value of $\beta$, we constrain NGR model with latest observations including type Ia supernovae (SNe Ia) from Union2 set (557 data), baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) galaxy sample and cosmic microwave background (CMB) observation from the 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) results. With Markov Chain Monte Carlo (MCMC) method, the constraint result is $\beta$=$0.08_{-0.21}^{+0.30}(1\sigma)_{-0.28}^{+0.43}(2\sigma)$, which manifests the observations prefer a XCDM universe rather than RDE model. It seems RDE model is ruled out in NGR scenario within $2\sigma$ regions. Furthermore, we compare it with some of successful cosmological models using AIC information criterion. NGR model seems to be a good choice for describing the universe.

We use the new gamma-ray bursts (GRBs) data, combined with the baryon acoustic oscillation(BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release, the newly obtained $A$ parameter at $z=0.6$ from the WiggleZ Dark Energy Survey, the cosmic microwave background (CMB) observations from the 7-Year Wilkinson Microwave Anisotropy Probe (WMAP7) results, and the type Ia supernovae (SNeIa) from Union2 set, to constrain a phenomenological model describing possible interactions between dark energy and dark matter, which was proposed to alleviate the coincidence problem of the standard $\Lambda$CDM model. By using the Markov Chain Monte Carlo (MCMC) method, we obtain the marginalized $1\sigma$ constraints $\Omega_{m}=0.2886\pm{0.0135}$, $r_m=-0.0047\pm{0.0046}$, and $w_X=-1.0658\pm{0.0564}$. We also consider other combinations of these data for comparison. These results show that: (1) the energy of dark matter is slightly transferring to that of dark energy; (2) even though the GRBs+BAO+CMB data present less stringent constraints than SNe+BAO+CMB data do, the GRBs can help eliminate the degeneracies among parameters.

We use the latest data to investigate observational constraints on the new generalized Chaplygin gas (NGCG) model. Using the Markov Chain Monte Carlo (MCMC) method, we constrain the NGCG model with the type Ia supernovae (SNe Ia) from Union2 set (557 data), the usual baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) galaxy sample, the cosmic microwave background (CMB) observation from the 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) results, the newly revised $H(z)$ data, as well as a value of $\theta_{BAO} (z=0.55) = (3.90 \pm 0.38)^{\circ}$ for the angular BAO scale. The constraint results for NGCG model are $\omega_X = -1.0510_{-0.1685}^{+0.1563}(1\sigma)_{-0.2398}^{+0.2226}(2\sigma)$, $\eta = 1.0117_{-0.0502}^{+0.0469}(1\sigma)_{-0.0716}^{+0.0693}(2\sigma)$, and $\Omega_X = 0.7297_{-0.0276}^{+0.0229}(1\sigma)_{-0.0402}^{+0.0329}(2\sigma)$, which give a rather stringent constraint. From the results, we can see a phantom model is slightly favored and the probability that energy transfers from dark matter to dark energy is a little larger than the inverse.

f(R) gravity is thought to be an alternative to dark energy which can explain the acceleration of the universe. It has been tested by different observations including type Ia supernovae (SNIa), the cosmic microwave background (CMB), the baryon acoustic oscillations (BAO) and so on. In this Letter, we use the Hubble constant independent ratio between two angular diameter distances $D=D_{ls}/D_s$ to constrain f(R) model in Palatini approach $f(R)=R-\alpha H^2_0(-\frac{R}{H^2_0})^\beta$. These data are from various large systematic lensing surveys and lensing by galaxy clusters combined with X-ray observations. We also combine the lensing data with CMB and BAO, which gives a stringent constraint. The best-fit results are $(\alpha,\beta)=(-1.50,0.696)$ or $(\Omega_m,\beta)=(0.0734,0.696)$ using lensing data only. When combined with CMB and BAO, the best-fit results are $(\alpha,\beta)=(-3.75,0.0651)$ or $(\Omega_m,\beta)=(0.286,0.0651)$. If we further fix $\beta=0$ (corresponding to $\Lambda$CDM), the best-fit value for $\alpha$ is $\alpha$=$-4.84_{-0.68}^{+0.91}(1\sigma)_{-0.98}^{+1.63}(2\sigma)$ for the lensing analysis and $\alpha$=$-4.35_{-0.16}^{+0.18}(1\sigma)_{-0.25}^{+0.3}(2\sigma)$ for the combined data, respectively. Our results show that $\Lambda$CDM model is within 1$\sigma$ range.

We constrain a unified dark matter (UDM) model from the latest observational data. This model assumes that the dark sector is degenerate. Dark energy and dark matter are the same component. It can be described by an affine equation of state $P_X= p_0 +\alpha \rho_X$. Our data set contains the newly revised $H(z)$ data, type Ia supernovae (SNe Ia) from Union2 set, baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) galaxy sample, as well as the cosmic microwave background (CMB) observation from the 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) results. By using the Markov Chain Monte Carlo (MCMC) method, we obtain the results in a flat universe: $\Omega_\Lambda$=$0.719_{-0.0305}^{+0.0264}(1\sigma)_{-0.0458}^{+0.0380}(2\sigma)$, $\alpha$=$1.72_{-4.79}^{+3.92}(1\sigma)_{-7.30}^{+5.47}(2\sigma)(\times10^{-3})$, $\Omega_bh^2$=$0.0226_{-0.0011}^{+0.0011}(1\sigma)_{-0.0015}^{+0.0016}(2\sigma)$. Moreover, when considering a non-flat universe, $\Omega_\Lambda$=$0.722_{-0.0447}^{+0.0362}(1\sigma)_{-0.0634}^{+0.0479}(2\sigma)$, $\alpha$=$0.242_{-0.775}^{+0.787}(1\sigma)_{-1.03}^{+1.10}(2\sigma)(\times10^{-2})$, $\Omega_bh^2$=$0.0227_{-0.0014}^{+0.0015}(1\sigma)_{-0.0018}^{+0.0021}(2\sigma)$, $\Omega_k$=$-0.194_{-1.85}^{+2.02}(1\sigma)_{-2.57}^{+2.75}(2\sigma)(\times10^{-2})$. These give a more stringent results than before. We also give the results from other combinations of these data for comparison. The observational Hubble parameter data can give a more stringent constraint than SNe Ia. From the constraint results, we can see the parameters $\alpha$ and $\Omega_k$ are very close to zero, which means a flat universe is strongly supported and the speed of sound of the dark sector seems to be zero.