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.

Sample selection is a necessary preparation for weak lensing measurement. It is well-known that selection itself may introduce bias in the measured shear signal. Using image simulation and the Fourier_Quad shear measurement pipeline, we quantify the selection bias in various commonly used selection function (signal-to-noise-ratio, magnitude, etc.). We proposed a new selection function defined in the power spectrum of the galaxy image. This new selection function has low selection bias, and it is particularly convenient for shear measurement pipelines based on Fourier transformation.

We use N-body simulation to study the structure formation in the Cubic Galileon Gravity model where along with the usual kinetic and potential term we also have a higher derivative self-interaction term. We find that the large scale structure provides a unique constraining power for this model. The matter power spectrum, halo mass function, galaxy-galaxy weak lensing signal, marked density power spectrum as well as count in cell are measured. The simulations show that there are less massive halos in the Cubic Galileon Gravity model than corresponding $\Lambda$CDM model and the marked density power spectrum in these two models are different by more than $10\%$. Furthermore, the Cubic Galileon model shows significant differences in voids compared to $\Lambda$CDM. The number of low density cells is far higher in the Cubic Galileon model than that in the $\Lambda$CDM model. Therefore, it would be interesting to put constraints on this model using future large scale structure observations, especially in void regions.

We test the Emergent Gravity(EG) theory using the galaxy-galaxy lensing technique based on SDSS DR7 data. In the EG scenario, we do not expect color dependence of the galaxy sample in the 'apparent dark matter' predicted by EG, which is exerted only by the baryonic mass. If the baryonic mass is similar, then the predicted lensing profiles from the baryonic mass should be similar according to the EG, regardless of the color of the galaxy sample. We use the stellar mass of the galaxy as a proxy of its baryonic mass. We divide our galaxy sample into 5 stellar mass bins, and further classify them as red and blue subsamples in each stellar mass bin. If we set halo mass and concentration as free parameters, $\Lambda$CDM is favored by our data in terms of the reduced $\chi^2$ while EG fails to explain the color dependence of ESDs from the galaxy-galaxy lensing measurement.

We use data from galaxy-galaxy weak lensing to perform a novel test on General Relativity (GR). In particular, we impose strong constraints using the torsional (teleparallel) formulation of gravity in which the deviation from GR is quantified by a single parameter $\alpha$, an approximation which is always valid at low-redshift Universe and weak gravitational fields. We calculate the difference in the deflection angle and eventually derive the modified Excess Surface Density profile, which is mainly affected at small scales. Hence, confronting the predictions with weak lensing data from the Sloan Digital Sky Survey Data Release 7 we obtain the upper bound on the deviation parameter, which, expressed via the dimensionless percentage in the universe energy content, reads as ${\rm{log_{10}}}\Omega_\alpha\le -18.52_{-0.42}^{+0.80}$({\it stat.})$_{-0.39}^{+1.50}$({\it sys.}) with systematics mainly arising from the modelling. To our knowledge, this is the first time that GR is verified at such an accuracy at the corresponding scales.

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

Forthcoming surveys such as the Large Synoptic Survey Telescope (LSST) and Euclid necessitate automatic and efficient identification methods of strong lensing systems. We present a strong lensing identification approach that utilizes a feature extraction method from computer vision, the Histogram of Oriented Gradients (HOG), to capture edge patterns of arcs. We train a supervised classifier model on the HOG of mock strong galaxy-galaxy lens images similar to observations from the Hubble Space Telescope (HST) and LSST. We assess model performance with the area under the curve (AUC) of a Receiver Operating Characteristic (ROC) curve. Models trained on 10,000 lens and non-lens containing images images exhibit an AUC of 0.975 for an HST-like sample, 0.625 for one exposure of LSST, and 0.809 for 10-year mock LSST observations. Performance appears to continually improve with the training set size. Models trained on fewer images perform better in absence of the lens galaxy light. However, with larger training data sets, information from the lens galaxy actually improves model performance, indicating that HOG captures much of the morphological complexity of the arc finding problem. We test our classifier on data from the Sloan Lens ACS Survey and find that small scale image features reduces the efficiency of our trained model. However, these preliminary tests indicate that some parameterizations of HOG can compensate for differences between observed mock data. One example best-case parameterization results in an AUC of 0.6 in the F814 filter image with other parameterization results equivalent to random performance.

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

