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).

We use HII starburst galaxy apparent magnitude measurements to constrain cosmological parameters in six cosmological models. A joint analysis of HII galaxy, quasar angular size, baryon acoustic oscillations peak length scale, and Hubble parameter measurements result in relatively model-independent and restrictive estimates of the current values of the non-relativistic matter density parameter $\Omega_{\rm m_0}$ and the Hubble constant $H_0$. These estimates favor a 2.0$\sigma$ to 3.4$\sigma$ (depending on cosmological model) lower $H_0$ than what is measured from the local expansion rate. The combined data are consistent with dark energy being a cosmological constant and with flat spatial hypersurfaces, but do not strongly rule out mild dark energy dynamics or slightly non-flat spatial geometries.

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.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

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.

Astrophysical observations provide a unique opportunity to test possible signatures of Lorentz Invariance Violation (LIV), due to the high energies and long distances involved. In quantum theory of gravity, one may expect the modification of the dispersion relation between energy and momentum for photons, which can be probed with the time-lag (the arrival time delay between light curves in different energy bands) of Gamma-ray bursts (GRBs). In this paper, by using the detailed time-delay measurements of GRB 160625B at different energy bands, as well as 23 time-delay GRBs covering the redshifts range of $z=0.168-2.5$ (which were measured at different energy channels from the light curves), we propose an improved model-independent method (based on the newly-compiled sample of $H(z)$ measurements) to probe the energy-dependent velocity due to the modified dispersion relation for photons. In the framework of a more complex and reasonable theoretical expression to describe the time delays, our results imply that the intrinsic time lags can be better described with more GRBs time delay data. More importantly, through direct fitting of the time-delay measurements of a sample of GRBs, our limit on the LIV energy scale is comparable to that with unknown constant for the intrinsic time lag, much lower than the Planck energy scale in both linear LIV and quadratic LIV cases.

Cosmological applications of HII galaxies (HIIGx) and giant extragalactic HII regions (GEHR) to construct the Hubble diagram at higher redshifts require knowledge of the "$L$--$\sigma$" relation of the standard candles used. In this paper, we study the properties of a large sample of 156 sources (25 high-$z$ HII galaxies, 107 local HII galaxies, and 24 giant extragalactic HII regions) compiled by Terlevich et al.(2015). Using the the cosmological distances reconstructed through two new cosmology-independent methods, we investigate the correlation between the H$\beta$ emission-line luminosity $L$ and ionized-gas velocity dispersion $\sigma$. The method is based on non-parametric reconstruction using the measurements of Hubble parameters from cosmic clocks, as well as the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard sirens. Assuming the emission-line luminosity versus ionized gas velocity dispersion relation, $\log L ($H$\beta) = \alpha \log \sigma($H$\beta)+\kappa$, we find the full sample provides a tight constraint on the correlation parameters. However, similar analysis done on three different sub-samples seems to support the scheme of treating HII galaxies and giant extragalactic HII regions with distinct strategies. Using the corrected "$L$--$\sigma$" relation for the HII observational sample beyond the current reach of Type Ia supernovae, we obtain a value of the matter density parameter, $\Omega_{m}=0.314\pm0.054$ (calibrated with standard clocks) and $\Omega_{m}=0.311\pm0.049$ (calibrated with standard sirens), in the spatially flat $\Lambda$CDM cosmology.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We use three different data sets, specifically $H(z)$ measurements from cosmic chronometers, the HII-galaxy Hubble diagram, and reconstructed quasar-core angular-size measurements, to perform a joint analysis of three flat cosmological models: the $R_{\rm h}=ct$ Universe, $\Lambda$CDM, and $w$CDM. For $R_{\rm h}=ct$, the 1$\sigma$ best-fit value of the Hubble constant $H_0$ is $62.336\pm1.464$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, which matches previous measurements ($\sim 63$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$) based on best fits to individual data sets. For $\Lambda$CDM, our inferred value of the Hubble constant, $H_0=67.013\pm2.578$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, is more consistent with the ${\it Planck}$ optimization than the locally measured value using $\mbox{Cepheid}$ variables, and the matter density $\Omega_{\rm m}=0.347\pm0.049$ similarly coincides with its ${\it Planck}$ value to within 1$\sigma$. For $w$CDM, the optimized parameters are $H_0=64.718\pm3.088$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, $\Omega_{\rm m}=0.247\pm0.108$ and $w=-0.693\pm0.276$, also consistent with ${\it Planck}$. A direct comparison of these three models using the Bayesian Information Criterion shows that the $R_{\rm h}=ct$ universe is favored by the joint analysis with a likelihood of $\sim 97\%$ versus $\lesssim 3\%$ for the other two cosmologies.

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.

