Source author record

Masamune Oguri

Masamune Oguri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

128works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

128 published item(s)

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2025arXiv

The SRG/eROSITA All-Sky Survey. The Weak-Lensing Mass Calibration and the Stellar Mass-to-Halo Mass Relation from the Hyper Suprime-Cam Subaru Strategic Program

We present the weak-lensing mass calibration and constrain the BCG (brightest cluster galaxy) stellar-mass-to-halo-mass-and-redshift ($M_{\star,\mathrm{BCG}}-M-z$) relation for a sample of $124$ galaxy clusters and groups at redshift $0.1<z<0.8$ from the first Data Release of the $eROSITA$ All-Sky Survey (eRASS1), using data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The cluster survey is conducted by the $eROSITA$ X-ray telescope aboard the Spectrum-Roentgen-Gamma (SRG) space observatory. The cluster sample is X-ray-selected and optically confirmed with a negligibly low contamination rate ($\approx5%$). On a basis of individual clusters, the shear profiles of $96$ clusters are derived using the HSC Three-Year (HSC-Y3) weak-lensing data, while the BCG stellar masses of $101$ clusters are estimated using the SED template fitting to the HSC five-band ($grizY$) photometry. The observed X-ray photon count rate is used as the mass proxy, based on which individual halo masses are obtained at the given count rate in a population modelling while accounting for systematic uncertainties in the weak-lensing modelling through a simulation-calibrated weak-lensing mass-to-halo-mass relation. The count rate and BCG stellar mass relations are simultaneously constrained in a forward and population modelling. In agreement with the results based on the weak-lensing data from the DES and KiDS surveys, we obtain a count rate relation with a self-similar redshift scaling and a mass trend that is steeper than the self-similar prediction. Our results suggest that the BCG stellar mass at a fixed halo mass has remained stable with a moderate increase at a level of $\left(20\pm8\right)%$ since redshift $z\approx0.8$. This finding supports the picture of the ``rapid-then-slow'' BCG formation, where the majority of the stellar mass must have been assembled at much earlier cosmic time.

preprint2022arXiv

ALMA Lensing Cluster Survey: ALMA-Herschel Joint Study of Lensed Dusty Star-Forming Galaxies across $z\simeq0.5-6$

We present an ALMA-Herschel joint analysis of sources detected by the ALMA Lensing Cluster Survey (ALCS) at 1.15 mm. Herschel/PACS and SPIRE data at 100-500 $μ$m are deblended for 180 ALMA sources in 33 lensing cluster fields that are either detected securely (141 sources; in our main sample) or tentatively at S/N$\geq$4 with cross-matched HST/Spitzer counterparts, down to a delensed 1.15-mm flux density of $\sim0.02$ mJy. We performed far-infrared spectral energy distribution modeling and derived the physical properties of dusty star formation for 125 sources (109 independently) that are detected at $>2σ$ in at least one Herschel band. 27 secure ALCS sources are not detected in any Herschel bands, including 17 optical/near-IR-dark sources that likely reside at $z=4.2\pm1.2$. The 16-50-84 percentiles of the redshift distribution are 1.15-2.08-3.59 for ALCS sources in the main sample, suggesting an increasing fraction of $z\simeq1-2$ galaxies among fainter millimeter sources ($f_{1150}\sim 0.1$ mJy). With a median lensing magnification factor of $μ= 2.6_{-0.8}^{+2.6}$, ALCS sources in the main sample exhibit a median intrinsic star-formation rate of $94_{-54}^{+84}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, lower than that of conventional submillimeter galaxies at similar redshifts by a factor of $\sim$3. Our study suggests weak or no redshift evolution of dust temperature with $L_\mathrm{IR}<10^{12}\,\mathrm{L}_\odot$ galaxies within our sample at $z \simeq 0 - 2$. At $L_\mathrm{IR}>10^{12}\,\mathrm{L}_\odot$, the dust temperatures show no evolution across $z \simeq 1 -4$ while being lower than those in the local Universe. For the highest-redshift source in our sample ($z=6.07$), we can rule out an extreme dust temperature ($>$80 K) that was reported for MACS0416 Y1 at $z=8.31$.

