The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

The H0LiCOW collaboration inferred via gravitational lensing time delays a Hubble constant $H_0=73.3^{+1.7}_{-1.8}$ km s$^{-1}{\rm Mpc}^{-1}$, describing deflector mass density profiles by either a power-law or stars plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in $H_0$. We quantify any potential effect of the MST with flexible mass models that are maximally degenerate with H0. Our calculation is based on a new hierarchical approach in which the MST is only constrained by stellar kinematics. The approach is validated on hydrodynamically simulated lenses. We apply the method to the TDCOSMO sample of 7 lenses (6 from H0LiCOW) and measure $H_0=74.5^{+5.6}_{-6.1}$ km s$^{-1}{\rm Mpc}^{-1}$. In order to further constrain the deflector mass profiles, we then add imaging and spectroscopy for 33 strong gravitational lenses from the SLACS sample. For 9 of the SLAC lenses we use resolved kinematics to constrain the stellar anisotropy. From the joint analysis of the TDCOSMO+SLACS sample, we measure $H_0=67.4^{+4.1}_{-3.2}$ km s$^{-1}{\rm Mpc}^{-1}$, assuming that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. While our new analysis does not statistically invalidate the mass profile assumptions by H0LiCOW, and thus their $H_0$ measurement relying on those, it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on $H_0$ derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data, chiefly spatially resolved kinematics of lens galaxies.

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We report the direct observation of neutrinos produced in the carbon-nitrogen-oxygen (CNO) fusion cycle in the Sun with the Borexino detector at the Laboratori Nazionali del Gran Sasso in Italy. This is the first experimental evidence of the existence of such reaction sequence in a star. The CNO solar neutrino interaction rate is $7.2^{+3.0}_{-1.7}$ counts per day per $100$ tonnes of target at 68% C.L., corresponding to a flux of neutrinos on Earth of $7.0^{+3.0}_{-2.0}\times10^8$ cm$^{-2}$ s$^{-1}$. The absence of CNO signal is disfavoured at 5$\sigma$.

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets presented in paper I. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

One of the main challenges in using high redshift active galactic nuclei to study the correlations between the mass of the supermassive Black Hole (MBH) and the properties of their active host galaxies is instrumental resolution. Strong lensing magnification effectively increases instrumental resolution and thus helps to address this challenge. In this work, we study eight strongly lensed active galactic nuclei (AGN) with deep Hubble Space Telescope imaging, using the lens modelling code Lenstronomy to reconstruct the image of the source. Using the reconstructed brightness of the host galaxy, we infer the host galaxy stellar mass based on stellar population models. MBH are estimated from broad emission lines using standard methods. Our results are in good agreement with recent work based on non-lensed AGN, providing additional evidence that the correlation evolves over cosmic time. At the moment, the sample size of lensed AGN is small and thus they provide mostly a consistency check on systematic errors related to resolution for the non-lensed AGN. However, the number of known lensed AGN is expected to increase dramatically in the next few years, through dedicated searches in ground and space based wide field surveys, and they may become a key diagnostic of black hole and galaxy co-evolution.

Neutrinos emitted in the carbon, nitrogen, oxygen (CNO) fusion cycle in the Sun are a sub-dominant, yet crucial component of solar neutrinos whose flux has not been measured yet. The Borexino experiment at the Laboratori Nazionali del Gran Sasso (Italy) has a unique opportunity to detect them directly thanks to the detector's radiopurity and the precise understanding of the detector backgrounds. We discuss the sensitivity of Borexino to CNO neutrinos, which is based on the strategies we adopted to constrain the rates of the two most relevant background sources, $pep$ neutrinos and $^{210}$Bi beta decays. Finally, the relevance of a detection of CNO neutrinos is extensively discussed in the context of the solar metallicity problem.

We examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analytic model (SAM). After applying an observational selection function to the simulations, we find that both MBII and SAM agree well with the data, in terms of the central distribution. However, the dispersion in the mass ratio between black hole mass and stellar mass is significantly more consistent with the MBII prediction ($\sim0.3~$dex), than with the SAM ($\sim0.7~$dex), even when accounting for observational uncertainties. Hence, our observations can distinguish between the different recipes adopted in the models. The mass relations in the MBII are highly dependent on AGN feedback while the relations in the SAM are more sensitive to galaxy merger events triggering nuclear activity. Moreover, the intrinsic scatter in the mass ratio of our high-$z$ sample is comparable to that observed in the local sample, all but ruling out the proposed scenario the correlations are purely stochastic in nature arising from some sort of cosmic central limit theorem. Our results support the hypothesis of AGN feedback being responsible for a causal link between the SMBH and its host galaxy, resulting in a tight correlation between their respective masses.

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.

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