In this paper, we present a scheme to investigate the opacity of the Universe in a cosmological-model-independent way, with the combination of current and future measurements of type Ia supernova sample and galactic-scale strong gravitational lensing systems with SNe Ia acting as background sources. The observational data include the current newly-compiled SNe Ia data (Pantheon sample) and simulated sample of SNe Ia observed by the forthcoming Large Synoptic Survey Telescope (LSST) survey, which are taken for luminosity distances ($D_L$) possibly affected by the cosmic opacity, as well as strongly lensed SNe Ia observed by the LSST, which are responsible for providing the observed time-delay distance ($D_{\Delta t}$) unaffected by the cosmic opacity. Two parameterizations, $\tau(z)=2\beta z$ and $\tau(z)=(1+z)^{2\beta}-1$ are adopted for the optical depth associated to the cosmic absorption. Focusing on only one specific type of standard cosmological probe, this provides an original method to measure cosmic opacity at high precision. Working on the simulated sample of strongly lensed SNe Ia observed by the LSST in 10 year $z$-band search, our results show that, with the combination of the current newly-compiled SNe Ia data (Pantheon sample), there is no significant deviation from the transparency of the Universe at the current observational data level. Moreover, strongly lensed SNe Ia in a 10 year LSST $z$-band search would produce more robust constraints on the validity of cosmic transparency (at the precision of $\Delta\beta=10^{-2}$), with a larger sample of unlensed SNe Ia detected in future LSST survey. We have also discussed the ways in which our methodology could be improved, with the combination of current and future available data in gravitational wave (GW) and electromagnetic (EM) domain.

Strongly gravitational lensing systems (SGL) encodes cosmology information in source/lens distance ratios $\mathcal{D}_{\rm obs}=\mathcal{D}_{\rm ls}/\mathcal{D}_{\rm s}$, which can be used to precisely constrain cosmological parameters. In this paper, based on future measurements of 390 strong lensing systems from the forthcoming LSST survey, we have successfully reconstructed the distance ratio $\mathcal{D}_{\rm obs}$ (with the source redshift $z_s\sim 4.0$), directly from the data without assuming any parametric form. A recently developed method based on model-independent reconstruction approach, Gaussian Processes (GP) is used in our study of these strong lensing systems. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. More specifically, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim0.015$ in the concordance $\Lambda$CDM model, which provides comparable constraint on $\Omega_m$ with Planck 2015 results. In the framework of the modified gravity theory (DGP), 390 detectable galactic lenses from future LSST survey would lead to stringent fits of $\Delta\Omega_m\sim0.030$. \textbf{Finally, we have discussed three possible sources of systematic errors (sample incompleteness, the determination of length of lens redshift bin, and the choice of lens redshift shells), and quantified their effects on the final cosmological constraints. Our results strongly indicate that future strong lensing surveys, with the accumulation of a larger and more accurate sample of detectable galactic lenses, will considerably benefit from the methodology described in this analysis. }

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

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.

In this paper, we present a scheme to investigate the opacity of the Universe in a cosmological-model-independent way, with the combination of current and future available data in gravitational wave (GW) and electromagnetic (EM) domain. In the FLRW metric, GWs propagate freely through a perfect fluid without any absorption and dissipation, which provides a distance measurement unaffected by the cosmic opacity. Focusing on the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), as well as the newly-compiled SNe Ia data (JLA and Pantheon sample), we find an almost transparent universe is strongly favored at much higher redshifts ($z\sim 2.26$). Our results suggest that, although the tests of cosmic opacity are not significantly sensitive to its parametrization, a strong degeneracy between the cosmic opacity parameter and the absolute \textit{B}-band magnitude of SNe Ia is revealed in this analysis. More importantly, we obtain that future measurements of the luminosity distances of gravitational waves sources will be much more competitive than the current analyses, which makes it expectable more vigorous and convincing constraints on the cosmic opacity (and consequently on background physical mechanisms) and a deeper understanding of the intrinsic properties of type Ia supernovae in a cosmological-model-independent way.

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.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