preprint2022arXiv

Amplitude and phase fluctuations of gravitational waves magnified by strong gravitational lensing

We discuss how small-scale density perturbations on the Fresnel scale affect amplitudes and phases of gravitational waves that are magnified by gravitational lensing in geometric optics. We derive equations that connect the small-scale density perturbations with the amplitude and phase fluctuations to show that such perturbative wave optics effects are significantly boosted in the presence of macro model magnifications such that amplitude and phase fluctuations can easily be observed for highly magnified gravitational waves. We discuss expected signals due to microlensing by stars, dark matter substructure, fuzzy dark matter, and primordial black holes.

preprint2022arXiv

Do blue galaxy-clusters have hot intracluster gas?

We present herein a systematic X-ray analysis of blue galaxy-clusters at $z=0.84$ discovered by the Subaru telescope. The sample consisted of 43 clusters identified by combining red-sequence and blue-cloud surveys, covering a wide range of emitter fractions (i.e., 0.3--0.8). The spatial extent of the over-density region of emitter galaxies was approximately 1~Mpc in radius. The average cluster mass was estimated as $0.6 (<1.5)\times10^{14}~{\rm M_\odot}$ from the stacked weak-lensing measurement. We analyzed the XMM-Newton archival data, and measured the X-ray luminosity of the hot intracluster medium. As a result, diffuse X-ray emission was marginally detected in 14 clusters, yielding an average luminosity of $5\times 10^{42}~{\rm erg\,s^{-1}}$. On the contrary, it was not significant in 29 clusters. The blue clusters were significantly fainter than the red-dominated clusters, and the X-ray luminosity did not show any meaningful correlation either with emitter fraction or richness. The X-ray surface brightness was low, but the amount of gas mass was estimated to be comparable to that observed in the $10^{13-14}~{\rm M_{\odot}}$ cluster. Based on the results, we suggest that the blue clusters are at the early formation stage, and the gas is yet to be compressed and heated up to produce appreciable X-rays. Follow-up spectroscopic measurements are essential to clarify the dynamical status and co-evolution of galaxies and hot gas in the blue clusters.

preprint2022arXiv

HSC-XXL : Baryon budget of the 136 XXL Groups and Clusters

We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution.

preprint2022arXiv

Line-of-sight Elongation and Hydrostatic Mass Bias of the Frontier Fields Galaxy Cluster Abell 370

We present a detailed weak-lensing and X-ray study of the Frontier Fields galaxy cluster Abell 370, one of the most massive known lenses on the sky, using wide-field BRz Subaru/Sprime-Cam and Chandra X-ray observations. By combining two-dimensional (2D) shear and azimuthally averaged magnification constraints derived from Subaru data, we perform a lensing mass reconstruction in a free-form manner, which allows us to determine both radial structure and 2D morphology of the cluster mass distribution. In a triaxial framework assuming a Navarro-Frenk-White density profile, we constrain the intrinsic structure and geometry of the cluster halo by forward modeling the reconstructed mass map. We obtain a halo mass $M_{200}=(1.54 \pm 0.29)\times 10^{15}h^{-1}M_\odot$, a halo concentration $c_{200}=5.27 \pm 1.28$, and a minor-major axis ratio $q_a=0.62 \pm 0.23$ with uninformative priors. Using a prior on the line-of-sight alignment of the halo major axis derived from binary merger simulations constrained by multi-probe observations, we find that the data favor a more prolate geometry with lower mass and lower concentration. From triaxial lens modeling with the line-of-sight prior, we find a spherically enclosed gas mass fraction of $f_\mathrm{gas}=(8.4 \pm 1.0)\%$ at $0.7h^{-1}$ Mpc $\sim 0.7r_{500}$. When compared to the hydrostatic mass estimate from Chandra observations, our triaxial weak-lensing analysis yields spherically enclosed mass ratios of $1-b \equiv M_\mathrm{HE}/M_\mathrm{WL} = 0.56 \pm 0.09$ and $0.51 \pm 0.09$ at $0.7h^{-1}$ Mpc with and without using the line-of-sight prior, respectively. Since the cluster is in a highly disturbed dynamical state, this represents the likely maximum level of hydrostatic bias in galaxy clusters.

