Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution.

We present improved modelling of the redshift-space distortions of galaxy clustering that arise from peculiar velocities. We create mock galaxy catalogues in the framework of the halo model, using data from the Bolshoi project. These mock galaxy populations are inserted into the haloes with additional degrees of freedom that govern spatial and kinematical biases of the galaxy populations relative to the dark matter. We explore this generalised halo model with an MCMC algorithm, comparing the predictions to data from the Galaxy And Mass Assembly (GAMA) survey, and thus derive one of the first constraints on the detailed kinematic degrees of freedom for satellite galaxies within haloes. With this approach, the distortions of the redshift-space galaxy autocorrelations can be accounted for down to spatial separations close to 10 kpc, opening the prospect of improved RSD measurements of the perturbation growth rate by the inclusion of data from nonlinear scales.

The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9<z<3.5$) Lyman-$\alpha$ line emitters. We demonstrate that the joint-analysis method yields a high signal-to-noise ratio measurement of the interloper fractions while only marginally increasing the uncertainties in the cosmological parameters relative to the case without interlopers. We also show the same is true for the high-latitude spectroscopic survey of Wide-Field Infrared Survey Telescope (WFIRST) mission where contamination occurs between the Balmer-$\alpha$ line emitters at lower redshifts ($1.1<z<1.9$) and Oxygen ([O III] $\lambda5007${\AA}) line emitters at higher redshifts ($1.7<z<2.8$).

We directly constrain the non-linear alignment (NLA) model of intrinsic galaxy alignments, analysing the most representative and complete flux-limited sample of spectroscopic galaxies available for cosmic shear surveys. We measure the projected galaxy position-intrinsic shear correlations and the projected galaxy clustering signal using high-resolution imaging from the Kilo Degree Survey (KiDS) overlapping with the GAMA spectroscopic survey, and data from the Sloan Digital Sky Survey. Separating samples by colour, we make no significant detection of blue galaxy alignments, constraining the blue galaxy NLA amplitude $A_{\textrm{IA}}^{\textrm{B}}=0.21^{+0.37}_{-0.36}$ to be consistent with zero. We make robust detections ($\sim9\sigma$) for red galaxies, with $A_{\textrm{IA}}^{\textrm{R}}=3.18^{+0.47}_{-0.46}$, corresponding to a net radial alignment with the galaxy density field, and we find no evidence for any scaling of alignments with galaxy luminosity. We provide informative priors for current and future weak lensing surveys, an improvement over de facto wide priors that allow for unrealistic levels of intrinsic alignment contamination. For a colour-split cosmic shear analysis of the final KiDS survey area, we forecast that our priors will improve the constraining power on $S_{8}$ and the dark energy equation of state $w_{0}$, by up to $62\%$ and $51\%$, respectively. Our results indicate, however, that the modelling of red/blue-split galaxy alignments may be insufficient to describe samples with variable central/satellite galaxy fractions.

This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

Statistical studies of galaxy-galaxy interactions often utilise net change in physical properties of progenitors as a function of the separation between their nuclei to trace both the strength and the observable timescale of their interaction. In this study, we use two-point auto, cross and mark correlation functions to investigate the extent to which small-scale clustering properties of star forming galaxies can be used to gain physical insight into galaxy-galaxy interactions between galaxies of similar optical brightness and stellar mass. The Halpha star formers, drawn from the highly spatially complete Galaxy And Mass Assembly (GAMA) survey, show an increase in clustering on small separations. Moreover, the clustering strength shows a strong dependence on optical brightness and stellar mass, where (1) the clustering amplitude of optically brighter galaxies at a given separation is larger than that of optically fainter systems, (2) the small scale clustering properties (e.g. the strength, the scale at which the signal relative to the fiducial power law plateaus) of star forming galaxies appear to differ as a function of increasing optical brightness of galaxies. According to cross and mark correlation analyses, the former result is largely driven by the increased dust content in optically bright star forming galaxies. The latter could be interpreted as evidence of a correlation between interaction-scale and optical brightness of galaxies, where physical evidence of interactions between optically bright star formers, likely hosted within relatively massive halos, persist over larger separations than those between optically faint star formers.

We measure the projected galaxy clustering and galaxy-galaxy lensing signals using the Galaxy And Mass Assembly (GAMA) survey and Kilo-Degree Survey (KiDS) to study galaxy bias. We use the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics to constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing. The bias function $\Gamma_{\text{gm}}(r_{\text{p}})$, where $r_{\text{p}}$ is the projected comoving separation, is evaluated using the analytical halo model from which the scale dependence of $\Gamma_{\text{gm}}(r_{\text{p}})$, and the origin of the non-linearity and stochasticity in halo occupation models can be inferred. Our observations unveil the physical reason for the non-linearity and stochasticity, further explored using hydrodynamical simulations, with the stochasticity mostly originating from the non-Poissonian behaviour of satellite galaxies in the dark matter haloes and their spatial distribution, which does not follow the spatial distribution of dark matter in the halo. The observed non-linearity is mostly due to the presence of the central galaxies, as was noted from previous theoretical work on the same topic. We also see that overall, more massive galaxies reveal a stronger scale dependence, and out to a larger radius. Our results show that a wealth of information about galaxy bias is hidden in halo occupation models. These models should therefore be used to determine the influence of galaxy bias in cosmological studies.