In this paper we analyse the implications of the latest cosmological data sets to test the Etherington's distance duality relation (DDR), which connects the luminosity distance $D_L$ and angular diameter distance $D_A$ at the same redshift. For $D_L$ we consider the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), which can be considered as standard candles (or standard siren), while the angular diameter distances $D_A$ are derived from the newly-compiled sample of compact radio quasars observed by very-long-baseline interferometry (VLBI), which represents a type of new cosmological standard ruler. Alleviating the absorption and scattering effects of dust in the Universe, this will create a valuable opportunity to directly test DDR at much higher precision with the combination of gravitational wave (GW) and electromagnetic (EM) signals. Our results show that, with the combination of the current radio quasar observations, the duality-distance relation can be verified at the precision of $10^{-2}$. Moreover, the Einstein Telescope ET would produce more robust constraints on the validity of such distance duality relation (at the precision of $10^{-3}$), with a larger sample of compact milliarcsecond radio quasars detected in future VLBI surveys.

In this paper, we assemble a well-defined sample of early-type gravitational lenses extracted from a large collection of 158 systems, and use the redshift distribution of galactic-scale lenses to test the standard cosmological model ($\Lambda$CDM) and the modified gravity theory (DGP). Two additional sub-samples are also included to account for possible selection effect introduced by the detectability of lens galaxies. Our results show that independent measurement of the matter density parameter ($\Omega_m$) could be expected from such strong lensing statistics. Based on future measurements of strong lensing systems from the forthcoming LSST survey, one can expect $\Omega_m$ to be estimated at the precision of $\Delta\Omega_m\sim 0.006$, which provides a better constraint on $\Omega_m$ than \textit{Planck} 2015 results. Moreover, use the lens redshift test is also used to constrain the characteristic velocity dispersion of the lensing galaxies, which is well consistent with that derived from the optical spectroscopic observations. A parameter $f_E$ is adopted to quantify the relation between the lensing-based velocity dispersion and the corresponding stellar value. Finally, the accumulation of detectable galactic lenses from future LSST survey would lead to more stringent fits of $\Delta f_E\sim10^{-3}$, which encourages us to test the global properties of early-type galaxies at much higher accuracy.

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.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

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.

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.

Aiming at exploring the nature of dark energy (DE), we use forty-three observational Hubble parameter data (OHD) in the redshift range $0 < z \leqslant 2.36$ to make a cosmological model-independent test of the $\Lambda$CDM model with two-point $Omh^2(z_{2};z_{1})$ diagnostic. In $\Lambda$CDM model, with equation of state (EoS) $w=-1$, two-point diagnostic relation $Omh^2 \equiv \Omega_m h^2$ is tenable, where $\Omega_m$ is the present matter density parameter, and $h$ is the Hubble parameter divided by 100 $\rm km s^{-1} Mpc^{-1}$. We utilize two methods: the weighted mean and median statistics to bin the OHD to increase the signal-to-noise ratio of the measurements. The binning methods turn out to be promising and considered to be robust. By applying the two-point diagnostic to the binned data, we find that although the best-fit values of $Omh^2$ fluctuate as the continuous redshift intervals change, on average, they are continuous with being constant within 1 $\sigma$ confidence interval. Therefore, we conclude that the $\Lambda$CDM model cannot be ruled out.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

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).

We focus on a series of $f(R)$ gravity theories in Palatini formalism to investigate the probabilities of producing the late-time acceleration for the flat Friedmann-Robertson-Walker (FRW) universe. We apply statefinder diagnostic to these cosmological models for chosen series of parameters to see if they distinguish from one another. The diagnostic involves the statefinder pair $\{r,s\}$, where $r$ is derived from the scale factor $a$ and its higher derivatives with respect to the cosmic time $t$, and $s$ is expressed by $r$ and the deceleration parameter $q$. In conclusion, we find that although two types of $f(R)$ theories: (i) $f(R) = R + \alpha R^m - \beta R^{-n}$ and (ii) $f(R) = R + \alpha \ln R - \beta$ can lead to late-time acceleration, their evolutionary trajectories in the $r-s$ and $r-q$ planes reveal different evolutionary properties, which certainly justify the merits of statefinder diagnostic. Additionally, we utilize the observational Hubble parameter data (OHD) to constrain these models of $f(R)$ gravity. As a result, except for $m=n=1/2$ of (i) case, $\alpha=0$ of (i) case and (ii) case allow $\Lambda$CDM model to exist in 1$\sigma$ confidence region. After adopting statefinder diagnostic to the best-fit models, we find that all the best-fit models are capable of going through deceleration/acceleration transition stage with late-time acceleration epoch, and all these models turn to de-Sitter point ($\{r,s\}=\{1,0\}$) in the future. Also, the evolutionary differences between these models are distinct, especially in $r-s$ plane, which makes the statefinder diagnostic more reliable in discriminating cosmological models.