preprint2022arXiv

Star Formation Properties of Sloan Digital Sky Survey BOSS Void Galaxies in the Hyper Suprime-Cam Survey

We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$σ$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias.

preprint2022arXiv

Survey of Gravitationally-lensed Objects in HSC Imaging (SuGOHI). VIII. New galaxy-scale lenses from the HSC SSP

We conduct a search for galaxy-scale strong gravitational lens systems in Data Release 4 of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP), consisting of data taken up to the S21A semester. We select 103191 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample that have deep multiband imaging from the HSC SSP and use the YattaLens algorithm to automatically identify lens candidates with blue arc-like features. The candidates are visually inspected and graded based on their likelihood of being a lens. We find 8 definite lenses, 28 probable lenses, and 138 possible lenses. The new lens candidates generally have lens redshifts in the range $0.3 \lesssim z_{\mathrm{L}} \lesssim 0.9$, a key intermediate redshift range to study the evolution of galaxy structure. Follow-up spectroscopy will confirm these new lenses and measure source redshifts to enable detailed lens modeling.

preprint2022arXiv

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements.

preprint2022arXiv

Three-Dimensional Reconstruction of Weak Lensing Mass Maps with a Sparsity Prior. I. Cluster Detection

We propose a novel method to reconstruct high-resolution three-dimensional mass maps using data from photometric weak-lensing surveys. We apply an adaptive LASSO algorithm to perform a sparsity-based reconstruction on the assumption that the underlying cosmic density field is represented by a sum of Navarro-Frenk-White halos. We generate realistic mock galaxy shape catalogues by considering the shear distortions from isolated halos for the configurations matched to Subaru Hyper Suprime-Cam Survey with its photometric redshift estimates. We show that the adaptive method significantly reduces line-of-sight smearing that is caused by the correlation between the lensing kernels at different redshifts. Lensing clusters with lower mass limits of $10^{14.0} h^{-1}M_{\odot}$, $10^{14.7} h^{-1}M_{\odot}$, $10^{15.0} h^{-1}M_{\odot}$ can be detected with 1.5-$σ$ confidence at the low ($z<0.3$), median ($0.3\leq z< 0.6$) and high ($0.6\leq z< 0.85$) redshifts, respectively, with an average false detection rate of 0.022 deg$^{-2}$. The estimated redshifts of the detected clusters are systematically lower than the true values by $Δz \sim 0.03$ for halos at $z\leq 0.4$, but the relative redshift bias is below $0.5\%$ for clusters at $0.4<z\leq 0.85$. The standard deviation of the redshift estimation is $0.092$. Our method enables direct three-dimensional cluster detection with accurate redshift estimates.

preprint2021arXiv

An analytic model for the sub-galactic matter power spectrum in fuzzy dark matter halos