We present cosmological parameter constraints from a joint analysis of three cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$ of data from the Kilo Degree Survey (KiDS), the galaxy-matter cross-correlation signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined with KiDS weak lensing, and the angular correlation function of the same GAMA galaxies. We use fast power spectrum estimators that are based on simple integrals over the real-space correlation functions, and show that they are practically unbiased over relevant angular frequency ranges. We test our full pipeline on numerical simulations that are tailored to KiDS and retrieve the input cosmology. By fitting different combinations of power spectra, we demonstrate that the three probes are internally consistent. For all probes combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.800_{-0.027}^{+0.029}$, consistent with Planck and the fiducial KiDS-450 cosmic shear correlation function results. Marginalising over wide priors on the mean of the tomographic redshift distributions yields consistent results for $S_8$ with an increase of $28\%$ in the error. The combination of probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the cosmic shear power spectra alone. The main gain from these additional probes comes through their constraining power on nuisance parameters, such as the galaxy intrinsic alignment amplitude or potential shifts in the redshift distributions, which are up to a factor of two better constrained compared to using cosmic shear alone, demonstrating the value of large-scale structure probe combination.

We describe data release 3 (DR3) of the Galaxy And Mass Assembly (GAMA) survey. The GAMA survey is a spectroscopic redshift and multi-wavelength photometric survey in three equatorial regions each of 60.0 deg^2 (G09, G12, G15), and two southern regions of 55.7 deg^2 (G02) and 50.6 deg^2 (G23). DR3 consists of: the first release of data covering the G02 region and of data on H-ATLAS sources in the equatorial regions; and updates to data on sources released in DR2. DR3 includes 154809 sources with secure redshifts across four regions. A subset of the G02 region is 95.5% redshift complete to r<19.8 over an area of 19.5 deg^2, with 20086 galaxy redshifts, that overlaps substantially with the XXL survey (X-ray) and VIPERS (redshift survey). In the equatorial regions, the main survey has even higher completeness (98.5%), and spectra for about 75% of H-ATLAS filler targets were also obtained. This filler sample extends spectroscopic redshifts, for probable optical counterparts to H-ATLAS sub-mm sources, to 0.8 mag deeper (r<20.6) than the GAMA main survey. There are 25814 galaxy redshifts for H-ATLAS sources from the GAMA main or filler surveys. GAMA DR3 is available at the survey website (www.gama-survey.org/dr3/).

We study present-day galaxy clustering in the EAGLE cosmological hydrodynamical simulation. EAGLE's galaxy formation parameters were calibrated to reproduce the redshift $z=0.1$ galaxy stellar mass function, and the simulation also reproduces galaxy colours well. The simulation volume is too small to correctly sample large-scale fluctuations and we therefore concentrate on scales smaller than a few mega parsecs. We find very good agreement with observed clustering measurements from the Galaxy And Mass Assembly (GAMA) survey, when galaxies are binned by stellar mass, color, or luminosity. However, low-mass red-galaxies are clustered too strongly, which is at least partly due to limited numerical resolution. Apart from this limitation, we conclude that EAGLE galaxies inhabit similar dark matter haloes as observed GAMA galaxies, and that the radial distribution of satellite galaxies as function of stellar mass and colour is similar to that observed as well.

We investigate possible signatures of halo assembly bias for spectroscopically selected galaxy groups from the GAMA survey using weak lensing measurements from the spatially overlapping regions of the deeper, high-imaging-quality photometric KiDS survey. We use GAMA groups with an apparent richness larger than 4 to identify samples with comparable mean host halo masses but with a different radial distribution of satellite galaxies, which is a proxy for the formation time of the haloes. We measure the weak lensing signal for groups with a steeper than average and with a shallower than average satellite distribution and find no sign of halo assembly bias, with the bias ratio of $0.85^{+0.37}_{-0.25}$, which is consistent with the $\Lambda$CDM prediction. Our galaxy groups have typical masses of $10^{13} M_{\odot}/h$, naturally complementing previous studies of halo assembly bias on galaxy cluster scales.

We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.

