We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping BOSS and 2dFLenS surveys. Cross-correlations between galaxy positions and ellipticities are incorporated into the analysis, developing a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and assessing signal contributions by astrophysical effects. All weak lensing signals are measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias is updated, and a more complete tally of residual calibration uncertainties is propagated into the likelihood. A dedicated suite of more than 20000 mocks is used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the $\chi^2$ goodness-of-fit statistic are validated, with proposed changes to calculating the effective number of degrees of freedom. Standard weak lensing point estimates on $S_8=\sigma_8\,(\Omega_{\rm m}/0.3)^{1/2}$ derived from its marginal posterior are easily misinterpreted to be biased low, and an alternative estimator and associated credible interval are proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias $S_8$ by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, $S_8$ constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis. [abridged]

We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $\sim 6.2$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1\sigma$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = \sigma_8\sqrt{\Omega_{\rm m}/0.3}$. Through a series of two-point null-tests we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is however shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.

In this work we study the shape of the projected surface mass density distribution of galaxy clusters using weak-lensing stacking techniques. In particular, we constrain the average aligned component of the projected ellipticity, $\epsilon$, for a sample of redMaPPer clusters ($0.1 \leq z < 0.4$). We consider six different proxies for the cluster orientation and measure $\epsilon$ for three ranges of projected distances from the cluster centres. The mass distribution in the inner region (up to $700\,$kpc) is better traced by the cluster galaxies with a higher membership probability, while the outer region (from $700\,$kpc up to $5\,$Mpc) is better traced by the inclusion of less probable galaxy cluster members. The fitted ellipticity in the inner region is $\epsilon = 0.21 \pm 0.04$, in agreement with previous estimates. We also study the relation between $\epsilon$ and the cluster mean redshift and richness. By splitting the sample in two redshift ranges according to the median redshift, we obtain larger $\epsilon$ values for clusters at higher redshifts, consistent with the expectation from simulations. In addition, we obtain higher ellipticity values in the outer region of clusters at low redshifts. We discuss several systematic effects that might affect the measured lensing ellipticities and their relation to the derived ellipticity of the mass distribution.

The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish General Relativity from its potential modifications. We use the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measure the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of General Relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics.

To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

Galaxy intrinsic alignment (IA) is both a source of systematic contamination to cosmic shear measurement and its cosmological applications, and a source of valuable information on the large scale structure of the universe and galaxy formation. The self-calibration (SC) method \citep{SC2008} was designed to separate IA from cosmic shear, free of IA modeling. It was first successfully applied to the KiDS450 and KV450 data \citep{Yao2019}. We apply the SC method to the DECaLS DR3 shear + photo-z catalog and significantly improve the IA detection to $\sim 20\sigma$. We find a strong dependence of IA on galaxy color, with strong IA signal for red galaxies, and weak/vanishing IA signal for blue galaxies. The detected IA for red galaxies is in reasonable agreement with the non-linear tidal alignment model and the inferred IA amplitude increases with redshift. We discuss various caveats and possible improvements in the measurement, theory and parameter fitting that will be addressed in future works.

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

We present the weak lensing analysis of 279 CODEX clusters using imaging data from 4200 $\text{deg}^{2}$ of the DECam Legacy Survey (DECaLS) Data Release 3. The cluster sample results from a joint selection in X-ray, optical richness in the range $20 \leq \lambda < 110$, and redshift in the range $0.1 \leq z \leq 0.2$. We model the cluster mass ($M_{\rm 200c}$) and the richness relation with the expression $\left\langle M_{\rm 200c} | \lambda \right\rangle \propto M_{0} \, (\lambda / 40)^{F_{\lambda}}$. By measuring the CODEX cluster sample as an individual cluster, we obtain the best-fit values, $M_{0} = 3.24^{+0.29}_{-0.27} \times 10^{14} \text{M}_{\odot}$, and $F_{\lambda} = 1.00 ^{+0.22}_{-0.22}$ for the richness scaling index, consistent with a power law relation. Moreover, we separate the cluster sample into three richness groups; $\lambda = 20 - 30, 30 - 50$ and $50 - 110$, and measure the stacked excess surface mass density profile in each group. The results show that both methods are consistent. In addition, we find an excellent agreement between our weak lensing based scaling relation and the relation obtained with dynamical masses estimated from cluster member velocity dispersions measured by the SDSS-IV/SPIDERS team. This suggests that the cluster dynamical equilibrium assumption involved in the dynamical mass estimates is statistically robust for a large sample of clusters.

Intrinsic galaxy alignments are a source of bias for weak lensing measurements as well as a tool for understanding galaxy formation and evolution. In this work, we measure the alignment of shapes of satellite galaxies, in galaxy groups, with respect to the brightest group galaxy (BGG), as well as alignments of the BGG shape with the satellite positions, using the highly complete Galaxy And Mass Assembly (GAMA) spectroscopic survey and deep imaging from the Kilo Degree Survey. We control systematic errors with dedicated image simulations and measure accurate shapes using the DEIMOS shape measurement method. We find a significant satellite radial alignment signal, which vanishes at large separations from the BGG. We do not identify any strong trends of the signal with galaxy absolute magnitude or group mass. The alignment signal is dominated by red satellites. We also find that the outer regions of galaxies are aligned more strongly than their inner regions, by varying the radial weight employed during the shape measurement process. This behaviour is evident for both red and blue satellites. BGGs are also found to be aligned with satellite positions, with this alignment being stronger when considering the innermost satellites, using red BGGs and the shape of the outer region of the BGG. Lastly, we measure the global intrinsic alignment signal in the GAMA sample for two different radial weight functions and find no significant difference.