In order to explore the generic properties of a backreaction model for explaining the accelerated expansion of the Universe, we exploit two metrics to describe the late time Universe. Since the standard FLRW metric cannot precisely describe the late time Universe on small scales, the template metric with an evolving curvature parameter $\kappa_{\mathcal{D}}(t)$ is employed. However, we doubt the validity of the prescription for $\kappa_{\mathcal{D}}$, which motivates us apply observational Hubble parameter data (OHD) to constrain parameters in dust cosmology. First, for FLRW metric, by getting best-fit constraints of $\Omega^{{\mathcal{D}}_0}_m = 0.25^{+0.03}_{-0.03}$, $n = 0.02^{+0.69}_{-0.66}$, and $H_{\mathcal{D}_0} = 70.54^{+4.24}_{-3.97}\ {\rm km \ s^{-1} \ Mpc^{-1}}$, the evolutions of parameters are explored. Second, in template metric context, by marginalizing over $H_{\mathcal{D}_0}$ as a prior of uniform distribution, we obtain the best-fit values of $n=-1.22^{+0.68}_{-0.41}$ and ${{\Omega}_{m}^{\mathcal{D}_{0}}}=0.12^{+0.04}_{-0.02}$. Moreover, we utilize three different Gaussian priors of $H_{\mathcal{D}_0}$, which result in different best-fits of $n$, but almost the same best-fit value of ${{\Omega}_{m}^{\mathcal{D}_{0}}}\sim0.12$. Also, the absolute constraints without marginalization of parameter are obtained: $n=-1.1^{+0.58}_{-0.50}$ and ${{\Omega}_{m}^{\mathcal{D}_{0}}}=0.13\pm0.03$. With these constraints, the evolutions of the effective deceleration parameter $q^{\mathcal{D}}$ indicate that the backreaction can account for the accelerated expansion of the Universe without involving extra dark energy component in the scaling solution context. Nevertheless, the results also verify that the prescription of $\kappa_{\mathcal{D}}$ is insufficient and should be improved.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

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.

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.

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.

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.

In this paper, by analyzing the dynamics of the inhomogeneous quintessence dark energy, we find that the gradient energy of dark energy will oscillate and gradually vanish, which indicates the gradient energy of the scalar field present in the early universe does not affect the current dynamics of the universe. Moreover, with the decaying of gradient energy, there exists a possible mutual transformation between kinetic energy and gradient energy. In the framework of interacting dark energy models, we argue that inhomogeneous dark energy may have a significant effect on the evolution of the cosmic background, the investigation of which still requires fully relativistic $N$-body numerical simulations in the future.

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.

In this paper we use a recently compiled data set, which comprises 118 galactic-scale strong gravitational lensing (SGL) systems to constrain the statistic property of SGL system, as well as the curvature of universe without assuming any fiducial cosmological model. Based on the singular isothermal ellipsoid (SIE) model of SGL system, we obtain that the constrained curvature parameter $\Omega_{\rm k}$ is close to zero from the SGL data, which is consistent with the latest result of planck measurement. More interestingly, we find that the parameter $f$ in the SIE model is strongly correlated with the curvature $\Omega_{\rm k}$. Neglecting this correlation in the analysis will significantly overestimate the constraining power of SGL data on the curvature. Furthermore, the obtained constraint on $f$ is different from previous results: $f=1.105\pm0.030$ ($68\%$ C.L.), which means that the standard singular isothermal sphere (SIS) model ($f=1$) is disfavored by the current SGL data at more than $3\sigma$ confidence level. We also divide the whole SGL data into two parts according to the centric stellar velocity dispersion $\sigma_{\rm c}$ and find that the larger value of $\sigma_{\rm c}$ the subsample has, the more favored the standard SIS model is. Finally, we extend the SIE model by assuming the 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}$, and obtain the constraints on the power-law indexes: $\alpha=1.95\pm0.04$ and $\delta=2.40\pm0.13$ at 68\% confidence level. When assuming the power-law index $\alpha=\delta=\gamma$, this scenario is totally disfavored by the current SGL data, $\chi^2_{\rm min,\gamma} - \chi^2_{\rm min,SIE} \simeq 53$.

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.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

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.

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.

