The non-detection of dark matter (DM) particles in increasingly stringent laboratory searches has encouraged alternative gravity theories where gravity is sourced only from visible matter. Here, we consider whether such theories can pass a two-dimensional test posed by gravitational lensing -- to reproduce a particularly detailed Einstein ring in the core of the galaxy cluster Abell 3827. We find that when we require the lensing mass distribution to strictly follow the shape (ellipticity and position angle) of the light distribution of cluster member galaxies, intracluster stars, and the X-ray emitting intracluster medium, we cannot reproduce the Einstein ring, despite allowing the mass-to-light ratios of these visible components to freely vary with radius to mimic alternative gravity theories. Alternatively, we show that the detailed features of the Einstein ring are accurately reproduced by allowing a smooth, freely oriented DM halo in the lens model, with relatively small contributions from the visible components at a level consistent with their observed brightnesses. This dominant DM component is constrained to have the same orientation as the light from the intracluster stars, indicating that the intracluster stars trace the gravitational potential of this component. The Einstein ring of Abell 3827 therefore presents a new challenge for alternative gravity theories: not only must such theories find agreement between the total lensing mass and visible mass, but they must also find agreement between the projected sky distribution of the lensing mass and that of the visible matter, a more stringent test than has hitherto been posed by lensing data.

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.

We explore the possibility of using amplitude and phase fluctuations of gravitational waves due to gravitational lensing as a probe of the small-scale matter power spectrum. The direct measurement of the small-scale matter power spectrum is made possible by making use of the frequency dependence of such gravitational lensing dispersions originating from the wave optics nature of the propagation of gravitational waves. We first study the small-scale behavior of the matter power spectrum in detail taking the so-called halo model approach including effects of baryons and subhalos. We find that the matter power spectrum at the wavenumber $k\sim 10^6h{\rm Mpc}^{-1}$ is mainly determined by the abundance of dark low mass halos with mass $1h^{-1}M_\odot \lesssim M \lesssim 10^4h^{-1}M_\odot$ and is relatively insensitive to baryonic effects. The matter power spectrum at this wavenumber is probed by gravitational lensing dispersions of gravitational waves at frequencies of $f\sim 0.1-1$~Hz with predicted signals of $\mathcal{O}(10^{-3})$. We also find that Primordial Black Holes (PBHs) with $M_{\rm PBH} \gtrsim 0.1~M_\odot$ can significantly enhance the matter power spectrum at $k \gtrsim 10^5h{\rm Mpc}^{-1}$ due to both the enhanced halo formation and the shot noise from PBHs. We find that gravitational lensing dispersions at $f\sim 10-100$~Hz are particularly sensitive to PBHs and can be enhanced by more than an order of magnitude depending on the mass and abundance of PBHs.

Gravitationally lensed quasars are useful for studying astrophysics and cosmology, and enlarging the sample size of lensed quasars is important for multiple studies. In this work, we develop a lens search algorithm for four-image (quad) lensed quasars based on their time variability. In the development of the lens search algorithm, we constructed a pipeline simulating multi-epoch images of lensed quasars in cadenced surveys, accounting for quasar variabilities, quasar hosts, lens galaxies, and the PSF variation. Applying the simulation pipeline to the Hyper Suprime-Cam (HSC) transient survey, we generated HSC-like difference images of the mock lensed quasars from Oguri & Marshall's lens catalog. We further developed a lens search algorithm that picks out variable objects as lensed quasar candidates based on their spatial extent in the difference images. We tested our lens search algorithm with the mock lensed quasars and variable objects from the HSC transient survey. Using difference images from multiple epochs, our lens search algorithm achieves a high true-positive rate (TPR) of 90.1% and a low false-positive rate (FPR) of 2.3% for the bright quads with wide separation. With a preselection of the number of blobs in the difference image, we obtain a TPR of 97.6% and a FPR of 2.6% for the bright quads with wide separation. Even when difference images are only available in one single epoch, our lens search algorithm can still detect the bright quads with wide separation at high TPR of 97.6% and low FPR of 2.4% in the optimal seeing scenario, and at TPR of $\sim94%$ and FPR of $\sim5%$ in typical scenarios. Therefore, our lens search algorithm is promising and is applicable to ongoing and upcoming cadenced surveys, particularly the HSC transient survey and the Rubin Observatory Legacy Survey of Space and Time, for finding new lensed quasar systems. [abridged]

We present spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27-022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s = 0.804$ and $z_l = 0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50+020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s = 2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_l = 0.439$. SDSS J135944.21+012809.8 has two point sources with quasar spectra at the same redshift $z_s = 1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that "observed" clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

We perform a self-calibration of the richness-to-mass ($N$-$M$) relation of CAMIRA galaxy clusters with richness $N\geq15$ at redshift $0.2\leq z<1.1$ by modeling redshift-space two-point correlation functions, namely, the auto-correlation function $\xi_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $\xi_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $\xi_{\mathrm{cg}}$ between these two samples. We focus on constraining the normalization $A_{\mathrm{N}}$ of the $N$-$M$ relation in a forward-modeling approach, carefully accounting for the redshift-space distortion, the Finger-of-God effect, and the uncertainty in photometric redshifts of CAMIRA clusters. The modeling also takes into account the projection effect on the halo bias of CAMIRA clusters. The parameter constraints are shown to be unbiased according to validation tests using a large set of mock catalogs constructed from N-body simulations. At the pivotal mass $M_{500}=10^{14}h^{-1}M_{\odot}$ and the pivotal redshift $z_{\mathrm{piv}} = 0.6$, the resulting normalization $A_{\mathrm{N}}$ is constrained as $13.8^{+5.8}_{-4.2}$, $13.2^{+3.4}_{-2.7}$, and $11.9^{+3.0}_{-1.9}$ by modeling $\xi_{\mathrm{cc}}$, $\xi_{\mathrm{cc}}+\xi_{\mathrm{cg}}$, and $\xi_{\mathrm{cc}} + \xi_{\mathrm{cg}} + \xi_{\mathrm{gg}}$, with average uncertainties at levels of $36\%$, $23\%$, and $21\%$, respectively. We find that the resulting $A_{\mathrm{N}}$ is statistically consistent with those independently obtained from weak-lensing magnification and from a joint analysis of shear and cluster abundance, with a preference for a lower value at a level of $\lesssim1.9\sigma$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

Green valley galaxies represent the population that is likely to transition from the star-forming to the quiescent phases. To investigate the role of the environment in quenching star formation, we use the wide-field data from the Hyper Suprime-Cam Strategic Subaru Proposal survey to quantify the frequency of green valley galaxies in different environments and their redshift evolution. We find that the green valley fraction, in general, is less than 20% in any redshift and environment. The green valley fraction, when normalized to the total population, is higher in the field than that in groups or clusters and decreases with a decreasing redshift and increasing mass. The lower fraction of transitional galaxies in denser environments could be a consequence of the lack of star-forming galaxies, which could be the progenitors of green valley galaxies. To assess the effect of the environment on star formation quenching, we define the effective green valley fraction as the ratio of the number of green valley galaxies to that of nonquiescent galaxies only. The effective green valley fraction for field galaxies is lower than that for group or cluster galaxies, which reveals a strong positive mass dependence and mild redshift evolution. Moreover, the specific star formation rate (sSFR) is reduced by 0.1-0.3 dex in groups or clusters. Our results thus imply that an ongoing slow quenching process has been acting in the dense environment since z~1.

We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 \texttt{CAMIRA} galaxy clusters with richness $N>15$ at redshift $0.2\leq z <1.1$ in the Subaru Hyper Suprime-Cam (HSC) survey. We use two distinct populations of color-selected, flux-limited background galaxies, namely the low-$z$ and high-$z$ samples at mean redshifts of $\approx1.1$ and $\approx1.4$, respectively, from which to measure the weak-lensing magnification signal by accounting for cluster contamination as well as masking effects. Our magnification bias measurements are found to be uncontaminated according to validation tests against the "null-test" samples for which the net magnification bias is expected to vanish. The magnification bias for the full \texttt{CAMIRA} sample is detected at a significance level of $9.51\sigma$, which is dominated by the high-$z$ background. We forward-model the observed magnification data to constrain the normalization of the richness-to-mass ($N$--$M$) relation for the \texttt{CAMIRA} sample with informative priors on other parameters. The resulting scaling relation is $N\propto {M_{500}}^{0.92\pm0.13} (1 + z)^{-0.48\pm0.69}$, with a characteristic richness of $N=\left(17.72\pm2.60\right)$ and intrinsic log-normal scatter of $0.15\pm0.07$ at $M_{500} = 10^{14}h^{-1}M_{\odot}$. With the derived $N$--$M$ relation, we provide magnification-calibrated mass estimates of individual \texttt{CAMIRA} clusters, with the typical uncertainty of $\approx39\%$ and $\approx32\%$ at richness$\approx20$ and $\approx40$, respectively. We further compare our magnification-inferred $N$--$M$ relation with those from the shear-based results in the literature, finding good agreement.

Recent constraints on the splashback radius around optically selected galaxy clusters from the redMaPPer cluster-finding algorithm in the literature have shown that the observed splashback radius is $\sim 20\%$ smaller than that predicted by N-body simulations. We present analyses on the splashback features around $\sim 3000$ optically selected galaxy clusters detected by the independent cluster-finding algorithm CAMIRA over a wide redshift range of $0.1<z_{\rm cl}<1.0$ from the second public data release of the Hyper Suprime-Cam (HSC) Subaru Strategic Program covering $\sim 427~{\rm deg}^2$ for the cluster catalog. We detect the splashback feature from the projected cross-correlation measurements between the clusters and photometric galaxies over the wide redshift range, including for high redshift clusters at $0.7<z_{\rm cl}<1.0$, thanks to deep HSC images. We find that constraints from red galaxy populations only are more precise than those without any color cut, leading to $1\sigma$ precisions of $\sim 15\%$ at $0.4<z_{\rm cl}<0.7$ and $0.7<z_{\rm cl}<1.0$. These constraints are more consistent with the model predictions ($\lesssim 1\sigma$) than their $20\%$ smaller values as suggested by the previous studies with the redMaPPer ($\sim 2\sigma$). We also investigate selection effects of the optical cluster-finding algorithms on the observed splashback features by creating mock galaxy catalogs from a halo occupation distribution model, and find that such effects to be sub-dominant for the CAMIRA cluster-finding algorithm. We also find that the redMaPPer-like cluster-finding algorithm induces a smaller inferred splashback radius in our mock catalog, especially at lower richness, which can well explain the smaller splashback radii in the literature. In contrast, these biases are significantly reduced when increasing its aperture size.

We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods.

We report the results of the Australia Telescope Compact Array (ATCA) 15 mm observation of the Phoenix galaxy cluster possessing an extreme star-burst brightest cluster galaxy (BCG) at the cluster center. We spatially resolved radio emission around the BCG, and found diffuse bipolar and bar-shape structures extending from the active galactic nucleus (AGN) of the BCG. They are likely radio jets/lobes, whose sizes are ~10-20 kpc and locations are aligned with X-ray cavities. If we assume that the radio jets/lobes expand with the sound velocity, their ages are estimated to be ~10 Myr. We also found compact radio emissions near the center and suggest that they are another young bipolar jets with ~1 Myr of age. Moreover, we found extended radio emission surrounding the AGN and discussed the possibility that the component is a product of the cooling flow, by considering synchrotron radiation partially absorbed by molecular clumps, free-free emission from the warm ionized gas, and the spinning dust emission from dusty circum-galactic medium.

Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

We report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We estimate the bias on weak lensing mass measurements of shear-selected galaxy cluster samples. The mass bias is expected to be significant because constructions of cluster samples from peaks in weak lensing mass maps and measurements of cluster masses from their tangential shear profiles share the same noise. We quantify this mass bias from large sets of mock cluster samples with analytical density profiles and realistic large-scale structure noise from ray-tracing simulations. We find that, even for peaks with signal-to-noise ratio larger than $4.0$ in weak lensing mass maps constructed in a deep survey with a high source galaxy number density of $30~\mathrm{arcmin}^{-2}$, derived weak lensing masses for these shear-selected clusters are still biased high by $\sim55\%$ on average. Such a large bias mainly originates from up-scattered low mass objects, which is an inevitable consequence of selecting clusters with a noisy observable directly linked to the mass measurement. We also investigate the dependence of the mass bias on different physical and observational parameters, finding that the mass bias strongly correlates with cluster redshifts, true halo masses, and selection signal-to-noise thresholds, but having moderate dependence on observed weak lensing masses and survey depths. This bias, albeit considerable, can still be modeled accurately in statistical studies of shear-selected clusters as the intrinsic scatter around the mean bias is found to be reasonable in size. We demonstrate that such a bias can explain the deviation in X-ray properties previously found on a shear-selected cluster sample. Our result will be useful for turning large samples of shear-selected clusters available in future surveys into potential probes of cosmology and cluster astrophysics.

Gravitational lensing sometimes dominates the observed properties of apparently very bright objects. We present morphological properties in the high-resolution (FWHM $\sim0.''15$) Atacama Large Millimeter/submillimeter Array (ALMA) 1-mm map for an ultra-luminous quasar (QSO) at $z=6.30$, SDSS J010013.02+280225.8 (hereafter J0100+2802), whose black hole mass $M_{\rm BH}$ is the most massive ($\sim$ 1.2$\times10^{10}M_{\odot}$) at $z>6$ ever known. We find that the continuum emission of J0100+2802 is resolved into a quadruple system within a radius of $0.''2$, which can be interpreted as either multiple dusty star-forming regions in the host galaxy or multiple images due to strong gravitational lensing. The Mg$_{\rm II}$ absorption and the potential Ly$\alpha$ line features have been identified at $z=2.33$ in the near-infrared spectroscopy towards J0100+2802, and a simple mass model fitting well reproduces the positions and flux densities of the quadruple system, both of which are consistent with the latter interpretation. Although a high-resolution map taken in the Advanced Camera for Survey (ACS) on board Hubble Space Telescope (HST) shows a morphology with an apparently single component, in our fiducial lens mass model it can simply be explained by a $\sim50$ pc scale offset between the ALMA and HST emission regions. In this case, the magnification factor for the observed HST emission is obtained to $\sim450$, reducing the intrinsic $M_{\rm BH}$ estimate to even below $10^{9}\,M_{\odot}$. The confirmation or the rejection of the gravitational lensing scenario is important for our understanding of the super-massive black holes in the early Universe.

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $\xi_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<\theta<56'$ for $\xi_+$ and $28'<\theta<178'$ for $\xi_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega_m=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints.

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16\AA to 50\AA and Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1\AA) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240\AA and CIV1549\AA with those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($\beta$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

A double source plane (DSP) system is a precious probe for the density profile of distant galaxies and cosmological parameters. However, these measurements could be affected by the surrounding environment of the lens galaxy. Thus, it is important to evaluate the cluster-scale mass for detailed mass modeling. We observed the {\it Eye of Horus}, a DSP system discovered by the Subaru HSC--SSP, with XMM--Newton. We detected two X-ray extended emissions, originating from two clusters, one centered at the {\it Eye of Horus}, and the other located $\sim100$ arcsec northeast to the {\it Eye of Horus}. We determined the dynamical mass assuming hydrostatic equilibrium, and evaluated their contributions to the lens mass interior of the Einstein radius. The contribution of the former cluster is $1.1^{+1.2}_{-0.5}\times10^{12}~M_{\odot}$, which is $21-76\%$ of the total mass within the Einstein radius. The discrepancy is likely due to the complex gravitational structure along the line of sight. On the other hand, the contribution of the latter cluster is only $\sim2\%$ on the {\it Eye of Horus}. Therefore, the influence associated with this cluster can be ignored.

It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. We use Chitah to hunt for new lens systems in the Hyper Suprime$-$Cam Subaru Strategic Program (HSC SSP) S16A. We present 46 lens candidates, of which 3 are previously known. Including 2 additional lenses found by YattaLens, we obtain X-shooter spectra of 6 promising candidates for lens confirmation and redshift measurements. We report new spectroscopic redshift measurements for both the lens and source galaxies in 4 lens systems. We apply the lens modeling software Glee to model our 6 X-shooter lenses uniformly. Through our analysis of the HSC images, we find that HSCJ022622$-$042522, HSCJ115252$+$004733, and HSCJ141136$-$010216 have point-like lensed images, and that the lens light distribution is well aligned with mass distribution within 6 deg. Thanks to the X-shooter spectra, we estimate fluxes on the Baldwin- Phillips-Terlevich (BPT) diagram, and find that HSCJ022622$-$042522 has a probable quasar source, based on the upper limit of the Nii flux intensity. We also measure the FWHM of Ly$\alpha$ emission of HSCJ141136$-$010216 to be $\sim$254 km/s, showing that it is a probable Lyman-$\alpha$ emitter.

We apply the Fourier Power Function Shapelets (FPFS) shear estimator to the first year data of the Hyper Suprime-Cam survey to construct a shape catalog. The FPFS shear estimator has been demonstrated to have multiplicative bias less than one percent in the absence of blending, regardless of complexities of galaxy shapes, smears of point spread functions (PSFs) and contamination from noise. The blending bias is calibrated with realistic image simulations, which include the impact of neighboring objects, using the COSMOS Hubble Space Telescope images. Here we carefully test the influence of PSF model residual on the FPFS shear estimation and the uncertainties in the shear calibration. Internal null tests are conducted to characterize potential systematics in the FPFS shape catalog and the results are compared with those measured using a catalog where the shapes were estimated using the re-Gaussianization algorithms. Furthermore, we compare various weak lensing measurements between the FPFS shape catalog and the re-Gaussianization shape catalog and conclude that the weak lensing measurements between these two shape catalogs are consistent with each other within the statistical uncertainty.

Recent rapid progress in time domain surveys makes it possible to detect various types of explosive transients in the Universe in large numbers, some of which will be gravitationally lensed into multiple images. Although a large number of strongly lensed distant galaxies and quasars have already been discovered, strong lensing of explosive transients opens up new applications, including improved measurements of cosmological parameters, powerful probes of small scale structure of the Universe, and new observational tests of dark matter scenarios, thanks to their rapidly evolving light curves as well as their compact sizes. In particular, compact sizes of emitting regions of these transient events indicate that wave optics effects play an important role in some cases, which can lead to totally new applications of these lensing events. Recently we have witnessed first discoveries of strongly lensed supernovae, and strong lensing events of other types of explosive transients such as gamma-ray bursts, fast radio bursts, and gravitational waves from compact binary mergers are expected to be observed soon. In this review article, we summarize the current state of research on strong gravitational lensing of explosive transients and discuss future prospects.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.

We present the first results of a pilot X-ray study of 37 rich galaxy clusters at $0.1<z<1.1$ in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) field. Diffuse X-ray emissions from these clusters were serendipitously detected in the XMM-Newton fields of view. We systematically analyze X-ray images of 37 clusters and emission spectra of a subsample of 17 clusters with high photon statistics by using the XMM-Newton archive data. The frequency distribution of the offset between the X-ray centroid or peak and the position of the brightest cluster galaxy was derived for the optical cluster sample. The fraction of relaxed clusters estimated from the X-ray peak offsets in 17 clusters is $29\pm11(\pm13)$\%, which is smaller than that of the X-ray cluster samples such as HIFLUGCS. Since the optical cluster search is immune to the physical state of X-ray-emitting gas, it is likely to cover a larger range of the cluster morphology. We also derived the luminosity-temperature relation and found that the slope is marginally shallower than those of X-ray-selected samples and consistent with the self-similar model prediction of 2. Accordingly, our results show that the X-ray properties of the optical clusters are marginally different from those observed in the X-ray samples. The implication of the results and future prospects are briefly discussed.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

Context: The determination of the stellar initial mass function (IMF) of massive galaxies is one of the open problems in cosmology. Strong gravitational lensing is one of the few methods that allow us to constrain the IMF outside of the Local Group. Aims: The goal of this study is to statistically constrain the distribution in the IMF mismatch parameter, defined as the ratio between the true stellar mass of a galaxy and that inferred assuming a reference IMF, of massive galaxies from the BOSS CMASS sample. Methods: We take 23 strong lenses drawn from the CMASS sample, measure their Einstein radii and stellar masses using multi-band photometry from the Hyper Suprime-Cam survey, then fit a model distribution for the IMF mismatch parameter and dark matter halo mass to the whole sample. We use a prior on halo mass from weak lensing measurements and account for strong lensing selection effects in our model. Results: Assuming an NFW density profile for the dark matter distribution, we infer a value $\mu_{\mathrm{IMF}} = -0.04\pm0.11$ for the average base-10 logarithm of the IMF mismatch parameter, defined with respect to a Chabrier IMF. A Salpeter IMF is in tension with our measurements. Conclusions: Our results are consistent with a scenario in which the region of massive galaxies where the IMF normalization is significantly heavier than that of the Milky Way is much smaller than the scales $5\sim10$~kpc probed by the Einstein radius of the lenses in our sample, as recent spatially resolved studies of the IMF in massive galaxies suggest. The MCMC chains describing the posterior probability distribution of the model are available online, together with the code used to obtain them.

Constraining the relation between the richness $N$ and the halo mass $M$ over a wide redshift range for optically-selected clusters is a key ingredient for cluster-related science in optical surveys, including the Subaru Hyper Suprime-Cam (HSC) survey. We measure stacked weak lensing profiles around 1747 HSC CAMIRA clusters over a redshift range of $0.1\leq z_{\rm cl}\leq 1.0$ with $N\geq 15$ using the HSC first-year shear catalog covering $\sim$$140$ ${\rm deg^2}$. The exquisite depth and image quality of the HSC survey allow us to measure lensing signals around the high-redshift clusters at $0.7\leq z_{\rm cl}\leq 1.0$ with a signal-to-noise ratio of 19 in the comoving radius range $0.5\lesssim R\lesssim 15 h^{-1}{\rm Mpc}$. We constrain richness-mass relations $P(\ln N|M,z)$ of the HSC CAMIRA clusters assuming a log-normal distribution without informative priors on model parameters, by jointly fitting to the lensing profiles and abundance measurements under both Planck and WMAP cosmological models. We show that our model gives acceptable $p$-values when we add redshift dependent terms which are proportional to $\ln (1+z)$ and $[\ln (1+z)]^{2}$ into the mean and scatter relations of $P(\ln N|M,z)$. Such terms presumably originate from the variation of photometric redshift errors as a function of the redshift. We show that the constraints on the mean relation $\langle M|N \rangle$ are consistent between the Planck and WMAP models, whereas the scatter values $\sigma_{\ln M|N}$ for the Planck model are systematically larger than those for the WMAP model. We also show that the scatter values for the Planck model increase toward lower richness values, whereas those for the WMAP model are consistent with constant values as a function of richness. This result highlights the importance of the scatter in the mass-richness relation for cluster cosmology.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$.

Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the Reionization Lensing Cluster Survey (RELICS). We find that for 6/8 \textit{HST}-selected $z\gtrsim8$ sources, the $z\gtrsim8$ solutions are still strongly preferred over $z\sim$1-2 solutions after the inclusion of \textit{Spitzer} fluxes, and two prefer a $z\sim 7$ solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses ($5\times10^6 M_{\odot}$ -- $4\times10^9$ $M_{\odot}$), star formation rates (0.2-14 $M_{\odot}\rm yr^{-1}$), and ages (30-600 Myr) among our sample. Of particular interest is Abell1763-1434, which shows evidence of an evolved stellar population at $z\sim8$, implying its first generation of star formation occurred just $< 100$ Myr after the Big Bang. SPT0615-JD, a spatially resolved $z\sim10$ candidate, remains at its high redshift, supported by deep \textit{Spitzer}/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these $z \gtrsim 8$ candidates (H = 26.1-27.8 AB mag), only the \textit{James Webb Space Telescope} will be able further confirm the presence of evolved stellar populations early in the universe.