Fuzzy dark matter (FDM), a scalar particle coupled to the gravitational field without self-interaction whose mass range is $m \sim 10^{-24} - 10^{-20}\ \rm{eV}$, is one of the promising alternative dark matter candidates to cold dark matter. The quantum interference pattern, which is a unique structure of FDM, can be seen in halos in cosmological FDM simulations. In this paper, we first provide an analytic model of the sub-galactic matter power spectrum originating from quantum clumps in FDM halos, in which the density distribution of the FDM is expressed by a superposition of quantum clumps whose size corresponds to the de Broglie wavelength of the FDM. These clumps are assumed to be distributed randomly such that the ensemble averaged density follows the halo profile such as the Navarro-Frenk-White profile. We then compare the convergence power spectrum projected along the line of sight around the Einstein radius, which is converted from the sub-galactic matter power spectrum, to that measured in the strong lens system SDSS J0252+0039. While we find that the current observation provides no useful constraint on the FDM mass, we show that future deep, high spatial resolution observations of strong lens systems can tightly constrain FDM with the mass around $10^{-22}\ \rm{eV}$.

preprint2021arXiv

FPFS Shear Estimator: Systematic Tests on the Hyper Suprime-Cam Survey First Year Data

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 $1\%$ 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.

preprint2021arXiv

HSC Year 1 cosmology results with the minimal bias method: HSC$\times$BOSS galaxy-galaxy weak lensing and BOSS galaxy clustering

We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three redshift bins, and use the HSC-Y1 galaxy catalog to define a secure sample of source galaxies at $z_\mathrm{ph}>0.75$ for the $Δ\!Σ$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $Λ$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $Δ\!Σ(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $σ_8=0.85^{+0.16}_{-0.11}$, and $Ω_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $Λ$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged)

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

preprint2021arXiv

Survey of Gravitationally lensed Objects in HSC Imaging (SuGOHI) -- VII. Discovery and Confirmation of Three Strongly Lensed Quasars

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_{\ell}=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_{\ell}=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.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): X-ray Observable-to-Mass-and-Redshift Relations of Galaxy Clusters and Groups with Weak-Lensing Mass Calibration from the Hyper Suprime-Cam Subaru Strategic Program Survey

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $η$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($η$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged)

preprint2021arXiv

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/.

preprint2020arXiv

