In this paper, we present a cosmological model-independent test for the cosmic opacity at high redshifts ($z\sim5$). We achieve this with the opacity-dependent luminosity distances derived from nonlinear relation between X-ray and UV emissions of quasars, combined with two types of opacity-independent luminosity distances derived from the Hubble parameter measurements and simulated gravitational wave (GW) events achievable with the Einstein Telescope (ET). In the framework of two phenomenological parameterizations adopted to describe cosmic opacity at high redshifts, our main results show that a transparent universe is supported by the current observational data at 2$\sigma$ confidence level. However, the derived value of the cosmic opacity is slightly sensitive to the parametrization of $\tau(z)$, which highlights the importance of choosing a reliable parametrization to describe the optical depth $\tau(z)$ in the early universe. Compared with the previous works, the combination of the quasar data and the $H(z)$/GW observations in similar redshift ranges provides a novel way to confirm a transparent universe ($\epsilon=0$ at higher redshifts $z\sim 5$), with an accuracy of $\Delta \epsilon\sim 10^{-2}$. More importantly, our findings indicate that a strong degeneracy between the cosmic opacity parameter and the parameters characterizing the $L_{UV}-L_X$ relation of quasars, which reinforces the necessity of proper calibration for such new type of high-redshift standard candle (in a cosmological model-independent way).

The cosmic curvature, a fundamental parameter for cosmology could hold deep clues to inflation and cosmic origins. We propose an improved model-independent method to constrain the cosmic curvature by combining the constructed Hubble diagram of high-redshift quasars with galactic-scale strong lensing systems expected to be seen by the forthcoming LSST survey. More specifically, the most recent quasar data are used as a new type of standard candles in the range $0.036<z<5.100$, whose luminosity distances can be directly derived from the non-linear relation between X-ray and UV luminosities. Compared with other methods, the proposed one involving the quasar data achieves constraints with higher precision ($\Delta \Omega_k\sim 10^{-2}$) at high redshifts ($z\sim 5.0$). We also investigate the influence of lens mass distribution in the framework of three types of lens models extensively used in strong lensing studies (SIS model, power-law spherical model, and extended power-law lens model), finding the strong correlation between the cosmic curvature and the lens model parameters. When the power-law mass density profile is assumed, the most stringent constraint on the cosmic curvature $\Omega_k$ can be obtained. Therefore, the issue of mass density profile in the early-type galaxies is still a critical one that needs to be investigated further.

Emergence of gravitational wave (GW) astronomy revived the interest in exploring the low frequency GW spectrum inaccessible from the ground. Satellite GW observatory DECIGO in its original configuration and the currently proposed smaller scale B-DECIGO are aimed to cover deci-Hertz part of the GW spectrum, which fills the gap between LISA mili-Hertz and deca- to kilo-Hertz range probed by ground-based detectors. In this paper we forecast the detection rates of inspiraling double compact objects (DCOs) and the unresolved confusion noise from these sources in DECIGO and B-DECIGO. In the context of DECIGO we use, for the first time, the population synthesis intrinsic inspiral rates of NS-NS, BH-NS and BH-BH systems. We also estimate the expected gravitational lensing rates of such sources for DECIGO and B-DECIGO. The result is that yearly detection of DCOs inspirals for the DECIGO is of order of $10^4 - 10^5$, while for a much smaller scale B-DECIGO they are about $10^2 - 10^4$ depending on the DCO population considered. Taking into account that considerable part of these events would be detectable by ground-based GW observatories the added value of DECIGO/B-DECIGO could be substantial. DECIGO will also be able to register $5-50$ lensed NS-NS inspirals and up to O(100) BH-BH inspirals. On the contrary, predictions for the B-DECIGO are pessimistic: only the optimistic common envelope scenario of DCO formation could yield a few lensed BH-BH inspirals per year. We also address the question of the magnification bias in the GW event catalogs of DECIGO and B-DECIGO.

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.

The assumptions of large-scale homogeneity and isotropy underly the familiar Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) metric that appears to be an accurate description of our Universe. In this paper, we propose a new strategy of testing the validity of the FLRW metric, based on the galactic-scale lensing systems where strongly lensed gravitational waves and their electromagnetic counterparts can be simultaneously detected. Each strong lensing system creates opportunity to infer the curvature parameter of the Universe. Consequently, combined analysis of many such systems will provide a model-independent tool to test the validity of the FLRW metric. Our study demonstrates that the third-generation ground based GW detectors, like the Einstein Telescope (ET) and space-based detectors, like the Big Bang Observer (BBO), are promising concerning determination of the curvature parameter or possible detection of deviation from the FLRW metric. Such accurate measurements of the FLRW metric can become a milestone in precision GW cosmology.

The gas depletion factor $\gamma(z)$, i.e., the average ratio of the gas mass fraction to the cosmic mean baryon fraction of galaxy clusters, plays a very important role in the cosmological application of the gas mass fraction measurements. In this paper, using the newest catalog of 182 galaxy clusters detected by the Atacama Cosmology Telescope (ACT) Polarization experiment, we investigate the possible redshift evolution of $\gamma(z)$ through a new cosmology-independent method. The method is based on non-parametric reconstruction using the measurements of Hubble parameters from cosmic chronometers. Unlike hydrodynamical simulations suggesting constant depletion factor, our results reveal the trend of $\gamma(z)$ decreasing with redshift. This result is supported by a parametric model fit as well as by calculations on the reduced ACTPol sample and on the alternative sample of 91 SZ clusters reported earlier in ACT compilation. Discussion of possible systematic effects leaves an open question about validity of the empirical relation $M_{tot}$-$f_{gas}$ obtained on very close clusters. These results might pave the way to explore the hot gas fraction within large radii of galaxy clusters as well as its possible evolution with redshift, which should be studied further on larger galaxy cluster samples in the upcoming X-ray/SZ cluster surveys.

We present a new idea of testing the validity of the Friedman-Lema\^{\i}tre-Robertson-Walker metric, through the multiple measurements of galactic-scale strong gravitational lensing systems with type Ia supernovae in the role of sources. Each individual lensing system will provide a model-independent measurement of the spatial curvature parameter referring only to geometrical optics independently of the matter content of the universe. This will create a valuable opportunity to test the FLRW metric directly. Our results show that with hundreds of strongly lensed SNe Ia observed by LSST, one would produce robust constraints on the spatial curvature with accuracy $\Delta \Omega_k=0.04$ comparable to the Planck 2015 results.

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.