We present very faint dropout galaxies at z~6-9 with a stellar mass M* down to M*~10^6Mo that are found in deep optical/near-infrared (NIR) images of the full data sets of the Hubble Frontier Fields (HFF) program in conjunction with deep ground-based and Spitzer images and gravitational lensing magnification effects. We investigate stellar populations of the HFF dropout galaxies with the optical/NIR photometry and BEAGLE models made of self-consistent stellar population synthesis and photoionization models, carefully including strong nebular emission impacting on the photometry. We identify 357 galaxies with M*~10^6-10^9Mo, and find that a stellar mass to UV luminosity L_UV ratio M*/L_UV is nearly constant at M*~10^6-10^9Mo. Our best-estimate M*/L_UV function is comparable to a model of star-formation duration time of 100 Myr, but 2-7 times higher than the one of 10 Myr assumed in a previous study (at the 5sigma level) that would probably underestimate M* of faint galaxies. We derive the galaxy stellar mass functions (GSMFs) at z~6-9 that agree with those obtained by previous studies with no M*/L_UV assumptions at M*>~10^8Mo, and that extends to M*~10^6Mo. Estimating the stellar mass densities rho* with the GSMFs, we find that rho* smoothly increases from log(rho*/[Mo Mpc^(-3)])=5.91 +0.75/-0.65 at z~9 to 6.21 +0.39/-0.37 at z~6-7, which is consistent with the one estimated from star-formation rate density measurements. In conjunction with the estimates of the galaxy effective radii R_e on the source plane, we have pinpointed two objects with low stellar masses (M*<=10^7Mo) and very compact morphologies (R_e<=40 physical pc) that are comparable with those of globular clusters (GCs) in the Milky Way today. These objects are candidates of star clusters that should be a part or a dominant component of high-redshift low-mass galaxy, some of which may be related to GCs today.

We present multi-wavelength studies of optically-defined merging clusters, based on the Hyper Suprime-Cam Subaru Strategic Program. Luminous red galaxies, tracing cluster mass distributions, enable to identify cluster subhalos at various merging stages, and thus make a homogeneous sample of cluster mergers, which is unbiased with respect to the merger boost of the intracluster medium (ICM). We define, using a peak-finding method, merging clusters with multiple-peaks and single clusters with single-peaks from the CAMIRA cluster catalog. Stacked weak-lensing analysis indicates that our sample of the merging clusters is categorized into major mergers. The average halo concentration for the merging clusters is $\sim70\%$ smaller than that of the single-peak clusters, which agrees well with predictions of numerical simulations. The spatial distribution of subhalos is less centrally concentrated than the mass distribution of the main halo. The fractions of red galaxies in the merging clusters are not higher than those of the single-peak clusters. We find a signature of the merger boost of the ICM from stacked Planck Sunyaev-Zeldovich effect and ROSAT X-ray luminosity, but not in optical richness. The stacked X-ray surface brightness distribution, aligned with the main-subhalo pairs of low redshift and massive clusters, shows that the central gas core is elongated along the merger axis and overall gas distribution is misaligned by $\sim60$ deg. The homogeneous, unbiased sample of cluster mergers and multi-wavelength follow-up studies provide a unique opportunity to make a complete picture of merger physics over the whole process.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

ULTIMATE-Subaru (Ultra-wide Laser Tomographic Imager and MOS with AO for Transcendent Exploration on Subaru) and WFIRST (Wide Field Infra-Red Survey Telescope) are the next generation near-infrared instruments that have a large field-of-view. They allow us to conduct deep and wide transient surveys in near-infrared. Such a near-infrared transient survey enables us to find very distant supernovae that are redshifted to the near-infrared wavelengths. We have performed the mock transient surveys with ULTIMATE-Subaru and WFIRST to investigate their ability to discover Population III pair-instability supernovae. We found that a 5-year 1 deg2 K-band transient survey with the point-source limiting magnitude of 26.5 mag with ULTIMATE-Subaru may find about 2 Population III pair-instability supernovae beyond the redshift of 6. A 5-year 10 deg2 survey with WFIRST reaching 26.5 mag in the F184 band may find about 7 Population III pair-instability supernovae beyond the redshift of 6. We also find that the expected numbers of the Population III pair-instability supernova detections increase about a factor of 2 if the near-infrared transient surveys are performed towards clusters of galaxies. Other supernovae such as Population II pair-instability supernovae would also be detected in the same survey. This study demonstrates that the future wide-field near-infrared instruments allow us to investigate the explosions of the first generation supernovae by performing the deep and wide near-infrared transient surveys.

Superluminous supernovae have been proposed to arise from Population III progenitors that explode as pair-instability supernovae. Pop III stars are the first generation of stars in the Universe, and are thought to form as late as $z \sim 6$. Future near-infrared imaging facilities such as ULTIMATE-Subaru can potentially detect and identify these PISNe with a dedicated survey. Gravitational lensing by intervening structure in the Universe can aid in the detection of these rare objects by magnifying the high-$z$ source population into detectability. We perform a mock survey with ULTIMATE-Subaru, taking into account lensing by line-of-sight structure to evaluate its impact on the predicted detection rate. We compare a LOS mass reconstruction using observational data from the Hyper Suprime Cam survey to results from cosmological simulations to test their consistency in calculating the magnification distribution in the Universe to high-$z$, but find that the data-based method is still limited by an inability to accurately characterize structure beyond $z \sim1.2$. We also evaluate a survey strategy of targeting massive galaxy clusters to take advantage of their large areas of high magnification. We find that targeting clusters can result in a gain of a factor of $\sim$two in the predicted number of detected PISNe at $z > 5$, and even higher gains with increasing redshift, given our assumed survey parameters. For the highest-redshift sources at $z \sim 7-9$, blank field surveys will not detect any sources, and lensing magnification by massive clusters will be necessary to observe this population.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

The cross-correlation study of the unresolved $\gamma$-ray background (UGRB) with galaxy clusters has a potential to reveal the nature of the UGRB. In this paper, we perform a cross-correlation analysis between $\gamma$-ray data by the Fermi Large Area Telescope (Fermi-LAT) and a galaxy cluster catalogue from the Subaru Hyper Suprime-Cam (HSC) survey. The Subaru HSC cluster catalogue provides a wide and homogeneous large-scale structure distribution out to the high redshift at $z=1.1$, which has not been accessible in previous cross-correlation studies. We conduct the cross-correlation analysis not only for clusters in the all redshift range ($0.1 < z < 1.1$) of the survey, but also for subsamples of clusters divided into redshift bins, the low redshift bin ($0.1 < z < 0.6$) and the high redshift bin ($0.6 < z < 1.1$), to utilize the wide redshift coverage of the cluster catalogue. We find the evidence of the cross-correlation signals with the significance of 2.0-2.3$\sigma$ for all redshift and low-redshift cluster samples. On the other hand, for high-redshift clusters, we find the signal with weaker significance level (1.6-1.9$\sigma$). We also compare the observed cross-correlation functions with predictions of a theoretical model in which the UGRB originates from $\gamma$-ray emitters such as blazars, star-forming galaxies and radio galaxies. We find that the detected signal is consistent with the model prediction.

Effect of gravitational magnification on the measurement of distance modulus of type Ia supernovae is presented. We investigate a correlation between magnification and Hubble residual to explore how the magnification affects the estimation of cosmological parameters. We estimate magnification of type Ia supernovae in two distinct methods: one is based on convergence mass reconstruction under the weak lensing limit and the other is based on the direct measurement from galaxies distribution. Both magnification measurements are measured from Subaru Hyper Suprime-Cam survey catalogue. For both measurements, we find no significant correlation between Hubble residual and magnification. Furthermore, we correct for the apparent supernovae fluxes obtained by Supernova Legacy Survey 3-year sample using direct measurement of the magnification. We find $\Omega_{\rm m0} = 0.287 ^{+0.104} _{-0.085}$ and $w = -1.161 ^{+0.595}_{-0.358}$ for supernovae samples corrected for lensing magnification when we use photometric redshift catalogue of Mizuki, while $\Omega_{\rm m0} = 0.253 ^{+0.113} _{-0.087}$ and $w = -1.078 ^{+0.498} _{-0.297}$ for DEmP photo-z catalogue. Therefore, we conclude that the effect of magnification on the supernova cosmology is negligibly small for the current surveys; however, it has to be considered for the future supernova survey like LSST.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

Using the two large cosmological hydrodynamical simulations, Horizon-AGN (H-AGN) and Horizon-noAGN (H-noAGN, no AGN feedback), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the total density profiles (dark matter + stars) at the effective radius of massive early-type galaxies (M*>10^11 Msun). We have studied the dependencies of the mass-weighted density slope gamma'_tot with the effective radius, the galaxy mass and the host halo mass at z~0.3 and found that the inclusion of AGN feedbackalways leads to a much better agreement with observational values and trends. Our analysis suggests also that the inclusion of AGN feedback favours a strong correlation between gamma'_tot and the density slope of the dark matter component while, in the absence of AGN activity, gamma'_tot is rather strongly correlated with the density slope of the stellar component. Finally, we find that gamma'_tot derived from our samples of galaxies increases from z=2 to z=0,in good agreement with the expected observational trend. The derived slopes are slightly lower than in the data when AGN is included because the simulated galaxies tend to be too extended, especially the least massive ones. However, the simulated compact galaxies without AGN feedback have gamma'_tot values that are significantly too high compared to observations.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

Primordial black holes (PBHs) have long been suggested as a viable candidate for the elusive dark matter (DM). The abundance of such PBHs has been constrained using a number of astrophysical observations, except for a hitherto unexplored mass window of $M_{\rm PBH}=[10^{-14},10^{-9}]M_\odot$. Here we carry out a dense-cadence (2~min sampling rate), 7 hour-long observation of the Andromeda galaxy (M31) with the Subaru Hyper Suprime-Cam to search for microlensing of stars in M31 by PBHs lying in the halo regions of the Milky Way (MW) and M31. Given our simultaneous monitoring of tens of millions of stars in M31, if such light PBHs make up a significant fraction of DM, we expect to find many microlensing events for the PBH DM scenario. However, we identify only a single candidate event, which translates into the most stringent upper bounds on the abundance of PBHs in the mass range $M_{\rm PBH}\simeq [10^{-11}, 10^{-6}]M_\odot$.

We reinterpret the shear estimator developed by Zhang & Komatsu (2011) within the framework of Shapelets and propose the Fourier Power Function Shapelets (FPFS) shear estimator. Four shapelet modes are calculated from the power function of every galaxy's Fourier transform after deconvolving the Point Spread Function (PSF) in Fourier space. We propose a novel normalization scheme to construct dimensionless ellipticity and its corresponding shear responsivity using these shapelet modes. Shear is measured in a conventional way by averaging the ellipticities and responsivities over a large ensemble of galaxies. With the introduction and tuning of a weighting parameter, noise bias is reduced below one percent of the shear signal. We also provide an iterative method to reduce selection bias. The FPFS estimator is developed without any assumption on galaxy morphology, nor any approximation for PSF correction. Moreover, our method does not rely on heavy image manipulations nor complicated statistical procedures. We test the FPFS shear estimator using several HSC-like image simulations and the main results are listed as follows. (i) For simulations which only contain isolated galaxies, the amplitude of multiplicative bias is below 1 % . (ii) For more realistic simulations which also contain blended galaxies, the blended galaxies are deblended by the first generation HSC deblender before shear measurement. Multiplicative bias of (-5.71 +- 0.31) % is found. The blending bias is calibrated by image simulations. Finally, we test the consistency and stability of this calibration.

Velocity dispersions have been employed as a method to measure masses of clusters. To complement this conventional method, we explore the possibility of constraining cluster masses from the stacked phase space distribution of galaxies at larger radii, where infall velocities are expected to have a sensitivity to cluster masses. First, we construct a two component model of the three-dimensional phase space distribution of haloes surrounding clusters up to 50 $h^{-1}$Mpc from cluster centres based on $N$-body simulations. We find that the three-dimensional phase space distribution shows a clear cluster mass dependence up to the largest scale examined. We then calculate the probability distribution function of pairwise line-of-sight velocities between clusters and haloes by projecting the three-dimensional phase space distribution along the line-of-sight with the effect of the Hubble flow. We find that this projected phase space distribution, which can directly be compared with observations, shows a complex mass dependence due to the interplay between infall velocities and the Hubble flow. Using this model, we estimate the accuracy of dynamical mass measurements from the projected phase space distribution at the transverse distance from cluster centres larger than $2h^{-1}$Mpc. We estimate that, by using $1.5\times 10^5$ spectroscopic galaxies, we can constrain the mean cluster masses with an accuracy of 14.5\% if we fully take account of the systematic error coming from the inaccuracy of our model. This can be improved down to 5.7\% by improving the accuracy of the model.

We investigate the local and line-of-sight overdensities of strong gravitational lens galaxies using wide-area multiband imaging from the Hyper Suprime-Cam Subaru Strategic Program. We present 41 new definite or probable lens candidates discovered in Data Release 2 of the survey. Using a combined sample of 87 galaxy-scale lenses out to a lens redshift of $z_{\mathrm{L}} \sim 0.8$, we compare galaxy number counts in lines of sight toward known and newly-discovered lenses in the survey to those of a control sample consisting of random lines of sight. We also compare the local overdensity of lens galaxies to a sample of "twin" galaxies that have a similar redshift and velocity dispersion to test whether lenses lie in different environments from similar non-lens galaxies. We find that lens fields contain higher number counts of galaxies compared to the control fields, but this effect arises from the local environment of the lens. Once galaxies in the lens plane are removed, the lens lines of sight are consistent with the control sample. The local environments of the lenses are overdense compared to the control sample, and are slightly overdense compared to those of the twin sample, although the significance is marginal. There is no significant evidence of the evolution of the local overdensity of lens environments with redshift.

RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5{\sigma}. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL).

The detailed observation of the distribution of redshifts and chirp masses of binary black hole mergers is expected to provide a clue to their origin. In this paper, we develop a hybrid model of the probability distribution function of gravitational lensing magnification taking account of both strong and weak gravitational lensing, and use it to study the effect of gravitational lensing magnification on the distribution of gravitational waves from distant binary black hole mergers detected in ongoing and future gravitational wave observations. We find that the effect of gravitational lensing magnification is significant at high ends of observed chirp mass and redshift distributions. While a high mass tail in the observed chirp mass distribution is produced by highly magnified gravitational lensing events, we find that highly demagnified images of strong lensing events produce a high redshift ($z_{\rm obs}> 15$) tail in the observed redshift distribution, which can easily be observed in the third-generation gravitational wave observatories. Such a demagnified, apparently high redshift event is expected to be accompanied by a magnified image that is observed typically $10-100$ days before the demagnified image. For highly magnified events that produce apparently very high chirp masses, we expect pairs of events with similar magnifications with time delays typically less than a day. This work suggests the critical importance of gravitational lensing (de-)magnification on the interpretation of apparently very high mass or redshift gravitational wave events.

While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

We study the dependence of surface mass density profiles, which can be directly measured by weak gravitational lensing, on the orientation of haloes with respect to the line-of-sight direction, using a suite of $N$-body simulations. We find that, when major axes of haloes are aligned with the line-of-sight direction, surface mass density profiles have higher amplitudes than those averaged over all halo orientations, over all scales from $0.1$ to $100\,\mathrm{Mpc}/h$ we studied. While the orientation dependence at small scales is ascribed to the halo triaxiality, our results indicate even stronger orientation dependence in the so-called two-halo regime, up to $100\,\mathrm{Mpc}/h$. The orientation dependence for the two-halo term is well approximated by a multiplicative shift of the amplitude and therefore a shift in the halo bias parameter value. The halo bias from the two-halo term can be overestimated or underestimated by up to $\sim 30 \%$ depending on the viewing angle, which translates into the bias in estimated halo masses by up to a factor of two from halo bias measurements. The orientation dependence at large scales originates from the anisotropic halo-matter correlation function, which has an elliptical shape with the axis ratio of $\sim 0.55$ up to $100\, \mathrm{Mpc}/h$. We discuss potential impacts of halo orientation bias on other observables such as optically selected cluster samples and a clustering analysis of large-scale structure tracers such as quasars.

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

A galaxy cluster acts as a cosmic telescope over background galaxies but also as a cosmic microscope of the lens imperfections. The diverging magnification of lensing caustics enhances the microlensing effect of substructure present within the lensing mass. Fine-scale structure can be accessed as a moving background source brightens and disappears when crossing these caustics. The recent recognition of a distant lensed star near the Einstein radius of the galaxy cluster MACSJ1149.5+2223 (Kelly et al. 2017) allows the rare opportunity to reach subsolar mass microlensing through a super-critical column of cluster matter. Here we compare these observations with high-resolution ray-tracing simulations that include stellar microlensing set by the observed intracluster starlight and also primordial black holes that may be responsible for the recently observed LIGO events. We explore different scenarios with microlenses from the intracluster medium and black holes, including primordial ones, and examine strategies to exploit these unique alignments. We find that the best constraints on the fraction of compact dark matter in the small-mass regime can be obtained in regions of the cluster where the intracluster medium plays a negligible role. This new lensing phenomenon should be widespread and can be detected within modest-redshift lensed galaxies so that the luminosity distance is not prohibitive for detecting individual magnified stars. Continuous {\it Hubble Space Telescope} monitoring of several such optimal arcs will be rewarded by an unprecedented mass spectrum of compact objects that can contribute to uncovering the nature of dark matter.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

We construct $z\sim6-7$, 8, and 9 faint Lyman break galaxy samples (334, 61, and 37 galaxies, respectively) with accurate size measurements with the software $\texttt{glafic}$ from the complete Hubble Frontier Fields (FF) cluster and parallel fields data. These are the largest samples hitherto and reach down to the faint ends of recently obtained deep luminosity functions. At faint magnitudes, however, these samples are highly incomplete for galaxies with large sizes, implying that derivation of the luminosity function (LF) sensitively depends on the intrinsic size--luminosity (RL) relation. We thus conduct simultaneous maximum-likelihood estimation of LF and RL relation parameters from the observed distribution of galaxies on the RL plane with help of a completeness map as a function of size and luminosity. At $z\sim6-7$, we find that the intrinsic RL relation expressed as $r_\textrm{e} \propto L^\beta$ has a notably steeper slope of $\beta=0.46^{+0.08}_{-0.09}$ than those at lower redshifts, which in turn implies that the LF has a relatively shallow faint-end slope of $\alpha=-1.86^{+0.17}_{-0.18}$. This steep $\beta$ can be reproduced by a simple analytical model in which smaller galaxies have lower specific angular momenta. The $\beta$ and $\alpha$ values for the $z\sim8$ and 9 samples are consistent with those for $z\sim6-7$ but with larger errors. For all three samples there is a large, positive covariance between $\beta$ and $\alpha$, implying that the simultaneous determination of these two parameters is important. We also provide new strong lens mass models of Abell S1063 and Abell 370 as well as updated mass models of Abell 2744 and MACS J0416.1$-$2403.

We examine the capability of pulsar timing arrays (PTAs) to detect very small-scale clumps of dark matter (DM), which are a natural outcome of the standard cold dark matter (CDM) paradigm. A clump streaming near the Earth or a pulsar induces an impulsive acceleration to encode residuals on pulsar timing data. We show that, assuming the standard abundance of DM clumps predicted by the CDM model, small-scale DM clumps with masses from $\sim 10^{-11} M_\odot$ to $\sim 10^{-8} \ M_\odot$ can be detectable by a PTA observation for a few decades with ${\cal O}(100)$ of pulsars with a timing noise of ${\cal O}(10)$ ns located at $\gtrsim 3$ kpc away from the Galactic center, as long as these mass scales are larger than the cutoff scale of the halo mass function that is determined by the particle nature of DM. Our result suggests that PTAs can provide a unique opportunity for testing one of the most fundamental predictions of the CDM paradigm. In addition, the detections and non-detections can constrain the cutoff mass scale inherent to the DM model.

We present UV luminosity functions of dropout galaxies at $z\sim6-10$ with the complete Hubble Frontier Fields data. We obtain a catalog of $\sim450$ dropout-galaxy candidates (350, 66, and 40 at $z\sim6-7$, 8, and 9, respectively), whose UV absolute magnitudes reach $\sim-14$ mag, $\sim2$ mag deeper than the Hubble Ultra Deep Field detection limits. We carefully evaluate number densities of the dropout galaxies by Monte-Carlo simulations, including all lensing effects such as magnification, distortion, and multiplication of images as well as detection completeness and contamination effects in a self-consistent manner. We find that UV luminosity functions at $z\sim6-8$ have steep faint-end slopes, $\alpha\sim-2$, and likely steeper slopes, $\alpha\lesssim-2$ at $z\sim9-10$. We also find that the evolution of UV luminosity densities shows a non-accelerated decline beyond $z\sim8$ in the case of $M_\mathrm{trunc}=-15$, while an accelerated in the case of $M_\mathrm{trunc}=-17$. We examine whether our results are consistent with the Thomson scattering optical depth from the Planck satellite and the ionized hydrogen fraction $Q_\mathrm{HII}$ at $z\lesssim7$ based on the standard analytic reionization model. We find that there exist reionizaiton scenarios that consistently explain all the observational measurements with the allowed parameters of $f_{\rm esc}=0.17^{+0.07}_{-0.03}$ and $M_\mathrm{trunc}>-14.0$ for $\log\xi_\mathrm{ion}/[\mathrm{erg}^{-1}\ \mathrm{Hz}]=25.34$, where $f_{\rm esc}$ is the escape fraction, $M_\mathrm{trunc}$ is the faint limit of the UV luminosity function, and $\xi_\mathrm{ion}$ is the conversion factor of the UV luminosity to the ionizing photon emission rate. The length of the reionization period is estimated to be $\Delta z=3.9^{+2.0}_{-1.6}$ (for $0.1<Q_\mathrm{HII}<0.99$), consistent with the recent estimate from Planck.