We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

We report on an improved measurement of the $^8$B solar neutrino interaction rate with the Borexino experiment at the Laboratori Nazionali del Gran Sasso. Neutrinos are detected via their elastic scattering on electrons in a large volume of liquid scintillator. The measured rate of scattered electrons above 3 MeV of energy is $0.223\substack{+0.015 \\ -0.016}\,(stat)\,\substack{+0.006 \\ -0.006}\,(syst)$ cpd/100 t, which corresponds to an observed solar neutrino flux assuming no neutrino flavor conversion of $\Phi\substack{\rm ES \\ ^8\rm B}=2.57\substack{+0.17 \\ -0.18}(stat)\substack{+0.07\\ -0.07}(syst)\times$10$^6$ cm$^{-2}\,$s$^{-1}$. This measurement exploits the active volume of the detector in almost its entirety for the first time, and takes advantage of a reduced radioactive background following the 2011 scintillator purification campaign and of novel analysis tools providing a more precise modeling of the background. Additionally, we set a new limit on the interaction rate of solar $hep$ neutrinos, searched via their elastic scattering on electrons as well as their neutral current-mediated inelastic scattering on carbon, $^{12}$C($\nu,\nu'$)$^{12}$C* ($E_{\gamma}$= 15.1 MeV).

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

The Borexino detector measures solar neutrino fluxes via neutrino-electron elastic scattering. Observed spectra are determined by the solar-$\nu_{e}$ survival probability $P_{ee}(E)$, and the chiral couplings of the neutrino and electron. Some theories of physics beyond the Standard Model postulate the existence of Non-Standard Interactions (NSI's) which modify the chiral couplings and $P_{ee}(E)$. In this paper, we search for such NSI's, in particular, flavor-diagonal neutral current interactions that modify the $\nu_e e$ and $\nu_\tau e$ couplings using Borexino Phase II data. Standard Solar Model predictions of the solar neutrino fluxes for both high- and low-metallicity assumptions are considered. No indication of new physics is found at the level of sensitivity of the detector and constraints on the parameters of the NSI's are placed. In addition, with the same dataset the value of $\sin^2\theta_W$ is obtained with a precision comparable to that achieved in reactor antineutrino experiments.

We present the first simultaneous measurement of the interaction rates of $pp$, $^7$Be, and $pep$ solar neutrinos performed with a global fit to the Borexino data in an extended energy range (0.19-2.93)$\,$MeV. This result was obtained by analyzing 1291.51$\,$days of Borexino Phase-II data, collected between December 2011 and May 2016 after an extensive scintillator purification campaign. We find: rate($pp$)$\,$=$\,$$134$$\,$$\pm$$\,$$10$$\,$($stat$)$\,$$^{\rm +6}_{\rm -10}$$\,$($sys$)$\,$cpd/100$\,$t, rate($^7$Be)$\,$=$\,$$48.3$$\,$$\pm$$\,$$1.1$$\,$($stat$)$\,$$^{\rm +0.4}_{\rm -0.7}$$\,$($sys$)$\,$cpd/100$\,$t, and rate($pep$)$\,$=$\,$$2.43$$\pm$$\,$$0.36$$\,$($stat$)$^{+0.15}_{-0.22}$$\,$($sys$)$\,$cpd/100$\,$t. These numbers are in agreement with and improve the precision of our previous measurements. In particular, the interaction rate of $^7$Be $\nu$'s is measured with an unprecedented precision of 2.7%, showing that discriminating between the high and low metallicity solar models is now largely dominated by theoretical uncertainties. The absence of $pep$ neutrinos is rejected for the first time at more than 5$\,$$\sigma$. An upper limit of $8.1$$\,$cpd/100$\,$t (95%$\,$C.L.) on the CNO neutrino rate is obtained by setting an additional constraint on the ratio between the $pp$ and $pep$ neutrino rates in the fit. This limit has the same significance as that obtained by the Borexino Phase-I (currently providing the tightest bound on this component), but is obtained by applying a less stringent constraint on the $pep$ $\nu$ flux.

Correlations between the mass of a supermassive black hole and the properties of its host galaxy (e.g., total stellar mass (M*), luminosity (Lhost)) suggest an evolutionary connection. A powerful test of a co-evolution scenario is to measure the relations MBH-Lhost and MBH-M* at high redshift and compare with local estimates. For this purpose, we acquired HST imaging with WFC3 of 32 X-ray-selected broad-line AGN at 1.2<z<1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the 2D model fitting. The black hole mass was determined using the broad Halpha line, detected in the near-infrared with Subaru/FMOS, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of MBH to total M* is 2.7 times larger at z~1.5 than in the local universe, while the scatter is equivalent between the two epochs. A non-evolving mass ratio is consistent with the data at the 2-3 sigma confidence level when accounting for selection effects and their uncertainties. The relationship between MBH-Lhost paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total M* and Lhost proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio MBH/M* is offset from the local value by a factor of ~7 which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum supported component of the galaxy to the pressure supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.

We report an experiment to test quantum interference, entanglement and nonlocality using two dissimilar photon sources, the Sun and a semiconductor quantum dot on the Earth, which are separated by 150 million kilometers. By making the otherwise vastly distinct photons indistinguishable all degrees of freedom, we observe time-resolved two-photon quantum interference with a raw visibility of 0.796(17), well above the 0.5 classical limit, providing the first evidence of quantum nature of thermal light. Further, using the photons with no common history, we demonstrate post-selected two-photon entanglement with a state fidelity of 0.826(24), and a violation of Bell's inequality by 2.20(6). The experiment can be further extended to a larger scale using photons from distant stars, and open a new route to quantum optics experiments at an astronomical scale.

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

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.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

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

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

Fast Radio bursts (FRBs), bright transients with millisecond durations at $\sim$ GHz and typical redshifts probably $>0.8$, are likely to be gravitationally lensed by intervening galaxies. Since the time delay between images of strongly lensed FRB can be measured to extremely high precision because of the large ratio $\sim10^9$ between the typical galaxy-lensing delay time $\sim\mathcal{O}$(10 days) and the width of bursts $\sim\mathcal{O}$(ms), we propose strongly lensed FRBs as precision probes of the universe. We show that, within the flat $\Lambda$CDM model, the Hubble constant $H_0$ can be constrained with a $\sim0.91\%$ uncertainty from 10 such systems probably observed with the Square Kilometer Array (SKA) in $<$ 30 years. More importantly, the cosmic curvature can be model-independently constrained to a precision of $\sim0.076$. This constraint can directly test the validity of the cosmological principle and break the intractable degeneracy between the cosmic curvature and dark energy.

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

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.

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

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

Strongly lensed active galactic nuclei (AGN) provide a unique opportunity to make progress in the study of the evolution of the correlation between the mass of supermassive black holes ($\mathcal M_{BH}$) and their host galaxy luminosity ($L_{host}$). We demonstrate the power of lensing by analyzing two systems for which state-of-the-art lens modelling techniques have been applied to Hubble Space Telescope imaging data. We use i) the reconstructed images to infer the total and bulge luminosity of the host and ii) published broad-line spectroscopy to estimate $\mathcal M_{BH}$ using the so-called virial method. We then enlarge our sample with new calibration of previously published measurements to study the evolution of the correlation out to z~4.5. Consistent with previous work, we find that without taking into account passive luminosity evolution, the data points lie on the local relation. Once passive luminosity evolution is taken into account, we find that BHs in the more distant Universe reside in less luminous galaxies than today. Fitting this offset as $\mathcal M_{BH}$/$L_{host}$ $\propto$ (1+z)$^{\gamma}$, and taking into account selection effects, we obtain $\gamma$ = 0.6 $\pm$ 0.1 and 0.8$\pm$ 0.1 for the case of $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively. To test for systematic uncertainties and selection effects we also consider a reduced sample that is homogeneous in data quality. We find consistent results but with considerably larger uncertainty due to the more limited sample size and redshift coverage ($\gamma$ = 0.7 $\pm$ 0.4 and 0.2$\pm$ 0.5 for $\mathcal M_{BH}$-$L_{bulge}$ and $\mathcal M_{BH}$-$L_{total}$, respectively), highlighting the need to gather more high-quality data for high-redshift lensed quasar hosts. Our result is consistent with a scenario where the growth of the black hole predates that of the host galaxy.