In this paper, by using the recently compiled set of 120 intermediate-luminosity quasars (ILQSO) observed in a single-frequency VLBI survey, we propose an improved model-independent method to probe cosmic curvature parameter $\Omega_k$ and make the first measurement of the cosmic curvature referring to a distant past, with redshifts up to $z\sim 3.0$. Compared with other methods, the proposed one involving the quasar data achieves constraints with higher precision in this redshift range. More importantly, our results indicate that the measured $\Omega_k$ is in good agreement with zero cosmic curvature, implying that there is no significant deviation from a flat Universe. Finally, we investigate the possibility of testing $\Omega_k$ with a much higher accuracy using quasars observed in the future VLBI surveys. It is shown that our method could provide a reliable and tight constraint on the prior $\Omega_k$ and one can expect the zero cosmic curvature to be established at the precision of $\Delta\Omega_k\sim 10^{-2}$ with 250 well-observed radio quasars.

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.

Testing the distance-sum-rule in strong lensing systems provides an interesting method to determine the curvature parameter $\Omega_k$ using more local objects. In this paper, we apply this method to a quite recent data set of strong lensing systems in combination with intermediate-luminosity quasars calibrated as standard rulers. In the framework of three types of lens models extensively used in strong lensing studies (SIS model, power-law spherical model, and extended power-law lens model), we show that the assumed lens model has a considerable impact on the cosmic curvature constraint, which is found to be compatible or marginally compatible with the flat case (depending on the lens model adopted). Analysis of low, intermediate and high-mass sub-samples defined according to the lens velocity dispersion demonstrates that, although it is not reasonable to characterize all lenses with a uniform model, such division has little impact on cosmic curvature inferred. Finally, thinking about future when massive surveys will provide their yields, we simulated a mock catalog of strong lensing systems expected to be seen by the LSST, together with a realistic catalog of quasars. We found that with about 16000 such systems, combined with the distance information provided by 500 compact milliarcsecond radio sources seen in future radio astronomical surveys, one would be able to constrain the cosmic curvature with an accuracy of $\Delta \Omega_k\simeq 10^{-3}$, which is comparable to the precision of \textit{Planck} 2015 results.

Probing the speed of light is as an important test of General Relativity but the measurements of $c$ using objects in the distant universe have been almost completely unexplored. In this letter, we propose an idea to use the multiple measurements of galactic-scale strong gravitational lensing systems with type Ia supernova acting as background sources to estimate the speed of light. This provides an original method to measure the speed of light using objects located at different redshifts which emitted their light in a distant past. Moreover, we predict that strongly lensed SNe Ia observed by the LSST would produce robust constraints on $\Delta c/c$ at the level of $10^{-3}$. We also discuss whether the future surveys such as LSST may succeed in detecting any hypothetical variation of $c$ predicted by theories in which fundamental constants have dynamical nature.

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.

In this paper, we use multi-frequency angular size measurements of 58 intermediate-luminosity quasars reaching the redshifts $z\sim 3$ and demonstrate that they can be used as standard rulers for cosmological inference. Our results indicate that, for the majority of radio-sources in our sample their angular sizes are inversely proportional to the observing frequency. From the physical point of view it means that opacity of the jet is governed by pure synchrotron self-absorption, i.e. external absorption does not play any significant role in the observed angular sizes at least up to 43 GHz. Therefore, we use the value of the intrinsic metric size of compact milliarcsecond radio quasars derived in a cosmology independent manner from survey conducted at 2 GHz and rescale it properly according to predictions of the conical jet model. This approach turns out to work well and produce quite stringent constraints on the matter density parameter $\Omega_m$ in the flat $\Lambda$CDM model and Dvali-Gabadadze-Porrati braneworld model. The results presented in this paper pave the way for the follow up engaging multi-frequency VLBI observations of more compact radio quasars with higher sensitivity and angular resolution.

As an important candidate gravity theory alternative to dark energy, a class of $f(R)$ modified gravity, which introduces a perturbation of the Ricci scalar $R$ in the Einstein-Hilbert action, has been extensively applied to cosmology to explain the acceleration of the universe. In this paper, we focus on the recently-released VLBI observations of the compact structure in intermediate-luminosity quasars combined with the angular-diameter-distance measurements from galaxy clusters, which consists of 145 data points performing as individual cosmological standard rulers in the redshift range $0.023\le z\le 2.80$, to investigate observational constraints on two viable models in $f(R)$ theories within the Palatini formalism: $f_1(R)=R-\frac{a}{R^b}$ and $f_2(R)=R-\frac{aR}{R+ab}$. We also combine the individual standard ruler data with the observations of CMB and BAO, which provides stringent constraints. Furthermore, two model diagnostics, $Om(z)$ and statefinder, are also applied to distinguish the two $f(R)$ models and $\Lambda$CDM model. Our results show that (1) The quasars sample performs very well to place constraints on the two $f(R)$ cosmologies, which indicates its potential to act as a powerful complementary probe to other cosmological standard rulers. (2) The $\Lambda$CDM model, which corresponds to $b=0$ in the two $f(R)$ cosmologies is still included within $1\sigma$ range. However, there still exists some possibility that $\Lambda$CDM may not the best cosmological model preferred by the current high-redshift observations. (3) The information criteria indicate that the cosmological constant model is still the best one, while the $f_1(R)$ model gets the smallest observational support. (4) The $f_2(R)$ model, which evolves quite different from $f_1(R)$ model at early times, still significantly deviates from both $f_1(R)$ and $\Lambda$CDM model at the present time.

Direct measurements of Hubble parameters $H(z)$ are very useful for cosmological model parameters inference. Based on them, Sahni, Shafieloo and Starobinski introduced a two-point diagnostic $Omh^2(z_i, z_j)$ as an interesting tool for testing the validity of the $\Lambda$CDM model. Applying this test they found a tension between observations and predictions of the $\Lambda$CDM model. We use the most comprehensive compilation $H(z)$ data from baryon acoustic oscillations (BAO) and differential ages (DA) of passively evolving galaxies to study cosmological models using the Hubble parameters itself and to distinguish whether $\Lambda$CDM model is consistent with the observational data with statistical analysis of the corresponding $Omh^2(z_i, z_j)$ two-point diagnostics. Our results show that presently available $H(z)$ data significantly improve the constraints on cosmological parameters. The corresponding statistical $Omh^2(z_i, z_j)$ two-point diagnostics seems to prefer the quintessence with $w>-1$ over the $\Lambda$CDM model. Better and more accurate prior knowledge of the Hubble constant, will considerably improve the performance of the statistical $Omh^2(z_i, z_j)$ method.