A clustering-based self-calibration of the richness-to-mass relation of CAMIRA galaxy clusters out to $z\approx1.1$ in the Hyper Suprime-Cam survey

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. These correlation functions are $ξ_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $ξ_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $ξ_{\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 $ξ_{\mathrm{cc}}$, $ξ_{\mathrm{cc}}+ξ_{\mathrm{cg}}$, and $ξ_{\mathrm{cc}} + ξ_{\mathrm{cg}} + ξ_{\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σ$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

preprint2020arXiv

Active gas features in three HSC-SSP CAMIRA clusters revealed by high angular resolution analysis of MUSTANG-2 SZE and XXL X-ray observations

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 main 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 the integrated Compton $y$ parameter and the temperature for the major merger are 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.

preprint2020arXiv

Cross-correlation of the thermal Sunyaev-Zel'dovich effect and weak gravitational lensing: Planck and Subaru Hyper Suprime-Cam first-year data

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.

preprint2020arXiv

Discovery of an unusually compact lensed Lyman Break Galaxy from the Hyper Suprime-Cam Survey

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 $σ_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$α$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $θ_{\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 $σ_{\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σ$ 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.

preprint2020arXiv

Discovery of Radio Jets in Phoenix Galaxy Cluster Center

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.

preprint2020arXiv

Dual supermassive black holes at close separation revealed by the Hyper Suprime-Cam Subaru Strategic Program

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.

preprint2020arXiv

Geometric Support for Dark Matter by an Unaligned Einstein Ring in Abell 3827

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.

preprint2020arXiv

Lensed quasar search via time variability with the HSC transient survey

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]

preprint2020arXiv

Mass Bias of Weak Lensing Shear-selected Galaxy Cluster Samples

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.

preprint2020arXiv

Redshift Evolution of Green Valley Galaxies in Different Environments from the Hyper Suprime-Cam Survey

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.

preprint2020arXiv

RELICS: A Very Large ($θ_{E}\sim40"$) Cluster Lens -- RXC J0032.1+1808

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 $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{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 $θ_{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.

preprint2020arXiv

Shapes and alignments of dark matter haloes and their brightest cluster galaxies in 39 strong lensing clusters

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.

preprint2020arXiv

Survey of Gravitationally-lensed Objects in HSC Imaging (SuGOHI). V. Group-to-cluster scale lens search from the HSC-SSP Survey

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.

preprint2020arXiv

The BUFFALO HST Survey

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.

preprint2020arXiv

The impact of projection effects on cluster observables: stacked lensing and projected clustering

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.

preprint2020arXiv

The Richness-to-Mass Relation of CAMIRA Galaxy Clusters from Weak-lensing Magnification in the Subaru Hyper Suprime-Cam Survey

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σ$, 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.

preprint2020arXiv

The splashback radius of optically selected clusters with Subaru HSC Second Public Data Release

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σ$ 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σ$) than their $20\%$ smaller values as suggested by the previous studies with the redMaPPer ($\sim 2σ$). 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.

preprint2020arXiv

Truth or delusion? A possible gravitational lensing interpretation of the ultra-luminous quasar SDSS J010013.02+280225.8 at $z=6.30$

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$α$ 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.

preprint2020arXiv

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

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 β= -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 σ(\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σ$ 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.

preprint2019arXiv

Cosmological evolution of orientations of cluster-sized dark matter haloes and their central galaxies in the Horizon-AGN simulation

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.

preprint2019arXiv

Early Low-Mass Galaxies and Star-Cluster Candidates at z~6-9 Identified by the Gravitational Lensing Technique and Deep Optical/Near-Infrared Imaging

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.

preprint2019arXiv

Stellar Properties of z ~ 8 Galaxies in the Reionization Lensing Cluster Survey

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.

preprint2019arXiv

Strong gravitational lensing of explosive transients

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.

preprint2019arXiv

Survey of Gravitationally-lensed Objects in HSC Imaging (SuGOHI). IV. Lensed quasar search in the HSC survey

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$α$ emission of HSCJ141136$-$010216 to be $\sim$254 km/s, showing that it is a probable Lyman-$α$ emitter.

preprint2019arXiv

X-ray properties of high-richness CAMIRA clusters in the Hyper Suprime-Cam Subaru Strategic Program field

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.

preprint2019arXiv

X-ray study of the double source plane gravitational lens system Eye of Horus observed with XMM-Newton

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.

preprint2016arXiv

A Spectroscopically Confirmed Double Source Plane Lens System in the Hyper Suprime-Cam Subaru Strategic Program

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.

preprint2016arXiv

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

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-$σ$ 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 $Λ$-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.

preprint2016arXiv

Measuring the distance-redshift relation with the cross-correlation of gravitational wave standard sirens and galaxies

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.

preprint2016arXiv

Precise Strong Lensing Mass Modeling of Four Hubble Frontier Fields Clusters and a Sample of Magnified High-Redshift Galaxies

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.

preprint2016arXiv

Relationship between the CMB, SZ Cluster Counts, and Local Hubble Parameter Measurements in a Simple Void Model

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.

preprint2016arXiv

Subaru Telescope adaptive optics observations of gravitationally lensed quasars in the Sloan Digital Sky Survey

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.

preprint2016arXiv

The BOSS Emission-Line Lens Survey. III. : Strong Lensing of Ly$α$ Emitters by Individual Galaxies

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ 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$α$ 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.

preprint2016arXiv

The BOSS Emission-Line Lens Survey. IV. : Smooth Lens Models for the BELLS GALLERY Sample

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$α$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ 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$α$ 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.

preprint2016arXiv

The Mass Distribution of the Unusual Merging Cluster Abell 2146 from Strong Lensing

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.

preprint2016arXiv

The SDSS-III BOSS quasar lens survey: discovery of thirteen gravitationally lensed quasars

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.

preprint2015arXiv

CHITAH: Strong-gravitational-lens hunter in imaging surveys

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.

preprint2015arXiv

Detection of a filament connected to CL0016 with weak gravitational lensing

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σ$ significance and the average surface mass density of the filament is $<Σ>=(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.

preprint2015arXiv

Euclid & SKA Synergies

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.

preprint2015arXiv

High-resolution ALMA observations of SDP.81. I. The innermost mass profile of the lensing elliptical galaxy probed by 30 milli-arcsecond images

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.

preprint2015arXiv

High-resolution ALMA Observations of SDP.81. II. Molecular Clump Properties of a Lensed Submillimeter Galaxy at z=3.042

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 ($Σ_{\rm H_2}$) and star-formation rate ($Σ_{\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 $Σ_{\rm H_2}$-$Σ_{\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.

preprint2015arXiv

Inferring Host Dark Matter Halo Masses of Individual Galaxies from Neighboring Galaxy Counts

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.

preprint2015arXiv

Overview of Complementarity and Synergy with Other Wavelengths in Cosmology in the SKA era

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.

preprint2015arXiv

Predicted properties of multiple images of the strongly lensed supernova SN Refsdal

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.

preprint2015arXiv

Properties of Weak Lensing Clusters Detected on Hyper Suprime-Cam 2.3 Square Degree Field

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.

preprint2015arXiv

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

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.

preprint2015arXiv

The sizes of $z\sim6-8$ lensed galaxies from the Hubble Frontier Fields Abell 2744 data

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érsic 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.

preprint2014arXiv

A cluster finding algorithm based on the multiband identification of red sequence galaxies

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.

preprint2014arXiv

Adaptive optics observations of the gravitationally lensed quasar SDSS J1405+0959

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.

preprint2014arXiv

Constraining primordial non-Gaussianity via a multitracer technique with surveys by Euclid and Square Kilometre Array

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 $σ(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 $σ(f_{\rm NL})= 0.61~(0.50)$.

preprint2014arXiv

Detection of the Gravitational Lens Magnifying a Type Ia Supernova

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.

preprint2014arXiv

Discovery of Four Doubly Imaged Quasar Lenses from the Sloan Digital Sky Survey

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.

preprint2014arXiv

Hubble Frontier Fields First Complete Cluster Data: Faint Galaxies at $z\sim 5-10$ for UV Luminosity Functions and Cosmic Reionization

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 ($μ\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, $α$, is as steep as $-2$ at $z\sim6-8$, and strengthen the evidence of the rapid decrease of UV luminosity densities, $ρ_\mathrm{UV}$, at $z>8$ from the large $z\sim9$ sample. We examine whether the rapid $ρ_\mathrm{UV}$ decrease trend can reconcile with the large Thomson scattering optical depth, $τ_\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 $ρ_\mathrm{UV}$ decrease and the large $τ_\mathrm{e}$. It is possible that the $ρ_\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.

preprint2014arXiv

Resolving the Clumpy Structure of the Outflow Winds in the Gravitationally Lensed Quasar SDSS J1029+2623

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, $θ$ $\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, $Δ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.

preprint2014arXiv

Statistical properties of filaments in weak gravitational lensing

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.

preprint2014arXiv

The Physical Conditions, Metallicity and Metal Abundance Ratios In a Highly Magnified Galaxy at z = 3.6252

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]$λ$$λ$3727 and H-$β$ 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]$λ$$λ$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$λ$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.

preprint2014arXiv

The stellar and dark matter distributions in elliptical galaxies from the ensemble of strong gravitational lenses

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 ρ(r) \propto r^γwith best-fit slope γ= -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.

preprint2013arXiv

A Two-Year Time Delay for the Lensed Quasar SDSS J1029+2623

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.

preprint2013arXiv

Baryon acoustic oscillations with the cross-correlation of spectroscopic and photometric samples

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 $Δ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, $σ_z/(1+z)$, are better than 10-20% level.

preprint2013arXiv

Extraordinary Magnification of the Ordinary Type Ia Supernova PS1-10afx

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.

preprint2013arXiv

Line-of-Sight Structure Toward Strong Lensing Galaxy Clusters

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 ($θ\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.

preprint2013arXiv

Measuring the mass distribution of voids with stacked weak lensing

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.

preprint2013arXiv

Reducing Systematic Error in Cluster Scale Weak Lensing

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 $κ$-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.

preprint2013arXiv

The quasar-galaxy cross SDSS J1320+1644: A probable large-separation lensed quasar

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.

preprint2012arXiv

Miscentring in Galaxy Clusters: Dark Matter to Brightest Cluster Galaxy Offsets in 10,000 SDSS 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, Δ, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(Δ[h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and σ=0.501\pm0.004 (95% confidence levels), or <log_{10}(Δ[\arcsec])>=0.564\pm0.005, and σ=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]

preprint2012arXiv

Revising the Halofit Model for the Nonlinear Matter Power Spectrum

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.

preprint2012arXiv

Scatter and bias in weak lensing selected clusters

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.

preprint2012arXiv

Spectroscopy along Multiple, Lensed Sightlines through Outflowing Winds in the Quasar SDSS J1029+2623

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.

preprint2012arXiv

Subaru weak-lensing measurement of a z = 0.81 cluster discovered by the Atacama Cosmology Telescope Survey

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ρ_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 Λ-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.

preprint2012arXiv

The Chandra view of the Largest Quasar Lens SDSS J1029+2623

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.

preprint2012arXiv

The Hidden Fortress: Structure and substructure of the complex strong lensing cluster SDSS J1029+2623

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 θ_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)

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. V. Final Catalog from the Seventh Data Release

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"<θ<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.

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. VI. Constraints on Dark Energy and the Evolution of Massive Galaxies

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 Ω_Λ=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 ν_n=d\lnϕ_*/d\ln(1+z) and the velocity dispersion evolution ν_σ=d\lnσ_*/d\ln(1+z) are constrained to ν_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and ν_σ=-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.

preprint2012arXiv

The Universal Einstein Radius Distribution from 10,000 SDSS Clusters

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(θ_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $σ=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $θ_{e}$ clusters than predicted by $Λ$CDM. We visually inspect each of the clusters with $θ_{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 $θ_{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.

preprint2012arXiv

Very small-scale clustering of quasars from a complete quasar lens survey

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.

preprint2011arXiv

Analyzing the Flux Anomalies of the Large-Separation Lensed Quasar SDSS J1029+2623

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.

preprint2011arXiv

Combined strong and weak lensing analysis of 28 clusters from the Sloan Giant Arcs Survey

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.

preprint2011arXiv

Detailed cluster lensing profiles at large radii and the impact on cluster weak lensing studies

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.

preprint2011arXiv

Gemini/GMOS Spectroscopy of 26 Strong Lensing Selected Galaxy Cluster Cores

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.

preprint2011arXiv

New lensed quasars from the MUSCLES survey

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.

preprint2011arXiv

Probability Distribution Functions of Cosmological Lensing: Convergence, Shear, and Magnification

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 $μ=(1-κ)^{-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.

preprint2011arXiv

SDSS J133401.39+331534.3: A New Subarcsecond Gravitationally Lensed Quasar

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.

preprint2011arXiv

The Redshift Distribution of Giant Arcs in the Sloan Giant Arcs Survey

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.

preprint2010arXiv

Combining cluster observables and stacked weak lensing to probe dark energy: Self-calibration of systematic uncertainties

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)

preprint2010arXiv

Detailed Cluster Mass and Light profiles of A1703, A370 and RXJ1347-11 from Deep Subaru Imaging

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, α~= -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.

preprint2010arXiv

Direct measurement of dark matter halo ellipticity from two-dimensional lensing shear maps of 25 massive clusters

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σsignificance level. The mean ellipticity is found to be e = 0.46 \pm 0.04 (1σ), 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).

preprint2010arXiv

Eight New Quasar Lenses from the Sloan Digital Sky Survey Quasar Lens Search

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 θ=1"13, source redshift z_s=1.826), SDSSJ1054+2733 (θ=1"27, z_s=1.452), SDSSJ1055+4628 (θ=1"15, z_s=1.249), SDSSJ1131+1915 (θ=1"46, z_s=2.915), SDSSJ1304+2001 (θ=1"87, z_s=2.175), SDSSJ1349+1227 (θ=3"00, z_s=1.722), SDSSJ1455+1447 (θ=1"73, z_s=1.424), and SDSSJ1620+1203 (θ=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.

preprint2010arXiv

Gravitationally lensed quasars and supernovae in future wide-field optical imaging surveys

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 σ(w_0)=0.15, σ(w_a)=0.41, and σ(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)

preprint2010arXiv

SDSS J094604.90+183541.8: A Gravitationally Lensed Quasar at z=4.8

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.

preprint2010arXiv

Sunyaev Zel'dovich Effect Observations of Strong Lensing Galaxy Clusters: Probing the Over-Concentration Problem

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.

preprint2010arXiv

The Mass Distribution of SDSS J1004+4112 Revisited

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 μ~ 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.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Lens Search. IV. Statistical Lens Sample from the Fifth Data Release

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"<θ<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 Ω_Λ=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.

preprint2009arXiv

A new gravitational lens from the MUSCLES survey: ULAS J082016.1+081216

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.

preprint2009arXiv

Discovery of a Giant Lya Emitter Near the Reionization Epoch

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.

preprint2009arXiv

LSST Science Book, Version 2.0

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.

preprint2009arXiv

Self-Calibrated Cluster Counts as a Probe of Primordial Non-Gaussianity

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 σ(f_NL) \sim 8, which is competitive with forecasted constraints from next-generation cosmic microwave background experiments.

preprint2009arXiv

Subaru Weak Lensing Measurements of Four Strong Lensing Clusters: Are Lensing Clusters Over-Concentrated?

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.

preprint2008arXiv

Comparison of Cluster Lensing Profiles with Lambda CDM Predictions

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.

preprint2008arXiv

Spectroscopic Confirmation of the Fifth Image of SDSS J1004+4112 and Implications for the M_BH-sigma_* Relation at z=0.68

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.

preprint2008arXiv

What is the largest Einstein radius in the universe?

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σlevel. (Abridged)

preprint2006arXiv

QSO Absorption Lines from QSOs

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.

preprint2006arXiv

Quasars Probing Quasars I: Optically Thick Absorbers Near Luminous Quasars

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.

preprint2006arXiv

Time Delays for the Giant Quadruple Lensed Quasar SDSS J1004+4112: Prospects for Determining the Density Profile of the Lensing Cluster

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 ρ\propto r^{-α}, 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 α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 α. 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 α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 (σ_α \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.

preprint2005arXiv

Binary Quasars in the Sloan Digital Sky Survey: Evidence for Excess Clustering on Small Scales

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.

preprint2005arXiv

Discovery of a Fifth Image of the Large Separation Gravitationally Lensed Quasar SDSS J1004+4112

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.

preprint2005arXiv

SDSS J0806+2006 and SDSS J1353+1138: Two New Gravitationally Lensed Quasars from the Sloan Digital Sky Survey

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.

preprint2005arXiv

SDSSJ102111.02+491330.4: A Newly Discovered Gravitationally Lensed Quasar

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.

preprint2004arXiv

Galaxy Structures and External Perturbations in Gravitational Lenses

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θ$ terms of the primary galaxy and/or of external perturbation significantly reduce the χ^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 χ^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.

preprint2004arXiv

Microlensing of the Broad Emission Line Region in the Quadruple Lens SDSS J1004+4112

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.

preprint2004arXiv

SDSS J1335+0118: A New Two-Image Gravitational Lens

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 Δt ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.