We performed spectroscopic observations of the small-separation lensed quasar SDSS J1001+5027, whose images have an angular separation $\theta \sim 2.^{\!\!\prime\prime}86$, and placed constraints on the physical properties of gas clouds in the vicinity of the quasar (i.e., in the outflowing wind launched from the accretion disk). The two cylinders of sight to the two lensed images go through the same region of the outflowing wind and they become fully separated with no overlap at a very large distance from the source ($\sim 330$ pc). We discovered a clear difference in the profile of the CIV broad absorption line (BAL) detected in the two lensed images in two observing epochs. Because the kinematic components in the BAL profile do not vary in concert, the observed variations cannot be reproduced by a simple change of ionization state. If the variability is due to gas motion around the background source (i.e., the continuum source), the corresponding rotational velocity is $v_{rot}\geq 18,000$ km/s, and their distance from the source is $r\leq 0.06$ pc assuming Keplerian motion. Among three MgII and three CIV NAL systems that we detected in the spectra, only the MgII system at $z_{abs} = 0.8716$ shows a hint of variability in its MgI profile on a rest-frame time scale of $\Delta t_{rest}$ $\leq 191$ days and an obvious velocity shear between the sightlines whose physical separation is $\sim 7$ kpc. We interpret this as the result of motion of a cosmologically intervening absorber, perhaps located in a foreground galaxy.

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

The recent discovery of fast transient events near critical curves of massive galaxy clusters, which are interpreted as highly magnified individual stars in giant arcs due to caustic crossing, opens up the possibility of using such microlensing events to constrain a range of dark matter models such as primordial black holes and scalar field dark matter. Based on a simple analytic model, we study lensing properties of a point mass lens embedded in a high magnification region, and derive the dependence of the peak brightness, microlensing time scales, and event rates on the mass of the point mass lens as well as the radius of a source star that is magnified. We find that the lens mass and source radius of the first event MACS J1149 Lensed Star 1 (LS1) are constrained, with the lens mass range of $0.1~M_\odot \lesssim M \lesssim 4\times 10^3M_\odot$ and the source radius range of $40~R_\odot \lesssim R \lesssim 260~R_\odot$. In the most plausible case with $M\approx 0.3~M_\odot$ and $R\approx 180~R_\odot$, the source star should have been magnified by a factor of $\approx 4300$ at the peak. The derived lens properties are fully consistent with the interpretation that MACS J1149 LS1 is a microlensing event produced by a star that contributes to the intra-cluster light. We argue that compact dark matter models with high fractional mass densities for the mass range $10^{-5}M_\odot \lesssim M\lesssim 10^2M_\odot$ are inconsistent with the observation of MACS J1149 LS1 because such models predict too low magnifications. Our work demonstrates a potential use of caustic crossing events in giant arcs to constrain compact dark matter.

We present a catalog of 182 galaxy clusters detected through the Sunyaev-Zel'dovich effect by the Atacama Cosmology Telescope in a contiguous 987.5 deg$^{2}$ field. The clusters were detected as SZ decrements by applying a matched filter to 148 GHz maps that combine the original ACT equatorial survey with data from the first two observing seasons using the ACTPol receiver. Optical/IR confirmation and redshift measurements come from a combination of large public surveys and our own follow-up observations. Where necessary, we measured photometric redshifts for clusters using a pipeline that achieves accuracy $\Delta z/(1 + z)=0.015$ when tested on SDSS data. Under the assumption that clusters can be described by the so-called Universal Pressure Profile and its associated mass-scaling law, the full signal-to-noise > 4 sample spans the mass range $1.6 < M^{\rm UPP}_{\rm 500c}/10^{14}{\rm M}_{\odot}<9.1$, with median $M^{\rm UPP}_{\rm 500c}=3.1 \times 10^{14}$ M$_{\odot}$. The sample covers the redshift range $0.1 < z < 1.4$ (median $z = 0.49$) and 28 clusters are new discoveries (median $z = 0.80$). We compare our catalog with other overlapping cluster samples selected using the SZ, optical,and X-ray wavelengths. We find the ratio of the UPP-based SZ mass to richness-based weak-lensing mass is $\langle M^{\rm UPP}_{\rm 500c} \rangle / \langle M^{\rm \lambda WL}_{\rm 500c} \rangle = 0.68 \pm 0.11$. After applying this calibration, the mass distribution for clusters with $M_{\rm 500c} > 4 \times 10^{14}$ M$_{\odot}$ is consistent with the number of such clusters found in the South Pole Telescope SZ survey.

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

The Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) is an excellent survey for the search for strong lenses, thanks to its area, image quality and depth. We use three different methods to look for lenses among 43,000 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample with photometry from the S16A internal data release of the HSC SSP. The first method is a newly developed algorithm, named YATTALENS, which looks for arc-like features around massive galaxies and then estimates the likelihood of an object being a lens by performing a lens model fit. The second method, CHITAH, is a modeling-based algorithm originally developed to look for lensed quasars. The third method makes use of spectroscopic data to look for emission lines from objects at a different redshift from that of the main galaxy. We find 15 definite lenses, 36 highly probable lenses and 282 possible lenses. Among the three methods, YATTALENS, which was developed specifically for this problem, performs best in terms of both completeness and purity. Nevertheless five highly probable lenses were missed by YATTALENS but found by the other two methods, indicating that the three methods are highly complementary. Based on these numbers we expect to find $\sim$300 definite or probable lenses by the end of the HSC SSP.

We study the internal structure of the Circum-Galactic Medium (CGM), using 29 spectra of 13 gravitationally lensed quasars with image separation angles of a few arcseconds, which correspond to 100 pc to 10 kpc in physical distances. After separating metal absorption lines detected in the spectra into high-ions with ionization parameter (IP) $>$ 40 eV and low-ions with IP $<$ 20 eV, we find that i) the fraction of absorption lines that are detected in only one of the lensed images is larger for low-ions ($\sim$16%) than high-ions ($\sim$2%), ii) the fractional difference of equivalent widths ($EW$s) between the lensed images is almost same (${\rm d}EW$ $\sim$ 0.2) for both groups although the low-ions have a slightly larger variation, and iii) weak low-ion absorbers tend to have larger ${\rm d}EW$ compared to weak high-ion absorbers. We construct simple models to reproduce these observed properties and investigate the distribution of physical quantities such as size and location of absorbers, using some free parameters. Our best models for absorbers with high-ions and low-ions suggest that i) an overall size of the CGM is at least $\sim$ 500 kpc, ii) a size of spherical clumpy cloud is $\sim$ 1 kpc or smaller, and iii) only high-ion absorbers can have diffusely distributed homogeneous component throughout the CGM. We infer that a high ionization absorber distributes almost homogeneously with a small-scale internal fluctuation, while a low ionization absorber consists of a large number of small-scale clouds in the diffusely distributed higher ionized region. This is the first result to investigate the internal small-scale structure of the CGM, based on the large number of gravitationally lensed quasar spectra.

The unprecedented depth and area surveyed by the Subaru Strategic Program with the Hyper Suprime-Cam (HSC-SSP) have enabled us to construct and publish the largest distant cluster sample out to $z\sim 1$ to date. In this exploratory study of cluster galaxy evolution from $z=1$ to $z=0.3$, we investigate the stellar mass assembly history of brightest cluster galaxies (BCGs), evolution of stellar mass and luminosity distributions, stellar mass surface density profile, as well as the population of radio galaxies. Our analysis is the first high redshift application of the top N richest cluster selection, which is shown to allow us to trace the cluster galaxy evolution faithfully. Over the 230 deg$^2$ area of the current HSC-SSP footprint, selecting the top 100 clusters in each of the 4 redshift bins allows us to observe the buildup of galaxy population in descendants of clusters whose $z\approx 1$ mass is about $2\times 10^{14}\,M_\odot$. Our stellar mass is derived from a machine-learning algorithm, which is found to be unbiased and accurate with respect to the COSMOS data. We find very mild stellar mass growth in BCGs (about 35% between $z=1$ and 0.3), and no evidence for evolution in both the total stellar mass-cluster mass correlation and the shape of the stellar mass surface density profile. We also present the first measurement of the radio luminosity distribution in clusters out to $z\sim 1$, and show hints of changes in the dominant accretion mode powering the cluster radio galaxies at $z\sim 0.8$.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

We present a morphological study of the 17 lensed Lyman-$\alpha$ emitter (LAE) galaxies of the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey (BELLS) for the GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) sample. This analysis combines the magnification effect of strong galaxy-galaxy lensing with the high resolution of the Hubble Space Telescope ($HST$) to achieve a physical resolution of $\sim$80 pc for this $2<z<3$ LAE sample, allowing a detailed characterization of the LAE rest-frame ultraviolet continuum surface brightness profiles and substructure. We use lens-model reconstructions of the LAEs to identify and model individual clumps, which we subsequently use to constrain the parameters of a generative statistical model of the LAE population. Since the BELLS GALLERY sample is selected primarily on the basis of Lyman-$\alpha$ emission, the LAEs that we study here are likely to be directly comparable to those selected in wide-field narrow-band LAE surveys, in contrast with the lensed LAEs identified in cluster lensing fields. We find an LAE clumpiness fraction of approximately 88%, significantly higher than found in previous (non-lensing) studies. We find a well-resolved characteristic clump half-light radii of $\sim$350 pc, a scale comparable to the largest H II regions seen in the local universe. This statistical characterization of LAE surface-brightness profiles will be incorporated into future lensing analyses using the BELLS GALLERY sample to constrain the incidence of dark-matter substructure in the foreground lensing galaxies.

We study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

We investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present the results of a blind millimeter line emitter search using ALMA Band 6 data with a single frequency tuning toward four gravitational lensing clusters (RXJ1347.5-1145, Abell S0592, MACS J0416.1-2403, and Abell 2744). We construct three-dimensional signal-to-noise ratio (S/N) cubes with 60 MHz and 100 MHz binning, and search for millimeter line emitters. We do not detect any line emitters with a peak S/N > 5, although we do find a line emitter candidate with a peak S/N ~ 4.5. These results provide upper limits to the CO(3-2), CO(4-3), CO(5-4), and [CII] luminosity functions at z ~ 0.3, 0.7, 1.2, and 6, respectively. Because of the magnification effect of gravitational lensing clusters, the new data provide the first constraints on the CO and [CII] luminosity functions at unprecedentedly low luminosity levels, i.e., down to $<~ 10^{-3}-10^{-1} $Mpc$^{-3}$ dex$^{-1}$ at $L'_{CO} ~ 10^8-10^{10}$ K km $s^{-1}$ pc$^2$ and $<~ 10^{-3}-10^{-2}$ Mpc$^{-3}$ dex$^{-1}$ at $L_[CII] ~ 10^8-10^{10} L_\odot$, respectively. Although the constraints to date are not stringent yet, we find that the evolution of the CO and [CII] luminosity functions are broadly consistent with the predictions of semi-analytical models. This study demonstrates that the wide observations with a single frequency tuning toward gravitational lensing clusters are promising for constraining the CO and [CII] luminosity functions.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

We present long-baseline ALMA observations of the strong gravitational lens H-ATLAS J090740.0-004200 (SDP.9), which consists of an elliptical galaxy at $z_{\mathrm{L}}=0.6129$ lensing a background submillimeter galaxy into two extended arcs. The data include Band 6 continuum observations, as well as CO $J$=6$-$5 molecular line observations, from which we measure an updated source redshift of $z_{\mathrm{S}}=1.5747$. The image morphology in the ALMA data is different from that of the HST data, indicating a spatial offset between the stellar, gas, and dust component of the source galaxy. We model the lens as an elliptical power law density profile with external shear using a combination of archival HST data and conjugate points identified in the ALMA data. Our best model has an Einstein radius of $\theta_{\mathrm{E}}=0.66\pm0.01$ and a slightly steeper than isothermal mass profile slope. We search for the central image of the lens, which can be used constrain the inner mass distribution of the lens galaxy including the central supermassive black hole, but do not detect it in the integrated CO image at a 3$\sigma$ rms level of 0.0471 Jy km s$^{-1}$.

We present wide-field (167 deg$^2$) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of $\rho=0.54\pm0.03$ (for a smoothing size of $8'$). This correlation is detected even with a smaller smoothing scale of $2'$ ($\rho=0.34\pm 0.01$). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25$ arcmin$^{-2}$). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses.

Accurate astronomical distance determination is crucial for all fields in astrophysics, from Galactic to cosmological scales. Despite, or perhaps because of, significant efforts to determine accurate distances, using a wide range of methods, tracers, and techniques, an internally consistent astronomical distance framework has not yet been established. We review current efforts to homogenize the Local Group's distance framework, with particular emphasis on the potential of RR Lyrae stars as distance indicators, and attempt to extend this in an internally consistent manner to cosmological distances. Calibration based on Type Ia supernovae and distance determinations based on gravitational lensing represent particularly promising approaches. We provide a positive outlook to improvements to the status quo expected from future surveys, missions, and facilities. Astronomical distance determination has clearly reached maturity and near-consistency.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We present the procedure to build and validate the bright-star masks for the Hyper-Suprime-Cam Strategic Subaru Proposal (HSC-SSP) survey. To identify and mask the saturated stars in the full HSC-SSP footprint, we rely on the Gaia and Tycho-2 star catalogues. We first assemble a pure star catalogue down to $G_{\rm Gaia} < 18$ after removing $\sim1.5\%$ of sources that appear extended in the Sloan Digital Sky Survey (SDSS). We perform visual inspection on the early data from the S16A internal release of HSC-SSP, finding that our star catalogue is $99.2\%$ pure down to $G_{\rm Gaia} < 18$. Second, we build the mask regions in an automated way using stacked detected source measurements around bright stars binned per $G_{\rm Gaia}$ magnitude. Finally, we validate those masks from visual inspection and comparison with the literature of galaxy number counts and angular two-point correlation functions. This version (Arcturus) supersedes the previous version (Sirius) used in the S16A internal and DR1 public releases. We publicly release the full masks and tools to flag objects in the entire footprint of the planned HSC-SSP observations at this address: ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/.

We present a list of galaxy-scale lens candidates including a highly probable interacting galaxy-scale lens in the Hyper Suprime-Cam (HSC) imaging survey. We combine HSC imaging with the blended-spectra catalog from the Galaxy And Mass Assembly (GAMA) survey to identify lens candidates, and use lens mass modeling to confirm the candidates. We find 46 matches between the HSC S14A_0b imaging data release and the GAMA catalog. Ten of them are probable lens systems according to their morphology and redshifts. There is one system with an interacting galaxy pair, HSC J084928+000949, that has a valid mass model. We predict the total mass enclosed by the Einstein radius of $\sim0.72$" ($\sim1.65$kpc) for this new expected lens system to be $\sim10^{10.59}M_{\odot}$. Using the photometry in the {\it grizy} bands of the HSC survey and stellar population synthesis modeling with a Salpeter stellar initial mass function, we estimate the stellar mass within the Einstein radius to be $\sim10^{10.46}\,M_{\odot}$. We thus find a dark matter mass fraction within the Einstein radius of $\sim25\%$. Further spectroscopy or high-resolution imaging would allow confirmation of the nature of these lens candidates. The particular system with the interacting galaxy pair, if confirmed, would provide an opportunity to study the interplay between dark matter and stars as galaxies build up through hierarchical mergers.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We utilize the HSC CAMIRA cluster catalog and the photo-$z$ galaxy catalog constructed in the HSC wide field (S16A), covering $\sim$ 174 deg$^{2}$, to study the star formation activity of galaxies in different environments over 0.2 $<$ $z$ $<$ 1.1. We probe galaxies down to $i \sim$ 26, corresponding to a stellar mass limit of log$_{10}$(M$_*$/M$_{\odot}$) $\sim$ 8.2 and $\sim$ 8.6 for star-forming and quiescent populations, respectively, at $z$ $\sim$ 0.2. The existence of the red sequence for low stellar mass galaxies in clusters suggests that the environmental quenching persists to halt the star formation in the low-mass regime. In addition, star-forming galaxies in groups or clusters are systematically biased toward lower values of specific star formation rate by 0.1 -- 0.3 dex with respect to those in the field and the offsets shows no strong redshift evolution over our redshift range, implying a universal slow quenching mechanism acting in the dense environments since $z$ $\sim$ 1.1. Moreover, the environmental quenching dominates the mass quenching in low mass galaxies, and the quenching dominance reverses in high mass ones. The transition mass is greater in clusters than in groups, indicating that the environmental quenching is more effective for massive galaxies in clusters compared to groups.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

We report the discovery of a very bright (r = 20.16), highly magnified, and yet intrinsically very luminous Ly{\alpha} emitter (LAE) at z = 2.82. This system comprises four images in the observer plane with a maximum separation of ~ 6" and it is lensed by a z = 0.55 massive early-type galaxy. It was initially identified in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey for GALaxy-Ly{\alpha} EmitteR sYstems (BELLS GALLERY) survey, and follow-up imaging and spectroscopic observations using the Gran Telescopio Canarias (GTC) and William Herschel Telescope (WHT) confirmed the lensing nature of this system. A lens model using a singular isothermal ellipsoid in an external shear field reproduces quite well the main features of the system, yielding an Einstein radius of 2.95" +/- 0.10", and a total magnification factor for the LAE of 8.8 +/- 0.4. This LAE is one of the brightest and most luminous galaxy-galaxy strong lenses known. We present initial imaging and spectroscopy showing the basic physical and morphological properties of this lensed system.

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

We present predictions for time delays between multiple images of the gravitationally lensed supernova, iPTF16geu, which was recently discovered from the intermediate Palomar Transient Factory (iPTF). As the supernova is of Type Ia where the intrinsic luminosity is usually well-known, accurately measured time delays of the multiple images could provide tight constraints on the Hubble constant. According to our lens mass models constrained by the {\it Hubble Space Telescope} F814W image, we expect the maximum relative time delay to be less than a day, which is consistent with the maximum of 100 hours reported by Goobar et al. but places a stringent upper limit. Furthermore, the fluxes of most of the supernova images depart from expected values suggesting that they are affected by microlensing. The microlensing timescales are small enough that they may pose significant problems to measure the time delays reliably. Our lensing rate calculation indicates that the occurrence of a lensed SN in iPTF is likely. However, the observed total magnification of iPTF16geu is larger than expected, given its redshift. This may be a further indication of ongoing microlensing in this system.

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$\alpha$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) survey. 17 systems are confirmed to be definite lenses with unambiguous evidence of multiple imaging. The lenses are primarily massive early-type galaxies (ETGs) at redshifts of approximately $0.55$, while the lensed sources are Ly$\alpha$ emitters (LAEs) at redshifts from 2 to 3. Although the \textsl{HST} imaging data are well fit by smooth lens models consisting of singular isothermal ellipsoids in an external shear field, a thorough exploration of dark substructures in the lens galaxies is required. The Einstein radii of the BELLS GALLERY lenses are on average $60\%$ larger than those of the BELLS lenses because of the much higher source redshifts which will allow a detailed investigation of the radius evolution of the mass profile in ETGs. With the aid of the average $\sim 13 \times$ lensing magnification, the LAEs are resolved to comprise individual star-forming knots of a wide range of properties with characteristic sizes from less than 100 pc to several kpc, rest-frame far UV apparent AB magnitudes from 29.6 to 24.2, and typical projected separations of 500 pc to 2 kpc.

Abell 2146 consists of two galaxy clusters that have recently collided close to the plane of the sky, and it is unique in showing two large shocks on $\textit{Chandra X-ray Observatory}$ images. With an early stage merger, shortly after first core passage, one would expect the cluster galaxies and the dark matter to be leading the X-ray emitting plasma. In this regard, the cluster Abell 2146-A is very unusual in that the X-ray cool core appears to lead, rather than lag, the Brightest Cluster Galaxy (BCG) in their trajectories. Here we present a strong lensing analysis of multiple image systems identified on $\textit{Hubble Space Telescope}$ images. In particular, we focus on the distribution of mass in Abell 2146-A in order to determine the centroid of the dark matter halo. We use object colours and morphologies to identify multiple image systems; very conservatively, four of these systems are used as constraints on a lens mass model. We find that the centroid of the dark matter halo, constrained using the strongly lensed features, is coincident with the BCG, with an offset of $\approx$ 2 kpc between the centres of the dark matter halo and the BCG. Thus from the strong lensing model, the X-ray cool core also leads the centroid of the dark matter in Abell 2146-A, with an offset of $\approx$ 30 kpc.

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) Survey, which is a Hubble Space Telescope program to image a sample of galaxy-scale strong gravitational lens candidate systems with high-redshift Ly$\alpha$ emitters (LAEs) as the background sources. The goal of the BELLS GALLERY Survey is to illuminate dark substructures in galaxy-scale halos by exploiting the small-scale clumpiness of rest-frame far-UV emission in lensed LAEs, and to thereby constrain the slope and normalization of the substructure-mass function. In this paper, we describe in detail the spectroscopic strong-lens selection technique, which is based on methods adopted in the previous Sloan Lens ACS (SLACS) Survey, BELLS, and SLACS for the Masses Survey. We present the BELLS GALLERY sample of the 21 highest-quality galaxy--LAE candidates selected from $\approx 1.4 \times 10^6$ galaxy spectra in the BOSS of the Sloan Digital Sky Survey III. These systems consist of massive galaxies at redshifts of approximately 0.5 strongly lensing LAEs at redshifts from 2--3. The compact nature of LAEs makes them an ideal probe of dark substructures, with a substructure-mass sensitivity that is unprecedented in other optical strong-lens samples. The magnification effect from lensing will also reveal the structure of LAEs below 100 pc scales, providing a detailed look at the sites of the most concentrated unobscured star formation in the universe. The source code used for candidate selection is available for download as a part of this release.

We report the serendipitous discovery of HSC J142449-005322, a double source plane lens system in the Hyper Suprime-Cam Subaru Strategic Program. We dub the system Eye of Horus. The lens galaxy is a very massive early-type galaxy with stellar mass of ~7x10^11 Msun located at z_L=0.795. The system exhibits two arcs/rings with clearly different colors, including several knots. We have performed spectroscopic follow-up observations of the system with FIRE on Magellan. The outer ring is confirmed at z_S2=1.988 with multiple emission lines, while the inner arc and counterimage is confirmed at z_S1=1.302. This makes it the first double source plane system with spectroscopic redshifts of both sources. Interestingly, redshifts of two of the knots embedded in the outer ring are found to be offset by delta_z=0.002 from the other knots, suggesting that the outer ring consists of at least two distinct components in the source plane. We perform lens modeling with two independent codes and successfully reproduce the main features of the system. However, two of the lensed sources separated by ~0.7 arcsec cannot be reproduced by a smooth potential, and the addition of substructure to the lens potential is required to reproduce them. Higher-resolution imaging of the system will help decipher the origin of this lensing feature and potentially detect the substructure.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

Gravitational waves from inspiraling compact binaries are known to be an excellent absolute distance indicator, yet it is unclear whether electromagnetic counterparts of these events are securely identified for measuring their redshifts, especially in the case of black hole-black hole mergers such as the one recently observed with the Advanced LIGO. We propose to use the cross-correlation between spatial distributions of gravitational wave sources and galaxies with known redshifts as an alternative means of constraining the distance-redshift relation from gravitational waves. In our analysis, we explicitly include the modulation of the distribution of gravitational wave sources due to weak gravitational lensing. We show that the cross-correlation analysis in next-generation observations will be able to tightly constrain the relation between the absolute distance and the redshift, and therefore constrain the Hubble constant as well as dark energy parameters.