In this paper, we propose a new parametrization of $Om(z)$ diagnostics and show how the most recent and significantly improved observations concerning the $H(z)$ and SN Ia measurements can be used to probe the consistency or tension between $\Lambda$CDM model and observations. Our results demonstrates that $H_0$ plays a very important role in the consistency test of $\Lambda$CDM with the $H(z)$ data. Adopting the Hubble constant priors from \textit{Planck} 2013 and Riess (2016), one finds a considerable tension between the current $H(z)$ data and $\Lambda$CDM model and confirms the conclusions obtained previously by the others. However, with the Hubble constant prior taken from WMAP9, the discrepancy between $H(z)$ data and $\Lambda$CDM disappears, i.e., the current $H(z)$ observations still support the cosmological constant scenario. This conclusion is also supported by the results derived from the JLA SNe Ia sample. The best-fit Hubble constant from the combination of $H(z)$+JLA ($H_0=68.81^{+1.50}_{-1.49}$ km/s/Mpc) is well consistent with the results derived both by Planck 2013 and WMAP9, which is significantly different from the recent local measurement by Riess (2016).

A new compilation of $120$ angular-size/redshift data for compact radio quasars from very-long-baseline interferometry (VLBI) surveys motivates us to revisit the interaction between dark energy and dark matter with these probes reaching high redshifts $z\sim 3.0$. In this paper, we investigate observational constraints on different phenomenological interacting dark energy (IDE) models with the intermediate-luminosity radio quasars acting as individual standard rulers, combined with the newest BAO and CMB observation from Planck results acting as statistical rulers. The results obtained from the MCMC method and other statistical methods including Figure of Merit and Information Criteria show that: (1) Compared with the current standard candle data and standard clock data, the intermediate-luminosity radio quasar standard rulers , probing much higher redshifts, could provide comparable constraints on different IDE scenarios. (2) The interaction between dark energy and dark matter seems to be vanishing or slightly smaller than zero. At the 68.3\% confidence level, the energy is seen transferred from dark matter to dark energy, which implies that those IDE models can not alleviate the coincidence problem or even more sever. However, the strong degeneracies between the interaction term and Hubble constant may contribute to alleviate the tension of $H_0$ between the recent Planck and HST measurements. (3) Concerning the ranking of competing dark energy models, IDE with more free parameters are substantially penalized by the BIC criterion, which agrees very well with the previous results derived from other cosmological probes.

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.

In this paper, we place constraints on four alternative cosmological models under the assumption of the spatial flatness of the Universe: CPL, EDE, GCG and MPC. A new compilation of 120 compact radio quasars observed by very-long-baseline interferometry, which represents a type of new cosmological standard rulers, are used to test these cosmological models. Our results show that the fits on CPL obtained from the quasar sample are well consistent with those obtained from BAO. For other cosmological models considered, quasars provide constraints in agreement with those derived with other standard probes at $1\sigma$ confidence level. Moreover, the results obtained from other statistical methods including Figure of Merit, $Om(z)$ and statefinder diagnostics indicate that: (1) Radio quasar standard ruler could provide better statistical constraints than BAO for all cosmological models considered, which suggests its potential to act as a powerful complementary probe to BAO and galaxy clusters. (2) Turning to $Om(z)$ diagnostics, CPL, GCG and EDE models can not be distinguished from each other at the present epoch. (3) In the framework of statefinder diagnostics, MPC and EDE will deviate from $\rm{\Lambda}$CDM model in the near future, while GCG model cannot be distinguished from $\rm{\Lambda}$CDM model unless much higher precision observations are available.

In this paper, we present a new compiled milliarcsecond compact radio data set of 120 intermediate-luminosity quasars in the redshift range $0.46< z <2.76$. These quasars show negligible dependence on redshifts and intrinsic luminosity, and thus represents, in the standard model of cosmology, a fixed comoving-length of standard ruler. We implement a new cosmology-independent technique to calibrate the linear size of of this standard ruler as $l_m= 11.03\pm0.25$ pc, which is the typical radius at which AGN jets become opaque at the observed frequency $\nu\sim 2$ GHz. In the framework of flat $\Lambda$CDM model, we find a high value of the matter density parameter, $\Omega_m=0.322^{+0.244}_{-0.141}$, and a low value of the Hubble constant, $H_0=67.6^{+7.8}_{-7.4}\; \rm{kms}^{-1}\rm{Mpc}^{-1}$, which is in excellent agreement with the CMB anisotropy measurements by \textit{Planck}. We obtain ${\Omega_m}=0.309^{+0.215}_{-0.151}$, $w=-0.970^{+0.500}_{-1.730}$ at 68.3% CL for the constant $w$ of a dynamical dark-energy model, which demonstrates no significant deviation from the concordance $\Lambda$CDM model. Consistent fitting results are also obtained for other cosmological models explaining the cosmic acceleration, like Ricci dark energy (RDE) or Dvali-Gabadadze-Porrati (DGP) brane-world scenario. While no significant change in $w$ with redshift is detected, there is still considerable room for evolution in $w$ and the transition redshift at which $w$ departing from -1 is located at $z\sim 2.0$. Our results demonstrate that the method extensively investigated in our work on observational radio quasar data can be used to effectively derive cosmological information. Finally, we find the combination of high-redshift quasars and low-redshift clusters may provide an important source of angular diameter distances, considering the redshift coverage of these two astrophysical probes.

Using a new recently compiled milliarcsecond compact radio data set of 120 intermediate-luminosity quasars in the redshift range $0.46< z <2.76$, whose statistical linear sizes show negligible dependence on redshifts and intrinsic luminosity and thus represent standard rulers in cosmology, we constrain three viable and most popular $f(T)$ gravity models, where $T$ is the torsion scalar in teleparallel gravity. Our analysis reveals that constraining power of the quasars data (N=120) is comparable to the Union2.1 SN Ia data (N=580) for all three $f(T)$ models. Together with other standard ruler probes such as Cosmic Microwave Background and Baryon Acoustic Oscillation distance measurements, the present value of the matter density parameter $\Omega_m$ obtained by quasars is much lager than that derived from other observations. For two of the models considered ($f_1$CDM and $f_2$CDM) a small but noticeable deviation from $\Lambda$CDM cosmology is present, while in the framework of $f_3$CDM the effective equation of state may cross the phantom divide line at lower redshifts. These results indicate that intermediate-luminosity quasars could provide an effective observational probe comparable to SN Ia at much higher redsifts, and $f(T)$ gravity is a reasonable candidate for the modified gravity theory.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

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.