In quantum theory of gravity, we expect the Lorentz Invariance Violation (LIV) and the modification of the dispersion relation between energy and momentum for photons. The effect of the energy-dependent velocity due to the modified dispersion relation for photons was studied in the standard cosmological context by using a sample of Gamma Ray Bursts (GRBs). In this paper we mainly discuss the possible LIV effect by using different cosmological models for the accelerating universe. Due to the degeneracies among model parameters, the GRBs' time delay data are combined with the cosmic microwave background data from the Planck first year release, the baryon acoustic oscillation data at six different redshifts, as well as Union2 type Ia supernovae data, to constrain both the model parameters and the LIV effect. We find no evidence of LIV.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

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.

We constrain the scalar field dark energy model with an inverse power-law potential, i.e., $V(\phi)\propto {\phi}^{-\alpha}$ ($\alpha>0$), from a set of recent cosmological observations by compiling an updated sample of Hubble parameter measurements including 30 independent data points. Our results show that the constraining power of the updated sample of $H(z)$ data with the HST prior on $H_0$ is stronger than those of the SCP Union2 and Union2.1 compilations. A recent sample of strong gravitational lensing systems is also adopted to confine the model even though the results are not significant. A joint analysis of the strong gravitational lensing data with the more restrictive updated Hubble parameter measurements and the Type Ia supernovae data from SCP Union2 indicates that the recent observations still can not distinguish whether dark energy is a time-independent cosmological constant or a time-varying dynamical component.

We present a new analysis of HST, Spitzer telescope imaging and VLT imaging and spectroscopic data of a bright lensed galaxy at $z$=1.0334 in the lensing cluster Abell~2667. Using this high-resolution imaging we present an updated lens model that allows us to fully understand the lensing geometry and reconstruct the lensed galaxy in the source plane. This giant arc gives a unique opportunity to peer into the structure of a high-redshift disk galaxy. We find that the lensed galaxy of Abell 2667 is a typical spiral galaxy with morphology similar to the structure of its counterparts at higher redshift $z\sim 2$. The surface brightness of the reconstructed source galaxy in the z$_{850}$ band reveals the central surface brightness $I(0)=20.28\pm0.22$ mag arcsec$^{-2}$ and the characteristic radius $r_s=2.01\pm0.16$ kpc at redshift $z \sim 1$. The morphological reconstruction in different bands shows obvious negative radial color gradients for this galaxy. Moreover, the redder central bulge tends to contain a metal-rich stellar population, rather than being heavily reddened by dust due to high and patchy obscuration. We analyze the VIMOS/IFU spectroscopic data and find that, in the given wavelength range ($\sim 1800-3200$ \AA), the combined arc spectrum of the source galaxy is characterized by a strong continuum emission with strong UV absorption lines (FeII and MgII) and shows the features of a typical starburst Wolf-Rayet galaxy NGC5253. More specifically, we have measured the EWs of FeII and MgII lines in the Abell 2667 spectrum, and obtained similar values for the same wavelength interval of the NGC5253 spectrum. Marginal evidence for CIII] 1909 emission at the edge of the grism range further confirms our expectation.