We present the application of the pseudo-spectrum method to galaxy-galaxy lensing. We derive explicit expressions for the pseudo-spectrum analysis of the galaxy-shear cross spectrum, which is the Fourier space counterpart of the stacked galaxy-galaxy lensing profile. The pseudo-spectrum method corrects observational issues such as the survey geometry, masks of bright stars and their spikes, and inhomogeneous noise, which distort the spectrum and also mix the E-mode and the B-mode signals. Using ray-tracing simulations in N-body simulations including realistic masks, we confirm that the pseudo-spectrum method successfully recovers the input galaxy-shear cross spectrum. We also investigate the covariance of the galaxy-shear cross spectrum using the ray-tracing simulations to show that there is an excess covariance relative to the Gaussian covariance at small scales where the shot noise is dominated in the Gaussian approximation. We find that the excess of the covariance is consistent with the expectation from the halo sample variance (HSV), which originates from the matter fluctuations at scales larger than the survey area. We apply the pseudo-spectrum method to the observational data of Canada-France-Hawaii Telescope Lensing survey (CFHTLenS) shear catalogue and three different spectroscopic samples of Sloan Digital Sky Survey Luminous Red Galaxy (SDSS LRG), and Baryon Oscillation Spectroscopic Survey (BOSS) CMASS and LOWZ galaxies. The galaxy-shear cross spectra are significantly detected at the level of 7-10$\sigma$ using the analytic covariance with the HSV contribution included. We also confirm that the observed spectra are consistent with the halo model predictions with the halo occupation distribution parameters estimated from previous work. This work demonstrates the viability of galaxy-galaxy lensing analysis in the Fourier space.