In this work, we use a test based on the differential ages of galaxies for distinguishing the dark energy models. As proposed by Jimenez and Loeb, relative ages of galaxies can be used to put constraints on various cosmological parameters. In the same vein, we reconstruct $H_0dt/dz$ and its derivative ($H_0d^2t/dz^2$) using a model independent technique called non-parametric smoothing. Basically, $dt/dz$ is the change in the age of the object as a function of redshift which is directly link with the Hubble parameter. Hence for reconstruction of this quantity, we use the most recent $H(z)$ data. Further, we calculate $H_0dt/dz$ and its derivative for several models like Phantom, Einstein de Sitter (EdS), $\Lambda$CDM, Chevallier-Polarski-Linder (CPL) parametrization, Jassal-Bagla-Padmanabhan (JBP) parametrization and Feng-Shen-Li-Li (FSLL) parametrization. We check the consistency of these models with the results of reconstruction obtained in model independent way from the data. It is observed that $H_0dt/dz$ as a tool is not able to distinguish between the $\Lambda$CDM, CPL, JBP and FSLL parametrizations but as expected EdS and Phantom models show noticeable deviation from the reconstructed results. Further, the derivative of $H_0dt/dz$ for various dark energy models is more sensitive at low redshift. It is observed that the FSLL model is not consistent with the reconstructed results at redshifts less than $0.5$, however, the $\Lambda$CDM model is in concordance with the 3$\sigma$ region of the reconstruction.

In this paper, based on a 2.29 GHz VLBI all-sky survey of 613 milliarcsecond ultra-compact radio sources with $0.0035<z<3.787$, we describe a method of identifying the sub-sample which can serve as individual standard rulers in cosmology. If the linear size of the compact structure is assumed to depend on source luminosity and redshift as $l_m=l L^\beta (1+z)^n$, only intermediate-luminosity quasars ($10^{27}$ W/Hz$<L<$ $10^{28}$ W/Hz) show negligible dependence ($|n|\simeq 10^{-3}$, $|\beta|\simeq 10^{-4}$), and thus represent a population of such rulers with fixed characteristic length $l=11.42$ pc. With a sample of 120 such sources covering the redshift range $0.46<z<2.80$, we confirm the existence of dark energy in the Universe with high significance under the assumption of a flat universe, and obtain stringent constraints on both the matter density $\Omega_m=0.323^{+0.245}_{-0.145}$ and the Hubble constant $H_0=66.30^{+7.00}_{-8.50}$ km sec$^{-1}$ Mpc$^{-1}$. Finally, with the angular diameter distances $D_A$ measured for quasars extending to high redshifts ($z\sim 3.0$), we reconstruct the $D_A(z)$ function using the technique of Gaussian processes. This allows us to identify the redshift corresponding to the maximum of the $D_A(z)$ function: $z_m=1.70$ and the corresponding angular diameter distance $D_A(z_m)=1719.01\pm43.46$ Mpc. Similar reconstruction of the expansion rate function $H(z)$ based on the data from cosmic chronometers and BAO gives us $H(z_m)=176.77\pm6.11$ km sec$^{-1}$ Mpc$^{-1}$. These measurements are used to estimate the speed of light: $c=3.039(\pm0.180)\times 10^5$ km/s. This is the first measurement of the speed of light in a cosmological setting referring to the distant past.

Based on a mass-selected sample of galaxy-scale strong gravitational lenses from the SLACS, BELLS, LSD and SL2S surveys and using a well-motivated fiducial set of lens-galaxy parameters we tested the weak-field metric on kiloparsec scales and found a constraint on the post-Newtonian parameter $\gamma = 0.995^{+0.037}_{-0.047}$ under the assumption of a flat $\Lambda$CDM universe with parameters taken from \textit{Planck} observations. General relativity (GR) predicts exactly $\gamma=1$. Uncertainties concerning the total mass density profile, anisotropy of the velocity dispersion and the shape of the light-profile combine to systematic uncertainties of $\sim 25\%$. By applying a cosmological model independent method to the simulated future LSST data, we found a significant degeneracy between the PPN $\gamma$ parameter and spatial curvature of the Universe. Setting a prior on the cosmic curvature parameter $-0.007< \Omega_k <0.006$, we obtained the following constraint on the PPN parameter: $\gamma=1.000^{+0.0023}_{-0.0025}$. We conclude that strong-lensing systems with measured stellar velocity dispersions may serve as another important probe to investigate validity of the GR, if the mass-dynamical structure of the lensing galaxies is accurately constrained in the future lens surveys.

We investigate two dark energy cosmological models (i.e., the $\Lambda$CDM and $\phi$CDM models) with massive neutrinos assuming two different neutrino mass hierarchies in both the spatially flat and non-flat scenarios, where in the $\phi$CDM model the scalar field possesses an inverse power-law potential, $V(\phi)\propto {\phi}^{-\alpha}$ ($\alpha>0$). Cosmic microwave background data from Planck 2015, baryon acoustic oscillations data from 6dFGS, SDSS-MGS, BOSS-LOWZ and BOSS CMASS-DR11, the JLA compilation of Type Ia supernova apparent magnitude observations, and the Hubble Space Telescope $H_0$ prior, are jointly employed to constrain the model parameters. We first determine constraints assuming three species of degenerate massive neutrinos. In the spatially flat (non-flat) $\Lambda$CDM model, the sum of neutrino masses is bounded as $\Sigma m_{\nu} < 0.165 (0.299)$ eV at 95% confidence level (CL). Correspondingly, in the flat (non-flat) $\phi$CDM model, we find $\Sigma m_{\nu} < 0.164 (0.301)$ eV at 95% CL. The inclusion of spatial curvature as a free parameter results in a significant broadening of confidence regions for $\Sigma m_{\nu}$ and other parameters. In the scenario where the total neutrino mass is dominated by the heaviest neutrino mass eigenstate, we can obtain the similar conclusions as those obtained in the degenerate neutrino mass scenario. In addition, the results show that the bounds on $\Sigma m_{\nu}$ based on two different neutrino mass hierarchies have insignificant differences in the spatially flat case for both the $\Lambda$CDM and $\phi$CDM models, however, the corresponding differences are larger in the non-flat case.

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.