Using a relatively complete observational data concerning four angular diameter distance (ADD) measurements and %synthetic combined SN+GRB observations representing current luminosity distance (LD) data, this paper investigates the %tension between compatibility of these two cosmological distances considering three classes of dark energy equation of state (EoS) reconstruction. In particular, we use strongly gravitationally lensed systems from various large systematic gravitational lens surveys and galaxy clusters, which yield the Hubble constant independent ratio between two angular diameter distances $D_{ls}/D_s$ data. Our results demonstrate that, with more general categories of standard ruler data, ADD and LD data are compatible at $1\sigma$ level. Secondly, we note that consistency between ADD and LD data %are blind is maintained irrespective of the EoS parameterizations: there is a good match between the universally explored CPL model and other formulations of cosmic equation of state. Especially for the truncated GEoS model with $\beta=-2$, the conclusions obtained with ADD and LD are almost the same. Finally, statistical analysis of generalized dark energy equation of state performed on four classes of ADD data provides stringent constraints on the EoS parameters $w_0$, $w_{\beta}$ and $\beta$, which suggest that dark energy was a subdominant component at early times. Moreover, the GEoS parametrization with $\beta\simeq 1$ seems to be a more favorable two-parameter model to characterize the cosmic equation of state, because the combined angular diameter distance data (SGL+CBF+BAO+WMAP9) provide the best-fit value $\beta=0.751^{+0.465}_{-0.480}$.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

In order to test if there is energy transfer between dark energy and dark matter, we investigate cosmological constraints on two forms of nontrivial interaction between the dark matter sector and the sector responsible for the acceleration of the universe, in light of the newly revised observations including OHD, CMB, BAO and SNe Ia. More precisely, we find the same tendencies for both phenomenological forms of the interaction term $Q=3\gamma H\rho$, i.e., the parameter $\gamma$ to be a small number, $|\gamma|\approx 10^{-2}$. However, concerning the sign of the interaction parameter, we observe that $\gamma>0$ when the interaction between dark sectors is proportional to the energy density of dust matter, whereas the negative coupling ($\gamma<0$) is preferred by observations when the interaction term is proportional to dark energy density. We further discuss two possible explanations to this incompatibility and apply a quantitative criteria to judge the severity of the coincidence problem. Results suggest that the $\gamma_m$IDE model with a positive coupling may alleviate the coincidence problem, since its coincidence index $C$ is smaller than that for the $\gamma_d$IDE model, the interacting quintessence and phantom models by four orders of magnitude.

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.

We present a multi-wavelength study of the gravitational lens COSMOS J095930+023427 (z=0.89), together with the associated galaxy group located at $z\sim0.7$ along the line of sight and the lensed background galaxy. The source redshift is currently unknown, but estimated to be at $z_s \sim 2$. The analysis is based on the available public HST, Subaru, Chandra imaging data, and VLT spectroscopy. The lensing system is an early-type galaxy showing a strong [OII] emission line, and produces 4 bright images of the distant background source. It has an Einstein radius of 0.79", about 4 times large than the effective radius. We perform a lensing analysis using both a Singular Isothermal Ellipsoid (SIE) and a Peudo-Isothermal Elliptical Mass Distribution (PIEMD) for the lensing galaxy, and find that the final results on the total mass, the dark matter (DM) fraction within the Einstein radius and the external shear due to a foreground galaxy group are robust with respect of the choice of the parametric model and the source redshift (yet unknown). We measure the luminous mass from the photometric data, and find the DM fraction within the Einstein radius $f_{\rm DM}$ to be between $0.71\pm 0.13$ and $0.79 \pm 0.15$, depending on the unknown source redshift. Meanwhile, the non-null external shear found in our lensing models supports the presence and structure of a galaxy group at $z\sim0.7$, and an independent measurement of the 0.5-2 keV X-ray luminosity within 20" around the X-ray centroid provides a group mass of $M=(3-10)\times 10^{13}$ M$_{\odot}$, in good agreement with the previous estimate derived through weak lensing analysis.

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 study the redshift distribution of two samples of early-type gravitational lenses, extracted from a larger collection of 122 systems, to constrain the cosmological constant in the LCDM model and the parameters of a set of alternative dark energy models (XCDM, Dvali-Gabadadze-Porrati and Ricci dark energy models), under a spatially flat universe. The likelihood is maximized for $\Omega_\Lambda= 0.70 \pm 0.09$ when considering the sample excluding the SLACS systems (known to be biased towards large image-separation lenses) and no-evolution, and $\Omega_\Lambda= 0.81\pm 0.05$ when limiting to gravitational lenses with image separation larger than 2" and no-evolution. In both cases, results accounting for galaxy evolution are consistent within 1$\sigma$. The present test supports the accelerated expansion, by excluding the null-hypothesis (i.e., $\Omega_\Lambda = 0 $) at more than 4$\sigma$, regardless of the chosen sample and assumptions on the galaxy evolution. A comparison between competitive world models is performed by means of the Bayesian information criterion. This shows that the simplest cosmological constant model - that has only one free parameter - is still preferred by the available data on the redshift distribution of gravitational lenses. We perform an analysis of the possible systematic effects, finding that the systematic errors due to sample incompleteness, galaxy evolution and model uncertainties approximately equal the statistical errors, with present-day data. We find that the largest sources of systemic errors are the dynamical normalization and the high-velocity cut-off factor, followed by the faint-end slope of the velocity dispersion function.

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.

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.