We report the discovery of a gravitationally lensed hyperluminous infrared galaxy (L_IR~10^13 L_sun) with strong radio emission (L_1.4GHz~10^25 W/Hz) at z=2.553. The source was identified in the citizen science project SpaceWarps through the visual inspection of tens of thousands of iJKs colour composite images of Luminous Red Galaxies (LRGs), groups and clusters of galaxies and quasars. Appearing as a partial Einstein ring (r_e~3") around an LRG at z=0.2, the galaxy is extremely bright in the sub-millimetre for a cosmological source, with the thermal dust emission approaching 1 Jy at peak. The redshift of the lensed galaxy is determined through the detection of the CO(3-2) molecular emission line with the Large Millimetre Telescope's Redshift Search Receiver and through [OIII] and H-alpha line detections in the near-infrared from Subaru/IRCS. We have resolved the radio emission with high resolution (300-400 mas) eMERLIN L-band and JVLA C-band imaging. These observations are used in combination with the near-infrared imaging to construct a lens model, which indicates a lensing magnification of ~10x. The source reconstruction appears to support a radio morphology comprised of a compact (<250 pc) core and more extended component, perhaps indicative of an active nucleus and jet or lobe.

We present the results of an imaging observation campaign conducted with the Subaru Telescope adaptive optics system (IRCS+AO188) on 28 gravitationally lensed quasars (23 doubles, 1 quad, and 1 possible triple, and 3 candidates) from the SDSS Quasar Lens Search. We develop a novel modelling technique that fits analytical and hybrid point spread functions (PSFs), while simultaneously measuring the relative astrometry, photometry, as well as the lens galaxy morphology. We account for systematics by simulating the observed systems using separately observed PSF stars. The measured relative astrometry is comparable with that typically achieved with the Hubble Space Telescope, even after marginalizing over the PSF uncertainty. We model for the first time the quasar host galaxies in 5 systems, without a-priory knowledge of the PSF, and show that their luminosities follow the known correlation with the mass of the supermassive black hole. For each system, we obtain mass models far more accurate than those previously published from low-resolution data, and we show that in our sample of lensing galaxies the observed light profile is more elliptical than the mass, for ellipticity > 0.25. We also identify eight doubles for which the sources of external and internal shear are more reliably separated, and should therefore be prioritized in monitoring campaigns aimed at measuring time-delays in order to infer the Hubble constant.

We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

The discrepancy between the amplitudes of matter fluctuations inferred from Sunyaev-Zel'dovich (SZ) cluster number counts, the primary temperature, and the polarization anisotropies of the cosmic microwave background (CMB) measured by the Planck satellite can be reconciled if the local universe is embedded in an under-dense region as shown by Lee, 2014. Here using a simple void model assuming the open Friedmann-Robertson-Walker geometry and a Markov Chain Monte Carlo technique, we investigate how deep the local under-dense region needs to be to resolve this discrepancy. Such local void, if exists, predicts the local Hubble parameter value that is different from the global Hubble constant. We derive the posterior distribution of the local Hubble parameter from a joint fitting of the Planck CMB data and SZ cluster number counts assuming the simple void model. We show that the predicted local Hubble parameter value of $H_{\rm loc}=70.1\pm0.34~{\rm km\,s^{-1}Mpc^{-1}}$ is in better agreement with direct local Hubble parameter measurements, indicating that the local void model may provide a consistent solution to the cluster number counts and Hubble parameter discrepancies.

We conduct precise strong lensing mass modeling of four ${\it Hubble}$ Frontier Fields (HFF) clusters, Abell$~$2744, MACS$~$J0416.1$-$2403, MACS$~$J0717.5$+$3745, and MACS$~$J1149.6$+$2223, for which HFF imaging observations are completed. We construct a refined sample of more than 100 multiple images for each cluster by taking advantage of the full depth HFF images, and conduct mass modeling using the ${\small \rm GLAFIC}$ software, which assumes simply parametrized mass distributions. Our mass modeling also exploits a magnification constraint from the lensed Type Ia supernova HFF14Tom for Abell$~$2744 and positional constraints from the multiple images S1$-$S4 of the lensed supernova SN Refsdal for MACS$~$J1149.6$+$2223. We find that our best-fitting mass models reproduce the observed image positions with RMS errors of $\sim 0.4$ arcsec, which are smaller than RMS errors in previous mass modeling that adopted similar numbers of multiple images. Our model predicts a new image of SN Refsdal with a relative time delay and magnification that are fully consistent with a recent detection of reappearance. We then construct catalogs of $z\sim 6-9$ dropout galaxies behind the four clusters and estimate magnification factors for these dropout galaxies with our best-fitting mass models. The dropout sample from the four cluster fields contains $\sim 120$ galaxies at $z\gtrsim 6$, about 20 of which are predicted to be magnified by a factor of more than 10. Some of the high-redshift galaxies detected in the HFF have lensing-corrected magnitudes of $M_{\rm UV}\sim -15$ to $-14$. Our analysis demonstrates that the HFF data indeed offer an ideal opportunity to study faint high-redshift galaxies. All lensing maps produced from our mass modeling will be made available on the STScI website.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. Current and upcoming imaging surveys will contain thousands of new lensed quasars, augmenting the existing sample by at least two orders of magnitude. To find such lens systems, we built a robot, CHITAH, that hunts for lensed quasars by modeling the configuration of the multiple quasar images. Specifically, given an image of an object that might be a lensed quasar, CHITAH first disentangles the light from the supposed lens galaxy and the light from the multiple quasar images based on color information. A simple rule is designed to categorize the given object as a potential four-image (quad) or two-image (double) lensed quasar system. The configuration of the identified quasar images is subsequently modeled to classify whether the object is a lensed quasar system. We test the performance of CHITAH using simulated lens systems based on the Canada-France-Hawaii Telescope Legacy Survey. For bright quads with large image separations (with Einstein radius $r_{\rm ein}>1.1"$) simulated using Gaussian point-spread functions, a high true-positive rate (TPR) of $\sim$90% and a low false-positive rate of $\sim$$3\%$ show that this is a promising approach to search for new lens systems. We obtain high TPR for lens systems with $r_{\rm ein}\gtrsim0.5"$, so the performance of CHITAH is set by the seeing. We further feed a known gravitational lens system, COSMOS 5921$+$0638, to CHITAH, and demonstrate that CHITAH is able to classify this real gravitational lens system successfully. Our newly built CHITAH is omnivorous and can hunt in any ground-based imaging surveys.

We present spatially-resolved properties of molecular gas and dust in a gravitationally-lensed submillimeter galaxy H-ATLAS J090311.6+003906 (SDP.81) at $z=3.042$ revealed by the Atacama Large Millimeter/submillimeter Array (ALMA). We identified 14 molecular clumps in the CO(5-4) line data, all with a spatial scale of $\sim$50-300 pc in the source plane. The surface density of molecular gas ($\Sigma_{\rm H_2}$) and star-formation rate ($\Sigma_{\rm SFR}$) of the clumps are more than three orders of magnitude higher than those found in local spiral galaxies. The clumps are placed in the `burst' sequence in the $\Sigma_{\rm H_2}$-$\Sigma_{\rm SFR}$ plane, suggesting that $z \sim 3$ molecular clumps follow the star-formation law derived for local starburst galaxies. With our gravitational lens model, the positions in the source plane are derived for the molecular clumps, dust clumps, and stellar components identified in the {\sl Hubble Space Telescope} image. The molecular and dust clumps coexist in a similar region over $\sim$2 kpc, while the stellar components are offset at most by $\sim$5 kpc. The molecular clumps have a systematic velocity gradient in the north-south direction, which may indicate a rotating gas disk. One possible scenario is that the components of molecular gas, dust, and stars are distributed in a several-kpc scale rotating disk, and the stellar emission is heavily obscured by dust in the central star-forming region. Alternatively, SDP.81 can be explained by a merging system, where dusty starbursts occur in the region where the two galaxies collide, surrounded by tidal features traced in the stellar components.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a S/N ratio of 4.5 are identified on the convergence S/N map of a 2.3 square degree field observed during the early commissioning phase of the camera. Multi-color photometric data is used to generate optically selected clusters using the CAMIRA algorithm. The optical cluster positions were correlated with the peak positions from the convergence map. All eight significant peaks have optical counterparts. The velocity dispersion of clusters are evaluated by adopting the Singular Isothemal Sphere (SIS) fit to the tangential shear profiles, yielding virial mass estimates, M500c, of the clusters which range from 2.7x10^13 to 4.4x10^14 solar mass. The number of peaks is considerably larger than the average number expected from LambdaCDM cosmology but this is not extremely unlikely if one takes the large sample variance in the small field into account. We could, however, safely argue that the peak count strongly favours the recent Planck result suggesting high sigma8$value of 0.83. The ratio of stellar mass to the dark matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, fg, in halos is universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible deficit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the cosmic microwave background.

We report a detailed modeling of a mass profile of a $z = 0.2999$ massive elliptical galaxy using 30 milli-arcsecond resolution 1-mm Atacama Large Millimeter/submillimeter Array (ALMA) images of the galaxy-galaxy lensing system SDP.81. The detailed morphology of the lensed multiple images of the $z = 3.042$ infrared-luminous galaxy, which is found to consist of tens of $\lesssim 100$-pc-sized star-forming clumps embedded in a $\sim 2$ kpc disk, are well reproduced by a lensing galaxy modeled by an isothermal ellipsoid with a 400 pc core. The core radius is consistent with that of the visible stellar light, and the mass-to-light ratio of $\sim 2\,M_{\odot}\,L_{\odot}^{-1}$ is comparable to the locally measured value, suggesting that the inner 1 kpc region is dominated by luminous matter. The position of the predicted mass centroid is consistent to within $\simeq 30$ mas with that of a non-thermal source detected with ALMA, which likely traces an active galactic nucleus of the foreground elliptical galaxy. A point source mass of $> 3 \times 10^8\,M_{\odot}$ mimicking a supermassive black hole is required to explain the non-detection of a central image of the background galaxy, although the black hole mass degenerates with the core radius of the elliptical galaxy. The required mass is consistent with that predicted from the well-known correlation between black hole mass and host velocity dispersion. Our analysis demonstrates the power of ALMA imaging of strong gravitational lensing for studying the innermost mass profiles and the central supermassive black hole of distant elliptical galaxies.

We report on the weak lensing detection of a filament between two galaxy clusters at $z=0.55$, CL0015.9+1609 and RX J0018.3+1618. We conduct weak lensing analysis of deep multi-band Subaru/Suprime-Cam images with $Lensfit$. The weak lensing signals from the filament are contaminated by signals from the adjacent massive clusters and we statistically subtract the cluster component using two different methods. Both methods yield consistent shear profiles on the filament with $\gtrsim2\sigma$ significance and the average surface mass density of the filament is $<\Sigma>=(3.20\pm0.10)\times10^{14}h$ M$_\odot$Mpc$^{-2}$, which is in broad agreement with previous studies. On-going surveys such as Hyper Suprime-Cam will identify more filaments, which will serve as a new probe of structure formation in the Universe.

We investigate sizes of $z\sim6-8$ dropout galaxies using the complete data of the Abell 2744 cluster and parallel fields in the Hubble Frontier Fields program. By directly fitting light profiles of observed galaxies with lensing-distorted S\'ersic profiles on the image plane with the \texttt{glafic} software, we accurately measure intrinsic sizes of 31 $z\sim6-7$ and eight $z\sim8$ galaxies, including those as faint as $M_{\mathrm{UV}}\simeq-16.6$. We find that half-light radii $r_\mathrm{e}$ positively correlates with UV luminosity at each redshift, although the correlation is not very tight. Largest ($r_\mathrm{e}>0.8$ kpc) galaxies are mostly red in UV color while smallest ($r_\mathrm{e} < 0.08$ kpc) ones tend to be blue. We also find that galaxies with multiple cores tend to be brighter. Combined with previous results at $2.5\lesssim z\lesssim12$, our result confirms that the average $r_{\mathrm{e}}$ of bright ($(0.3-1)L^\star_{z=3}$) galaxies scales as $r_{\mathrm{e}}\propto(1+z)^{-m}$ with $m=1.24\pm0.1$. We find that the ratio of $r_\mathrm{e}$ to virial radius is virtually constant at $3.3\pm0.1\%$ over a wide redshift range, where the virial radii of hosting dark matter halos are derived based on the abundance matching. This constant ratio is consistent with the disk formation model by Mo et al. (1998) with $j_\mathrm{d}\sim m_\mathrm{d}$, where $j_\mathrm{d}$ and $m_\mathrm{d}$ are the fractions of the angular momentum and mass within halos confined in the disks. A comparison with various types of local galaxies indicates that our galaxies are most similar to circumnuclear star-forming regions of barred galaxies in the sense that a sizable amount of stars are forming in a very small area.

We construct a mass model of the cluster MACS J1149.6+2223 to study the expected properties of multiple images of SN Refsdal, the first example of a gravitationally lensed supernova with resolved multiple images recently reported by Kelly et al. We find that the best-fit model predicts six supernova images in total, i.e., two extra images in addition to the observed four Einstein cross supernova images S1--S4. One extra image is predicted to have appeared about 17 years ago, whereas the other extra image is predicted to appear in about one year from the appearance of S1--S4, which is a testable prediction with near future observations. The predicted magnification factors of individual supernova images range from $\sim 18$ for the brightest image to $\sim 4$ for the faint extra images. Confronting these predictions with future observations should provide an unprecedented opportunity to improve our understanding of cluster mass distributions.

Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.

We give an overview of complementarity and synergy in cosmology between the Square Kilometre Array and future survey projects in other wavelengths. In the SKA era, precision cosmology will be limited by systematic errors and cosmic variance, rather than statistical errors. However, combining and/or cross-correlating multi-wavelength data, from the SKA to the cosmic microwave background, optical/infrared and X-ray, substantially reduce these limiting factors. In this chapter, we summarize future survey projects and show highlights of complementarity and synergy, which can be very powerful to probe major cosmological problems such as dark energy, modified gravity and primordial non-Gaussianity.

How well can we infer host dark matter halo masses of individual galaxies? Based on the halo occupation distribution (HOD) framework, we analytically compute the number of neighboring galaxies within a cylinder of some redshift interval and radius in transverse comoving distance. The result is used to derive the conditional probability distribution function (PDF) of the host halo mass of a galaxy, given the neighboring galaxy counts. We compare our analytic results with those obtained using a realistic mock galaxy catalog, finding reasonable agreements. We find the optimal cylinder radius to be $\sim 0.5-1h^{-1}{\rm Mpc}$ for the inference of halo masses. The PDF is generally broad, and sometimes has two peaks at low- and high-mass regimes, because of the effect of chance projection along the line-of-sight. Potential applications and extensions of the new theoretical framework developed herein are also discussed.

We present the comprehensive analyses of faint dropout galaxies up to $z\sim10$ with the first full-depth data set of Abell 2744 lensing cluster and parallel fields observed by the Hubble Frontier Fields (HFF) program. We identify $54$ dropouts at $z\sim5-10$ in the HFF fields, and enlarge the size of $z\sim9$ galaxy sample obtained to date. Although the number of highly magnified ($\mu\sim10$) galaxies is small due to the tiny survey volume of strong lensing, our study reaches the galaxies' intrinsic luminosities comparable to the deepest-field HUDF studies. We derive UV luminosity functions with these faint dropouts, carefully evaluating the combination of observational incompleteness and lensing effects in the image plane by intensive simulations including magnification, distortion, and multiplication of images, with the evaluations of mass model dependences. Our results confirm that the faint-end slope, $\alpha$, is as steep as $-2$ at $z\sim6-8$, and strengthen the evidence of the rapid decrease of UV luminosity densities, $\rho_\mathrm{UV}$, at $z>8$ from the large $z\sim9$ sample. We examine whether the rapid $\rho_\mathrm{UV}$ decrease trend can reconcile with the large Thomson scattering optical depth, $\tau_\mathrm{e}$, measured by CMB experiments allowing a large space of free parameters such as average ionizing photon escape fraction and stellar-population dependent conversion factor. No parameter set can reproduce both the rapid $\rho_\mathrm{UV}$ decrease and the large $\tau_\mathrm{e}$. It is possible that the $\rho_\mathrm{UV}$ decrease moderates at $z\gtrsim11$, that the free parameters significantly evolve towards high-$z$, or that there exist additional sources of reionization such as X-ray binaries and faint AGNs.

We forecast future constraints on local-type primordial non-Gaussianity parameter $f_{\rm NL}$ with a photometric galaxy survey by Euclid, a continuum galaxy survey by Square Kilometre Array (SKA), and their combination. We derive a general expression for the covariance matrix of the power spectrum estimates of multiple tracers to show how the so-called multitracer technique improves constraints on $f_{\rm NL}$. In particular we clarify the role of the overlap fraction of multiple tracers and the division method of the tracers. Our Fisher matrix analysis indicates that stringent constraints of $\sigma (f_{\rm NL})\lesssim 1$ can be obtained even with a single survey, assuming five mass bins. When Euclid and SKA phase 1 (2) are combined, constraints on $f_{\rm NL}$ are improved to $\sigma (f_{\rm NL})= 0.61~(0.50)$.

We study the geometry and the internal structure of the outflowing wind from the accretion disk of a quasar by observing multiple sightlines with the aid of strong gravitational lensing. Using Subaru/HDS, we performed high-resolution ($R$ $\sim$ 36,000) spectroscopic observations of images A and B of the gravitationally lensed quasar SDSS J1029+2623 (at $z_{em}$ $\sim$ 2.197) whose image separation angle, $\theta$ $\sim$ 22$^{\prime\prime}\!\!$.5, is the largest among those discovered so far. We confirm that the difference in absorption profiles in the images A and B discovered by Misawa et al. (2013) remains unchanged since 2010, implying the difference is not due to time variability of the absorption profiles over the delay between the images, $\Delta t$ $\sim$ 744 days, but rather due to differences along the sightlines. We also discovered time variation of C IV absorption strength in both images A and B, due to change of ionization condition. If a typical absorber's size is smaller than its distance from the flux source by more than five orders of magnitude, it should be possible to detect sightline variations among images of other smaller separation, galaxy-scale gravitationally lensed quasars.

We present the result of Subaru Telescope multi-band adaptive optics observations of the complex gravitationally lensed quasar SDSS J1405+0959, which is produced by two lensing galaxies. These observations reveal dramatically enhanced morphological detail, leading to the discovery of an additional object 0. 26'' from the secondary lensing galaxy, as well as three collinear clumps located in between the two lensing galaxies. The new object is likely to be the third quasar image, although the possibility that it is a galaxy cannot be entirely excluded. If confirmed via future observations, it would be the first three image lensed quasar produced by two galaxy lenses. In either case, we show based on gravitational lensing models and photometric redshift that the collinear clumps represent merging images of a portion of the quasar host galaxy, with a magnification factor of 15 - 20, depending on the model.

We present a new algorithm, CAMIRA, to identify clusters of galaxies in wide-field imaging survey data. We base our algorithm on the stellar population synthesis model to predict colours of red-sequence galaxies at a given redshift for an arbitrary set of bandpass filters, with additional calibration using a sample of spectroscopic galaxies to improve the accuracy of the model prediction. We run the algorithm on ~11960 deg^2 of imaging data from the Sloan Digital Sky Survey (SDSS) Data Release 8 to construct a catalogue of 71743 clusters in the redshift range 0.1<z<0.6 with richness after correcting for the incompleteness of the richness estimate greater than 20. We cross-match the cluster catalogue with external cluster catalogues to find that our photometric cluster redshift estimates are accurate with low bias and scatter, and that the corrected richness correlates well with X-ray luminosities and temperatures. We use the publicly available Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) shear catalogue to calibrate the mass-richness relation from stacked weak lensing analysis. Stacked weak lensing signals are detected significantly for 8 subsamples of the SDSS clusters divided by redshift and richness bins, which are then compared with model predictions including miscentring effects to constrain mean halo masses of individual bins. We find the richness correlates well with the halo mass, such that the corrected richness limit of 20 corresponds to the cluster virial mass limit of about 1 \times 10^14 M_Sun/h for the SDSS DR8 cluster sample.

We present optical and near-IR imaging and spectroscopy of SGAS J105039.6$+$001730, a strongly lensed galaxy at z $=$ 3.6252 magnified by $>$30$\times$, and derive its physical properties. We measure a stellar mass of log(M$_{*}$/M$_{\odot}$) $=$ 9.5 $\pm$ 0.35, star formation rates from [O II]$\lambda$$\lambda$3727 and H-$\beta$ of 55 $\pm$ 20 and 84 $\pm$ 17 M$_{\odot}$ yr$^{-1}$, respectively, an electron density of n$_{e} \leq$ 10$^{3}$ cm$^{-2}$, an electron temperature of T$_{e} \leq$ 14000 K, and a metallicity of 12+log(O/H) $=$ 8.3 $\pm$ 0.1. The strong C III]$\lambda$$\lambda$1907,1909 emission and abundance ratios of C, N, O and Si are consistent with well-studied starbursts at z $\sim$ 0 with similar metallicities. Strong P Cygni lines and He II$\lambda$1640 emission indicate a significant population of Wolf-Rayet stars, but synthetic spectra of individual populations of young, hot stars do not reproduce the observed integrated P Cygni absorption features. The rest-frame UV spectral features are indicative of a young starburst with high ionization, implying either 1) an ionization parameter significantly higher than suggest by rest-frame optical nebular lines, or 2) differences in one or both of the initial mass function and the properties of ionizing spectra of massive stars. We argue that the observed features are likely the result of a superposition of star forming regions with different physical properties. These results demonstrate the complexity of star formation on scales smaller than individual galaxies, and highlight the importance of systematic effects that result from smearing together the signatures of individual star forming regions within galaxies.

Objects of known brightness, like Type Ia supernovae (SNIa), can be used to measure distances. If a massive object warps spacetime to form multiple images of a background SNIa, a direct test of cosmic expansion is also possible. However, these lensing events must first be distinguished from other rare phenomena. Recently, a supernova was found to shine much brighter than normal for its distance, which resulted in a debate: was it a new type of superluminous supernova or a normal SNIa magnified by a hidden gravitational lens? Here we report that a spectrum obtained after the supernova faded away shows the presence of a foreground galaxy--the first found to strongly magnify a SNIa. We discuss how more lensed SNIa may be found than previously predicted.

We study weak lensing properties of filaments that connect clusters of galaxies through large cosmological $N$-body simulations. We select 4639 halo pairs with masses higher than $10^{14}h^{-1}\mathrm{M}_\odot$ from the simulations and investigate dark matter distributions between two haloes with ray-tracing simulations. In order to classify filament candidates, we estimate convergence profiles and perform profile fitting. We find that matter distributions between haloes can be classified in a plane of fitting parameters, which allow us to select straight filaments from the ray-tracing simulations. We also investigate statistical properties of these filaments, finding them to be consistent with previous studies. We find that $35\%$ of halo pairs possess straight filaments, $4\%$ of which can directly be detected at $S/N\geq2$ with weak lensing. Furthermore, we study statistical properties of haloes at the edges of filaments. We find that haloes are preferentially elongated along filamentary structures and are less massive with increasing filament masses. However, the dependence of these halo properties on masses of straight filaments is very weak.

We derive the average mass profile of elliptical galaxies from the ensemble of 161 strong gravitational lens systems selected from several surveys, assuming that the mass profile scales with the stellar mass and effective radius of each lensing galaxy. The total mass profile is well fitted by a power-law \rho(r) \propto r^\gamma with best-fit slope \gamma = -2.11+/-0.05. The decomposition of the total mass profile into stellar and dark matter distributions is difficult due to a fundamental degeneracy between the stellar initial mass function (IMF) and the dark matter fraction f_DM. We demonstrate that this IMF-f_DM degeneracy can be broken by adding direct stellar mass fraction measurements by quasar microlensing observations. Our best-fit model prefers the Salpeter IMF over the Chabrier IMF, and a smaller central dark matter fraction than that predicted by adiabatic contraction models.

At low redshift, a handful of gamma-ray bursts (GRBs) have been discovered with peak luminosities ($L_{\rm iso} < 10^{48.5}~\rm{erg\,s}^{-1}$) substantially lower than the average of the more distant ones ($L_{\rm iso} > 10^{49.5}~\rm{erg\,s}^{-1}$). The properties of several low-luminosity (low-$L$) GRBs indicate that they can be due to shock break-out, as opposed to the emission from ultrarelativistic jets. Owing to this, it is highly debated how both populations are connected, and whether there is a continuum between them. The burst at redshift $z=0.283$ from 2012 April 22 is one of the very few examples of intermediate-$L$ GRBs with a $\gamma$-ray luminosity of $L\sim10^{48.9}~\rm{erg\,s}^{-1}$ that have been detected up to now. Together with the robust detection of its accompanying supernova SN 2012bz, it has the potential to answer important questions on the origin of low- and high-$L$ GRBs and the GRB-SN connection. We carried out a spectroscopy campaign using medium- and low-resolution spectrographs at 6--10-m class telescopes, covering the time span of 37.3 days, and a multi-wavelength imaging campaign from radio to X-ray energies over a duration of $\sim270$ days. Furthermore, we used a tuneable filter centred at H$\alpha$ to map star formation in the host galaxy and the surrounding galaxies. We used these data to extract and model the properties of different radiation components and incorporate spectral-energy-distribution fitting techniques to extract the properties of the host galaxy. Modelling the light curve and spectral energy distribution from the radio to the X-rays revealed the blast-wave to expand with an initial Lorentz factor of $\Gamma_0\sim60$, low for a high-$L$ GRB, and that the afterglow had an exceptional low peak luminosity-density of $\lesssim2\times10^{30}~\rm{erg\,s}^{-1}\,\rm{Hz}^{-1}$ in the sub-mm. [Abridged]

We report the discovery of four doubly imaged quasar lenses. All the four systems are selected as lensed quasar candidates from the Sloan Digital Sky Survey data. We confirm their lensing hypothesis with additional imaging and spectroscopic follow-up observations. The discovered lenses are SDSS J0743+2457 with the source redshift z_s=2.165, the lens redshift z_l=0.381, and the image separation theta=1.034", SDSS J1128+2402 with z_s=1.608 and theta=0.844", SDSS J1405+0959 with z_s=1.810, z_l~0.66, and theta=1.978", and SDSS J1515+1511 with z_s=2.054, z_l=0.742, and theta=1.989". It is difficult to estimate the lens redshift of SDSS J1128+2402 from the current data. Two of the four systems (SDSS J1405+0959 and SDSS J1515+1511) are included in our final statistical lens sample to derive constraints on dark energy and the evolution of massive galaxies.

We present an analysis of the line-of-sight structure toward a sample of ten strong lensing cluster cores. Structure is traced by groups that are identified spectroscopically in the redshift range, 0.1 $\leq$ z $\leq$ 0.9, and we measure the projected angular and comoving separations between each group and the primary strong lensing clusters in each corresponding line of sight. From these data we measure the distribution of projected angular separations between the primary strong lensing clusters and uncorrelated large scale structure as traced by groups. We then compare the observed distribution of angular separations for our strong lensing selected lines of sight against the distribution of groups that is predicted for clusters lying along random lines of sight. There is clear evidence for an excess of structure along the line of sight at small angular separations ($\theta \leq 6'$) along the strong lensing selected lines of sight, indicating that uncorrelated structure is a significant systematic that contributes to producing galaxy clusters with large cross sections for strong lensing. The prevalence of line-of-sight structure is one of several biases in strong lensing clusters that can potentially be folded into cosmological measurements using galaxy cluster samples. These results also have implications for current and future studies -- such as the Hubble Space Telescope Frontier Fields -- that make use of massive galaxy cluster lenses as precision cosmological telescopes; it is essential that the contribution of line-of-sight structure be carefully accounted for in the strong lens modeling of the cluster lenses.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

The baryon acoustic oscillation (BAO) experiment requires a sufficiently dense sampling of large-scale structure tracers with spectroscopic redshift, which is observationally expensive especially at high redshifts $z\simgt 1$. Here we present an alternative route of the BAO analysis that uses the cross-correlation of sparse spectroscopic tracers with a much denser photometric sample, where the spectroscopic tracers can be quasars or bright, rare galaxies that are easier to access spectroscopically. We show that measurements of the cross-correlation as a function of the transverse comoving separation rather than the angular separation avoid a smearing of the BAO feature without mixing the different scales at different redshifts in the projection, even for a wide redshift slice $\Delta z\simeq 1$. The bias, scatter, and catastrophic redshift errors of the photometric sample affect only the overall normalization of the cross-correlation which can be marginalized over when constraining the angular diameter distance. As a specific example, we forecast an expected accuracy of the BAO geometrical test via the cross-correlation of the SDSS and BOSS spectroscopic quasar sample with a dense photometric galaxy sample that is assumed to have a full overlap with the SDSS/BOSS survey region. We show that this cross-correlation BAO analysis allows us to measure the angular diameter distances to a fractional accuracy of about 10% at each redshift bin over $1\simlt z\simlt 3$, if the photometric redshift errors of the galaxies, $\sigma_z/(1+z)$, are better than 10-20% level.

Weak lensing provides an important route toward collecting samples of clusters of galaxies selected by mass. Subtle systematic errors in image reduction can compromise the power of this technique. We use the B-mode signal to quantify this systematic error and to test methods for reducing this error. We show that two procedures are efficient in suppressing systematic error in the B-mode: (1) refinement of the mosaic CCD warping procedure to conform to absolute celestial coordinates and (2) truncation of the smoothing procedure on a scale of 10$^{\prime}$. Application of these procedures reduces the systematic error to 20% of its original amplitude. We provide an analytic expression for the distribution of the highest peaks in noise maps that can be used to estimate the fraction of false peaks in the weak lensing $\kappa$-S/N maps as a function of the detection threshold. Based on this analysis we select a threshold S/N = 4.56 for identifying an uncontaminated set of weak lensing peaks in two test fields covering a total area of $\sim 3$deg$^2$. Taken together these fields contain seven peaks above the threshold. Among these, six are probable systems of galaxies and one is a superposition. We confirm the reliability of these peaks with dense redshift surveys, x-ray and imaging observations. The systematic error reduction procedures we apply are general and can be applied to future large-area weak lensing surveys. Our high peak analysis suggests that with a S/N threshold of 4.5, there should be only 2.7 spurious weak lensing peaks even in an area of 1000 deg$^2$ where we expect $\sim$ 2000 peaks based on our Subaru fields.

Recently, Chornock and co-workers announced the Pan-STARRS discovery of a transient source reaching an apparent peak luminosity of ~4x10^44 erg s^-1. We show that the spectra of this transient source are well fit by normal Type Ia supernova (SNIa) templates. The multi-band colors and light-curve shapes are also consistent with normal SNeIa at the spectroscopically determined redshift of z=1.3883; however, the observed flux is a constant factor of ~30 times too bright in each band over time as compared to the templates. At minimum, this shows that the peak luminosities inferred from the light-curve widths of some SNeIa will deviate significantly from the established, empirical relation used by cosmologists. We argue on physical grounds that the observed fluxes do not reflect an intrinsically luminous SNIa, but rather PS1-10afx is a normal SNIa whose flux has been magnified by an external source. The only known astrophysical source capable of such magnification is a gravitational lens. Given the lack of obvious lens candidates, such as galaxy clusters, in the vicinity, we further argue that the lens is a supermassive black hole or a comparatively low-mass dark matter halo. In this case, the lens continues to magnify the underlying host galaxy light. If confirmed, this discovery could impact a broad range of topics including cosmology, gamma-ray bursts, and dark matter halos.

We study the prospects for measuring the dark matter distribution of voids with stacked weak lensing. We select voids from a large set of $N$-body simulations, and explore their lensing signals with the full ray-tracing simulations including the effect of the large-scale structure along the line-of-sight. The lensing signals are compared with simple void model predictions to infer the three-dimensional mass distribution of voids. We show that the stacked weak lensing signals are detected at significant level (S/N$\geq5$) for a 5000 degree$^2$ survey area, for a wide range of void radii up to $\sim50$ Mpc. The error from the galaxy shape noise little affects lensing signals at large scale. It is also found that dense ridges around voids have a great impact on the weak lensing signals, suggesting that proper modeling of the void density profile including surrounding ridges is essential for extracting the average total underdens mass of voids.

We report the discovery of a pair of quasars at $z=1.487$, with a separation of $8\farcs585\pm0\farcs002$. Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located almost symmetrically between the quasars, creating a cross-like configuration. Based on absorption lines in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at redshift $z=0.899$. This, as well as the similarity of the quasar spectra, suggests that the system is a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing hypothesis implies these galaxies are not isolated, but must be embedded in a dark matter halo of virial mass $\sim 4 \times 10^{14}\ h_{70}^{-1}\ {M}_\odot$ assuming an NFW model with a concentration parameter of $c_{vir}=6$, or a singular isothermal sphere profile with a velocity dispersion of $\sim 670$ km s$^{-1}$. We place constraints on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar spectra and broadband photometry.

We present 279 epochs of optical monitoring data spanning 5.4 years from 2007 January to 2012 June for the largest image separation (22.6 arcsec) gravitationally lensed quasar, SDSS J1029+2623. We find that image A leads the images B and C by dt_AB = (744+-10) days (90% confidence); the uncertainty includes both statistical uncertainties and systematic differences due to the choice of models. With only a ~1% fractional error, the interpretation of the delay is limited primarily by cosmic variance due to fluctuations in the mean line-of-sight density. We cannot separate the fainter image C from image B, but since image C trails image B by only 2-3 days in all models, the estimate of the time delay between image A and B is little affected by combining the fluxes of images B and C. There is weak evidence for a low level of microlensing, perhaps created by the small galaxy responsible for the flux ratio anomaly in this system. Interpreting the delay depends on better constraining the shape of the gravitational potential using the lensed host galaxy, other lensed arcs and the structure of the X-ray emission.

We study the origin of absorption features on the blue side of the C IV broad emission line of the large-separation lensed quasar SDSS J1029+2623 at z_em ~ 2.197. The quasar images, produced by a foreground cluster of galaxies, have a maximum separation angle of ~ 22".5. The large angular separation suggests that the sight-lines to the quasar central source can go through different regions of outflowing winds from the accretion disk of the quasar, providing a unique opportunity to study the structure of outflows from the accretion disk, a key ingredient for the evolution of quasars as well as for galaxy formation and evolution. Based on medium- and high-resolution spectroscopy of the two brightest images conducted at the Subaru telescope, we find that each image has different intrinsic levels of absorptions, which can be attributed either to variability of absorption features over the time delay between the lensed images, ~ 774 days, or to the fine structure of quasar outflows probed by the multiple sight-lines toward the quasar. While both these scenarios are consistent with the current data, we argue that they can be distinguished with additional spectroscopic monitoring observations.

Based on a suite of state-of-the-art high-resolution $N$-body simulations, we revisit the so-called halofit model (Smith et al. 2003) as an accurate fitting formula for the nonlinear matter power spectrum. While the halofit model has been frequently used as a standard cosmological tool to predict the nonlinear matter power spectrum in a universe dominated by cold dark matter, its precision has been limited by the low-resolution of $N$-body simulations used to determine the fitting parameters, suggesting the necessity of improved fitting formula at small scales for future cosmological studies. We run high-resolution $N$-body simulations for 16 cosmological models around the Wilkinson Microwave Anisotropy Probe (WMAP) best-fit cosmological parameters (1, 3, 5, and 7 year results), including dark energy models with a constant equation of state. The simulation results are used to re-calibrate the fitting parameters of the halofit model so as to reproduce small-scale power spectra of the $N$-body simulations, while keeping the precision at large scales. The revised fitting formula provides an accurate prediction of the nonlinear matter power spectrum in a wide range of wavenumber ($k \leq 30h$\,Mpc$^{-1}$) at redshifts $0 \leq z \leq 10$, with 5% precision for $k\leq1 h$ Mpc$^{-1}$ at $0 \leq z \leq 10$ and 10% for $1 \leq k\leq 10 h$ Mpc$^{-1} $ at $0 \leq z \leq 3$. We discuss the impact of the improved halofit model on weak lensing power spectra and correlation functions, and show that the improved model better reproduces ray-tracing simulation results.

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be \theta_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

We characterise the typical offset between the Dark Matter (DM) projected centre and the Brightest Cluster Galaxy (BCG) in 10,000 SDSS clusters. To place constraints on the centre of DM, we use an automated strong-lensing analysis, mass-modelling technique which is based on the well-tested assumption that light traces mass. The cluster galaxies are modelled with a steep power-law, and the DM component is obtained by smoothing the galaxy distribution fitting a low-order 2D polynomial (via spline interpolation), while probing a whole range of polynomial degrees and galaxy power laws. We find that the offsets between the BCG and the peak of the smoothed light map representing the DM, \Delta, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(\Delta [h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and \sigma=0.501\pm0.004 (95% confidence levels), or <log_{10}(\Delta [\arcsec])>=0.564\pm0.005, and \sigma=0.475\pm0.007. Some of the offsets originate in prior misidentifications of the BCG or other bright cluster members by the cluster finding algorithm, whose level we make an additional effort to assess, finding that ~10% of the clusters in the probed catalogue are likely to be misidentified, contributing to higher-end offsets in general agreement with previous studies. Our results constitute the first statistically-significant high-resolution distributions of DM-to-BCG offsets obtained in an observational analysis, and importantly show that there exists such a typical non-zero offset in the probed catalogue. The offsets show a weak positive correlation with redshift, so that higher separations are generally found for higher-z clusters in agreement with the hierarchical growth of structure, which in turn could help characterise the merger, relaxation and evolution history of clusters, in future studies. [ABRIDGED]

We examine scatter and bias in weak lensing selected clusters, employing both an analytic model of dark matter haloes and numerical mock data of weak lensing cluster surveys. We pay special attention to effects of the diversity of dark matter distributions within clusters. We find that peak heights of the lensing convergence map correlates rather poorly with the virial mass of haloes. The correlation is tighter for the spherical overdensity mass with a higher mean interior density. We examine the dependence of the halo shape on the peak heights, and find that the rms scatter caused by the halo diversity scales linearly with the peak heights with the proportionality factor of 0.1-0.2. The noise originated from the halo shape is found to be comparable to the source galaxy shape noise and the cosmic shear noise. We find the significant halo orientation bias, i.e., weak lensing selected clusters on average have their major axes aligned with the line-of-sight direction. We compute the orientation bias using an analytic triaxial halo model and obtain results quite consistent with the ray-tracing results. We develop a prescription to analytically compute the number count of weak lensing peaks taking into account all the main sources of scatters in peak heights. We find that the improved analytic predictions agree well with the simulation results for high S/N peaks. We also compare the expected number count with our weak lensing analysis results for 4 sq deg of Subaru/Suprime-Cam observations and find a good agreement.

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

We measure the small-scale (comoving separation 10 kpc/h < r_p < 200 kpc/h) two-point correlation function of quasars using a sample of 26 spectroscopically confirmed binary quasars at 0.6<z<2.2 from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). Thanks to careful candidate selections and extensive follow-up observations of the SQLS, which is aimed at constructing a complete quasar lens sample, our sample of binary quasars is also expected to be nearly complete within a specified range of angular separations and redshifts. The measured small-scale correlation function rises steeply toward smaller scales, which is consistent with earlier studies based on incomplete or smaller binary quasar samples. We find that the quasar correlation function can be fitted by a power-law reasonably well over 4 order of magnitudes, with the best-fit slope of xi(r)\propto r^{-1.92}. We interpret the measured correlation function within the framework of the Halo Occupation Distribution (HOD). We propose a simple model which assumes a constant fraction of quasars that appear as satellites in dark matter haloes, and find that measured small-scale clustering signals constrain the satellite fraction to f_sat=0.054_{-0.016}^{+0.017} for a singular isothermal sphere number density profile of satellites. We note that the HOD modelling appears to underpredict clustering signals at the smallest separations of r_p ~ 10 kpc/h unless we assume very steep number density profiles (such as an NFW profile with the concentration parameter c_vir > 30), which may be suggestive of enhanced quasar activities by direct interactions.

We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<\theta<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.

We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to \Omega_\Lambda=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as \nu_n=d\ln\phi_*/d\ln(1+z) and the velocity dispersion evolution \nu_\sigma=d\ln\sigma_*/d\ln(1+z) are constrained to \nu_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and \nu_\sigma=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.

We present results from strong-lens modelling of 10,000 SDSS clusters, to establish the universal distribution of Einstein radii. Detailed lensing analyses have shown that the inner mass distribution of clusters can be accurately modelled by assuming light traces mass, successfully uncovering large numbers of multiple-images. Approximate critical curves and the effective Einstein radius of each cluster can therefore be readily calculated, from the distribution of member galaxies and scaled by their luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS and HST to calibrate and test this method, and show that an accurate determination of the Einstein radius and mass can be achieved by this approach "blindly", in an automated way, and without requiring multiple images as input. We present the results of the first 10,000 clusters analysed in the range $0.1<z<0.55$, and compare them to theoretical expectations. We find that for this all-sky representative sample the Einstein radius distribution is log-normal in shape, with $< Log(\theta_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $\sigma=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $\theta_{e}$ clusters than predicted by $\Lambda$CDM. We visually inspect each of the clusters with $\theta_{e}>40 \arcsec$ ($z_{s}=2$) and find that $\sim20%$ are boosted by various projection effects detailed here, remaining with $\sim40$ real giant-lens candidates, with a maximum of $\theta_{e}=69\pm12 \arcsec$ ($z_{s}=2$) for the most massive candidate, in agreement with semi-analytic calculations. The results of this work should be verified further when an extended calibration sample is available.

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

Gravitational lens systems containing lensed quasars are important as cosmological probes, as diagnostics of structural properties of the lensing galaxies and as tools to study the quasars themselves. The largest lensed quasar sample is the SDSS Quasar Lens Search, drawn from the Sloan Digital Sky Survey (SDSS). We are attempting to extend this survey using observations of lens candidates selected from a combination of the quasar sample from the SDSS and the UKIRT Infrared Deep Sky Survey (UKIDSS). This adds somewhat higher image quality together with a wider range of wavelength for the selection process. In previous pilot surveys we observed 5 objects, finding 2 lenses; here we present further observations of 20 objects in which we find 4 lenses, of which 2 are independently discovered in SQLS (in preparation). Following earlier work on the combination of these two surveys, we have refined our method and find that use of a colour-separation diagnostic, where we select for separations between components which appear to decrease in wavelength, is an efficient method to find lensed quasars and may be useful in ongoing and future large-scale strong lensing surveys with instruments such as Pan-STARRS and LSST. The new lenses have mostly high flux ratios, with faint secondaries buried in the lensing galaxy and typically 6-10 times less bright than the primary. Our survey brings the total number of lenses discovered in the SDSS quasar sample to 46, plus 13 lenses already known. This is likely to be up to 60-70% of the total number of lensed quasars; we briefly discuss strategies by which the rest might be found.

We perform high resolution ray-tracing simulations to investigate probability distribution functions (PDFs) of lensing convergence, shear, and magnification on distant sources up to the redshift of $z_s=20$. We pay particular attention to the shot noise effect in $N$-body simulations by explicitly showing how it affects the variance of the convergence. We show that the convergence and magnification PDFs are closely related with each other via the approximate relation $\mu=(1-\kappa)^{-2}$, which can reproduce the behavior of PDFs surprisingly well up to the high magnification tail. The mean convergence measured in the source plane is found to be systematically negative, rather than zero as often assumed, and is correlated with the convergence variance. We provide simple analytical formulae for the PDFs, which reproduce simulated PDFs reasonably well for a wide range of redshifts and smoothing sizes. As explicit applications of our ray-tracing simulations, we examine the strong lensing probability and the magnification effects on the luminosity functions of distant galaxies and quasars.

The quasar SDSS J133401.39+331534.3 at z = 2.426 is found to be a two-image gravitationally lensed quasar with the image separation of 0.833. The object is first identified as a lensed quasar candidate in the Sloan Digital Sky Survey Quasar Lens Search, and then confirmed as a lensed system from follow-up observations at the Subaru and University of Hawaii 2.2-meter telescopes. We estimate the redshift of the lensing galaxy to be 0.557 based on absorption lines in the quasar spectra as well as the color of the galaxy. In particular, we observe the system with the Subaru Telescope AO188 adaptive optics with laser guide star, in order to derive accurate astrometry, which well demonstrates the usefulness of the laser guide star adaptive optics imaging for studying strong lens systems. Our mass modeling with improved astrometry implies that a nearby bright galaxy $\sim 4"$ apart from the lensing galaxy is likely to affect the lens potential.

Using a large set of ray-tracing in N-body simulations, we examine lensing profiles around massive dark haloes in detail, with a particular emphasis on the profile at around the virial radii. We compare radial convergence profiles, which are measured accurately in the ray-tracing simulations by stacking many dark haloes, with our simple analytic model predictions. Our analytic models consist of a main halo, which is modelled by the Navarro-Frenk-White (NFW) density profile with three different forms of the truncation, plus the correlated matter (2-halo term) around the main halo. We find that the smoothly truncated NFW profile best reproduces the simulated lensing profiles, out to more than 10 times the virial radius. We then use this analytic model to investigate potential biases in cluster weak lensing studies in which a single, untruncated NFW component is usually assumed in interpreting observed signals. We find that cluster masses, inferred by fitting reduced tangential shear profiles with the NFW profile, tend to be underestimated by ~5-10% if fitting is performed out to ~10'-30'. In contrast, the concentration parameter is overestimated typically by ~20% for the same fitting range. We also investigate biases in computing the signal-to-noise ratio of weak lensing mass peaks, finding them to be <4% for significant mass peaks. In the Appendices, we provide useful formulae for the smoothly truncated NFW profile.

We present results from a spectroscopic program targeting 26 strong lensing cluster cores that were visually identified in the Sloan Digital Sky Survey (SDSS) and the Second Red-Sequence Cluster Survey (RCS-2). The 26 galaxy cluster lenses span a redshift range of 0.2 < z < 0.65, and our spectroscopy reveals 69 unique background sources with redshifts as high as z=5.200. We also identify redshifts for 262 cluster member galaxies and measure the velocity dispersions and dynamical masses for 18 clusters where we have redshifts for N \geq 10 cluster member galaxies. We include an accounting for the expected biases in dynamical masses of strong lensing selected clusters as predicted by results from numerical simulations and discuss possible sources of bias in our observations. The median dynamical mass of the 18 clusters with N \geq 10 spectroscopic cluster members is M_{Vir} = 7.84 x 10^14 M_sun h_{0.7}^{-1}, which is somewhat higher than predictions for strong lensing selected clusters in simulations. The disagreement is not significant considering the large uncertainty in our dynamical data, systematic uncertainties in the velocity dispersion calibration, and limitations of the theoretical modeling. Nevertheless our study represents an important first step toward characterizing large samples of clusters that are identified in a systematic way as systems exhibiting dramatic strong lensing features.

Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing substructure close to one of the components can account for a flux anomaly with negligible changes in the component positions. Our best fit model has a substructure mass of ~10^8 solar masses up to the mass-sheet degeneracy, located roughly 0.1 arcsecs West and 0.1 arcsecs North of component B. We demonstrate that a positional offset between the centers of the source components can explain the differences between the optical and radio flux ratios.

We measure the redshift distribution of a sample of 28 giant arcs discovered as a part of the Sloan Giant Arcs Survey (SGAS). Gemini/GMOS-North spectroscopy provides precise redshifts for 24 arcs, and "redshift desert" constraints for the remaining four. This is a direct measurement of the redshift distribution of a uniformly selected sample of bright giant arcs, which is an observable that can be used to inform efforts to predict giant arc statistics. Our primary giant arc sample has a median redshift z=1.821 and nearly two thirds of the arcs - 64% - are sources at z \gtrsim 1.4, indicating that the population of background sources that are strongly lensed into bright giant arcs resides primarily at high redshift. We also analyze the distribution of redshifts for 19 secondary strongly lensed background sources that are not visually apparent in SDSS imaging, but were identified in deeper follow-up imaging of the lensing cluster fields. Our redshift sample for the secondary sources is not spectroscopically complete, but combining it with our primary giant arc sample suggests that a large fraction of all background galaxies which are strongly lensed by foreground clusters reside at z \gtrsim 1.4. Kolmogorov-Smirnov (KS) tests indicate that our well-selected, spectroscopically complete primary giant arc redshift sample can be reproduced with a model distribution that is constructed from a combination of results from studies of strong lensing clusters in numerical simulations, and observational constraints on the galaxy luminosity function.

We develop a new method of combining cluster observables (number counts and cluster-cluster correlation functions) and stacked weak lensing signals of background galaxy shapes, both of which are available in a wide-field optical imaging survey. Assuming that the clusters have secure redshift estimates, we show that the joint experiment enables a self-calibration of important systematic errors including the source redshift uncertainty and the cluster mass-observable relation, by adopting a single population of background source galaxies for the lensing analysis. It allows us to use the relative strengths of stacked lensing signals at different cluster redshifts for calibrating the source redshift uncertainty, which in turn leads to accurate measurements of the mean cluster mass in each bin. In addition, our formulation of stacked lensing signals in Fourier space simplifies the Fisher matrix calculations, as well as the marginalization over the cluster off-centering effect, the most significant uncertainty in stacked lensing. We show that upcoming wide-field surveys yield stringent constraints on cosmological parameters including dark energy parameters, without any priors on nuisance parameters that model systematic uncertainties. Specifically, the stacked lensing information improves the dark energy FoM by a factor of 4, compared to that from the cluster observables alone. The primordial non-Gaussianity parameter can also be constrained with a level of f_NL~10. In this method, the mean source redshift is well calibrated to an accuracy of 0.1 in redshift, and the mean cluster mass in each bin to 5-10% accuracies, which demonstrates the success of the self-calibration of systematic uncertainties from the joint experiment. (Abridged)

We have measured the Sunyaev Zel'dovich (SZ) effect for a sample of ten strong lensing selected galaxy clusters using the Sunyaev Zel'dovich Array (SZA). The SZA is sensitive to structures on spatial scales of a few arcminutes, while the strong lensing mass modeling constrains the mass at small scales (typically < 30"). Combining the two provides information about the projected concentrations of the strong lensing clusters. The Einstein radii we measure are twice as large as expected given the masses inferred from SZ scaling relations. A Monte Carlo simulation indicates that a sample randomly drawn from the expected distribution would have a larger median Einstein radius than the observed clusters about 3% of the time. The implied overconcentration has been noted in previous studies with smaller samples of lensing clusters. It persists for this sample, with the caveat that this could result from a systematic effect such as if the gas fractions of the strong lensing clusters are substantially below what is expected.

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<\theta<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be \Omega_\Lambda=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

We report the discovery of a gravitationally lensed quasar identified serendipitously in the Sloan Digital Sky Survey (SDSS). The object, SDSS J094604.90+183541.8, was initially targeted for spectroscopy as a luminous red galaxy, but the SDSS spectrum has the features of both a z=0.388 galaxy and a z=4.8 quasar. We have obtained additional imaging that resolves the system into two quasar images separated by 3.06 arcsec and a bright galaxy that is strongly blended with one of the quasar images. We confirm spectroscopically that the two quasar images represent a single lensed source at z=4.8 with a total magnification of 3.2, and we derive a model for the lensing galaxy. This is the highest redshift lensed quasar currently known. We examine the issues surrounding the selection of such an unusual object from existing data and briefly discuss implications for lensed quasar surveys.

We present a strong lens analysis of SDSS J1004+4112, a unique quasar lens produced by a massive cluster of galaxies at z=0.68, using a newly developed software for gravitational lensing. We find that our parametric mass model well reproduces all observations including the positions of quasar images as well as those of multiply imaged galaxies with measured spectroscopic redshifts, time delays between quasar images, and the positions of faint central images. The predicted large total magnification of \mu ~ 70 suggests that the lens system is indeed a useful site for studying the fine structure of a distant quasar and its host galaxy. The dark halo component is found to be unimodal centered on the brightest cluster galaxy and the Chandra X-ray surface brightness profile. In addition, the orientation of the halo component is quite consistent with those of the brightest cluster galaxy and member galaxy distribution, implying that the lensing cluster is a relaxed system. The radial profile of the best-fit mass model is in good agreement with a mass profile inferred from the X-ray observation. While the inner radial slope of the dark halo component is consistent with being -1, a clear dependence of the predicted A-D time delay on the slope indicates that an additional time delay measurement will improve constraints on the mass model.

We present new measurements of dark matter distributions in 25 X-ray luminous clusters by making a full use of the two-dimensional (2D) weak lensing signals obtained from high-quality Subaru/Suprime-Cam imaging data. Our approach to directly compare the measured lensing shear pattern with elliptical model predictions allows us to extract new information on the mass distributions of individual clusters, such as the halo ellipticity and mass centroid. We find that these parameters on the cluster shape are little degenerate with cluster mass and concentration parameters. By combining the 2D fitting results for a subsample of 18 clusters, the elliptical shape of dark matter haloes is detected at 7\sigma significance level. The mean ellipticity is found to be e = 0.46 \pm 0.04 (1\sigma), which is in excellent agreement with the standard collisionless CDM model prediction. The mass centroid can be constrained with a typical accuracy of ~20" (~50 kpc/h) in radius for each cluster with some significant outliers, enabling to assess one of the most important systematic errors inherent in the stacked cluster weak lensing technique, the mass centroid uncertainty. In addition, the shape of the dark mass distribution is found to be only weakly correlated with that of the member galaxy distribution. We carefully examine possible sources of systematic errors in our measurements, finding none of them to be significant. Our results demonstrate the power of high-quality imaging data for exploring the detailed spatial distribution of dark matter (Abridged).

We report the discovery and confirmation of eight new two-image lensed quasars by the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. The lenses are SDSSJ0904+1512 (image separation \theta=1"13, source redshift z_s=1.826), SDSSJ1054+2733 (\theta=1"27, z_s=1.452), SDSSJ1055+4628 (\theta=1"15, z_s=1.249), SDSSJ1131+1915 (\theta=1"46, z_s=2.915), SDSSJ1304+2001 (\theta=1"87, z_s=2.175), SDSSJ1349+1227 (\theta=3"00, z_s=1.722), SDSSJ1455+1447 (\theta=1"73, z_s=1.424), and SDSSJ1620+1203 (\theta=2"77, z_s=1.158). Three of them, SDSSJ1055+4628, SDSSJ1455+1447, and SDSSJ1620+1203, satisfy the criteria for constructing our statistical sample for studying the cosmological model. Based on galactic absorption lines of the lens galaxies, we also derive lens redshifts of z_l=0.398 and z_l=0.513 for SDSSJ1620+1203 and the previously discovered lens SDSSJ0746+4403, respectively.

Cadenced optical imaging surveys in the next decade will be capable of detecting time-varying galaxy-scale strong gravitational lenses in large numbers, increasing the size of the statistically well-defined samples of multiply-imaged quasars by two orders of magnitude, and discovering the first strongly-lensed supernovae. We carry out a detailed calculation of the likely yields of several planned surveys, using realistic distributions for the lens and source properties and taking magnification bias and image configuration detectability into account. We find that upcoming wide-field synoptic surveys should detect several thousand lensed quasars. In particular, the LSST should find 8000 lensed quasars, 3000 of which will have well-measured time delays, and also ~130 lensed supernovae, which is compared with ~15 lensed supernovae predicted to be found by the JDEM. We predict the quad fraction to be ~15% for the lensed quasars and ~30% for the lensed supernovae. Generating a mock catalogue of around 1500 well-observed double-image lenses, we compute the available precision on the Hubble constant and the dark energy equation parameters for the time delay distance experiment (assuming priors from Planck): the predicted marginalised 68% confidence intervals are \sigma(w_0)=0.15, \sigma(w_a)=0.41, and \sigma(h)=0.017. While this is encouraging in the sense that these uncertainties are only 50% larger than those predicted for a space-based type-Ia supernova sample, we show how the dark energy figure of merit degrades with decreasing knowledge of the the lens mass distribution. (Abridged)

Weak lensing work can be badly compromised by unlensed foreground and cluster members which dilute the true lensing signal. We show how the lensing amplitude in multi-colour space can be harnessed to securely separate cluster members from the foreground and background populations for three massive clusters, A1703 (z=0.258), A370 (z=0.375) and RXJ1347-11 (z=0.451) imaged with Subaru. The luminosity functions of these clusters when corrected for dilution, show similar faint-end slopes, \alpha ~= -1.0, with no marked faint-end upturn to our limit of M_R ~= -15.0, and only a mild radial gradient. In each case, the radial profile of the M/L ratio peaks at intermediate radius, ~=0.2r_{vir}, at a level of 300-500(M/L_R)_\odot, and then falls steadily towards ~100(M/L_R)_{\odot} at the virial radius, similar to the mean field level. This behaviour is likely due to the relative paucity of central late-type galaxies, whereas for the E/S0-sequence only a mild radial decline in M/L is found for each cluster. We discuss this behaviour in the context of detailed simulations where predictions for tidal stripping may now be tested accurately with observations.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present observations of a new double-image gravitational lens system, ULAS J082016.1+081216, of image separation 2.3" and high (~6) flux ratio. The system is selected from the Sloan Digital Sky Survey spectroscopic quasar list using new high-quality images from the UKIRT Deep Sky Survey (UKIDSS). The lensed quasar has a source redshift of 2.024, and we identify the lens galaxy as a faint red object of redshift 0.803+/-0.001. Three other objects from the UKIDSS survey, selected in the same way, were found not to be lens systems. Together with the earlier lens found using this method, the SDSS-UKIDSS lenses have the potential to significantly increase the number of quasar lenses found in SDSS, to extend the survey to higher flux ratios and lower separations, and to give greater completeness which is important for statistical purposes.

We show that the ability to probe primordial non-Gaussianity with cluster counts is drastically improved by adding the excess variance of counts which contains information on the clustering. The conflicting dependences of changing the mass threshold and including primordial non-Gaussianity on the mass function and biasing indicate that the self-calibrated cluster counts well break the degeneracy between primordial non-Gaussianity and the observable-mass relation. Based on the Fisher matrix analysis, we show that the count variance improves constraints on f_NL by more than an order of magnitude. It exhibits little degeneracy with dark energy equation of state. We forecast that upcoming Hyper Suprime-cam cluster surveys and Dark Energy Survey will constrain primordial non-Gaussianity at the level \sigma(f_NL) \sim 8, which is competitive with forecasted constraints from next-generation cosmic microwave background experiments.

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

We report the discovery of five gravitationally lensed quasars from the Sloan Digital Sky Survey (SDSS). All five systems are selected as two-image lensed quasar candidates from a sample of high-redshift (z>2.2) SDSS quasars. We confirmed their lensing nature with additional imaging and spectroscopic observations. The new systems are SDSS J0819+5356 (source redshift z_s=2.237, lens redshift z_l=0.294, and image separation \theta=4.04"), SDSS J1254+2235 (z_s=3.626, \theta=1.56"), SDSS J1258+1657 (z_s=2.702, \theta=1.28"), SDSS J1339+1310 (z_s=2.243, \theta=1.69"), and SDSS J1400+3134 (z_s=3.317, \theta=1.74"). We estimate the lens redshifts of the latter four systems to be z_l=0.2-0.8 from the colors and magnitudes of the lensing galaxies. We find that the image configurations of all systems are well reproduced by standard mass models. Although these lenses will not be included in our statistical sample of z_s<2.2 lenses, they expand the number of lensed quasars which can be used for high-redshift galaxy and quasar studies.

Gravitational lenses on galaxy scales are plausibly modelled as having ellipsoidal symmetry and a universal dark matter density profile, with a Sersic profile to describe the distribution of baryonic matter. Predicting all lensing effects requires knowledge of the total lens potential: in this work we give analytic forms for that of the above hybrid model. Emphasising that complex lens potentials can be constructed from simpler components in linear combination, we provide a recipe for attaining elliptical symmetry in either projected mass or lens potential. We also provide analytic formulae for the lens potentials of Sersic profiles for integer and half-integer index. We then present formulae describing the gravitational lensing effects due to smoothly-truncated universal density profiles in cold dark matter model. For our isolated haloes the density profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, functional forms that allow all orders of lens potential derivatives to be calculated analytically, while ensuring a non-divergent total mass. We show how the observables predicted by this profile differ from that of the original infinite-mass NFW profile. Expressions for the gravitational flexion are highlighted. We show how decreasing the tidal radius allows stripped haloes to be modelled, providing a framework for a fuller investigation of dark matter substructure in galaxies and clusters. Finally we remark on the need for finite mass halo profiles when doing cosmological ray-tracing simulations, and the need for readily-calculable higher order derivatives of the lens potential when studying catastrophes in strong lenses.

The Einstein radius plays a central role in lens studies as it characterises the strength of gravitational lensing. The distribution of Einstein radii near the upper cutoff should probe the largest mass concentrations in the universe. Adopting a triaxial halo model, we compute expected distributions of large Einstein radii. To assess the cosmic variance, we generate a number of all-sky Monte-Carlo realisations. We find that the expected largest Einstein radius in the universe is sensitive to the cosmological model: for a source redshift z=1, they are 42^{+9}_{-7}, 35^{+8}_{-6}, and 54^{+12}_{-7} arcseconds, assuming best-fit parameters of the WMAP5, WMAP3 and WMAP1 data, respectively. These values are broadly consistent with current observations given their incompleteness. For the same source redshift, we expect in all-sky 35 (WMAP5), 15 (WMAP3), and 150 (WMAP1) clusters that have Einstein radii larger than 20". Whilst the values of the largest Einstein radii are almost unaffected by the primordial non-Gaussianity currently of interest, the abundance of large lens clusters should probe non-Gaussianity competitively with CMB, but only if other cosmological parameters are well-measured. We also find that these "superlens" clusters constitute a highly biased population. For instance, a substantial fraction of these superlens clusters have major axes preferentially aligned with the line-of-sight. As a consequence, the projected mass distributions of the clusters are rounder by an ellipticity of 0.2 and have 40%-60% larger concentrations compared with typical clusters with similar redshifts and masses. We argue that the large concentration measured in A1689 is consistent with our model prediction at the 1.2\sigma level. (Abridged)

We present the discovery of a new quadruply lensed quasar. The lens system, SDSS J1330+1810 at z_s=1.393, was identified as a lens candidate from the spectroscopic sample of the Sloan Digital Sky Survey. Optical and near-infrared images clearly show four quasar images with a maximum image separation of 1.76", as well as a bright lensing galaxy. We measure a redshift of the lensing galaxy of z_l=0.373 from absorption features in the spectrum. We find a foreground group of galaxies at z=0.31 centred ~120" southwest of the lens system. Simple mass models fit the data quite well, including the flux ratios between images, although the lens galaxy appears to be ~1 mag brighter than expected by the Faber-Jackson relation. Our mass modelling suggests that shear from nearby structure is affecting the lens potential.

We report on the discovery of a very bright z = 2.00 star-forming galaxy that is strongly lensed by a foreground z=0.422 luminous red galaxy (LRG). This system was found in a systematic search for bright arcs lensed by LRGs and brightest cluster galaxies in the Sloan Digital Sky Survey Data Release 5 sample. Follow-up observations on the Subaru 8.2m telescope on Mauna Kea and the Astrophysical Research Consortium 3.5m telescope at Apache Point Observatory confirmed the lensing nature of this system. A simple lens model for the system, assuming a singular isothermal ellipsoid mass distribution, yields an Einstein radius of 3.82 +/- 0.03 arcsec or 14.8 +/- 0.1 kpc/h at the lens redshift. The total projected mass enclosed within the Einstein radius is 2.10 +/- 0.03 x 10^12 M_sun/h, and the magnification factor for the source galaxy is 27 +/- 1. Combining the lens model with our gVriz photometry, we find an (unlensed) star formation rate for the source galaxy of 32 M_sun/h / yr, adopting a fiducial constant star formation rate model with an age of 100 Myr and E(B-V) = 0.25. With an apparent magnitude of r = 19.9, this system is among the very brightest lensed z >= 2 galaxies, and provides an excellent opportunity to pursue detailed studies of the physical properties of an individual high-redshift star-forming galaxy.

We derive lens distortion and magnification profiles of four well known clusters observed with Subaru. Each cluster is very well fitted by the general form predicted for Cold Dark Matter (CDM) dominated halos, with good consistency found between the independent distortion and magnification measurements. The inferred level of mass concentration is surprisingly high, 8<c_{vir}<15 (<c_{vir}> = 10.4 \pm 0.9), compared to the relatively shallow profiles predicted by the Lambda CDM model, c_{vir}=5.1 \pm 1.1 (for <M_{vir}> =1.25\times 10^{15}M_{\odot}/h). This represents a 4sigma discrepancy, and includes the relatively modest effects of projection bias and profile evolution derived from N-body simulations, which oppose each other with little residual effect. In the context of CDM based cosmologies, this discrepancy implies clusters collapse earlier (z\geq 1) than predicted (z<0.5), when the Universe was correspondingly denser.

We present the results of deep spectroscopy for the central region of the cluster lens SDSS J1004+4112 with the Subaru telescope. A secure detection of an emission line of the faint blue stellar object (component E) near the center of the brightest cluster galaxy (G1) confirms that it is the central fifth image of the lensed quasar system. In addition, we measure the stellar velocity dispersion of G1 to be sigma_* = 352+-13 km/s. We combine these results to obtain constraints on the mass M_BH of the putative black hole (BH) at the center of the inactive galaxy G1, and hence on the M_BH-sigma_* relation at the lens redshift z_l=0.68. From detailed mass modeling, we place an upper limit on the black hole mass, M_BH < 2.1x10^{10}M_sun at 1-sigma level (<3.1x10^{10}M_sun at 3-sigma), which is consistent with black hole masses expected from the local and redshift-evolved M_BH-sigma_* relations, M_BH~10^{9}-10^{10}M_sun.

We present the first result of a survey for strong galaxy-galaxy lenses in Sloan Digital Sky Survey (SDSS) images. SDSS J082728.70+223256.4 was selected as a lensing candidate using selection criteria based on the color and positions of objects in the SDSS photometric catalog. Follow-up imaging and spectroscopy showed this object to be a lensing system. The lensing galaxy is an elliptical at z = 0.349 in a galaxy cluster. The lensed galaxy has the spectrum of a post-starburst galaxy at z = 0.766. The lensing galaxy has an estimated mass of $\sim 1.2 \times 10^{12} M_{\odot}$ and the corresponding mass to light ratio in the B-band is $\sim 26 M_{\odot}/L_{\odot}$ inside 1.1 effective radii of the lensing galaxy. Our study shows how catalogs drawn from multi-band surveys can be used to find strong galaxy-galaxy lenses having multiple lens images. Our strong lensing candidate selection based on photometry-only catalogs will be useful in future multi-band imaging surveys such as SNAP and LSST.

We report the discovery of a new gravitational lens system. This object, ULAS J234311.93-005034.0, is the first to be selected by using the new UKIRT Infrared Deep Sky Survey (UKIDSS), together with the Sloan Digital Sky Survey (SDSS). The ULAS J234311.93-005034.0 system contains a quasar at redshift 0.788 which is doubly imaged, with separation 1.4". The two quasar images have the same redshift and similar, though not identical, spectra. The lensing galaxy is detected by subtracting point-spread functions from R-band images taken with the Keck telescope. The lensing galaxy can also be detected by subtracting the spectra of the A and B images, since more of the galaxy light is likely to be present in the latter. No redshift is determined from the galaxy, although the shape of its spectrum suggests a redshift of about 0.3. The object's lens status is secure, due to the identification of two objects with the same redshift together with a lensing galaxy. Our imaging suggests that the lens is found in a cluster environment, in which candidate arc-like structures, that require confirmation, are visible in the vicinity. Further discoveries of lenses from the UKIDSS survey are likely as part of this programme, due to the depth of UKIDSS and its generally good seeing conditions.

We identify a third image in the unique quasar lens SDSS J1029+2623, the second known quasar lens produced by a massive cluster of galaxies. The spectrum of the third image shows similar emission and absorption features, but has a redder continuum than the other two images which can be explained by differential extinction or microlensing. We also identify several lensed arcs. Our observations suggest a complicated structure of the lens cluster at z~0.6. We argue that the three lensed images are produced by a naked cusp on the basis of successful mass models, the distribution of cluster member galaxies, and the shapes and locations of the lensed arcs. Lensing by a naked cusp is quite rare among galaxy-scale lenses but is predicted to be common among large-separation lensed quasars. Thus the discovery can be viewed as support for an important theoretical prediction of the standard cold dark matter model.

Difference imaging provides a new way to discover gravitationally lensed quasars because few non-lensed sources will show spatially extended, time variable flux. We test the method on lens candidates in the Sloan Digital Sky Survey (SDSS) Supernova Survey region from the SDSS Quasar Lens Search (SQLS) and their surrounding fields. Starting from 20768 sources, including 49 SDSS quasars and 36 candidate lenses/lensed images, we find that 21 sources including 15 SDSS QSOs and 7 candidate lenses/lensed images are non-periodic variable sources. We can measure the spatial structure of the variable flux for 18 of these sources and identify only one as a non-point source. This source does not display the compelling spatial structure of the variable flux of known lensed quasars, so we reject it as a lens candidate. None of the lens candidates from the SQLS survive our cuts. Given our effective survey area of order 0.71 square degrees, this indicates a false positive rate of order one per square degree for themethod. The fraction of quasars not found to be variable and the false positive rate should both fall if we analyze the full, later data releases for the SDSS fields. While application of the method to the SDSS is limited by the resolution, depth, and sampling of the survey, several future surveys such as Pan-STARRS, LSST, and SNAP will avoid these limitations.

We present the discovery of four gravitationally lensed quasars selected from the spectroscopic quasar catalog of the Sloan Digital Sky Survey. We describe imaging and spectroscopic follow-up observations that support the lensing interpretation of the following four quasars: SDSS J0832+0404 (image separation \theta=1.98", source redshift z_s=1.115, lens redshift z_l=0.659); SDSS J1216+3529 (\theta=1.49", z_s=2.012); SDSS J1322+1052 (\theta=2.00", z_s=1.716); and SDSS J1524+4409 (\theta=1.67", z_s=1.210, z_l=0.320). Each system has two lensed images. We find that the fainter image component of SDSS J0832+0404 is significantly redder than the brighter component, perhaps because of differential reddening by the lensing galaxy. The lens potential of SDSS J1216+3529 might be complicated by the presence of a secondary galaxy near the main lensing galaxy.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We report the discovery of a new gravitationally lensed quasar, SDSS J131339.98+515128.3, at a redshift of 1.875 with an image separation of 1.24". The lensing galaxy is clearly detected in visible-light follow-up observations. We also identify three absorption-line doublets in the spectra of the lensed quasar images, from which we measure the lens redshift to be 0.194. Like several other known lenses, the lensed quasar images have different continuum slopes. This difference is probably the result of reddening and microlensing in the lensing galaxy. The lensed quasar was selected by correlating Sloan Digital Sky Survey (SDSS) spectroscopic quasars with Two Micron All Sky Survey (2MASS) sources and choosing quasars that show near-infrared (IR) excess. The near-IR excess can originate, for example, from the contribution of the lensing galaxy at near-IR wavelengths. We show that the near-IR excess technique is indeed an efficient method to identify lensed systems from a large sample of quasars.

We report the discovery of a quadruply imaged quasar, SDSSJ125107.57+293540.5, selected from the Sloan Digital Sky Survey. Follow-up imaging reveals that the system consists of four blue point-like components in a typical cusp lens geometry surrounding a central red galaxy. The source redshift is 0.802 and the lens redshift is 0.410. The maximum image separation between the lensed components is 1"79. While the image configuration is well reproduced by standard mass models with reasonable parameter values, the flux ratios predicted by these models differ from the observed ratios in all bands. This is suggestive of small-scale structures in this lens, although the definitive identification of the anomaly requires more accurate photometry and astrometry.

We use high-resolution N-body numerical simulations to study the number of predicted large-separation multiply-imaged systems produced by clusters of galaxies in the SDSS photometric and spectroscopic quasar samples. We incorporate the condensation of baryons at the centre of clusters by (artificially) adding a brightest central galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cosmological models: a LCDM model with a matter density $\Omega_0=0.3$, and $\sigma_8=0.9$ (LCDM0), and a model favoured by the WMAP three-year data with $\Omega_0=0.238$, and $\sigma_8=0.74$ (WMAP3). We found that the predicted multiply-imaged quasars with separation >10" is about 6.2 and 2.6 for the SDSS photometric (with an effective area 8000 deg$^2$) and spectroscopic (with an effective area 5000 deg$^2$) quasar samples respectively in the LCDM0 model; the predicted numbers of large-separation lensed quasars agree well with the observations. These numbers are reduced by a factor of 7 or more in the WMAP3 model, and are consistent with data at < 8% level. The predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.3 and 1. We find that the BCG creates a central circular region, comparable to the Einstein ring of the BCG, where the central image disappears in the usual three-image and five-image configurations. If we include four image systems as an extreme case of five-image systems (with an infinitely demagnified central image), we find that 68% of the central images are fainter by a factor of 100 than the brightest image, and about 80% are within 1.5"of the BCG.

Most of optical gravitational lenses recently discovered in the Sloan Digital Sky Survey Quasar Lens Search (SQLS) have two-images rather than four-images, in marked contrast to radio lenses for which the fraction of four-image lenses (quad fraction) is quite high. We revisit the quad fraction among optical lenses by taking the selection function of the SQLS into account. We find that the current observed quad fraction in the SQLS is indeed lower than, but consistent with, the prediction of our theoretical model. The low quad fraction among optical lenses, together with the high quad fraction among radio lenses, implies that the quasar optical luminosity function has a relatively shallow faint end slope.

(Abridged) We study the two-point correlation function of a uniformly selected sample of 4,426 luminous optical quasars with redshift $2.9 \le z\le 5.4$ selected over 4041 deg$^2$ from the Fifth Data Release of the Sloan Digital Sky Survey. For a real-space correlation function of the form $\xi(r)=(r/r_0)^{-\gamma}$, the fitted parameters in comoving coordinates are $r_0 = 15.2 \pm 2.7 h^{-1}$ Mpc and $\gamma = 2.0 \pm 0.3$, over a scale range $4\le r_p\le 150 h^{-1}$ Mpc. Thus high-redshift quasars are appreciably more strongly clustered than their $z \approx 1.5$ counterparts, which have a comoving clustering length $r_0 \approx 6.5 h^{-1}$ Mpc. Dividing our sample into two redshift bins: $2.9\le z\le 3.5$ and $z\ge 3.5$, and assuming a power-law index $\gamma=2.0$, we find a correlation length of $r_0 = 16.9 \pm 1.7 h^{-1}$ Mpc for the former, and $r_0 = 24.3 \pm 2.4 h^{-1}$ Mpc for the latter. Following Martini & Weinberg, we relate the clustering strength and quasar number density to the quasar lifetimes and duty cycle. Using the Sheth & Tormen halo mass function, the quasar lifetime is estimated to lie in the range $4\sim 50$ Myr for quasars with $2.9\le z\le 3.5$; and $30\sim 600$ Myr for quasars with $z\ge 3.5$. The corresponding duty cycles are $0.004\sim 0.05$ for the lower redshift bin and $0.03\sim 0.6$ for the higher redshift bin. The minimum mass of halos in which these quasars reside is $2-3\times 10^{12}\ h^{-1}M_\odot$ for quasars with $2.9\le z\le 3.5$ and $4-6\times 10^{12}\ h^{-1}M_\odot$ for quasars with $z\ge 3.5$.

Time delays between lensed multiple images have been known to provide an interesting probe of the Hubble constant, but such application is often limited by degeneracies with the shape of lens potentials. We propose a new statistical approach to examine the dependence of time delays on the complexity of lens potentials, such as higher-order perturbations, non-isothermality, and substructures. Specifically, we introduce a reduced time delay and explore its behavior as a function of the image configuration that is characterized by the asymmetry and opening angle of the image pair. In particular we derive a realistic conditional probability distribution. We find that the probability distribution is sensitive to the image configuration such that more asymmetric and/or smaller opening angle image pairs are more easily affected by perturbations on the primary lens potential. On average time delays of double lenses are less scattered than those of quadruple lenses. Furthermore, the realistic conditional distribution allows a new statistical method to constrain the Hubble constant. We find that 16 published time delay quasars constrain the Hubble constant to be H_0=70+/-6 km/s/Mpc, where the value and its error are estimated using jackknife resampling. After including rough estimates of the sizes of important systematic errors, we find H_0=68+/-6(stat.)+/-8(syst.) km/s/Mpc. The reasonable agreement with other estimates indicates the usefulness of our new approach as a cosmological and astrophysical probe, particularly in the era of large-scale synoptic surveys. (Abridged)

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z_s=2.197. The image separation of 22.5" makes it the largest separation lensed quasar discovered to date. The similarity of the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. Images of the field show a cluster of galaxies at z_l~0.55 that is responsible for the large image separation. The lensed images and the cluster light center are not collinear, which implies that the lensing cluster has a complex structure.

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

We report the discovery of the two-image gravitationally lensed quasar SDSS J133222.62+034739.9 (SDSS J1332+0347) with an image separation of Delta_theta=1.14". This system consists of a source quasar at z_s=1.445 and a lens galaxy at z_l=0.191. The agreement of the luminosity, ellipticity and position angle of the lens galaxy with those expected from lens model confirms the lensing hypothesis.

We report the discoveries of the two-image gravitationally lensed quasars, SDSS J0746+4403 and SDSS J1406+6126, selected from the Sloan Digital Sky Survey (SDSS). SDSS J0746+4403, which will be included in our lens sample for statistics and cosmology, has a source redshift of z_s=2.00, an estimated lens redshift of z_l~0.3, and an image separation of 1.08". SDSS J1406+6126 has a source redshift of z_s=2.13, a spectroscopically measured lens redshift of z_l=0.27, and an image separation of 1.98". We find that the two quasar images of SDSS J1406+6126 have different intervening MgII absorption strengths, which are suggestive of large variations of absorbers on kpc scales. The positions and fluxes of both the lensed quasar systems are easily reproduced by simple mass models with reasonable parameter values. These objects bring to 18 the number of lensed quasars that have been discovered from the SDSS data.

We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong lensing events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, at constraining cosmological parameters from statistical lensing, or at determining the inner structure of galaxy clusters through gravitational arcs. We do so by creating two-dimensional gradually smoothed, differently elliptical and asymmetric versions of some numerical models. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~40%, ~10% and ~30% of the total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and asymmetric clusters, even if located at large distances from the cluster centers (~1Mpc/h). Conversely, their effect is less important in highly asymmetric lenses. The morphology, position and flux of individual arcs are strongly affected by the presence of substructures in the clusters. Removing substructures on spatial scales <~50kpc/h, roughly corresponding to mass scales <~5 10^{10}M_\odot/h, alters the image multiplicity of ~35% of the sources used in the simulations and causes position shifts larger than 5'' for ~40% of the arcs longer than 5''. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong gravitational lensing is potentially a powerfull tool to measure the level of substructures and asymmetries in clusters.

Gamma-Ray Bursts (GRBs) offer a potential way to extend the Hubble diagram to very high redshifts and to constrain the nature of dark energy in a way complementary to distant type Ia supernovae. However, gravitational lensing systematically brightens distant GRBs through the magnification bias, in addition to increasing the dispersions of distance measurements. We investigate how the magnification bias limits the cosmological usage of GRBs. We perform Monte-Carlo simulations of Swift GRBs assuming a cosmological constant dominated universe and then constrain the dark energy equation of state neglecting gravitational lens effects. The originally assumed model is recovered with 68% confidence limit even when the dispersion of inferred luminosities is comparable to that of type Ia supernovae. This implies that the bias is not so drastic for Swift GRBs as to change constraints on dark energy and its evolution. However, the precise degree of the bias in cosmological parameter determinations depends strongly on the shape of the luminosity function of GRBs. Therefore, an accurate determination of the shape of the luminosity function is required to remove the effect of gravitational lensing and to obtain an unbiased Hubble diagram.

We present the results of a search for metal absorption lines in the spectra of background QSOs whose sightlines pass close to foreground QSOs. We detect MgII(2796,2803) absorption in Sloan Digital Sky Survey (SDSS) spectra of four z>1.5 QSOs whose lines of sight pass within 26-98 kpc of lower redshift (z~0.5-1.5) QSOs. The 100% [4/4 pairs] detection of MgII in the background QSOs is clearly at odds with the incidence of associated (z_abs ~ z_em) systems -- absorbers which exist towards only a few percent of QSOs. Although the quality of our foreground QSO spectra is not as high as the SDSS data, absorption seen towards one of the background QSOs clearly does not show up at the same strength in the spectrum of the corresponding foreground QSO. This implies that the absorbing gas is distributed inhomogeneously around the QSO, presumably as a direct consequence of the anisotropic emission from the central AGN. We discuss possible origins for the MgII lines, including: absorption by gas from the foreground QSO host galaxy; companion galaxies fuelling the QSO through gravitational interactions; and tidal debris left by galaxy mergers or interactions which initiated the QSO activity. No single explanation is entirely satisfactory, and we may well be seeing a mixture of phenomena.

We present an algorithm for selecting an uniform sample of gravitationally lensed quasar candidates from low-redshift (0.6<z<2.2) quasars brighter than i=19.1 that have been spectroscopically identified in the SDSS. Our algorithm uses morphological and color selections that are intended to identify small- and large-separation lenses, respectively. Our selection algorithm only relies on parameters that the SDSS standard image processing pipeline generates, allowing easy and fast selection of lens candidates. The algorithm has been tested against simulated SDSS images, which adopt distributions of field and quasar parameters taken from the real SDSS data as input. Furthermore, we take differential reddening into account. We find that our selection algorithm is almost complete down to separations of 1'' and flux ratios of 10^-0.5. The algorithm selects both double and quadruple lenses. At a separation of 2'', doubles and quads are selected with similar completeness, and above (below) 2'' the selection of quads is better (worse) than for doubles. Our morphological selection identifies a non-negligible fraction of single quasars: To remove these we fit images of candidates with a model of two point sources and reject those with unusually small image separations and/or large magnitude differences between the two point sources. We estimate the efficiency of our selection algorithm to be at least 8% at image separations smaller than 2'', comparable to that of radio surveys. The efficiency declines as the image separation increases, because of larger contamination from stars. We also present the magnification factor of lensed images as a function of the image separation, which is needed for accurate computation of magnification bias.

We present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system with a maximum image separation of 15". All four bright images of the quasar, as well as resolved X-ray emission originating from the lensing cluster, are clearly detected. The emission from the lensing cluster extends out to approximately 1.5 arcmin. We measure the bolometric X-ray luminosity and temperature of the lensing cluster to be 4.7e44 erg s^-1 and 6.4 keV, consistent with the luminosity-temperature relation for distant clusters. The mass estimated from the X-ray observation shows excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the images possesses an unusually strong Fe Kalpha emission line, both of which are indicative of microlensing.

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We search 149 moderate resolution background quasar spectra, from Gemini, Keck, the MMT, and the SDSS to survey Lyman Limit Systems (LLSs) and Damped Ly-alpha systems (DLAs) in the vicinity of 1.8 < z < 4.0 luminous foreground quasars. A sample of 27 new quasar-absorber pairs is uncovered with column densities, 17.2 < log (N_HI/cm^2) < 20.9, and transverse (proper) distances of 22 kpc/h < R < 1.7 Mpc/h, from the foreground quasars. If they emit isotropically, the implied ionizing photon fluxes are a factor of ~ 5-8000 times larger than the ambient extragalactic UV background over this range of distances. The observed probability of intercepting an absorber is very high for small separations: six out of eight projected sightlines with transverse separations R < 150 kpc/h have an absorber coincident with the foreground quasar, of which four have log N_HI > 10^19. The covering factor of log N_HI > 10^19 absorbers is thus ~ 50 % (4/8) on these small scales, whereas < 2% would have been expected at random. There are many cosmological applications of these new sightlines: they provide laboratories for studying fluorescent Ly-alpha recombination radiation from LLSs, constrain the environments, emission geometry, and radiative histories of quasars, and shed light on the physical nature of LLSs and DLAs.

We have studied the generation of magnetic fields by the Biermann mechanism in the pair-instability supernovae explosions of the first stars. The Biermann mechanism produces magnetic fields in the shocked region between the bubble and interstellar medium (ISM), even if magnetic fields are absent initially. We have performed a series of two-dimensional magnetohydrodynamic simulations with the Biermann term and estimate the amplitude and total energy of the produced magnetic fields. We find that magnetic fields with amplitude $10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of magnetic fields generated depend on specific values of initial conditions. This corresponds to magnetic fields of $10^{28}-10^{31}$ ergs per each supernova remnant, which is strong enough to be the seed magnetic field for a galactic and/or interstellar dynamo.

We investigate the dependence of the time delays for the large-separation gravitationally lensed quasar SDSS J1004+4112 on the inner mass profile of the lensing cluster. Adopting the mass model whose innermost density profile is parameterized as \rho \propto r^{-\alpha}, we derive a series of mass models which can fit observational data and then compute the probability distribution functions of time delays. We find that larger \alpha has longer time delays, longer tails at the higher end of the probability distribution, and larger model uncertainties. The ratios of time delays slightly depend on the slope \alpha. Among others, time delays between images C and A (or B) have little dependence on the inner slope, particularly when the time delays are short. The dependence of time delays on \alpha is well fitted by a linear form, which reflects well-known degeneracy between the mass profile and time delays. We perform a Monte-Carlo simulation to illustrate how well the inner slope can be constrained from measurements of time delays. We find that measurements of more than one time delays result in reasonably tight constraints on the inner slope (\sigma_{\alpha} \lesssim 0.25), while only one time delay cannot determine the inner slope very well. Our result indicates that time delays indeed serve as a powerful tool to determine the mass profile, despite the complexity of the lensing cluster.

We present a halo model prediction of the image separation distribution of strong lenses. Our model takes account of the subhalo population, which has been ignored in previous studies, as well as the conventional halo population. Halos and subhalos are linked to central and satellite galaxies by adopting an universal scaling relation between masses of (sub-)halos and luminosities of galaxies. Our model predicts that 10%-20% of lenses should be caused by the subhalo population. The fraction of lensing by satellite galaxies (subhalos) peaks at ~1'' and decreases rapidly with increasing image separations. We compute fractions of lenses which lie in groups and clusters, and find them to be ~14% and ~4%, respectively: Nearly half of such lenses are expected to be produced by satellite galaxies, rather than central parts of halos. We also study mass distributions of lensing halos and find that even at image separations of ~3'' the deviation of lens mass distributions from isothermal profiles is large: At or beyond ~3'' image separations are enhanced significantly by surrounding halos. Our model prediction agrees reasonably well with observed image separation distributions from galaxy to cluster scales.

We report the discoveries of two, two-image gravitationally lensed quasars selected from the Sloan Digital Sky Survey: SDSS J0806+2006 at z_s=1.540 and SDSS J1353+1138 at z_s=1.629 with image separations of 1.40" and 1.41" respectively. Spectroscopic and optical/near-infrared imaging follow-up observations show that the quasar images have identical redshifts and possess extended objects between the images that are likely to be lens galaxies at z_l~0.6 in SDSS J0806+2006 and z_l~0.3 in SDSS J1353+1138. The field of SDSS J0806+2006 contains several nearby galaxies that may significantly perturb the system, and SDSS J1353+1138 has an extra component near its Einstein ring that is probably a foreground star. Simple mass models with reasonable parameters reproduce the quasar positions and fluxes of both systems.

The deep imaging of the field of an observed lensing event by a cosmic string reveals many additional lensing events. We study the statistics of such string lensing. We derive explicit expressions for the distributions of image separations of lensing by a cosmic string and point out that they are quite sensitive to parameters which characterize the cosmic string, such as the redshift and tension of the cosmic string. Thus the statistics of string lensing events add new important information on the cosmic string which cannot be obtained from the detailed investigation of one lensing event.

We study the impact of lens galaxy environments on the image separation distribution of lensed quasars. We account for both environmental convergence and shear, using a joint distribution derived from galaxy formation models calibrated by galaxy-galaxy lensing data and number counts of massive elliptical galaxies. We find that the external field enhances lensing probabilities, particularly at large image separations; the increase is ~30% at \theta=3'' and ~200% at \theta=5'', when we adopt a power-law source luminosity function \Phi(L) \propto L^-2.1. The enhancement is mainly driven by convergence, which boosts both the image separation and magnification bias (for a fixed lens galaxy mass). These effects have been neglected in previous studies of lens statistics. Turning the problem around, we derive the posterior convergence and shear distributions and point out that they are strong functions of image separation; lens systems with larger image separations are more likely to lie in dense environments.

We report follow-up observations of two gravitational lens candidates identified in the Sloan Digital Sky Survey (SDSS) dataset. We have confirmed that SDSS J102111.02+491330.4 is a previously unknown gravitationally lensed quasar. This lens system exhibits two images of a $z = 1.72$ quasar, with an image separation of $1{\farcs}14 \pm 0.04$. Optical and near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Observations of SDSS J112012.12+671116.0 indicate that it is more likely a binary quasar than a gravitational lens. This system has two quasars at a redshift of $z = 1.49$, with an angular separation of $1{\farcs}49 \pm 0.02$. However, the two quasars have markedly different SEDs and no lens galaxy is apparent in optical and near-IR images of this system. We also present a list of 31 SDSS lens candidates which follow-up observations have confirmed are \textit{not} gravitational lenses.

We have identified three multiply imaged galaxies in Hubble Space Telescope images of the redshift z=0.68 cluster responsible for the large-separation quadruply lensed quasar, SDSS J1004+4112. Spectroscopic redshifts have been secured for two of these systems using the Keck I 10m telescope. The most distant lensed galaxy, at z=3.332, forms at least four images, and an Einstein ring encompassing 3.1 times more area than the Einstein ring of the lensed QSO images at z=1.74, due to the greater source distance. For a second multiply imaged galaxy, we identify Ly_alpha emission at a redshift of z=2.74. The cluster mass profile can be constrained from near the center of the brightest cluster galaxy, where we observe both a radial arc and the fifth image of the lensed quasar, to the Einstein radius of the highest redshift galaxy, ~110 kpc. Our preliminary modeling indicates that the mass approximates an elliptical body, with an average projected logarithmic gradient of ~-0.5. The system is potentially useful for a direct measurement of world models in a previously untested redshift range.

The steep mass profile of A1689 derived from recent detailed lensing observations is not readily reconciled with the low concentration halos predicted by the standard CDM model. However, halo triaxiality may act to bias the profile constraints derived assuming a spherically symmetric mass distribution, since lensing relates only to the projected mass distribution. A degree of halo triaxiality is inherent to the CDM structure formation, arising from the collision-less nature of the dark matter. Here we compare the CDM-based model predictions of triaxial halo with the precise lensing measurements of A1689 based on the ACS/HST and Subaru data, over a wide range of 10kpc<r<2Mpc. The model lensing profiles cover the intrinsic spread of halo mass and shape, and are projected over all inclinations when comparing with the data. We show that the model parameters are only weakly constrained and strongly degenerate mainly because of the lack of information along the line of sight. In particular, the limits on the concentration parameter become less restrictive with increasing triaxiality. Yet, by comparing the obtained constraints with expected probability distributions for the axis ratio and concentration parameters computed from numerical simulations, we find that ~6% of cluster-size halos in the CDM model can match the A1689 lensing observations at the 2-sigma level, corresponding to cases where the major-axis of the halo is aligned with the line of sight. Thus halo triaxiality could reduce the apparent discrepancy between theory and observation. This hypothesis needs to be further explored by a statistical lensing study for other clusters as well as by complementary three-dimensional information derived using X-ray, kinematics, and SZ effect observations.

We study generation of magnetic fields by the Biermann mechanism in the pair-instability supernovae explosions of first stars. The Biermann mechanism produces magnetic fields in the shocked region between the bubble and interstellar medium (ISM), even if magnetic fields are absent initially. We perform a series of two-dimensional magnetohydrodynamic simulations with the Biermann term and estimate the amplitude and total energy of the produced magnetic fields. We find that magnetic fields with amplitude $10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of magnetic fields generated depend on specific values of initial conditions. This corresponds to magnetic fields of $10^{28}-10^{31}$ erg per each supernova remnant, which is strong enough to be the seed magnetic field for galactic and/or interstellar dynamo.

We report the discovery of a new two-image gravitationally lensed quasar, SDSS J024634.11-082536.2 (SDSS J0246-0825). This object was selected as a lensed quasar candidate from the Sloan Digital Sky Survey (SDSS) by the same algorithm that was used to discover other SDSS lensed quasars (e.g., SDSS J0924+0219). Multicolor imaging with the Magellan Consortium's Walter Baade 6.5-m telescope and the spectroscopic observations using the W. M. Keck Observatory's Keck II telescope confirm that SDSS J0246-0825 consists of two lensed images ($\Delta{\theta}=$1\farcs04) of a source quasar at z=1.68. Imaging observations with the Keck telescope and the Hubble Space Telescope reveal an extended object between the two quasar components, which is likely to be a lensing galaxy of this system. From the absorption lines in the spectra of quasar components and the apparent magnitude of the galaxy, combined with the expected absolute magnitude from the Faber-Jackson relation, we estimate the redshift of the lensing galaxy to be z=0.724. A highly distorted ring is visible in the Hubble Space Telescope images, which is likely to be the lensed host galaxy of the source quasar. Simple mass modeling predicts the possibility that there is a small (faint) lensing object near the primary lensing galaxy.

The fraction of substructures required to account for anomalous flux ratios in gravitational lens systems appears to be higher than that predicted in the standard cold dark matter cosmology. We present a possible alternative route to anomalous flux ratios from lens galaxy environments. We consider the compound lens system such that a lens galaxy lie in a group or cluster, and estimate the contribution of substructures in the group/cluster to the fraction using an analytic model of substructures. We find that the contribution becomes dominant when the impact parameter of the lens is less than ~30% of the virial radius of the group/cluster. This indicates that the environmental effect can partly explain the high incidence of anomalous flux ratios.

We present a sample of 218 new quasar pairs with proper transverse separations R_prop < 1 Mpc/h over the redshift range 0.5 < z < 3.0, discovered from an extensive follow up campaign to find companions around the Sloan Digital Sky Survey and 2dF Quasar Redshift Survey quasars. This sample includes 26 new binary quasars with separations R_prop < 50 kpc/h (theta < 10 arcseconds), more than doubling the number of such systems known. We define a statistical sample of binaries selected with homogeneous criteria and compute its selection function, taking into account sources of incompleteness. The first measurement of the quasar correlation function on scales 10 kpc/h < R_prop < 400 kpc/h is presented. For R_prop < 40 kpc/h, we detect an order of magnitude excess clustering over the expectation from the large scale R_prop > 3 Mpc/h quasar correlation function, extrapolated down as a power law to the separations probed by our binaries. The excess grows to ~ 30 at R_prop ~ 10 kpc/h, and provides compelling evidence that the quasar autocorrelation function gets progressively steeper on sub-Mpc scales. This small scale excess can likely be attributed to dissipative interaction events which trigger quasar activity in rich environments. Recent small scale measurements of galaxy clustering and quasar-galaxy clustering are reviewed and discussed in relation to our measurement of small scale quasar clustering.

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

We report the discovery of two doubly-imaged quasars, SDSS J100128.61+502756.9 and SDSS J120629.65+433217.6, at redshifts of 1.838 and 1.789 and with image separations of 2.86'' and 2.90'', respectively. The objects were selected as lens candidates from the Sloan Digital Sky Survey (SDSS). Based on the identical nature of the spectra of the two quasars in each pair and the identification of the lens galaxies, we conclude that the objects are gravitational lenses. The lenses are complicated; in both systems there are several galaxies in the fields very close to the quasars, in addition to the lens galaxies themselves. The lens modeling implies that these nearby galaxies contribute significantly to the lens potentials. On larger scales, we have detected an enhancement in the galaxy density near SDSS J100128.61+502756.9. The number of lenses with image separation of ~3'' in the SDSS already exceeds the prediction of simple theoretical models based on the standard Lambda-dominated cosmology and observed velocity function of galaxies.

We construct an analytic model for the mass distributions of superclusters and clusters in each supercluster. Our model is a modification of the Press-Schechter theory, and defines superclusters as the regions that have some overdensity smaller than those for usual virialized objects. We compare the mass functions with a catalog of superclusters in the Sloan Digital Sky Survey Early Data Release and found that they are in reasonable agreement with each other.

We present theoretical and observational studies of strong gravitational lenses produced by clusters of galaxies. Our purpose is to test the Cold Dark Matter (CDM) model at small and highly non-linear scales where it has been claimed that the CDM model may confront several difficulties. We concentrate our attention on the statistics of strong gravitational lenses. We use two complementary statistics, lensed arcs and quasars, to probe the mass distributions. First, we construct a triaxial lens model, and develop a new method to include triaxiality of dark halos in the lens statistics. We find that the effect of triaxiality is significant; it enhances lensing probabilities by a factors of a few to ten, assuming the degree of triaxiality predicted in the CDM model. In particular, we argue that both central concentration and large triaxiality of dark halos are required to reproduce the observed number of arcs in clusters; thus the result can be interpreted as a strong evidence for the cold and collisionless dark matter. One of the most notable advantages of the triaxial modeling is that the triaxial modeling allows us to predict image multiplicities. We find that the CDM halos predict significant fraction (more than 20%) of naked cusp lenses. In addition, we point out the image multiplicities depend strongly on the central concentration of dark halos. Therefore we propose image multiplicities as a new powerful test of the CDM model. Finally we searched for large-separation lensed quasars from the data of the Sloan Digital Sky Survey, and succeeded in discovering the first large-separation lensed quasar SDSS J1004+4112. We find that the discovery of the large-separation quadruple lens SDSS J1004+4112 is quite consistent with the theoretical predictions based on the CDM model. (Abridged)

We construct a realistic model for the dark halo substructures, and derive analytically their spatial, mass, and velocity distributions. Basically, our model is a modification of the Press-Schechter theory to account for the dominant dynamical processes that mark the evolution of dark halo substructures such as the tidal stripping and dynamical friction. Our analytic model successfully reproduces all the well known behaviors of the substructure distributions that have been found in recent numerical simulations: the weak dependence of the mass distributions on the host halo mass; the anti-bias of the spatial distribution relative to the dark matter particle components; the nearly power-law shapes of the mass and velocity distributions. We compare our analytic results with recent high-resolution N-body simulation data and find that they are in excellent agreement with each other.

We re-formulate cosmological perturbations in the decaying cold dark matter model, and calculate cosmological microwave background (CMB) anisotropies. By comparing our theoretical predictions with recent observational data from the Wilkinson Microwave Anisotropy Probe(WMAP), we derive a new bound on the abundance and lifetime of decaying dark matter particles. We show that the data of WMAP alone do not prefer the decaying cold dark matter model: the lifetime is constrained to \Gamma^{-1}> 123 Gyr at 68% C.L. (52 Gyr at 95.4% C.L.) when cold dark matter consists only of such decaying particles. We also consider a more general case in which cold dark matter consists of both stable and decaying particles, and show that the constraint generalizes down to \Omega_{DDM}h^2 < -0.5 (Gamma^{-1}/1 Gyr)^{-1}+0.12 for \Gamma^{-1} > 5 Gyr at 95.4% C.L.. These bounds are robust and widely applicable, because they are derived only from gravitational effects and therefore do not rely on the details of decay channels or decay products.

We study the statistics of large-separation gravitational lens systems produced by non-spherical halos in the Cold Dark Matter (CDM) model. Specifically, we examine how the triaxiality of CDM halos affects the overall lensing probabilities and the relative numbers of different image configurations (double, quadruple, and naked cusp lenses). We find that triaxiality significantly enhances lensing probabilities by a factor of ~2-4, so it cannot be ignored. If CDM halos have central density slopes alpha<~1.5, we predict that a significant fraction (>~20%) of large-separation lenses should have naked cusp image configurations; this contrasts with lensing by isothermal (\alpha~2) galaxies where naked cusp configurations are rare. The image multiplicities depend strongly on the inner density slope \alpha: for \alpha=1, the naked cusp fraction is >~60%; while for \alpha=1.5, quadruple lenses are actually the most probable. Thus, the image multiplicities in large-separation lenses offer a simple new probe of the inner density profiles of dark matter halos. We also compute the expected probabilities and image multiplicities for lensed quasars in the Sloan Digital Sky Survey, and argue that the recent discovery of the large-separation quadruple lens SDSS J1004+4112 is consistent with expectations for CDM.

In modeling strong gravitational lens systems, one often adopts simple models, such as singular isothermal elliptical density plus lowest-order external perturbation. However, such simple models may mislead us if the real mass distribution is more complicated than that in the assumed models. In particular, assumptions on mass models are crucial in studying flux ratio anomalies that have been suggested as evidence for a cold dark-matter substructure. We reinvestigated four quadruple lens systems using power-law Fourier models, which have advantages of clear physical meanings and the applicability of a linear method, as well as a simple singular isothermal elliptical density model. We also investigated the effect of external perturbations, including a singular isothermal sphere, lowest-order expansion, and next-order expansion. We have found that the $\cos3\theta$ terms of the primary galaxy and/or of external perturbation significantly reduce the \chi^2 in PG 1115+080 and B1422+231. In particular, we could reproduce the flux ratios of B1422+231 with next-order external perturbation assuming 5% flux uncertainties, suggesting that external perturbation cannot be described by a simple singular isothermal sphere approximation. On the other hand, we could not fit B0712+472 and B2045+265 very well even with our models, although the \chi^2 were reduced compared with the case of using simple models. Our results clearly demonstrate that both the primary lens galaxy and the external perturbation are often more complicated than we usually assume.

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

Statistics of lensed arcs in clusters of galaxies serve as a powerful probe of both the non-sphericity and the inner slope of dark matter halos. We develop a semi-analytic method to compute the number of arcs in triaxial dark matter halos. This combines the lensing cross section from the Monte Carlo ray-tracing simulations, and the probability distribution function (PDF) of the axis ratios evaluated from cosmological N-body simulations. This approach enables one to incorporate both asymmetries in the projected mass density and elongations along the line-of-sight analytically, for the first time in cosmological lensed arc statistics. As expected, triaxial dark matter halos significantly increase the number of arcs relative to spherical models; the difference amounts to more than one order of magnitude while the value of enhancement depends on the specific properties of density profiles. Then we compare our theoretical predictions with the observed number of arcs from 38 X-ray selected clusters. In contrast to the previous claims, our triaxial dark matter halos with inner density profile \rho \propto r^{-1.5} in a Lambda-dominated cold dark matter (CDM) universe reproduces well the observation. Since both the central mass concentration and triaxial axis ratios (minor to major axis ratio ~0.5) required to account for the observed data are consistent with cosmological N-body simulations, our result may be interpreted to lend strong support for the collisionless CDM paradigm at the mass scale of clusters.

We report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow-up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z=1.57. The image separation is 1.56''. We also detect an extended component between the two quasar images, likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay \Delta t ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

Recent studies of structures of galaxies and clusters imply that dark matter might be unstable and decay with lifetime $\Gamma^{-1}$ about the age of universe. We study the effects of the decay of cold dark matter on the Sunyaev-Zel'dovich (SZ) power spectrum. We analytically calculate the SZ power spectrum taking finite lifetime of cold dark matter into account. We find the finite lifetime of dark matter decreases the power at large scale ($l < 4000$) and increases at small scale ($l > 4000$). This is in marked contrast with the dependence of other cosmological parameters such as the amplitude of mass fluctuations $\sigma_{8}$ and the cosmological constant $\Omega_{\lambda}$ (under the assumption of a flat universe) which mainly change the normalization of the angular power spectrum. This difference allows one to determine the lifetime and other cosmological parameters rather separately. We also investigate sensitivity of a future SZ survey to the cosmological parameters including the life time, assuming a fiducial model $\Gamma^{-1} = 10 h^{-1} {\rm Gyr}$, $\sigma_{8} = 1.0$, and $\Omega_{\lambda} = 0.7$. We show that future SZ surveys such as ACT, AMIBA, and BOLOCAM can determine the lifetime within factor of two even if $\sigma_{8}$ and $\Omega_{\lambda}$ are marginalized.

The cluster abundance and its redshift evolution are known to be a powerful tool to constrain the amplitude of mass fluctuations \sigma_8 and the mass density parameter \Omega_m0. We study the impact of the finite decay rate of cold dark matter particles on the cluster abundances. On the basis of spherical model in decaying cold dark matter universe, we calculate the mass function of clusters and compare it with observed cluster abundance. We find the decay of cold dark matter particles significantly changes the evolution of the cluster abundance. In particular, we point out that the lifetime of dark matter particles comparable to the age of the universe lowers the ratio of the local cluster abundance to the high-redshift cluster abundance and can account for the observed evolution of the cluster abundance quite well. The strong dependence of the cluster abundance on the decay rate of dark matter suggests that distant cluster surveys may offer clues to the nature of dark matter.

We propose a new method to use gamma-ray bursts (GRBs) as an alternative probe of the dark energy. By calibrating luminosity-variability and luminosity-lag time relations at low redshift where distance-redshift relations have been already determined from type Ia supernova (SNIa), GRBs at high redshift can be used as a distance indicator. We investigate the potential impact of future GRB data on determining the current matter density $\Omega_{m}$ and the dark energy equation of state $w$ which is assumed to be constant. We show that a combined analysis of a data set expected from the {\it Swift} and the current SNIa data results in excellent determination of both $\Omega_{m}$ and $w$.

Wide separation lensing statistics offer information about the density profile and abundance of dark halos. Recently a possible discovery of six quasar pairs, which may be lensed multiple images, was reported by Miller et al. (2003, astro-ph/0210644). These pairs are selected from a catalog of the 2dF quasar (QSO) survey comprising 22163 quasars. We calculate expected numbers of lensed quasars taking account of the redshift and magnitude distributions of the quasar catalog. Given some of the six pairs are genuine lensed systems, we put interesting constraints on the inner slope of dark halos, \Omega_0, and \sigma_8. We show that the detection of even one lens with separation >30'' is marginally consistent with models that have cuspy inner density profile, \rho \propto r^{-1.5}, and very large \sigma_8, \sigma_8>~1.2 for \Omega_0=0.3. To reconcile with constraints from X-ray clusters or cosmic shear, much lower \Omega_0 and much higher \sigma_8 are needed, although such high \sigma_8 seems too extreme. The shallower inner density profile \rho \propto r^{-1} is hardly acceptable. In particular, the expected number of lenses with separation >200'' is too small to explain the discovery of such anomalously wide separation lens systems. These results imply that we miss some important systematic effects, there is a problem in the cold dark matter scenario, or none of these six quasar pairs is likely to be lensed images.

An attempt to measure the Hubble constant with gravitational lens time delays is often limited by the strong degeneracy between radial mass profiles of lens galaxies and the Hubble constant. We show that strong gravitational lensing of type Ia supernovae breaks this degeneracy; the standard candle nature of type Ia supernova luminosity function allows one to measure the magnification factor directly, and this information is essential to constrain radial mass profiles and the Hubble constant separately. Our numerical simulation demonstrates that the Hubble constant can be determined with \sim 5% accuracy from only several lens events if magnification factors are used as constraints. Therefore, distant supernova survey is a promising way to measure the global Hubble constant independently with the local estimates.

Strong gravitational lensing of distant supernovae (SNe), particularly Type Ia's, has some exploitable properties not available when other sorts of cosmologically distant sources are lensed. One such property is that the ``standard candle'' nature of the SN at peak brightness allows a direct determination of the lensing magnification factor for each well observed image. Another is that the duration of an SN event is of the same order as the differential time delays between the various lens images. The comparable time scales of the photometric event and the time delay invite a variety of applications, including high precision measurements of the delay and the targeting of especially interesting phases of the explosion for intensive observation when they appear in trailing images. As an initial exploration of these possibilities we present calculations of SN lensing statistics in a ``concordance cosmology''. We emphasize magnification and time delay effects. Plausible SN surveys, such as the proposed SNAP space mission, would discover several to some tens of strongly lensed SNe Ia per year, and at least a few of these will be at redshifts well beyond those that would be accessible via unlensed events. The total number of such anomalously high redshift SNe Ia will be a useful test of high redshift star formation models. SN surveys of finite duration will, of course, miss the appearance of some images; we quantify this selection bias. Finally, we investigate how well the appearance of trailing images can be predicted based on various amounts of available information on the lensing event. Knowledge of the magnification factor for the leading image makes it possible to predict the appearance of a trailing image relatively accurately if the lens redshift is also known. (Abridged)

We predict the fraction of dark halo lenses, that is, the fraction of lens systems produced by the gravitational potential of dark halos, on the basis of a simple parametric model of baryonic compression. The fraction of dark halo lenses primarily contains information on the effect of baryonic compression and the density profile of dark halos, and is expected to be insensitive to cosmological parameters and source population. The model we adopt comprises the galaxy formation probability p_g(M) which describes the global efficiency of baryonic compression and the ratio of circular velocities of galaxies to virial velocities of dark halos gamma_v=v_c/v_{vir} which means how the inner structure of dark halos is modified due to baryonic compression. The model parameters are constrained from the velocity function of galaxies and the distribution of image separations in gravitational lensing, although the degeneracy between model parameters still remains. We show that the fraction of dark halo lenses depends strongly on gamma_v and the density profile of dark halos such as inner slope alpha. This means that the observation of the fraction of dark halos can break the degeneracy between model parameters if the density profile of dark halo lenses is fully settled. On the other hand, by restricting gamma_v to physically plausible range we can predict the lower limit of the fraction of dark halo lenses on the basis of our model. Our result indicates that steeper inner cusps of dark halos (alpha >~ 1.5) or too centrally concentrated dark halos are inconsistent with the lack of dark halo lenses in observations.

It has been recognized that the arc statistics of gravitational lensing is a useful probe of the density profile of clusters of galaxies. We examine several systematic effects which are important in predicting the number of arcs, with particular attention to the difference between tangential and radial arcs. First we derive an analytic expression of the cross section for radial arcs taking account of the source size and find that the moderate source size enhances the cross section for radial arcs while larger source size (\gtrsim 1'' in our example) suppresses the number of radial arcs. On the other hand, tangential arcs are much less sensitive to the source size. Next we numerically calculate the cross section for arcs considering the lens and source ellipticities. We find that the numbers of both tangential and radial arcs are highly enhanced by both ellipticities, by one or two orders of magnitude. The number ratio of radial to tangential arcs is, however, not so affected if the threshold axis ratio of arcs is large (\gtrsim 7). The number ratio therefore still remains good statistics which probe the density profile of the lens objects, if the source size effect is correctly taken into account.

The distribution of differential time delays \Delta t between images produced by strong gravitational lensing contains information on the mass distributions in the lensing objects as well as on cosmological parameters such as H_0. We derive an explicit expression for the conditional probability distribution function of time delays P(\Delta t|\theta), given an image separation between multiple images \theta, and related statistics. We consider lensing halos described by the singular isothermal sphere (SIS) approximation and by its generalization as proposed by Navarro, Frenk, & White (NFW) which has a density profile \rho \propto r^{-\alpha} in the innermost region. The time delay distribution is very sensitive to these profiles; steeper inner slopes tend to produce larger time delays. For example, if H_0=70 km/s/Mpc, a \Lambda-dominated cosmology and a source redshift z_S=1.27 are assumed, lenses with \theta=5'' produce a time delay of \Delta t[yr]=1.5^{+1.7}_{-0.9}, 0.39^{+0.37}_{-0.22}, 0.15^{+0.11}_{-0.09}, and 0.071^{+0.054}_{-0.038} (50% confidence interval), for SIS, generalized NFW with \alpha=1.5, \alpha=1.0, and \alpha=0.5, respectively. At a fixed image separation, the time delay is determined by the difference in the lensing potential between the position of the two images, which typically occur at different impact parameters. Although the values of \Delta t are proportional to the inverse of H_0, P(\Delta t|\theta) is rather insensitive to all other cosmological model parameters, source redshifts, magnification biases and so on. A knowledge of P(\Delta t|\theta) will also be useful in designing the observing program of future large scale synoptic variability surveys and for evaluating possible selection biases operating against large splitting lens systems.

We study the arc statistics of gravitational lensing generated by dark matter halos in order to probe their density profile. We characterize the halo profile by two parameters, the inner slope of the central cusp $\alpha$, and the median amplitude of the concentration parameter, $c_{\rm norm}$, for a halo of mass $10^{14}h^{-1}M_\odot$ at $z=0$, and compute the numbers of tangential and radial arcs produced by gravitational lensing of galaxy clusters. We find that the number of arcs divided by the number of halos is a good statistic which is sensitive to both $c_{\rm norm}$ and $\alpha$ with very weak dependence on the cosmological parameters. If the arc samples with well-defined selection criteria for the clusters become available, one can strongly constrain both $c_{\rm norm}$ and $\alpha$. While our tentative comparison with the existing observational data indicates that the inner density profile of dark halos is indeed as steep as predicted by recent simulations ($\alpha \sim 1.5$), the homogeneous samples of tangential and radial arcs are required for more quantitative discussions.