We investigate the abundance, small-scale clustering and galaxy-galaxy lensing signal of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). To this end, we present new measurements of the redshift and stellar mass dependence of the lensing properties of the galaxy sample. We analyse to what extent models assuming the Planck18 cosmology fit to the number density and clustering can accurately predict the small-scale lensing signal. In qualitative agreement with previous BOSS studies at redshift $z \sim 0.5$ and with results from the Sloan Digital Sky Survey, we find that the expected signal at small scales ($0.1 < r_{\rm p} < 3 \, h^{-1} \mathrm{Mpc}$) is higher by $\sim 25\%$ than what is measured. Here, we show that this result is persistent over the redshift range $0.1 < z < 0.7$ and for galaxies of different stellar masses. If interpreted as evidence for cosmological parameters different from the Planck CMB findings, our results imply $S_8 = \sigma_8 \sqrt{\Omega_{\rm m} / 0.3} = 0.744 \pm 0.015$, whereas $S_8 = 0.832 \pm 0.013$ for Planck18. However, in addition to being in tension with CMB results, such a change in cosmology alone does not accurately predict the lensing amplitude at larger scales. Instead, other often neglected systematics like baryonic feedback or assembly bias are likely contributing to the small-scale lensing discrepancy. We show that either effect alone, though, is unlikely to completely resolve the tension. Ultimately, a combination of the two effects in combination with a moderate change in cosmological parameters might be needed.

In the CDM paradigm, the halo mass function is a sensitive probe of the cosmic structure. In observations, halo mass is typically estimated from its relation with other observables. The resulting halo mass function is subject to systematic bias, such as the Eddington bias, due to the scatter or uncertainty in the observable - mass relation. Exact correction for the bias is not easy, as predictions for the observables are typically model-dependent in simulations. In this paper, we point out an interesting feature in the halo mass function of the concordence $\Lambda$CDM model: the total halo mass within each evenly-spaced logarithmic mass bin is approximately the same over a large mass range. We show that this property allows us to construct an almost bias-free halo mass function using only an observable (as a halo mass estimator) and stacked weak lensing measurements as long as the scatter between the true halo mass and the observable-inferred mass has a stable form in logarithmic units. The method is not sensitive to the form of the mass-observable relation. We test the idea using cosmological simulations, and show that the method performs very well for realistic observables.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

We combine constraints from linear and nonlinear scales, for the first time, to study the interaction between dark matter and dark energy. We devise a novel N-body simulation pipeline for cosmological models beyond $\Lambda$CDM. This pipeline is fully self-consistent and opens a new window to study the nonlinear structure formation in general phenomenological interacting dark energy models. By comparing our simulation results with the SDSS galaxy-galaxy weak lensing measurements, we are able to constrain the strength of interaction between dark energy and dark matter. Compared with the previous studies using linear examinations, we point to plausible improvements on the constraints of interaction strength by using small scale information from weak lensing. This improvement is mostly due to the sensitivity of weak lensing measurements on nonlinear structure formation at low redshift. With this new pipeline, it is possible to look for smoking gun signatures of dark matter-dark energy interaction.

Low density regions are less affected by the nonlinear structure formation and baryonic physics. They are ideal places for probing the nature of dark energy, a possible explanation for the cosmic acceleration. Unlike void lensing, which requires identifications of individual voids, we study the stacked lensing signals around the low-density-positions (LDP), defined as places that are devoid of foreground bright galaxies in projection. The method allows a direct comparison with numerical results by drawing correspondence between the bright galaxies with halos. It leads to lensing signals that are significant enough for differentiating several dark energy models. In this work, we use the CFHTLenS catalogue to define LDPs, as well as measuring their background lensing signals. We consider several different definitions of the foreground bright galaxies (redshift range \& magnitude cut). Regarding the cosmological model, we run six simulations: the first set of simulations have the same initial conditions, with $\rm{w_{de}=-1,-0.5,-0.8,-1.2}$; the second set of simulations include a slightly different $\Lambda$CDM model and a w(z) model from \cite{2017NatAs...1..627Z}. The lensing results indicate that the models with $\rm{w_{de}=-0.5,-0.8}$ are not favored, and the other four models all achieve comparable agreement with the data.

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline.

In our previous work, we propose a cross spectrum based method to extract single pulse signals from RFI contaminated data, which is originated from geodetic VLBI postprocessing. This method fully utilizes fringe phase information of the cross spectrum and hence maximizes signal power, however the localization was not discussed in that work yet. As the continuation of that work, in this paper, we further study how to localize single pulses using astrometric solving method. Assuming that the burst is a point source, we derive the burst position by solving a set of linear equations given the relation between residual delay and offset to a priori position. We find that the single pulse localization results given by both astrometric solving and radio imaging are consistent within 3 sigma level. Therefore we claim that it is possible to derive the position of a single pulse with reasonable precision based on only 3 or even 2 baselines with 4 milliseconds integration. The combination of cross spectrum based detection and the localization proposed in this work then provide a thorough solution for searching single pulse in VLBI observation. According to our calculation, our pipeline gives comparable accuracy as radio imaging pipeline. Moreover, the computational cost of our pipeline is much smaller, which makes it more practical for FRB search in regular VLBI observation. The pipeline is now publicly available and we name it as "VOLKS", which is the acronym of "VLBI Observation for frb Localization Keen Searcher".