Two-point diagnostics $Om(z_i,z_j)$ and $Omh^2(z_i,z_j)$ have been introduced as an interesting tool for testing the validity of the $\Lambda$CDM model. Quite recently, Sahni, Shafieloo $\&$ Starobinsky (2014) combined two independent measurements of $H(z)$ from BAO data with the value of the Hubble constant $H_0$, and used the second of these diagnostics to test the $\Lambda$CDM model. Their result indicated a considerable tension between observations and predictions of the $\Lambda$CDM model. Since reliable data concerning expansion rates of the Universe at different redshifts $H(z)$ are crucial for the successful application of this method, we investigate both two-point diagnostics on the most comprehensive set of $N=36$ measurements of $H(z)$ coming from the BAO and differential ages (DA) of passively evolving galaxies. We discuss the uncertainties of two-point diagnostics and find that they are strongly non-Gaussian and follow the patterns deeply rooted in their very construction. Therefore we propose that non-parametric median statistics is the most appropriate way of treating this problem. Our results support the claims that $\Lambda$CDM is in tension with $H(z)$ data according to the two-point diagnostics developed by Shafieloo, Sahni and Starobinsky. However, other alternatives to the $\Lambda$CDM, such as wCDM or CPL models perform even worse. We also notice that there are serious systematic differences between BAO and DA methods which ought to be better understood before $H(z)$ measurements can become competitive to the other probes.

We use 118 strong gravitational lenses observed by the SLACS, BELLS, LSD and SL2S surveys to constrain the total mass profile and the profile of luminosity density of stars (light-tracers) in elliptical galaxies up to redshift $z \sim 1$. Assuming power-law density profiles for the total mass density, $\rho=\rho_0(r/r_0)^{-\alpha}$, and luminosity density, $\nu=\nu_0(r/r_0)^{-\delta}$, we investigate the power law index and its first derivative with respect to the redshift. Using Monte Carlo simulations of the posterior likelihood taking the Planck's best-fitted cosmology as a prior, we find $\gamma= 2.132\pm0.055$ with a mild trend $\partial \gamma/\partial z_l= -0.067\pm0.119$ when $\alpha=\delta=\gamma$, suggesting that the total density profile of massive galaxies could have become slightly steeper over cosmic time. Furthermore, similar analyses performed on sub-samples defined by different lens redshifts and velocity dispersions, indicate the need of treating low, intermediate and high-mass galaxies separately. Allowing $\delta$ to be a free parameter, we obtain $\alpha=2.070\pm0.031$, $\partial \alpha/\partial z_l= -0.121\pm0.078$, and $\delta= 2.710\pm0.143$. The model in which mass traces light is rejected at $>95\%$ confidence and our analysis robustly indicates the presence of dark matter in the form of a mass component that is differently spatially extended than the light. In this case, intermediate-mass elliptical galaxies ($200$ km/s $ < \sigma_{ap} \leq 300$ km/s) show the best consistency with the singular isothermal sphere as an effective model of galactic lenses.

In this paper, assuming the validity of distance duality relation, $\eta=D_L(z)(1+z)^{-2}/D_A(z)=1$, where $D_A(z)$ and $D_L(z)$ are the angular and the luminosity distance respectively, we explore two kinds of gas mass density profiles of clusters: the isothermal $\beta$ model and the non-isothermal double-$\beta$ model. In our analysis, performed on 38 massive galaxy clusters observed by \textit{Chandra} (within the redshift range of $0.14<z<0.89$), we use two types of cluster gas mass fraction data corresponding to different mass density profiles fitted to the X-ray data. Using two general parameterizations of $\eta(z)$ (phenomenologically allowing for distance duality violation), we find that the non-isothermal double-$\beta$ model agrees better with the distance duality relation, while the isothermal $\beta$ model tends to be marginally incompatible with the Etherington theorem at 68.3% CL. However, current accuracy of the data does not allow to distinguish between the two models for the gas-density distribution at a significant level.

In this paper, we used standard rulers and standard candles (separately and jointly) to explore five popular dark energy models under assumption of spatial flatness of the Universe. As standard rulers, we used a data set comprising 118 galactic-scale strong lensing systems (individual standard rulers if properly calibrated for the mass density profile) combined with BAO diagnostics (statistical standard ruler). Supernovae Ia served asstandard candles. Unlike in the most of previous statistical studies involving strong lensing systems, we relaxed the assumption of singular isothermal sphere (SIS) in favor of its generalization: the power-law mass density profile. Therefore, along with cosmological model parameters we fitted the power law index and its first derivative with respect to the redshift (thus allowing for mass density profile evolution). It turned out that the best fitted $\gamma$ parameters are in agreement with each other irrespective of the cosmological model considered. This demonstrates that galactic strong lensing systems may provide a complementary probe to test the properties of dark energy. Fits for cosmological model parameters which we obtained are in agreement with alternative studies performed by the others. Because standard rulers and standard candles have different parameter degeneracies, combination of standard rulers and standard candles gives much more restrictive results for cosmological parameters. At last, we attempted at model selection based on information theoretic criteria (AIC and BIC). Our results support the claim, that cosmological constant model is still the best one and there is no (at least statistical) reason to prefer any other more complex model.

With a fantastic sensitivity improving significantly over the advanced GW detectors, Einstein Telescope (ET) will be able to observe hundreds of thousand inspiralling double compact objects per year. By virtue of gravitational lensing effect, intrinsically unobservable faint sources can be observed by ET due to the magnification by intervening galaxies. We explore the possibility of observing such faint sources amplified by strong gravitational lensing. Following our previous work, we use the merger rates of DCO (NS-NS,BH-NS,BH-BH systems) as calculated by Dominik et al.(2013). It turns out that tens to hundreds of such (lensed) extra events will be registered by ET. This will strongly broaden the ET's distance reach for signals from such coalescences to the redshift range z=2 - 8. However, with respect to the full inspiral event catalog this magnification bias is at the level of 0.001 and should not affect much cosmological inferences.