In this paper, we propose an accurate test of the distance-duality (DD) relation, $\eta=D_{L}(z)(1+z)^{-2}/D_{A}(z)=1$ (where $D_{L}$ and $D_{A}$ are the luminosity distances and angular diameter distances, respectively), with a combination of cosmological observational data of Type Ia Supernave (SNe Ia) from Union2 set and the galaxy cluster sample under an assumption of spherical model. In order to avoid bias brought by redshift incoincidence between observational data and to consider redshift error bars of both clusters and SNe Ia in analysis, we carefully choose the SNe Ia points which have the minimum acceptable redshift difference of the galaxy cluster sample ($|\Delta z|_{\rm min} =\sigma_{z, \rm SN}+\sigma_{z, \rm cluster}$). By assuming $\eta$ a constant and functions of the redshift parameterized by six different expressions, we find that there exists no conceivable evidence for variations in the DD relation concerning with observational data, since it is well satisfied within $1\sigma$ confidence level for most cases. Further considering different values of $\Delta z$ in constraining, we also find that the choosing of $\Delta z$ may play an important role in this model-independent test of the distance-duality relation for the spherical sample of galaxy clusters.

In order to test the possible interaction between dark energy and dark matter, we investigate observational constraints on a phenomenological scenario, in which the ratio between the dark energy and matter densities is proportional to the power law case of the scale factor, $r\equiv (\rho_X/\rho_m)\propto a^{\xi}$. By using the Markov chain Monte Carlo method, we constrain the phenomenological interacting dark energy model with the newly revised $H(z)$ data, as well as the cosmic microwave background (CMB) observation from the 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) results, the baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) galaxy sample and the type Ia supernovae (SNe Ia) from Union2 set. The best-fit values of the model parameters are $\Omega_{m0}=0.27_{-0.02}^{+0.02}(1\sigma)_{-0.03}^{+0.04}(2\sigma)$, $\xi=3.15_{-0.50}^{+0.48}(1\sigma)_{-0.71}^{+0.72}(2\sigma)$, and $w_X=-1.05_{-0.14}^{+0.15}(1\sigma)_{-0.21}^{+0.21}(2\sigma)$, which are more stringent than previous results. These results show that the standard $\Lambda$CDM model without any interaction remains a good fit to the recent observational data; however, the interaction that the energy transferring from dark matter to dark energy is slightly favored over the interaction from dark energy to dark matter. It is also shown that the $H(z)$ data can give more stringent constraints on the phenomenological interacting scenario when combined to CMB and BAO observations, and the confidence regions of $H(z)$+BAO+CMB, SNe+BAO+CMB, and $H(z)$+SNe+BAO+CMB combinations are consistent with each other.

Strong lensing has developed into an important astrophysical tool for probing both cosmology and galaxies (their structures, formations, and evolutions). Now several hundreds of strong lens systems produced by massive galaxies have been discovered, which may form well-defined samples useful for statistical analyses. To collect a relatively complete lens redshift data from various large systematic surveys of gravitationally lensed quasars and check the possibility to use it as a future complementarity to other cosmological probes. We use the distribution of gravitationally-lensed image separations observed in the Cosmic Lens All-Sky Survey (CLASS), the PMN-NVSS Extragalactic Lens Survey (PANELS), the Sloan Digital Sky Survey (SDSS) and other surveys, considering a singular isothermal ellipsoid (SIE) model for galactic potentials as well as improved new measurements of the velocity dispersion function of galaxies based on the SDSS DR5 data and recent semi-analytical modeling of galaxy formation, to constrain two dark energy models ($\Lambda$CDM and constant $w$) under a flat universe assumption. We find that the current lens redshift data give a relatively weak constraint on the model parameters. However, by combing the redshift data with the baryonic acoustic oscillation peak and the comic macrowave background data, we obtain more stringent results, which show that the flat $\Lambda$ CDM model is still included at 1$\sigma$.