We present the results of a low-energy neutrino search using the Borexino detector in coincidence with the gravitational wave (GW) events GW150914, GW151226 and GW170104. We searched for correlated neutrino events with energies greater than 250 keV within a time window of $\pm500$ s centered around the GW detection time. A total of five candidates were found for all three GW150914, GW151226 and GW170104. This is consistent with the number of expected solar neutrino and background events. As a result, we have obtained the best current upper limits on the GW event neutrino fluence of all flavors ($\nu_e, \nu_{\mu}, \nu_{\tau}$) in the energy range $(0.5 - 5.0)$ MeV.

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

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

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.

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

GRB 080503 is a short gamma-ray burst (GRB) detected by \emph{Swift} and has been classified as a compact-star-merger-origin GRB. The soft extended emission and the simultaneous late re-brightening in both the X-ray and optical afterglow lightcurves raise interesting questions regarding its physical origin. We show that the broad-band data of GRB 080503 can be well explained within the framework of the double neutron star merger model, provided that the merger remnant is a rapidly-rotating massive neutron star with an extremely high magnetic field (i.e. a millisecond magnetar). We show that the late optical re-brightening is consistent with the emission from a magnetar-powered "merger-nova". This adds one more case to the growing sample of merger-novae associated with short GRBs. The soft extended emission and the late X-ray excess emission are well connected through a magnetar dipole spin-down luminosity evolution function, suggesting that direct magnetic dissipation is the mechanism to produce these X-rays. The X-ray emission initially leaks from a hole in the merger ejecta pierced by the short GRB jet. The hole subsequently closes after the magnetar spins down and the magnetic pressure drops below ram pressure. The X-ray photons are then trapped behind the mergernova ejecta until the ejecta becomes optically thin at a later time. This explains the essentially simultaneous re-brightening in both the optical and X-ray lightcurves. Within this model, future gravitational wave sources could be associated with a bright X-ray counterpart along with the mergernova, even if the short GRB jet beams away from Earth.