In this paper, we assemble a catalog of 118 strong gravitational lensing systems from SLACS, BELLS, LSD and SL2S surveys and use them to constrain the cosmic equation of state. In particular we consider two cases of dark energy phenomenology: $XCDM$ model where dark energy is modeled by a fluid with constant $w$ equation of state parameter and in Chevalier - Polarski - Linder (CPL) parametrization where $w$ is allowed to evolve with redshift: $w(z) = w_0 + w_1 \frac{z}{1+z}$. We assume spherically symmetric mass distribution in lensing galaxies, but relax the rigid assumption of SIS model in favor to more general power-law index $\gamma$, also allowing it to evolve with redshifts $\gamma(z)$. Our results for the $XCDM$ cosmology show the agreement with values (concerning both $w$ and $\gamma$ parameters) obtained by other authors. We go further and constrain the CPL parameters jointly with $\gamma(z)$. The resulting confidence regions for the parameters are much better than those obtained with a similar method in the past. They are also showing a trend of being complementary to the supernova Ia data. Our analysis demonstrates that strong gravitational lensing systems can be used to probe cosmological parameters like the cosmic equation of state for dark energy. Moreover, they have a potential to judge whether the cosmic equation of state evolved with time or not.

Recently, Sahni, Shafielo o & Starobinsky (2014) combined two independent measurements of $H(z)$ from BAO data with the value of the Hubble constant $H_0 = H(z=0)$, in order to test the cosmological constant hypothesis by means of an improved version of the $Om$ diagnostic. Their result indicated a considerable tension between observations and predictions of the $\Lambda$CDM model. However, such strong conclusion was based only on three measurements of $H(z)$. This motivated us to repeat similar work on a larger sample. By using a comprehensive data set of 29 $H(z)$, we find that discrepancy indeed exists. Even though the value of $\Omega_{m,0} h^2$ inferred from $Omh^2$ diagnostic depends on the way one chooses to make a summary statistics (weighted mean or the median), the persisting discrepancy supports the claims of Sahni, Shafielo o & Starobinsky (2014) that $\Lambda$CDM model may not be the best description of our Universe.

Cosmological applications of the "redshift - angular size" test require knowledge of the linear size of the "standard rod" used. In this paper, we study the properties of a large sample of 140 milliarcsecond compact radio sources with flux densities measured at 6 cm and 20 cm, compiled by Gurvits et al.(1999). Using the best-fitted cosmological parameters given by Planck/WMAP9 observations, we investigate the characteristic length $l_m$ as well as its dependence on the source luminosity $L$ and redshift $l_m=l L^\beta (1+z)^n$. For the full sample, measurements of the angular size $\theta$ provide a tight constraint on the linear size parameters. We find that cosmological evolution of the linear size is small ($|n|\simeq 10^{-2}$) and consistent with previous analysis. However, a substantial evolution of linear sizes with luminosity is still required ($\beta\simeq 0.17$). Furthermore, similar analysis done on sub-samples defined by different source optical counterparts and different redshift ranges, seems to support the scheme of treating radio galaxies and quasars with distinct strategies. Finally, a cosmological-model-independent method is discussed to probe the properties of angular size of milliarcsecond radio quasars. Using the corrected redshift - angular size relation for quasar sample, we obtained a value of the matter density parameter, $\Omega_m=0.292^{+0.065}_{-0.090}$, in the spatially flat $\Lambda$CDM cosmology.

Gravitational wave (GW) experiments are entering their advanced stage which should soon open a new observational window on the Universe. Looking into this future, the Einstein Telescope (ET) was designed to have a fantastic sensitivity improving significantly over the advanced GW detectors. One of the most important astrophysical GW sources supposed to be detected by the ET in large numbers are double compact objects (DCO) and some of such events should be gravitationally lensed by intervening galaxies. We explore the prospects of observing gravitationally lensed inspiral DCO events in the ET. This analysis is a significant extension of our previous paper. We are using the intrinsic merger rates of the whole class of DCO (NS-NS,BH-NS,BH-BH)located at different redshifts as calculated by Dominik et al. by using {\tt StarTrack} population synthesis evolutionary code. We discuss in details predictions from each evolutionary scenario. Our general conclusion is that ET would register about $50-100$ strongly lensed inspiral events per year. Only the scenario in which nascent BHs receive strong kick gives the predictions of a few events per year. Such lensed events would be dominated by the BH-BH merging binary systems. Our results suggest that during a few years of successful operation ET will provide a considerable catalog of strongly lensed events.

Gravitational wave experiments have entered a new stage which gets us closer to the opening a new observational window on the Universe. In particular, the Einstein Telescope (ET) is designed to have a fantastic sensitivity that will provide with tens or hundreds of thousand NS-NS inspiral events per year up to the redshift z = 2. Some of such events should be gravitationally lensed by intervening galaxies. We explore the prospects of observing gravitationally lensed inspiral NS-NS events in the Einstein telescope. Being conservative we consider the lens population of elliptical galaxies. It turns out that depending on the local insipral rate ET should detect from one per decade detection in the pessimistic case to a tens of detections per year for the most optimistic case. The detection of gravitationally lensed source in gravitational wave detectors would be an invaluable source of information con- cerning cosmography, complementary to standard ones (like supernovae or BAO) independent of the local cosmic distance ladder calibrations.

Classical Cepheids form one of the foundations of modern cosmology and the extragalactic distance scale, however, cosmic microwave background observations measure cosmological parameters and indirectly the Hubble Constant, H0, to unparalleled precision. The coming decade will provide opportunities to measure H0 to 2% uncertainty thanks to the GAIA satellite, JWST, ELTs and other telescopes using Cepheids and other standard candles. In this work, we discuss the upcoming role for variable stars and asteroseismology in calibrating the distance scale and measuring H0 and what problems exist in understanding these stars that will feedback on these measurements.

Strong lensing has developed into an important astrophysical tool for probing both cosmology and galaxies (their structure, formation, and evolution). Using the gravitational lensing theory and cluster mass distribution model, we try to collect a relatively complete observational data concerning the Hubble constant independent ratio between two angular diameter distances $D_{ds}/D_s$ from various large systematic gravitational lens surveys and lensing by galaxy clusters combined with X-ray observations, and check the possibility to use it in the future as complementary to other cosmological probes. On one hand, strongly gravitationally lensed quasar-galaxy systems create such a new opportunity by combining stellar kinematics (central velocity dispersion measurements) with lensing geometry (Einstein radius determination from position of images). We apply such a method to a combined gravitational lens data set including 70 data points from Sloan Lens ACS (SLACS) and Lens Structure and Dynamics survey (LSD). On the other hand, a new sample of 10 lensing galaxy clusters with redshifts ranging from 0.1 to 0.6 carefully selected from strong gravitational lensing systems with both X-ray satellite observations and optical giant luminous arcs, is also used to constrain three dark energy models ($\Lambda$CDM, constant $w$ and CPL) under a flat universe assumption. For the full sample ($n=80$) and the restricted sample ($n=46$) including 36 two-image lenses and 10 strong lensing arcs, we obtain relatively good fitting values of basic cosmological parameters, which generally agree with the results already known in the literature. This results encourages further development of this method and its use on larger samples obtained in the future.