We use the WISE-2MASS infrared galaxy catalog matched with Pan-STARRS1 (PS1) galaxies to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. Our imaging catalog has median redshift $z\simeq 0.14$, and we obtain photometric redshifts from PS1 optical colours to create a tomographic map of the galaxy distribution. The radial profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. The counts in photometric redshift bins show significantly low densities at high detection significance, $\gtrsim 5 \sigma$ and $\gtrsim 6 \sigma$, respectively, for the two fiducial radii. The line-of-sight position of the deepest region of the void is $z\simeq 0.15-0.25$. Our data, combined with an earlier measurement by Granett et al. 2010, are consistent with a large $R_{\rm void}=(220 \pm 50) h^{-1}Mpc $ supervoid with $\delta_{m} \simeq -0.14 \pm 0.04$ centered at $z=0.22\pm0.03$. Such a supervoid, constituting at least a $\simeq 3.3\sigma$ fluctuation in a Gaussian distribution of the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

The Pan-STARRS1 (PS1) survey has obtained imaging in 5 bands (grizy_P1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering new superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of M_AB ~ 23.5, of which 76 were classified as Type Ia SNe. There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of Type Ic at redshifts between 0.5-1.4. This illustrates that the discovery rate of Type Ic SLSNe can be maximised by concentrating on hostless transients and removing normal SNe Ia. We present data for two new possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of SLSNe-Ic to be between 3^{+3}_{-2} * 10^{-5} and 8^{+2}_{-1} * 10^{-5} of the CCSNe rate within 0.3 <= z <= 1.4 by applying a Monte-Carlo technique. The rate of slowly evolving, SN2007bi-like explosions is estimated as a factor of 10 lower than this range.

We investigate the dependence of the galaxy luminosity function on geometric environment within the Galaxy And Mass Assembly (GAMA) survey. The tidal tensor prescription, based on the Hessian of the pseudo-gravitational potential, is used to classify the cosmic web and define the geometric environments: for a given smoothing scale, we classify every position of the surveyed region, $0.04<{z}<0.26$, as either a void, a sheet, a filament or a knot. We consider how to choose appropriate thresholds in the eigenvalues of the Hessian in order to partition the galaxies approximately evenly between environments. We find a significant variation in the luminosity function of galaxies between different geometric environments; the normalisation, characterised by $\phi^{*}$ in a Schechter function fit, increases by an order of magnitude from voids to knots. The turnover magnitude, characterised by $M^*$, brightens by approximately $0.5$ mag from voids to knots. However, we show that the observed modulation can be entirely attributed to the indirect local-density dependence. We therefore find no evidence of a direct influence of the cosmic web on the galaxy luminosity function.

We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc $ $(2\sigma)$ supervoid with $\delta \simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 \sigma$ fluctuation in the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and the leading explanation of it posits a large void that imprints this extremely cold area via the linear Integrated Sachs-Wolfe (ISW) effect due to the decay of gravitational potentials over cosmic time, or via the Rees-Sciama (RS) effect due to late-time non-linear evolution. Despite several observational campaigns targeting the Cold Spot region, to date no suitably large void was found at higher redshifts $z > 0.3$. Here we report the detection of an $R =(192 \pm 15) h^{-1}Mpc$ size supervoid of depth $\delta = -0.13 \pm 0.03$, and centred at redshift $z = 0.22$. This supervoid, possibly the largest ever found, is large enough to significantly affect the CMB via the non-linear RS effect, as shown in our Lemaitre-Tolman-Bondi framework. This discovery presents the first plausible explanation for any of the physical CMB anomalies, and raises the possibility that local large-scale structure could be responsible for other anomalies as well.

The Pan-STARRS1 survey is acquiring multi-epoch imaging in 5 bands (grizy) over the entire sky north of declination -30deg (the $3\pi$ survey). In July 2011 a test area of about 70 sq.deg. was observed to the expected final depth of the main survey. In this, the first of a series of papers targetting the galaxy count and clustering properties of the combined multi-epoch test area data, we present a detailed investigation into the depth of the survey and the reliability of the Pan-STARRS1 analysis software. We show that the Pan-STARRS1 reduction software can recover the properties of fake sources, and show good agreement between the magnitudes measured by Pan-STARRS1 and those from Sloan Digital Sky Survey. We also examine the number of false detections apparent in the Pan-STARRS1 data. Our comparisons show that the test area survey is somewhat deeper than the Sloan Digital Sky Survey in all bands, and, in particular, the z band approaches the depth of the stacked Sloan Stripe 82 data.

The Pan-STARRS1 survey is currently obtaining imaging in 5 bands (grizy) for the $3\pi$ steradian survey, one of the largest optical surveys ever conducted. The finished survey will have spatially varying depth, due to the survey strategy. This paper presents a method to correct galaxy number counts and galaxy clustering for this potential systematic based on a simplified signal to noise measurement. A star and galaxy separation method calibrated using realistic synthetic images is also presented, along with an approach to mask bright stars. By using our techniques on a ~69 sq. degree region of science verification data this paper shows PS1 measurements of the two point angular correlation function as a function of apparent magnitude agree with measurements from deeper, smaller surveys. Clustering measurements appear reliable down to a magnitude limit of rps<22.5. Additionally, stellar contamination and false detection issues are discussed and quantified. This work is the second of two papers which pave the way for the exploitation of the full $3\pi$ survey for studies of large scale structure.