We perform unbiased tests for the clumpiness of universe by confronting the Zel'dovich-Kantowski-Dyer-Roeder luminosity distance which describes the effect of local inhomogeneities on the propagation of light with the observational one estimated from measurements of standard candles, i.e., type Ia supernovae (SNe Ia) and gamma-ray bursts (GRBs). Methodologically, we first determine the light-curve fitting parameters which account for distance estimation in SNe Ia observations and luminosity/energy relations which are responsible for distance estimation of GRBs in the global fit to reconstruct the Hubble diagrams in the context of a clumpy universe. Subsequently, these Hubble diagrams allow us to achieve unbiased constraints on the matter density parameter $\Omega_m$ as well as clumpiness parameter $\eta$ which quantifies the fraction of homogeneously distributed matter within a given light cone. At 1$\sigma$ confidence level, the constraints are $\Omega_m=0.34\pm0.02$ and $\eta=1.00^{+0.00}_{-0.02}$ from the joint analysis. The results suggest that the universe full of Friedman-Lema\^{i}tre-Robertson-Walker fluid is favored by observations of standard candles with very high statistical significance. On the other hand, they may also indicate that the Zel'dovich-Kantowski-Dyer-Roeder approximation is a sufficient accurate form to describe the effects of local homogeneity in the expanding universe.

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.

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

The Palomar Transient Factory (PTF) team recently reported the discovery of a rapidly fading optical transient source, PTF11agg. A long-lived scintillating radio counterpart was identified, but the search for a high energy counterpart showed negative results. The PTF team speculated that PTF11agg may represent a new class of relativistic outbursts. Here we suggest that a neutron star (NS)-NS merger system with a supra-massive magnetar central engine could be a possible source to power such a transient, if our line of sight is not on the jet axis direction of the system. These systems are also top candidates for gravitational wave sources to be detected in the advanced LIGO/Virgo era. We find that the PTF11agg data could be explained well with such a model, suggesting that at least some gravitational wave bursts due to NS-NS mergers may be associated with such a bright electromagnetic counterpart without a \gamma-ray trigger.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

If double neutron star mergers leave behind a massive magnetar rather than a black hole, a bright early afterglow can follow the gravitational wave burst (GWB) even if there is no short gamma-ray burst (SGRB) - GWB association or there is an association but the SGRB does not beam towards earth. Besides directly dissipating the proto-magnetar wind as suggested by Zhang, we here suggest that the magnetar wind could push the ejecta launched during the merger process, and under certain conditions, would reach a relativistic speed. Such a magnetar-powered ejecta, when interacting with the ambient medium, would develop a bright broad-band afterglow due to synchrotron radiation. We study this physical scenario in detail, and present the predicted X-ray, optical and radio light curves for a range of magnetar and ejecta parameters. We show that the X-ray and optical lightcurves usually peak around the magnetar spindown time scale (10^3-10^5s), reaching brightness readily detectable by wide-field X-ray and optical telescopes, and remain detectable for an extended period. The radio afterglow peaks later, but is much brighter than the case without a magnetar energy injection. Therefore, such bright broad-band afterglows, if detected and combined with GWBs in the future, would be a probe of massive millisecond magnetars and stiff equation-of-state for nuclear matter.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.