Accelerating expansion of the Universe is a great challenge for both physics and cosmology. In light of lacking the convincing theoretical explanation, an effective description of this phenomenon in terms of cosmic equation of state turns out useful. The strength of modern cosmology lies in consistency across independent, often unrelated pieces of evidence. Therefore, every alternative method of restricting cosmic equation of state is important. Strongly gravitationally lensed quasar-galaxy systems create such new opportunity by combining stellar kinematics (central velocity dispersion measurements) with lensing geometry (Einstein radius determination form position of images). In this paper we apply such method to a combined data sets from SLACS and LSD surveys of gravitational lenses. In result we obtain the cosmic equation of state parameters, which generally agree with results already known in the literature. This demonstrates that the method can be further used on larger samples obtained in the future. Independently noticed systematic deviation between fits done on standard candles and standard rulers is revealed in our findings. We also identify an important selection effect crucial to our method associated with geometric configuration of the lensing system along line of sight, which may have consequences for sample construction from the future lensing surveys.

Lorentz Invariance Violation (LIV) manifesting itself by energy dependent modification of standard relativistic dispersion relation has recently attracted a considerable attention. Ellis et al. previously investigated the energy dependent time offsets in different energy bands on a sample of gamma ray bursts and, assuming standard cosmological model, they found a weak indication for redshift dependence of time delays suggestive of LIV. Going beyond the $\Lambda$CDM cosmology we extend this analysis considering also four alternative models of dark energy (quintessence with constant and variable equation of state, Chaplygin gas and brane-world cosmology). It turns out that the effect noticed by Ellis et al. is also present in those models and is the strongest for quintessence with variable equation of state.

Despite the fact that quantum gravity theory still remains elusive, it is generally expected that it will bring the picture of a space-time foam at short distances leading to Lorenz Invariance Violation (LIV) manifested e.g. by energy dependent modification of standard relativistic dispersion relation. One direction of research, pursued intensively is to measure energy dependent time of arrival delays in photons emitted by astrophysical sources located at cosmological distances. This is tempered however by our ignorance of either intrinsic delay (at source frame) in different energy channels or as recently shown, the correct model of the background cosmology. In this paper we propose a new test based on gravitational lensing. Monitoring time delays between images performed in different energy channels (e.g. optical - low energy and TeV photons) may reveal extra delays due to distorted dispersion relation typical in LIV theories - a test which is free from the systematics inherent in other settings.

Modern ideas in quantum gravity predict the possibility of Lorenz Invariance Violation (LIV) manifested e.g. by energy dependent modification of standard relativistic dispersion relation. In a recent paper Jacob and Piran proposed that time of flight delays in high energy neutrinos emitted by gamma ray bursts (GRBs) located at cosmological distances can become a valuable tool for setting limits on LIV theories. However, current advances in observational cosmology suggest that our Universe is dominated by dark energy with relatively little guidance on its nature thus leading to several cosmological scenarios compatible with observations. In this paper we raise the issue of how important, in the context of testing LIV theories, is our knowledge of background cosmological model. Specifically we calculate expected time lags for high-energy (100 TeV) neutrinos in different cosmological models. Out of many particular models of dark energy we focus on five: $\Lambda$CDM, quintessence, quintessence with time varying equation of state, brane-world and generalized Chaplygin gas model as representative for various competing approaches. The result is that uncertainty introduced by our ignorance concerning the right phenomenological model describing dark energy dominated universe is considerable and may obscure bounds derived from studying time delays from cosmological sources.

There are several different theoretical ideas invoked to explain the dark energy with relatively little guidance of which one of them might be right. Therefore the emphasis of ongoing and forthcoming research in this field shifts from estimating specific parameters of cosmological model to the model selection. In this paper we apply information-theoretic model selection approach based on Akaike criterion as an estimator of Kullback-Leibler entropy. In particular, we present the proper way of ranking the competing models based on Akaike weights (in Bayesian language - posterior probabilities of the models). Out of many particular models of dark energy we focus on four: quintessence, quintessence with time varying equation of state, brane-world and generalized Chaplygin gas model and test them on Riess' Gold sample. As a result we obtain that the best model - in terms of Akaike Criterion - is the quintessence model. The odds suggest that although there exist differences in the support given to specific scenarios by supernova data most of the models considered receive similar support. One can also notice that models similar in structure i.e. $\Lambda$CDM, quintessence and quintessence with variable equation of state are closer to each other in terms of Kullback-Leibler entropy. Models having different structure i.e. Chaplygin gas or brane-world scenario are more distant (in Kullback-Leibler sense) from the best one.

Could cosmic topology imply dark energy? We use a weak field (Newtonian) approximation of gravity and consider the gravitational effect from distant, multiple copies of a large, collapsed (virialised) object today (i.e. a massive galaxy cluster), taking into account the finite propagation speed of gravity, in a flat, multiply connected universe, and assume that due to a prior epoch of fast expansion (e.g. inflation), the gravitational effect of the distant copies is felt locally, from beyond the naively calculated horizon. We find that for a universe with a $T^1xR^2$ spatial section, the residual Newtonian gravitational force (to first order) provides an anisotropic effect that repels test particles from the cluster in the compact direction, in a way algebraically similar to that of dark energy. For a typical test object at comoving distance $\chi$ from the nearest dense nodes of the cosmic web of density perturbations, the pressure-to-density ratio $w$ of the equation of state in an FLRW universe, is w \sim - (\chi/L)^3, where $L$ is the size of the fundamental domain, i.e. of the universe. Clearly, |w|<<1. For a T^3 spatial section of exactly equal fundamental lengths, the effect cancels to zero. For a T^3 spatial section of unequal fundamental lengths, the acceleration effect is anisotropic in the sense that it will *tend to equalise the three fundamental lengths*. Provided that at least a modest amount of inflation occurred in the early Universe, and given some other conditions, multiple connectedness does generate an effect similar to that of dark energy, but the amplitude of the effect at the present epoch is too small to explain the observed dark energy density and its anisotropy makes it an unrealistic candidate for the observed dark energy.