We devise a fully self-consistent simulation pipeline for the first time to study the interaction between dark matter and dark energy. We perform convergence tests and show that our code is accurate on different scales. Using the parameters constrained by Planck, Type Ia Supernovae, Baryon Acoustic Oscillations (BAO) and Hubble constant observations, we perform cosmological N-body simulations. We calculate the resulting matter power spectra and halo mass functions for four different interacting dark energy models. In addition to the dark matter density distribution, we also show the inhomogeneous density distribution of dark energy. With this new simulation pipeline, we can further refine and constrain interacting dark energy models.

Photodisintegration reaction rates involving charged particles are of relevance to the p-process nucleosynthesis that aims at explaining the production of the stable neutron-deficient nuclides heavier than iron. In this study, the cross sections and astrophysical rates of (g,p) and (g,a) reactions for about 3000 target nuclei with 10<Z<100 ranging from stable to proton dripline nuclei are computed. To study the sensitivity of the calculations to the optical model potentials (OMPs), both the phenomenological Woods-Saxon and the microscopic folding OMPs are taken into account. The systematic comparisons show that the reaction rates, especially for the (g,a) reaction, are dramatically influenced by the OMPs. Thus the better determination of the OMP is crucial to reduce the uncertainties of the photodisintegration reaction rates involving charged particles. Meanwhile, a gamma-beam facility at ELI-NP is being developed, which will open new opportunities to experimentally study the photodisintegration reactions of astrophysics interest. Considering both the important reactions identified by the nucleosynthesis studies and the purpose of complementing the experimental results for the reactions involving p-nuclei, the measurements of six (g,p) and eight (g,a) reactions based on the gamma-beam facility at ELI-NP and the ELISSA detector for the charged particles detection are proposed, and the GEANT4 simulations are correspondingly performed. The minimum required energies of the gamma-beam to measure these reactions are estimated. It is shown that the direct measurements of these photonuclear reactions within the Gamow windows at T_9=2.5 for p-process are fairly feasible and promising at ELI-NP. The expected experimental results will be used to constrain the OMPs of the charged particles, which can eventually reduce the uncertainties of the reaction rates for the p-process nucleosynthesis.

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

Weak lensing statistics is typically measured as weighted sum of shear estimators or their products (shear-shear correlation). The weighting schemes are designed in the hope of minimizing the statistical error without introducing systematic errors. It would be ideal to approach the Cramer-Rao bound (the lower bound of the statistical uncertainty) in shear statistics, though it is generally difficult to do so in practice. The reasons may include: difficulties in galaxy shape measurement, inaccurate knowledge of the probability-distribution-function (PDF) of the shear estimator, misidentification of point sources as galaxies, etc.. Using the shear estimators defined in Zhang et al. (2015), we show that one can overcome all these problems, and allow shear measurement accuracy to approach the Cramer-Rao bound. This can be achieved by symmetrizing the PDF of the shear estimator, or the joint PDF of shear estimator pairs (for shear-shear correlation), without any prior knowledge of the PDF. Using simulated galaxy images, we demonstrate that under general observing conditions, this idea works as expected: it minimizes the statistical uncertainty without introducing systematic error.

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

For cosmic shear to become an accurate cosmological probe, systematic errors in the shear measurement method must be unambiguously identified and corrected for. Previous work of this series has demonstrated that cosmic shears can be measured accurately in Fourier space in the presence of background noise and finite pixel size, without assumptions on the morphologies of galaxy and PSF. The remaining major source of error is source Poisson noise, due to the finiteness of source photon number. This problem is particularly important for faint galaxies in space-based weak lensing measurements, and for ground-based images of short exposure times. In this work, we propose a simple and rigorous way of removing the shear bias from the source Poisson noise. Our noise treatment can be generalized for images made of multiple exposures through MultiDrizzle. This is demonstrated with the SDSS and COSMOS/ACS data. With a large ensemble of mock galaxy images of unrestricted morphologies, we show that our shear measurement method can achieve sub-percent level accuracy even for images of signal-to-noise ratio less than 5 in general, making it the most promising technique for cosmic shear measurement in the ongoing and upcoming large scale galaxy surveys.