The KiDS Strongly lensed QUAsar Detection project (KiDS-SQuaD) aims at finding as many previously undiscovered gravitational lensed quasars as possible in the Kilo Degree Survey. This is the second paper of this series where we present a new, automatic object classification method based on machine learning technique. The main goal of this paper is to build a catalogue of bright extragalactic objects (galaxies and quasars), from the KiDS Data Release 4, with a minimum stellar contamination, preserving the completeness as much as possible, to then apply morphological methods to select reliable gravitationally lensed (GL) quasar candidates. After testing some of the most used machine learning algorithms, decision trees based classifiers, we decided to use CatBoost, that was specifically trained with the aim of creating a sample of extragalactic sources as clean as possible from stars. We discuss the input data, define the training sample for the classifier, give quantitative estimates of its performances, and finally describe the validation results with Gaia DR2, AllWISE, and GAMA catalogues. We have built and make available to the scientific community the KiDS Bright EXtraGalactic Objects catalogue (KiDS-BEXGO), specifically created to find gravitational lenses. This is made of $\approx6$ millions of sources classified as quasars ($\approx 200\,000$) and galaxies ($\approx 5.7$M), up to $r<22^m$. From this catalog we selected 'Multiplets': close pairs of quasars or galaxies surrounded by at least one quasar, presenting the 12 most reliable gravitationally lensed quasar candidates, to demonstrate the potential of the catalogue, which will be further explored in a forthcoming paper. We compared our search to the previous one, presented in the first paper from this series, showing that employing a machine learning method decreases the stars-contaminators within the GL candidates.

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope, specifically designed for measuring weak gravitational lensing by galaxies and large-scale structure. When completed it will consist of 1350 square degrees imaged in four filters (ugri). Here we present the fourth public data release which more than doubles the area of sky covered by data release 3. We also include aperture-matched ZYJHKs photometry from our partner VIKING survey on the VISTA telescope in the photometry catalogue. We illustrate the data quality and describe the catalogue content. Two dedicated pipelines are used for the production of the optical data. The Astro-WISE information system is used for the production of co-added images in the four survey bands, while a separate reduction of the r-band images using the theli pipeline is used to provide a source catalogue suitable for the core weak lensing science case. All data have been re-reduced for this data release using the latest versions of the pipelines. The VIKING photometry is obtained as forced photometry on the theli sources, using a re-reduction of the VIKING data that starts from the VISTA pawprints. Modifications to the pipelines with respect to earlier releases are described in detail. The photometry is calibrated to the Gaia DR2 G band using stellar locus regression. In this data release a total of 1006 square-degree survey tiles with stacked ugri images are made available, accompanied by weight maps, masks, and single-band source lists. We also provide a multi-band catalogue based on r-band detections, including homogenized photometry and photometric redshifts, for the whole dataset. Mean limiting magnitudes (5 sigma in a 2" aperture) are 24.23, 25.12, 25.02, 23.68 in ugri, respectively, and the mean r-band seeing is 0.70".

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The combination of Galaxy-Galaxy Lensing (GGL) and Redshift Space Distortion of galaxy clustering (RSD) is a privileged technique to test General Relativity predictions, and break degeneracies between the growth rate of structure parameter $f$ and the amplitude of the linear power-spectrum $\sigma_8$. We perform a joint GGL and RSD analysis on 250 sq. degrees using shape catalogues from CFHTLenS and CFHT-Stripe 82, and spectroscopic redshifts from the BOSS CMASS sample. We adjust a model that includes non-linear biasing, RSD and Alcock-Paczynski effects. We find $f(z=0.57) =0.95\pm0.23$, $\sigma_8(z=0.57)=0.55\pm0.07$ and $\Omega_{\rm m} = 0.31\pm0.08$, in agreement with Planck cosmological results 2018. We also estimate the probe of gravity $E_{\rm G} = 0.43\pm0.10$ in agreement with $\Lambda$CDM-GR predictions of $E_{\rm G} = 0.40$. This analysis reveals that RSD efficiently decreases the GGL uncertainty on $\Omega_{\rm m}$ by a factor of 4, and by 30\% on $\sigma_8$. We use an N-body simulation supplemented by an abundance matching prescription for CMASS to build a set of overlapping lensing and clustering mocks. Together with additional spectroscopic data, this helps us to quantify and correct several systematic errors, such as photometric redshifts. We make our mock catalogues available on the Skies and Universe database.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention.

The Square Kilometre Array (SKA) will be the first low-frequency instrument with the capability to directly image the structures of the Epoch of Reionization (EoR). Indeed, deep imaging of the EoR over 5 targeted fields of 20 square degrees each has been selected as the highest priority science objective for SKA1, probably the most exciting scientific outcome to be expected for the SKA1-Low. Aiming at preparing for this highly challenging observation, we perform an extensive pre-selection of the 'quietest' and 'cleanest' candidate fields in the southern sky to be suited for deep imaging of the EoR using existing catalogs and observations over a broad band. The candidate fields should meet a number of strict criteria to avoid contamination from foreground structures and sources such as the Galactic plane, LMC/SMC, clusters of galaxies and bright radio sources. The candidate fields should also exhibit both the lowest average surface brightness and smallest variance to ensure uniformity and high quality deep imaging over the fields. We use two sets of field of view in our selection: the exactly same size of 20 square degrees as defined in current SKA1 documentation and 40 square degrees - an area of twice big as the EoR imaging fields. The choice of a larger area allows us to maximally reduce the environmental effect from radio structures/sources in the vicinity of the candidate field. Our selection yields a sample of 29 (9) `ideal' fields of 20 (40) square degrees in the southern sky that could be targeted for deep imaging of the EoR. We have launched a campaign to observe these candidate fields with some of the SKA pathfinders to set further constraints on these candidates. These pre-selected candidate fields can provide a valuable guidance for eventual selection of the five targeted fields for imaging the EoR structures with the future SKA1-low.

In this work we analyze a sample of close galaxy pairs (relative projected separation $<25 h^{-1}$kpc and relative radial velocities $< 350$ km s$^{-1}$) using a weak lensing analysis based on the CFHT Stripe 82 Survey. We determine halo masses for the Total sample of pairs as well as for Interacting, Red and Higher luminosity pair subsamples with $\sim 3\sigma$ confidence. The derived lensing signal for the total sample can be fitted either by a singular isothermal sphere with $\sigma_V = 223 \pm 24$ km s$^{-1}$ or a NFW profile with $R_{200} = 0.30 \pm 0.03\,h^{-1}$ Mpc. The pair total masses and total $r$ band luminosities imply an average mass-to-light ratio of $\sim 200\,h\,M_\odot/L_\odot$. On the other hand, Red pairs which include a larger fraction of elliptical galaxies, show a larger mass-to-light ratio of $\sim 345\,h\,M_\odot/L_\odot$. Derived lensing masses were compared to a proxy of the dynamical mass, obtaining a good correlation. However, there is a large discrepancy between lensing masses and the dynamical mass estimates, which could be accounted by astrophysical processes such as dynamical friction, by the inclusion of unbound pairs, and by significant deviations of the density distribution from a SIS and NFW profiles in the inner regions. We also compared lensing masses with group mass estimates obtained from the Yang et al. galaxy group catalog, finding a very good agreement with the sample of groups with 2 members. Red and Blue pairs show large differences between group and lensing masses, which is likely due to the single mass-to-light ratio adopted to compute the group masses.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

In this work we present a lensing study of Compact Groups (CGs) using data obtained from the high quality Canada-France-Hawaii Telescope Stripe 82 Survey. Using stacking techniques we obtain the average density contrast profile. We analyse the lensing signal dependence on the groups surface brightness and morphological content, for CGs in the redshift range $z = 0.2 - 0.4$. We obtain a larger lensing signal for CGs with higher surface brightness, probably due to their lower contamination by interlopers. Also, we find a strong dependence of the lensing signal on the group concentration parameter, with the most concentrated quintile showing a significant lensing signal, consistent with an isothermal sphere with $\sigma_V =336 \pm 28$ km/s and a NFW profile with $R_{200}=0.60\pm0.05$ $h_{70}^{-1}$Mpc. We also compare lensing results with dynamical estimates finding a good agreement with lensing determinations for CGs with higher surface brightness and higher concentration indexes. On the other hand, CGs that are more contaminated by interlopers show larger dynamical dispersions, since interlopers bias dynamical estimates to larger values, although the lensing signal is weakened.

We present the first weak lensing calibration of $\mu_{\star}$, a new galaxy cluster mass proxy corresponding to the total stellar mass of red and blue members, in two cluster samples selected from the SDSS Stripe 82 data: 230 redMaPPer clusters at redshift $0.1\leq z<0.33$ and 136 Voronoi Tessellation (VT) clusters at $0.1 \leq z < 0.6$. We use the CS82 shear catalog and stack the clusters in $\mu_{\star}$ bins to measure a mass-observable power law relation. For redMaPPer clusters we obtain $M_0 = (1.77 \pm 0.36) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 1.74 \pm 0.62$. For VT clusters, we find $M_0 = (4.31 \pm 0.89) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.59 \pm 0.54$ and $M_0 = (3.67 \pm 0.56) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.68 \pm 0.49$ for a low and a high redshift bin, respectively. Our results are consistent, internally and with the literature, indicating that our method can be applied to any cluster finding algorithm. In particular, we recommend that $\mu_{\star}$ be used as the mass proxy for VT clusters. Catalogs including $\mu_{\star}$ measurements will enable its use in studies of galaxy evolution in clusters and cluster cosmology.

In Montero-Dorta et al. 2017, we show that luminous red galaxies (LRGs) from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at $z\sim0.55$ can be divided into two groups based on their star formation histories. So-called fast-growing LRGs assemble $80\%$ of their stellar mass at $z\sim5$, whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. We further demonstrate that these two subpopulations present significantly different clustering properties on scales of $\sim1 - 30 \mathrm{Mpc}$. Here, we measure the mean halo mass of each subsample using the galaxy-galaxy lensing technique, in the $\sim190\deg^2$ overlap of the LRG catalogue and the CS82 and CFHTLenS shear catalogues. We show that fast- and slow-growing LRGs have similar lensing profiles, which implies that they live in haloes of similar mass: $\log\left(M_{\rm halo}^{\rm fast}/h^{-1}\mathrm{M}_{\odot}\right) = 12.85^{+0.16}_{-0.26}$ and $\log\left(M_{\rm halo}^{\rm slow}/h^{-1}\mathrm{M}_{\odot}\right) =12.92^{+0.16}_{-0.22}$. This result, combined with the clustering difference, suggests the existence of galaxy assembly bias, although the effect is too subtle to be definitively proven given the errors on our current weak-lensing measurement. We show that this can soon be achieved with upcoming surveys like DES.

We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.

The spatial distribution of satellite galaxies encodes rich information of the structure and assembly history of galaxy clusters. In this paper, we select a redMaPPer cluster sample in SDSS Stripe 82 region with $0.1 \le z \le 0.33$, $20<\lambda<100$ and $P_{\rm cen}>0.7$. Using the high-quality weak lensing data from CS82 Survey, we constrain the mass profile of this sample. Then we compare directly the mass density profile with the satellite number density profile. We find that the total mass and number density profiles have the same shape, both well fitted by an NFW profile. The scale radii agree with each other within 1$\sigma$ error ($r_{\rm s,gal}=0.34_{-0.03}^{+0.04}$Mpc vs $r_{\rm s}=0.37_{-0.10}^{+0.15}$Mpc ).

We report the result of a faint quasar survey in a one square degree field. The aim is to test the Y-K/g-z and J-K/i-Y color selection criteria for quasars at faint magnitude, to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection, and to measure the faint end of quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y-K/g-z and J-K/i-Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one square degree field. The survey significantly increases the number of faint quasars in this field, especially at z ~ 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z ~ 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars, and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

In this work, we reinvestigate the electron fraction $Y_{e}$ and electron Fermi energy $E_{F}(e)$ of neutron stars, based on our previous work of Li et al.(2016), in which we firstly deduced a special solution to $E_{F}(e)$, and then obtained several useful analytical formulae for $Y_{\rm e}$ and matter density $\rho$ within classical models and the relativistic mean field(RMF) theory using numerically fitting. The advantages of this work include the following aspects:(1) The linear functions are substituted for the nonlinear exponential functions used in the previous work. This method may be more simple, and closer to realistic equation of state\,(EoS) of a neutron star(NS), because there are linear or quasi-linear relationships between number fractions of leptons and matter density, which can be seen by solving NS EoS; (2)we introduce a dimensionless variable $\varrho$\,($\varrho=\rho/\rho_0$, $\rho_{0}$ is the standard saturated nuclear density), which greatly reduces the scope of the fitting coefficients;(3)we present numerical errors including absolute and relative deviations between the data and fit. By numerically simulating, we have obtained several analytical formulae for $Y_{e}$ and $\rho$ for both APR98 and RMF models. Combining these analytical formulae with the special solution, we can calculate the value of $E_{\rm F}(e)$ for any given matter density. Since $Y_e$ and $E_{ F}(e)$ are important in assessing cooling rate of a NS and the possibility of kaon/pion condensation in the NS interior, this study could be useful in the future study on the thermal evolution of a NS.

The low braking-index pulsar PSR J1734$-$3333 could be born with superhigh internal magnetic fields $B_{\rm in}\sim10^{15}-10^{16}$ G, and undergo a supercritical accretion soon after its formation in a supernova explosion. The buried multipole magnetic fields will merger into a dipole magnetic field. Since the magnetic flow transfers from the core to the crust of the pulsar, its surface dipole field grows quickly at a power-law form assumed until it saturates at the level of internal dipole field. The increase in surface dipole magnetic field results in the observed low braking index of $n=0.9(2)$. Keeping an average field growth index $\varepsilon=1.34(6)$, this pulsar will become a magnetar with surface dipole magnetic field at the equator $B_{\rm d}\sim 2.6(1)\times 10^{14}$\,G and $\sim 5.3(2)\times 10^{14}$\,G after next 50\,kyrs and 100\,kys, respectively.

In this work, we interpreted the high braking index of PSR J1640$-$4631 with a combination of the magneto-dipole radiation and dipole magnetic field decay models. By introducing a mean rotation energy conversion coefficient $\overline{\zeta}$, the ratio of the total high-energy photon energy to the total rotation energy loss in the whole life of the pulsar, and combining the pulsar's high-energy and timing observations with reliable nuclear equation of state, we estimate the pulsar's initial spin period, $P_{0}\sim (17-44)$ ms, corresponding to the moment of inertia $I\sim (0.8-2.1)\times 10^{45}$ g cm$^{2}$. Assuming that PSR J1640$-$4631 has experienced a long-term exponential decay of the dipole magnetic field, we calculate the true age $t_{\rm age}$, the effective magnetic field decay timescale $\tau_{D}$, and the initial surface dipole magnetic field at the pole $B_{p}(0)$ of the pulsar to be $(2900-3100)$ yrs, $1.07(2)\times10^{5}$ yrs, and $(1.84-4.20)\times10^{13}$ G, respectively. The measured braking index of $n=3.15(3)$ for PSR J1640$-$4631 is attributed to its long-term dipole magnetic field decay and a low magnetic field decay rate, $dB_{\rm p}/dt\sim -(1.66-3.85)\times10^{8}$ G yr$^{-1}$. Our model can be applied to both the high braking index ($n>3$) and low braking index ($n<3$) pulsars, tested by the future polarization, timing, and high-energy observations of PSR J1640$-$4631.

This paper is the first of a series of papers constraining cosmological parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of imaging data from the Kilo Degree Survey (KiDS-450). We measure high signal-to-noise ratio (SNR: $\nu$) weak lensing convergence peaks in the range of $3<\nu<5$, and employ theoretical models to derive expected values. These models are validated using a suite of simulations. We take into account two major systematic effects, the boost factor and the effect of baryons on the mass-concentration relation of dark matter haloes. In addition, we investigate the impacts of other potential astrophysical systematics including the projection effects of large scale structures, intrinsic galaxy alignments, as well as residual measurement uncertainties in the shear and redshift calibration. Assuming a flat $\Lambda$CDM model, we find constraints for $S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{-0.107}$ according to the degeneracy direction of the cosmic shear analysis and $\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.38}=0.696^{+0.048}_{-0.050}$ based on the derived degeneracy direction of our high-SNR peak statistics. The difference between the power index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and $\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic correlation analysis of the same dataset and $\sim 2\sigma$ lower than the Planck 2016 results.

We study the statistics of peaks in a weak lensing reconstructed mass map of the first 450 square degrees of the Kilo Degree Survey. The map is computed with aperture masses directly applied to the shear field with an NFW-like compensated filter. We compare the peak statistics in the observations with that of simulations for various cosmologies to constrain the cosmological parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$, which probes the ($\Omega_{\rm m}, \sigma_8$) plane perpendicularly to its main degeneracy. We estimate $S_8=0.750\pm0.059$, using peaks in the signal-to-noise range $0 \leq {\rm S/N} \leq 4$, and accounting for various systematics, such as multiplicative shear bias, mean redshift bias, baryon feedback, intrinsic alignment, and shear-position coupling. These constraints are $\sim25\%$ tighter than the constraints from the high significance peaks alone ($3 \leq {\rm S/N} \leq 4$) which typically trace single-massive halos. This demonstrates the gain of information from low-S/N peaks. However we find that including ${\rm S/N} < 0$ peaks does not add further information. Our results are in good agreement with the tomographic shear two-point correlation function measurement in KiDS-450. Combining shear peaks with non-tomographic measurements of the shear two-point correlation functions yields a $\sim20\%$ improvement in the uncertainty on $S_8$ compared to the shear two-point correlation functions alone, highlighting the great potential of peaks as a cosmological probe.

We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques.

In the hierarchical formation model, galaxy clusters grow by accretion of smaller groups or isolated galaxies. During the infall into the centre of a cluster, the properties of accreted galaxies change. In particular, both observations and numerical simulations suggest that its dark matter halo is stripped by the tidal forces of the host. We use galaxy-galaxy weak lensing to measure the average mass of dark matter haloes of satellite galaxies as a function of projected distance to the centre of the host, for different stellar mass bins. Assuming that the stellar component of the galaxy is less disrupted by tidal stripping, stellar mass can be used as a proxy of the infall mass. We study the stellar to halo mass relation of satellites as a function of the cluster-centric distance to measure tidal stripping. We use the shear catalogues of the DES science verification archive, the CFHTLenS and the CFHT Stripe 82 (CS82) surveys, and we select satellites from the redMaPPer catalogue of clusters. For galaxies located in the outskirts of clusters, we find a stellar to halo mass relation in good agreement with the theoretical expectations from \citet{moster2013} for central galaxies. In the centre of the cluster, we find that this relation is shifted to smaller halo mass for a given stellar mass. We interpret this finding as further evidence for tidal stripping of dark matter haloes in high density environments.

We set out to quantify the number density of quiescent massive compact galaxies at intermediate redshifts. We determine structural parameters based on i-band imaging using the CFHT equatorial SDSS Stripe 82 (CS82) survey (~170 sq. degrees) taking advantage of an exquisite median seeing of ~0.6''. We select compact massive (M > 5x10^10 M_sun) galaxies within the redshift range of 0.2<z<0.6. The large volume sampled allows to decrease the effect of cosmic variance that has hampered the calculation of the number density for this enigmatic population in many previous studies. We undertake an exhaustive analysis in an effort to untangle the various findings inherent to the diverse definition of compactness present in the literature. We find that the absolute number of compact galaxies is very dependent on the adopted definition and can change up to a factor of >10. We systematically measure a factor of ~5 more compacts at the same redshift than what was previously reported on smaller fields with HST imaging, which are more affected by cosmic variance. This means that the decrease in number density from z ~ 1.5 to z ~ 0.2 might be only of a factor of ~2-5, significantly smaller than what previously reported. This supports progenitor bias as the main contributor to the size evolution. This milder decrease is roughly compatible with the predictions from recent numerical simulations. Only the most extreme compact galaxies, with Reff < 1.5x( M/10^11 M_sun)^0.75 and M > 10^10.7 M_sun, appear to drop in number by a factor of ~20 and hence likely experience a noticeable size evolution.

We present a stacked weak lensing analysis of 27 richness selected galaxy clusters at $0.40 \leqslant z \leqslant 0.62$ in the CODEX survey. The fields were observed in 5 bands with the CFHT. We measure the stacked surface mass density profile with a $14\sigma$ significance in the radial range $0.1 < R\ Mpc\ h^{-1} < 2.5$. The profile is well described by the halo model, with the main halo term following an NFW profile and including the off-centring effect. We select the background sample using a conservative colour-magnitude method to reduce the potential systematic errors and contamination by cluster member galaxies. We perform a Bayesian analysis for the stacked profile and constrain the best-fit NFW parameters $M_{200c} = 6.6^{+1.0}_{-0.8} \times 10^{14} h^{-1} M_{\odot}$ and $c_{200c} = 3.7^{+0.7}_{-0.6}$. The off-centring effect was modelled based on previous observational results found for redMaPPer SDSS clusters. Our constraints on $M_{200c}$ and $c_{200c}$ allow us to investigate the consistency with numerical predictions and select a concentration-mass relation to describe the high richness CODEX sample. Comparing our best-fit values for $M_{200c}$ and $c_{200c}$ with other observational surveys at different redshifts, we find no evidence for evolution in the concentration-mass relation, though it could be mitigated by particular selection functions. Similar to previous studies investigating the X-ray luminosity-mass relation, our data suggests a lower evolution than expected from self-similarity.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions.

Mass calibration uncertainty is the largest systematic effect for using clusters of galaxies to constrain cosmological parameters. We present weak lensing mass measurements from the Canada-France-Hawaii Telescope Stripe 82 Survey for galaxy clusters selected through their high signal-to-noise thermal Sunyaev-Zeldovich (tSZ) signal measured with the Atacama Cosmology Telescope (ACT). For a sample of 9 ACT clusters with a tSZ signal-to-noise greater than five the average weak lensing mass is $\left(4.8\pm0.8\right)\,\times10^{14}\,\mathrm{M}_\odot$, consistent with the tSZ mass estimate of $\left(4.70\pm1.0\right)\,\times10^{14}\,\mathrm{M}_\odot$ which assumes a universal pressure profile for the cluster gas. Our results are consistent with previous weak-lensing measurements of tSZ-detected clusters from the Planck satellite. When comparing our results, we estimate the Eddington bias correction for the sample intersection of Planck and weak-lensing clusters which was previously excluded.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass Function (IMF) and a truncated NFW model for the subhalo mass distribution, we find that the sub-halo mass to galaxy stellar mass ratio increases as a function of projected halo-centric radius $r_p$, from $M_{\rm sub}/M_{\rm star}=4.43^{+ 6.63}_{- 2.23}$ at $r_p \in [0.1,0.3]$ $h^{-1}Mpc$ to $M_{\rm sub}/M_{\rm star}=75.40^{+ 19.73}_{- 19.09}$ at $r_p \in [0.6,0.9]$ $h^{-1}Mpc$. We also investigate the dependence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar mass bins: $10<\log(M_{\rm star}/M_{\rm sun})<10.5$ and $11<\log(M_{\rm star}/M_{\rm sun})<12$. The best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of the less massive satellites: $\log(M_{\rm sub}/M_{\rm sun})=11.14 ^{+ 0.66 }_{- 0.73}$ ($M_{\rm sub}/M_{\rm star}=19.5^{+19.8}_{-17.9}$) versus $\log(M_{\rm sub}/M_{\rm sun})=12.38 ^{+ 0.16 }_{- 0.16}$ ($M_{\rm sub}/M_{\rm star}=21.1^{+7.4}_{-7.7}$).

In this paper we present a large database of weak lensing light cones constructed using different snapshots from the Big MultiDark simulation (BigMDPL). The ray-tracing through different multiple planes has been performed with the GLAMER code accounting both for single source redshifts and for sources distributed along the cosmic time. This first paper presents weak lensing forecasts and results according to the geometry of the VIPERS-W1 and VIPERS-W4 field of view. Additional fields will be available on our database and new ones can be run upon request. Our database also contains some tools for lensing analysis. In this paper we present results for convergence power spectra, one point and high order weak lensing statistics useful for forecasts and for cosmological studies. Covariance matrices have also been computed for the different realisations of the W1 and W4 fields. In addition we compute also galaxy-shear and projected density contrasts for different halo masses at two lens redshifts according to the CFHTLS source redshift distribution both using stacking and cross-correlation techniques, finding very good agreement.

Observations of the cosmic microwave background indicate that baryons account for 5% of the Universe's total energy content. In the local Universe, the census of all observed baryons falls short of this estimate by a factor of two. Cosmological simulations indicate that the missing baryons might not have condensed into virialized haloes, but reside throughout the filaments of the cosmic web (where matter density is larger than average) as a low-density plasma at temperatures of $10^5-10^7$ kelvin, known as the warm-hot intergalactic medium. There have been previous claims of the detection of warm baryons along the line of sight to distant blazars and of hot gas between interacting clusters. These observations were, however, unable to trace the large-scale filamentary structure, or to estimate the total amount of warm baryons in a representative volume of the Universe. Here we report X-ray observations of filamentary structures of gas at $10^7$ kelvin associated with the galaxy cluster Abell 2744. Previous observations of this cluster were unable to resolve and remove coincidental X-ray point sources. After subtracting these, we reveal hot gas structures that are coherent over scales of 8 mergaparsecs. The filaments coincide with over-densities of galaxies and dark matter, with 5-10% of their mass in baryonic gas. This gas has been heated up by the cluster's gravitational pull and is now feeding its core. Our findings strengthen evidence for a picture of the Universe in which a large fraction of the missing baryons reside in the filaments of the cosmic web.

Based on weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), in this paper we study the mass-concentration ($M$-$c$) relation for $\sim 200$ redMaPPer clusters in the fields. We extract the $M$-$c$ relation by measuring the density profiles of individual clusters instead of using stacked weak lensing signals. By performing Monte Carlo simulations, we demonstrate that although the signal-to-noise ratio for each individual cluster is low, the unbiased $M$-$c$ relation can still be reliably derived from a large sample of clusters by carefully taking into account the impacts of shape noise, cluster center offset, dilution effect from member or foreground galaxies, and the projection effect. Our results show that within error bars the derived $M$-$c$ relation for redMaPPer clusters is in agreement with simulation predictions. There is a weak deviation in that the halo concentrations calibrated by Monte Carlo simulations are somewhat higher than that predicted from ${\it Planck}$ cosmology.

By introducing Dirac ?-function in superhigh magnetic field, we deduce a general formula for pressure of degenerate and relativistic electrons, Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of electrons, and consider the quantum electrodynamic(QED) effects on the equations of states (EOSs) for different matter systems. The main conclusions are as follows: the stronger the magnetic field strength, the higher the electron pressure becomes; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because of the positive contribution of the magnetic field energy to the EOS of the star. Since this is an original work in which some uncertainties could exist, to further modify and perfect our theory model should be considered in our future studies.

Observations show that there is a proton spectral "break" with E$_{break}$ at 1-10MeV in some large CME-driven shocks. Theoretical model usually attribute this phenomenon to a diffusive shock acceleration. However, the underlying physics of the shock acceleration still remains uncertain. Although previous numerical models can hardly predict this "break" due to either high computational expense or shortcomings of current models, the present paper focuses on simulating this energy spectrum in converged two shocks by Monte Carlo numerical method. Considering the Dec 13 2006 CME-driven shock interaction with an Earth bow shock, we examine whether the energy spectral "break" could occur on an interaction between two shocks. As result, we indeed obtain the maximum proton energy up to 10MeV, which is the premise to investigate the existence of the energy spectral "break". Unexpectedly, we further find a proton spectral "break" appears distinctly at the energy $\sim$5MeV.

Observations from multiple spacecraft show that there are energy spectral "breaks" at 1-10MeV in some large CME-driven shocks. However, numerical models can hardly simulate this property due to high computational expense. The present paper focuses on analyzing these energy spectral "breaks" by Monte Carlo particle simulations of an isolated CME-driven shock. Taking the Dec 14 2006 CME-driven shock as an example, we investigate the formation of this energy spectral property. For this purpose, we apply different values for the scattering time in our isolated shock model to obtain the highest energy "tails", which can potentially exceed the "break" energy range. However, we have not found the highest energy "tails" beyond the "break" energy range, but instead find that the highest energy "tails" reach saturation near the range of energy at 5MeV. So, we believe that there exists an energy spectral "cut off" in an isolated shock. If there is no interaction with another shock, there would not be formation of the energy spectral "break" property.

We derived constraints on cosmological parameters using weak lensing peak statistics measured on the $\sim130~{\rm deg}^2$ of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). This analysis demonstrates the feasibility of using peak statistics in cosmological studies. For our measurements, we considered peaks with signal-to-noise ratio in the range of $\nu=[3,6]$. For a flat $\Lambda$CDM model with only $(\Omega_{\rm m}, \sigma_8)$ as free parameters, we constrained the parameters of the following relation $\Sigma_8=\sigma_8(\Omega_{\rm m}/0.27)^{\alpha}$ to be: $\Sigma_8=0.82 \pm 0.03 $ and $\alpha=0.43\pm 0.02$. The $\alpha$ value found is considerably smaller than the one measured in two-point and three-point cosmic shear correlation analyses, showing a significant complement of peak statistics to standard weak lensing cosmological studies. The derived constraints on $(\Omega_{\rm m}, \sigma_8)$ are fully consistent with the ones from either WMAP9 or Planck. From the weak lensing peak abundances alone, we obtained marginalised mean values of $\Omega_{\rm m}=0.38^{+0.27}_{-0.24}$ and $\sigma_8=0.81\pm 0.26$. Finally, we also explored the potential of using weak lensing peak statistics to constrain the mass-concentration relation of dark matter halos simultaneously with cosmological parameters.

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.

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.

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2{\sigma}, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements.

We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections.

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.

The baryon acoustic oscillation (BAO) feature in the two-point correlation function of galaxies supplies a standard ruler to probe the expansion history of the Universe. We study here several galaxy selection schemes, aiming at building an emission-line galaxy (ELG) sample in the redshift range $0.6<z<1.7$, that would be suitable for future BAO studies, providing a highly biased galaxy sample. We analyse the angular galaxy clustering of galaxy selections at the redshifts 0.5, 0.7, 0.8, 1 and 1.2 and we combine this analysis with a halo occupation distribution (HOD) model to derive the properties of the haloes these galaxies inhabit, in particular the galaxy bias on large scales. We also perform a weak lensing analysis (aperture statistics) to extract the galaxy bias and the cross-correlation coefficient and compare to the HOD prediction. We apply this analysis on a data set composed of the photometry of the deep co-addition on Sloan Digital Sky Survey (SDSS) Stripe 82 (225 deg$^2$), of Canda-France-Hawai Telescope/Stripe 82 deep \emph{i}-band weak lensing survey and of the {\it Wide-Field Infrared Survey Explorer}infrared photometric band W1. The analysis on the SDSS-III/constant mass galaxies selection at $z=0.5$ is in agreement with previous studies on the tracer, moreover we measure its cross-correlation coefficient $r=1.16\pm0.35$. For the higher redshift bins, we confirm the trends that the brightest galaxy populations selected are strongly biased ($b>1.5$), but we are limited by current data sets depth to derive precise values of the galaxy bias. A survey using such tracers of the mass field will guarantee a high significance detection of the BAO.

Using several cosmological observations, i.e. the cosmic microwave background anisotropies (WMAP), the weak gravitational lensing (CFHTLS), the measurements of baryon acoustic oscillations (SDSS+WiggleZ), the most recent observational Hubble parameter data, the Union2.1 compilation of type Ia supernovae, and the HST prior, we impose constraints on the sum of neutrino masses ($\mnu$), the effective number of neutrino species ($\neff$) and dark energy equation of state ($w$), individually and collectively. We find that a tight upper limit on $\mnu$ can be extracted from the full data combination, if $\neff$ and $w$ are fixed. However this upper bound is severely weakened if $\neff$ and $w$ are allowed to vary. This result naturally raises questions on the robustness of previous strict upper bounds on $\mnu$, ever reported in the literature. The best-fit values from our most generalized constraint read $\mnu=0.556^{+0.231}_{-0.288}\rm eV$, $\neff=3.839\pm0.452$, and $w=-1.058\pm0.088$ at 68% confidence level, which shows a firm lower limit on total neutrino mass, favors an extra light degree of freedom, and supports the cosmological constant model. The current weak lensing data are already helpful in constraining cosmological model parameters for fixed $w$. The dataset of Hubble parameter gains numerous advantages over supernovae when $w=-1$, particularly its illuminating power in constraining $\neff$. As long as $w$ is included as a free parameter, it is still the standardizable candles of type Ia supernovae that play the most dominant role in the parameter constraints.

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We study the 64 square degrees W1 field, the largest of the CFHTLS-Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet made. 2.66 million galaxy shapes are measured, using a KSB pipeline verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio {\nu}>3.5, consistent with predictions of a {\Lambda}CDM model. Of these peaks, 126 lie within 3.0' of a BCG identified from multicolor optical imaging in an independent, red sequence survey. We also identify 7 counterparts for massive clusters previously seen in X-ray emission within 6 square degrees XMM-LSS survey. With photometric redshift estimates for the source galaxies, we use a tomographic lensing method to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we confirm 85 groups/clusters with \chi^2 reduced < 3.0, at a mean redshift <z_c> = 0.36 and velocity dispersion <\sigma_c> = 658.8 km/s. Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to apply these same techniques to map clusters in much larger volumes and thus tightly constrain cosmological models.

Taking into account the noise from intrinsic ellipticities of source galaxies, we study the efficiency and completeness of halo detections from weak lensing convergence maps. Particularly, with numerical simulations, we compare the Gaussian filter with the so called MRLens treatment based on the modification of the Maximum Entropy Method. For a pure noise field without lensing signals, a Gaussian smoothing results a residual noise field that is approximately Gaussian in statistics if a large enough number of galaxies are included in the smoothing window. On the other hand, the noise field after the MRLens treatment is significantly non-Gaussian, resulting complications in characterizing the noise effects. Considering weak-lensing cluster detections, although the MRLens treatment effectively deletes false peaks arising from noise, it removes the real peaks heavily due to its inability to distinguish real signals with relatively low amplitudes from noise in its restoration process. The higher the noise level is, the larger the removal effects are for the real peaks. For a survey with a source density n_g~30 arcmin^(2), the number of peaks found in an area of 3x3 sq.deg after MRLens filtering is only ~50 for the detection threshold kappa=0.02, while the number of halos with M>5x10^{13} M_{\odot} and with redshift z<=2 in the same area is expected to be ~530. For the Gaussian smoothing treatment, the number of detections is ~260, much larger than that of the MRLens. The Gaussianity of the noise statistics in the Gaussian smoothing case adds further advantages for this method to circumvent the problem of the relatively low efficiency in weak-lensing cluster detections. Therefore, in studies aiming to construct large cluster samples from weak-lensing surveys, the Gaussian smoothing method performs significantly better than the MRLens.

Taking into account noise from intrinsic ellipticities of source galaxies, in this paper, we study the peak statistics in weak-lensing convergence maps around clusters of galaxies and beyond. We emphasize how the noise peak statistics is affected by the density distribution of nearby clusters, and also how cluster-peak signals are changed by the existence of noise. These are the important aspects to be understood thoroughly in weak-lensing analyses for individual clusters as well as in cosmological applications of weak-lensing cluster statistics. We adopt Gaussian smoothing with the smoothing scale $\theta_G=0.5\hbox{ arcmin}$ in our analyses. It is found that the noise peak distribution near a cluster of galaxies depends sensitively on the density profile of the cluster. For a cored isothermal cluster with the core radius $R_c$, the inner region with $R\le R_c$ appears noisy containing on average $\sim 2.4$ peaks with $\nu\ge 5$ for $R_c= 1.7\hbox{ arcmin}$ and the true peak height of the cluster $\nu=5.6$, where $\nu$ denotes the convergence signal to noise ratio. For a NFW cluster of the same mass and the same central $\nu$, the average number of peaks with $\nu\ge 5$ within $R\le R_c$ is $\sim 1.6$. Thus a high peak corresponding to the main cluster can be identified more cleanly in the NFW case. In the outer region with $R_c<R\le 5R_c$, the number of high noise peaks is considerably enhanced in comparison with that of the pure noise case without the nearby cluster. (abridged)

Using the Hyperz code and a template spectral library which consists of 4 observed galaxy spectra from Coleman, Wu & Weedman (CWW, 1980) and 8 spectral families built with evolutionary population synthesis models, we present photometric redshift estimates (photo-z) for a spectroscopic sample of 6,531 galaxies, and morphologies for a morphological sample of 1,502 bright galaxies. All galaxies are matched with the SDSS DR7 and GALEX DR4. The inclusion of Fuv or Nuv or both photometry decreases the number of catastrophic identifications (CIs, |z_phot -z_spec| > 1.0). If CIs are removed, the inclusion of both Fuv and Nuv photometry mainly increases the number of non-CIs in the low redshift, g-r < 0.8 and fainter r-magnitude regions. The inclusion of binary interactions (BIs) mainly increases the number of non-CIs and decreases the deviations in the 0.3 < g-r < 0.8 region in the case of only using optical photometry. The inclusion of UV photometry would decrease and increase the probability that early types are classified as Burst and E types, respectively, and increase that late types are classified as CWW-Sbc and CWW-Scd types. If CIs are excluded, the inclusion of UV data mainly raises the identifications of late types in all redshift, bluer g-r and r > 14 regions. Moreover, BIs mainly affect the determinations of E and S0 types. Nuv -u = 1.94 and 5.77-1.47(u-r) = Fuv discriminators can be used as morphology selection indicators. These two criteria have comparable reliability and completeness for selecting early- and late-type galaxies to C=2.6 criterion and higher completeness for early-type selection than u-r=2.22 criterion.

Recent studies of lensing clusters reveal that it might be fairly common for a galaxy cluster that the X-ray center has an obvious offset from its gravitational center which is measured by strong lensing. We argue that if these offsets exist, then X-rays and lensing are indeed measuring different regions of a cluster, and may thus naturally result in a discrepancy in the measured gravitational masses by the two different methods. Here we investigate theoretically the dynamical effects of such lensing-X-ray offsets, and compare with observational data. We find that for typical values, the offset alone can give rise to a factor of two difference between the lensing and X-ray determined masses for the core regions of a cluster, suggesting that such "offset effect" may play an important role and should not be ignored in our dynamical measurements of clusters.

We compile a sample of 38 galaxy clusters which have both X-ray and strong lensing observations, and study for each cluster the projected offset between the dominant component of baryonic matter center (measured by X-rays) and the gravitational center (measured by strong lensing). Among the total sample, 45% clusters have offsets >10". The >10" separations are significant, considering the arcsecond precision in the measurement of the lensing/X-ray centers. This suggests that it might be a common phenomenon in unrelaxed galaxy clusters that gravitational field is separated spatially from the dominant component of baryonic matter. It also has consequences for lensing models of unrelaxed clusters since the gas mass distribution may differ from the dark matter distribution and give perturbations to the modeling. Such offsets can be used as a statistical tool for comparison with the results of Lambda-CDM simulations and to test the modified dynamics.

Quasar outflows may play important role in the evolution of its host galaxy and central black hole. In this paper, we present a detailed study of multiple outflows in the obscured ultra-luminous infrared quasar Q1321+058. The outflows reveal themselves in the complex optical and UV emission line spectrum, with a broad component blueshifted by 1650 km/s and a narrow component by 360 km/s, respectively.The higher velocity component shows ever strong N III] and strong Si III], in addition to strong [O III]5007 and [Ne III]3869 emission lines, suggesting an overabundance of N and Si relative to C. The abundance pattern is consistent a fast chemical enriching process associated with a recent starburst. The outflow extends to several tens to hundred parsecs from the quasar, and covers only a very small sky. We find that the outflow with line emitting gas is energetically insufficient to remove the ISM of the host galaxy. The velocity range and the column density suggest that the outflow might be part of the low ionization broad absorption line region as seen in a small class of quasars. The optical and UV continuum is starlight-dominated and can be modelled with a young-aged (1 Myr) plus an intermediate-aged (~0.5-1 Gyr) stellar populations, suggesting a fast building of the stellar mass in the host galaxy, consistent with the starburst-type metal abundances inferred from the high velocity outflow spectrum. The broad band spectral energy distribution shows that it is an obscured quasar with its bulk emission in the middle infrared. The star formation rate, independently estimated from UV, far-infrared, and emission line luminosity, is much lower than that is required for the co-evolution of the black hole and its host spheroid.

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), $\gamma_{N}\sim 10^{-4}-10^{-3}$. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, a single filament inclined by $\approx 45^\circ$ from the line of sight can cause substantial distortion of background sources pointing towards the filament's axis ($\kappa=\gamma=(1-A^{-1})/2\sim 0.01$); this is rigorous for infinitely long uniform filaments, but also qualitatively true for short filaments ($\sim 30$Mpc), and even in regions where the projected matter density of the filament is equal to zero. Since galaxies and galaxy clusters are generally embedded in filaments or are projected on such structures, this contribution complicates the interpretation of the weak lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous weak lensing signals, e.g. the "bullet cluster" of Clowe et al., the "cosmic train wreck" of Abell 520 from Mahdavi et al., and the "dark clusters" of Erben et al., filamentary structures might contribute in a significant and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing signal could, in principle, be used to falsify MOND/TeVeS and its variations.

Strong gravitational lensing by galaxies in MOdified Newtonian Dynamics (MOND) has until now been restricted to spherically symmetric models. These models were able to account for the size of the Einstein ring of observed lenses, but were unable to account for double-imaged systems with collinear images, as well as four-image lenses. Non-spherical models are generally cumbersome to compute numerically in MOND, but we present here a class of analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged systems. We use them to obtain a reasonable MOND fit to ten double-imaged systems, as well as to the quadruple-imaged system Q2237+030 which is an isolated bulge-disc lens producing an Einstein cross. However, we also find five double-imaged systems and three quadruple-imaged systems for which no reasonable MOND fit can be obtained with our models. We argue that this is mostly due to the intrinsic limitation of the analytic models, even though the presence of small amounts of additional dark mass on galaxy scales in MOND is also plausible.

In a galaxy cluster, galaxies are mostly collisionless particles in recent epoches. They resemble collisionless cold dark matter particles in some way. Therefore, the spatial distributions of dark matter and cluster galaxies might be expected to possess similar features in the gravitational potential of a cluster. Here we use the galaxy distribution in cluster Cl0024+17 to probe for the ringlike dark matter structure recently discovered by means of strong and weak lensing observations. The galaxies are taken from the catalog of Czoske et al., which contains 650 objects with measured redshifts, of which ~300 galaxies have redshifts in the range 0.37<z<0.41 (and are therefore probable cluster members). We find that, at about the 3-sigma level, the ringlike structure seen in the dark matter measurement is not observed in the projected two-dimensional galaxy distribution.

We report the discovery of superstrong, fading, high-ionization iron line emission in the galaxy SDSSJ095209.56+214313.3 (SDSSJ0952+2143 hereafter), which must have been caused by an X-ray outburst of large amplitude. SDSSJ0952+2143 is unique in its strong multiwavelength variability; such a broadband emission-line and continuum response has not been observed before. The strong iron line emission is accompanied by unusual Balmer line emission with a broad base, narrow core and double-peaked narrow horns, and strong HeII emission. These lines, while strong in the SDSS spectrum taken in 2005, have faded away significantly in new spectra taken in December 2007. Comparison of SDSS, 2MASS, GALEX and follow-up GROND photometry reveals variability in the NUV, optical and NIR band. Taken together, these unusual observations can be explained by a giant outburst in the EUV--X-ray band, detected even in the optical and NIR. The intense and variable iron, Helium and Balmer lines represent the ``light echo'' of the flare, as it traveled through circumnuclear material. The outburst may have been caused by the tidal disruption of a star by a supermassive black hole. Spectroscopic surveys such as SDSS are well suited to detect emission-line light echoes of such rare flare events. Reverberation-mapping of these light echoes can then be used as a new and efficient probe of the physical conditions in the circumnuclear material in non-active or active galaxies.

We summarize the status of Modified Newtonian Dynamics (MOND) in galaxy clusters. The observed acceleration is typically larger than the acceleration threshold of MOND in the central regions, implying that some dark matter is necessary to explain the mass discrepancy there. A plausible resolution of this issue is that the unseen mass in MOND is in the form of ordinary neutrinos with masses just below the experimentally detectable limit. In particular, we show that the lensing mass reconstructions of the clusters 1E0657-56 (the bullet cluster) and Cl0024+17 (the ring) do not pose a new challenge to this scenario. However, the mass discrepancy for cool X-ray emitting groups, in which neutrinos cannot cluster, pose a more serious problem, meaning that dark baryons could present a more satisfactory solution to the problem of unseen mass in MOND clusters.

We present an analysis of the emission line property and the broad band spectral energy distribution of the ultra-luminous infrared Type II quasar Q132+058. The optical and ultraviolet emission lines show four distinct components: a LINER-like component at the systematic velocity, a heavily reddened HII-like component blueshifted by about 400 km/s, and two broad components blueshifted by about 400 km/s and 1900 km/s, respectively. The emission line ratios suggest that broad components are produced in dense and alpha-enriched outflows with a metalicity of 3-5 Zsun and a density of n_H ~ 10^7 cm^-3. The optical--UV continuum is dominated by stellar light and can be modeled with a young (<1 Myr) plus an intermediate age (0.5-0.8Gyr) stellar populations. The near to mid-infrared light is dominated by hot and warm dust heated by the hidden quasar. We derive a star formation rate (SFR) of 300-500 M\sun yr$^{-1}$ from the UV spectrum and far-infrared luminosity, which is two orders of magnitude larger than that indicated by reddening uncorrected [OII] luminosity.

We report the identification of 2MASX J032441.19+341045.9 (hereafter 2MASX J0324+3410) with an appealing object which shows the dual properties of both a narrow line Seyfert 1 galaxy (NLS1) and a blazar. Its optical spectrum, which has a H\beta line width about 1600 km s^-1 (FWHM), an [OIII] to H\beta line ratio ~0.12, and strong FeII emission, clearly fulfills the conventional definition of NLS1s. On the other hand, 2MASX J0324+3410 also exhibits some behavior which is characteristic of blazars, including a flat radio spectrum above 1 GHz, a compact core plus a one-sided jet structure on mas-scale at 8.4 GHz, highly variable fluxes in the radio, optical, and X-ray bands, and a possible detection of TeV gamma-ray emission. On its optical image, obtained with the HST WFPC2, the active nucleus is displaced from the center of the host galaxy, which exhibits an apparent one-armed spiral structure extended to 16 kpc. The remarkable hybrid behavior of this object presents a challenge to current models of NLS1 galaxies and $\gamma$-ray blazars.

We have conducted a systematic search of AGNs with IMBHs from the SDSS DR4. As results we found 245 candidates of broad-line AGN with M_{BH}<10^6 \Msun estimated from the luminosity and width of the broad Halpha component. Compared to the pioneer Greene & Ho (2004) sample of 19 IMBH AGNs, our sample has improved in covering a larger range of the Eddington ratio, as well as black hole mass and redshift, taking the advantage of our AGN-galaxy spectral decomposition algorithm. Among these, thirty-six have L_{bol}/L_{Edd} < 0.1, hinting that a significant fraction of IMBHs might exist with weak or no nuclear activity.

We report the discovery of a dwarf Seyfert 1 active galactic nucleus (AGN) with a candidate intermediate-mass black hole hosted by the dwarf galaxy SDSS J160531.84+174826.1 at z=0.032. A broad component of the H-alpha line with FWHM=781 km/s is detected in its optical spectrum, and a bright, point-like nucleus is evident from a HST imaging observation. Non-thermal X-ray emission is also detected from the nucleus. The black hole mass, as estimated from the luminosity and width of the broad H-alpha component, is about 7x10^4 \msun. The host galaxy appears to be a disk galaxy with a boxy bulge or nuclear bar; with an absolute magnitude of M_R = -17.8 (M_B = -16.4), it is among the least luminous host galaxies ever identified for a Seyfert 1.

We develop a new method to predict the density associated with weak lensing maps of (un)relaxed clusters in a range of theories interpolating between GR and MOND (General Relativity and Modified Newtonian Dynamics). We apply it to fit the lensing map of the bullet merging cluster 1E0657-56, in order to constrain more robustly the nature and amount of collisionless matter in clusters {\it beyond} the usual assumption of spherical equilibrium (Pointecouteau & Silk 2005) and the validity of GR on cluster scales (Clowe et al. 2006). Strengthening the proposal of previous authors we show that the bullet cluster is dominated by a collisionless -- most probably non-baryonic -- component in GR as well as in MOND, a result consistent with the dynamics of many X-ray clusters. Our findings add to the number of known pathologies for a purely baryonic MOND, including its inability to fit the latest data from the Wilkinson Microwave Anisotropy Probe. A plausible resolution of all these issues and standard issues of Cold Dark Matter with galaxy rotation curves is the "marriage" of MOND with ordinary hot neutrinos of 2eV. This prediction is just within the GR-independent maximum of neutrino mass from current $\beta$-decay experiments, and is falsifiable by the Karlsruhe Tritium Neutrino (KATRIN) experiment by 2009. Issues of consistency with strong lensing arcs and the large relative velocity of the two clusters comprising the bullet cluster are also addressed.