The so called Generalized Chaplygin Gas (GCG) with the equation of state $p = - \frac{A}{{\rho}^{\alpha}}$ was recently proposed as a candidate for dark energy in the Universe. In this paper we confront the GCG with SNIa data. Specifically we have tested the GCG cosmology in three different classes of models with (1) $\Omega_m= 0.3$, $\Omega_{Ch}= 0.7$; (2) $\Omega_m= 0.05$, $\Omega_{Ch}= 0.95$ and (3) $\Omega_m = 0$, $\Omega_{Ch} = 1$, as well as the model withouth any assumption on $\Omega_m$. The best fitted models are obtained by minimalizing the $\chi^2$ function and $\chi^2$ levels in the $(A_0, \alpha)$ plane. We supplemented our analysis with confidence intervals in the $(A_0, \alpha)$ plane, as well as one-dimensional probability distribution functions for models parameter. The general conclusion is that SNIa data strongly support the Chaplygin gas (with $\alpha = 1$). Extending our analysisby relaxing the flat prior lead to the result that even though the best fitted values of $\Omega_k$ are formally non-zero, still they are close to flat case. It should be viewed as an advantage of the GCG model since in similar analysisof $\Lambda$CDM model high negative value of $\Omega_{k}$ were found to be bestfitted to the data and independent inspiration from CMBR and extragalactic astronomy has been invoked to fix the curvature problem. Our results show clearly that in Generalized Chaplygin Gas cosmology distant $z >1$ supernovae should be brighter than in $\Lambda$CDM model.This prediction seems to be confirmed with new Riess high redshift SNIa sample. Moreover, we argue that with the future SNAP data it would be possible to differentiate between models with various value of $\alpha$ parameter and/or discriminated between GCG, Cardassian and $\Lambda$CDM models

In this paper we derive a bound on the rate of change of the gravitational constant G coming from the pulsating white dwarf G117-B15A. This star is a ZZ Ceti pulsator extensively studied with astroseismological techniques for last three decades. The most recent determination of {\dot P} = (2.3 \pm 1.4) * 10^{-15} s/s^{-1} for the 215.2s fundamental mode agrees very well with predictions of the best fit theoretical model. The rate of change of the oscillation period can be explained by two effects: the cooling (dominant factor) and change of gravitational binding energy (residual gravitational contraction). Since the white dwarfs are pulsating in g-modes whose frequencies are related to the Brunt-Vaisala frequency (explicitly dependent on G) observational determination of the change of the period (more precisely the difference between observed and calculated \dot P) can be used to set the upper bound on the rate of change of G. In the light of the current data concerning G117-B15A we derive the following bound: |{\frac {\dot G}{G}}| \leq 4.10 \times 10^{-10} yr^{-1}. Our result is model independent in the sense that it does not need to invoke a concrete physical theory (such like Brans-Dicke theory)underlying the temporal variability of G. We also demonstrate that varying gravitational constant G does not modify cooling of white dwarfs in a significant way. This result implies that some earlier claims present in the literature that varying G can be reflected in the WD luminosity function are not correct.

Theories of fundamental interactions with large extra dimensions have recently become very popular. Astrophysical bounds from the Sun, red-giants and SN1987a have already been derived by other authors for the theory proposed by Arkani-Hamed, Dimopoulos and Dvali. In this paper we consider G117-B15A pulsating white dwarf (ZZ Ceti star) for which the secular rate at which the period of its fundamental mode increases has been accurately measured and claimed that this mode of G117-B15A is perhaps the most stable oscillation ever recorded in the optical band. Because an additional channel of energy loss (Kaluza-Klein gravitons) would speed up the cooling rate, one is able to use the aforementioned stability to derive a bound on theories with large extra dimensions. Within the framework of the theory with large extra dimensions proposed by Arkani-Hamed, Dimopoulos and Dvali we find the lower bound on string comapctification scale which is more stringent than solar or red-giant bounds.

Future generation of interferometric gravitational wave detectors is hoped to provide accurate measurements of the final stages of binary inspirals. The sources probed by such experiments are of extragalactic origin and the observed chirp mass is the intrinsic chirp mass multiplied by $(1+z)$ where $z$ is the redshift of the source. Moreover the luminosity distance is a direct observable is such experiments. This creates the possibility to establish a new kind of cosmological tests, supplementary to more standard ones. Recent observations of distant type Ia supernovae light-curves suggest that the expansion of the universe has recently begun to accelerate. A popular explanation of present accelerating expansion of the universe is to assume that some part $\Omega_Q$ of the matter-energy density is in the form of dark component called ``the quintessence'' with the equation of state $p_Q = w \rho_Q$ with $w \geq -1.$ In this paper we consider the predictions concerning observations of binary inspirals in future LIGO type interferometric experiments assuming a ``quintessence cosmology''. In particular we compute the expected redshift distributions of observed events in the a priori admissible range of parameters describing the equation of state for the quintessence. We find that this distribution has a robust dependence on the cosmic equation of state.

In this paper we discuss the QCD phase-transitions in the nontopological soliton model of quark confinement and explore possible astrophysical consequences. Our key idea is to look at quark stars (which are believed to exist since the quark matter is energetically preferred over the ordinary matter) from the point of view of soliton model. We propose that the phase transition taking place during the core collapse of massive evolved star may provide a new physical effect not taken into account in modeling the supernova explosions. We also point out the possibility that merging quark stars may produce gamma-ray bursts energetic enough to be at cosmological distances. Our idea based on the finite-temperature nontopologiocal soliton model overcomes major difficulties present in neutron star merger scenario --- the baryon loading problem and nonthermal spectra of the bursts.

In this letter we propose a physical explanation for recently reported correlations between pairs of close and antipodal gamma-ray bursts from publicly available BATSE catalogue. Our model is based on the cosmological scenario in which bursters are located at cosmological distances of order of 0.5--2~Gpc. Observed distribution of gamma-ray bursts strongly suports this assumption. If so gamma-ray bursts may provide a very good probe for investigating the topological structure of the Universe. We notice that correlation between antipodal events may in fact indicate that we live in the so called Ellis' small universe which has Friedman-Roberston-Walker metric structure and nontrivial topology.