Using a sample of spiral galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Galaxy Zoo 2 (GZ2), we investigate the alignment of spin axes of spiral galaxies with their surrounding large scale structure, which is characterized by the large-scale tidal field reconstructed from the data using galaxy groups above a certain mass threshold. We find that the spin axes of only have weak tendency to be aligned with (or perpendicular to) the intermediate (or minor) axis of the local tidal tensor. The signal is the strongest in a \cluster environment where all the three eigenvalues of the local tidal tensor are positive. Compared to the alignments between halo spins and local tidal field obtained in N-body simulations, the above observational results are in best agreement with those for the spins of inner regions of halos, suggesting that the disk material traces the angular momentum of dark matter halos in the inner regions.

Major mergers and interactions between gas-rich galaxies with comparable masses are thought to be the main triggers of starburst. In this work, we study, for a large stellar mass range, the interaction rate of the starburst galaxies in the local universe. We focus independently on central and satellite star forming galaxies extracted from the Sloan Digital Sky Survey. Here the starburst galaxies are selected in the star formation rate (SFR) stellar mass plane with SFR five times larger than the median value found for "star forming" galaxies of the same stellar mass. Through visual inspection of their images together with close companions determined using spectroscopic redshifts, we find that ~50% of the "starburst" populations show evident merger features, i.e., tidal tails, bridges between galaxies, double cores and close companions. In contrast, in the control sample we selected from the normal star forming galaxies, only ~19% of galaxies are associated with evident mergers. The interaction rates may increase by ~5% for the starburst sample and 2% for the control sample if close companions determined using photometric redshifts are considered. The contrast of the merger rate between the two samples strengthens the hypothesis that mergers and interactions are indeed the main causes of starburst.

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

We select satellite galaxies from the galaxy group catalog constructed with the SDSS spectroscopic galaxies and measure the tangential shear around these galaxies with source catalog extracted from CFHT/MegaCam Stripe-82 Survey to constrain the mass of subhalos associated with them. The lensing signal is measured around satellites in groups with masses in the range [10^{13}, 5x10^{14}]h^{-1}M_{sun}, and is found to agree well with theoretical expectation. Fitting the data with a truncated NFW profile, we obtain an average subhalo mass of log M_{sub}= 11.68 \pm 0.67 for satellites whose projected distances to central galaxies are in the range [0.1, 0.3] h^{-1}Mpc, and log M_{sub}= 11.68 \pm 0.76 for satellites with projected halo-centric distance in [0.3, 0.5] h^{-1}Mpc. The best-fit subhalo masses are comparable to the truncated subhalo masses assigned to satellite galaxies using abundance matching and about 5 to 10 times higher than the average stellar mass of the lensing satellite galaxies.

We use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$\sigma$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

Aims and Methods: We present the results of VLBI observations of nineteen GHz-Peaked-Spectrum (GPS) radio sources at 1.6 GHz. Of them, 15 sources are selected from the Parkes Half Jansky (PHJ) sample (Snellen 2002), 4 others are from our previous observation list. We aimed at imaging the structure of GPS sources, searching for Compact Symmetric Objects (CSOs) and studying the absorption for the convex radio spectra of GPS sources. Results: We obtained total intensity 1.6 GHz VLBI images of 17 sources for the first time. Of them, 80% show mini-double-lobe radio structure, indicating that they are CSOs or candidates, and their host AGNs could be edge-on to us. This result suggests that there is a high incidence of mini double-lobe sources (or CSOs) in the PHJ sample. The sources J0323+0534, J1135-0021, J1352+0232, J2058+0540, J2123-0112 and J2325-0344 with measured redshift, showing double-lobe structure with sizes of <1 kpc, are classified as CSOs. Three sources J1057+0012, J1600-0037 and J1753+2750 are considered as core-jet sources according to their morphologies and flux variability.

We present images of quasi-simultaneous VLBI observations of the GHz-Peaked-Spectrum radio source OQ 208 with the Very Long Baseline Array at 1.4, 1.7, 2.3, 5.0, 8.4, 15.4 GHz and the European VLBI Network at 6.7 GHz. The low frequency (1.4, 1.7 and 2.3 GHz) observations reveal a weak and extended steep-spectrum component at about 30 mas away in the position angle of $- 110^\circ$ which may be a remnant emission. The radio structure of OQ 208 consists of two mini-lobes at 5.0, 6.7, 8.4 and 15.4 GHz. Our spectral analysis further confirms that the southwest lobe undergoes free-free absorption and finds that the free-free absorption is stronger in the inner region. By fitting the 8.4 GHz images from 1994 to 2005, we obtain a separation speed of 0.031 $\pm$ 0.006 mas yr$^{-1}$ between the two mini-lobes. This indicates a jet proper motion of 0.105 $\pm$ 0.020 $c$ and a kinematic age of 219 $\pm$ 42 yr for the radio source.