The interacting dark matter (IDM) scenario allows for the acceleration of the Universe without Dark Energy. We constrain the IDM model by using the newly revised observational data including $H(z)$ data and Union2 SNe Ia via the Markov chain Monte Carlo method. When mimicking the $\Lambda$CDM model, we obtain a more stringent upper limit to the effective annihilation term at $\kappa C_1\approx 10^{-3.4}\rm{Gyr}^{-1}$, and a tighter lower limit to the relevant mass of Dark Matter particles at $M_x\approx 10^{-8.6}\rm{Gev}$. When mimicking the $w$CDM model, we find that the effective equation of state of IDM is consistent with the concordance $\Lambda$CDM model and appears to be most consistent with the effective phantom model with a constant EoS for which $w<-1$.

The validity of distance duality relation, $\eta=D_L(z)(1+z)^{-2}/D_A(z)=1$, an exact result required by the Etherington reciprocity theorem, where $D_A(z)$ and $D_L(z)$ are the angular and luminosity distances, plays an essential part in cosmological observations and model constraints. In this paper, we investigate some consequences of such a relation by assuming $\eta$ a constant or a function of the redshift. In order to constrain the parameters concerning $\eta$, we consider two groups of cluster gas mass fraction data including 52 X-ray luminous galaxy clusters observed by Chandra in the redshift range $0.3\sim 1.273$ and temperature range $T_{\rm gas}> 4$ keV, under the assumptions of two different temperature profiles [1]. We find that the constant temperature profile is in relatively good agreement with no violation of the distance duality relation for both parameterizations of $\eta$, while the one with temperature gradient (the Vikhlinin et al. temperature profile) seems to be incompatible even at 99% CL.

The activity of active galaxy may be triggered by the merge of galaxies and present-day galaxies are probably the product of successive minor mergers. The frequent galactic merges at high redshift imply that active galaxy harbors supermassive unequal-mass binary black holes in its center at least once during its life time. In this paper, we showed that the recently discovered double-lobed FR II radio galaxies are the remnants of such supermassive binary black holes. The inspiraling secondary black hole opens a gap in the accretion disk and removes the inner accretion disk when it merges into the primary black hole, leaving a big hole of about several hundreds of Schwarzschild radius in the vicinity of the post-merged supermassive black hole and leading to an interruption of jet formation. When the outer accretion disk slowly refills the big hole on a viscous time scale, the jet formation restarts and the interaction of the recurrent jets and the inter-galactic medium forms a secondary pair of lobes. We applied the model to a particular double-lobed radio source B1834+620 and showed that the orbit of the secondary is elliptical with a typical eccentricity $e \simeq 0.68$ and the mass ratio $q$ of the secondary and the primary is $0.01 \la q \la 0.4$. The accretion disk is a standard $\alpha$-disk with $0.01 \la \alpha \la 0.04$ and the ratio of disk half height $H$ and radius $r$ is $\delta \simeq 0.01$. The model predicates that double-lobed radio structure forms only in FR II radio galaxies.

In this paper, the correlation between the gamma-ray and the radio bands is investigated. The results show that there is a closer correlation between the gamma-ray emission and the high frequency (1.3mm, 230GHz) radio emission for maximum data than between the gamma-ray and the lower frequency (5GHz) radio emissions, which means that the gamma-ray is associated with the radio emission from the jet.

In this paper, the long-term historical optical (UBVRI) and near-infrared (JHK) data are presented with some new observations in the optical (February 1994-January 1995) and near-infrared (November 1995) bands included for BL Lac object OJ287. The new optical data in V-band are in agreement with the results reported by other authors (Sillanpaa et al. 1996a; Arimoto et al. 1997), a close correlation between the color index of B-V and the magnitude V has been obtained from our new observations; The new infrared observations presented here indicate that the source was at a high level in the infrared band during the observation period; From the available literature, we have got that the largest variations for UBVRIJHK bands are respectively: 4.72mag, 5.93mag, 5.18mag, 4.45mag, 4.07mag, 3.87mag, 3.78mag, and 3.54mag. A strong correlation is found between the optical and near-infrared bands when the DCF method is used, which suggests that these two bands have the same emission mechanism.