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.

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 simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the 4th data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 \pm 0.39$ dex (68% confidence, $3\sigma$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $\log(M_{\star}/M_{\odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness, and does not provide any improvement in precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection.

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.

We cross-correlate maps of the thermal Sunyaev-Zeldovich (tSZ) Compton-$y$ parameter published by Planck with the projected distribution of galaxies in a set of low-redshift tomographic bins. We use the nearly full-sky 2MASS Photometric Redshift and WISE $\times$ SuperCOSMOS public catalogues, covering the redshift range $z\lesssim0.4$. Our measurements allow us to place constraints on the redshift dependence of the mass-observable relation for tSZ cluster count analyses in terms of the so-called 'hydrostatic mass bias' parameter $1-b_{\rm H}$. These results can also be interpreted as measurements of the bias-weighted average gas pressure $\langle bP_e\rangle$ as a function of redshift, a quantity that can be related to the thermodynamics of gas inside haloes and used to constrain energy injection processes. We measure $1-b_{\rm H}$ with $\sim13\%$ precision in 6 equispaced redshift bins, and find no evidence for a redshift-dependent mass bias parameter, in agreement with previous analyses. Our mean value of $1-b_{\rm H} = 0.59\pm0.03$ is also in good agreement with the one estimated by the joint analysis of Planck cluster counts and CMB anisotropies. Our measurements of $\langle bP_e\rangle$, at the level of $\sim10\%$ in each bin, are the most stringent constraints on the redshift dependence of this parameter to date, and agree well both with previous measurements and with theoretical expectations from shock-heating models.

Dust lanes bisect the plane of a typical edge-on spiral galaxy as a dark optical absorption feature. Their appearance is linked to the gravitational stability of spiral disks; the fraction of edge-on galaxies that displays a dust lane is a direct indicator of the typical vertical balance between gravity and turbulence; a balance struck between the energy input from star-formation and the gravitational pull into the plane of the disk. Based on morphological classifications by the Galaxy~Zoo project on the Kilo-Degree Survey (KiDS) imaging data in the Galaxy and Mass Assembly (GAMA) fields, we explore the relation of dust lanes to the galaxy characteristics, most of which were determined using the {\sc magphys} spectral energy distribution fitting tool: stellar mass, total and specific star-formation rates, and several parameters describing the cold dust component. We find that the fraction of dust lanes does depend on the stellar mass of the galaxy; they start to appear at $M^* \sim 10^9 M_\odot$. A dust lane also implies strongly a dust mass of at least $10^5 M_\odot$, but otherwise does not correlate with cold dust mass parameters of the {\sc magphys} spectral energy distribution analysis, nor is there a link with star-formation rate, specific or total. Dust lane identification does not depend on disk ellipticity (disk thickness) or Sersic profile but correlates with bulge morphology; a round bulge favors dust lane votes. The central component along the line of sight that produces the dust lane is not associated with either one of the components fit by {\sc magphys}, the cold diffuse component or the localized, heated component in HII regions, but a mix of these two.

Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statistically large samples of galaxies, totalling over 200\,000, over a large redshift range, 0.0 to 4.0. We train and use convolutional neural networks to create binary merger identifications (merger or non-merger) in the SDSS, KiDS and CANDELS imaging surveys. We then compare the galaxy main sequence subtracted SFR of the merging and non-merging galaxies to determine what effect, if any, a galaxy merger has on SFR. We find that the SFR of merging galaxies are not significantly different from the SFR of non-merging systems. The changes in the average SFR seen in the star forming population when a galaxy is merging are small, of the order of a factor of 1.2. However, the higher the SFR above the galaxy main sequence, the higher the fraction of galaxy mergers. Galaxy mergers have little effect on the SFR of the majority of merging galaxies compared to the non-merging galaxies. The typical change in SFR is less than 0.1~dex in either direction. Larger changes in SFR can be seen but are less common. The increase in merger fraction as the distance above the galaxy main sequence increases demonstrates that galaxy mergers can induce starbursts.

Aiming to extract cosmological information from linear scales of the WISE$\times$SuperCOSMOS photometric redshift catalog, we perform a characterization of the systematic effects associated with stellar content, evidencing the presence of contamination and obscuration. We create an integrated model for these effects (which together we call `usurper contamination'), devise a method to remove both of them simultaneously and show its functionality by applying it to a set of mock catalogs. When administered to WISE$\times$SuperCOSMOS data, our method shows to improve the measurements of angular power spectra on scales $\ell\lesssim15$ and the extraction of cosmological parameters therefrom, even though a significant excess of power remains at these scales. When ignoring scales $\ell<15$, we still find strong indications of systematics, albeit these can be localized in the southern equatorial hemisphere. An independent analysis of the northern hemisphere at $\ell\geq 15$ agrees with a $\Lambda$CDM model with parameters from the Planck satellite and gives $\Omega_{\mathrm{c}}=0.254\pm0.020$ and $\Omega_{\mathrm{b}}<0.065$ at 95% confidence limit when combined with priors on $H_0$, $A_s$ and $n_s$.

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.

We have observed the G23 field of the Galaxy And Mass Assembly (GAMA) survey using the Australian Square Kilometre Array Pathfinder (ASKAP) in its commissioning phase, to validate the performance of the telescope and to characterize the detected galaxy populations. This observation covers $\sim$48 deg$^2$ with synthesized beam of 32.7$^{\prime\prime}$ by 17.8$^{\prime\prime}$ at 936 MHz, and $\sim$39 deg$^2$ with synthesized beam of 15.8$^{\prime\prime}$ by 12.0$^{\prime\prime}$ at 1320 MHz. At both frequencies, the r.m.s. (root-mean-square) noise is $\sim$0.1 mJy/beam. We combine these radio observations with the GAMA galaxy data, which includes spectroscopy of galaxies that are i-band selected with a magnitude limit of 19.2. Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry is used to determine which galaxies host an active galactic nucleus (AGN). In properties including source counts, mass distributions, and IR vs. radio luminosity relation, the ASKAP detected radio sources behave as expected. Radio galaxies have higher stellar mass and luminosity in IR, optical and UV than other galaxies. We apply optical and IR AGN diagnostics and find that they disagree for $\sim$30% of the galaxies in our sample. We suggest possible causes for the disagreement. Some cases can be explained by optical extinction of the AGN, but for more than half of the cases we do not find a clear explanation. Radio sources are more likely ($\sim$6%) to have an AGN than radio quiet galaxies ($\sim$1%), but the majority of AGN are not detected in radio at this sensitivity.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present a catalog of quasars selected from broad-band photometric ugri data of the Kilo-Degree Survey Data Release 3 (KiDS DR3). The QSOs are identified by the random forest (RF) supervised machine learning model, trained on SDSS DR14 spectroscopic data. We first cleaned the input KiDS data from entries with excessively noisy, missing or otherwise problematic measurements. Applying a feature importance analysis, we then tune the algorithm and identify in the KiDS multiband catalog the 17 most useful features for the classification, namely magnitudes, colors, magnitude ratios, and the stellarity index. We used the t-SNE algorithm to map the multi-dimensional photometric data onto 2D planes and compare the coverage of the training and inference sets. We limited the inference set to r<22 to avoid extrapolation beyond the feature space covered by training, as the SDSS spectroscopic sample is considerably shallower than KiDS. This gives 3.4 million objects in the final inference sample, from which the random forest identified 190,000 quasar candidates. Accuracy of 97%, purity of 91%, and completeness of 87%, as derived from a test set extracted from SDSS and not used in the training, are confirmed by comparison with external spectroscopic and photometric QSO catalogs overlapping with the KiDS footprint. The robustness of our results is strengthened by number counts of the quasar candidates in the r band, as well as by their mid-infrared colors available from WISE. An analysis of parallaxes and proper motions of our QSO candidates found also in Gaia DR2 suggests that a probability cut of p(QSO)>0.8 is optimal for purity, whereas p(QSO)>0.7 is preferable for better completeness. Our study presents the first comprehensive quasar selection from deep high-quality KiDS data and will serve as the basis for versatile studies of the QSO population detected by this survey.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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

We investigate the clustering properties of radio sources in the Alternative Data Release 1 of the TIFR GMRT Sky Survey (TGSS), focusing on large angular scales, where previous analyses have detected a large clustering signal. After appropriate data selection, the TGSS sample we use contains ~110,000 sources selected at 150 MHz over ~70% of the sky. The survey footprint is largely superimposed on that of the NRAO VLA Sky Survey (NVSS) with the majority of TGSS sources having a counterpart in the NVSS sample. These characteristics make TGSS suitable for large-scale clustering analyses and facilitate the comparison with the results of previous studies. In this analysis we focus on the angular power spectrum, although the angular correlation function is also computed to quantify the contribution of multiple-component radio sources. We find that on large angular scales, corresponding to multipoles $2 \leq \ell \leq 30$, the amplitude of the TGSS angular power spectrum is significantly larger than that of the NVSS. We do not identify any observational systematic effects that may explain this mismatch. We have produced a number of physically motivated models for the TGSS angular power spectrum and found that all of them fail to match observations, even when taking into account observational and theoretical uncertainties. The same models provide a good fit to the angular spectrum of the NVSS sources. These results confirm the anomalous nature of the TGSS large-scale power, which has no obvious physical origin and seems to indicate that unknown systematic errors are present in the TGSS dataset.

This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-large-scale density fluctuations required to explore primordial non-Gaussianity and therefore inflation. We expect a measurement of the local parameter $f_{NL}$ to a precision $\sigma(f_{NL}) = 1.8$. Combining the MSE High-z Cosmology Survey data with data from a next generation CMB stage 4 experiment and existing DESI data will provide the first $5\sigma$ confirmation of the neutrino mass hierarchy from astronomical observations. In addition, the Baryonic Acoustic Oscillations (BAO) observed within the sample will provide measurements of the distance-redshift relationship in six different redshift bins between $z=1.6$ and 4.0, each with an accuracy of $\sim0.6\%$. The simultaneous measurements of Redshift Space Distortions (RSD) constrain the amplitude of the fluctuations, at a level ranging from $1.9\%$ to $3.6\%$. The proposed survey covers 10,000 ${\rm deg}^2$, measuring redshifts for three classes of target objects: Emission Line Galaxies (ELGs) with $1.6<z<2.4$, Lyman Break Galaxies (LBGs) with $2.4<z<4.0$, and quasars $2.1<z<3.5$. The ELGs and LBGs will be used as direct tracers of the underlying density field, while the Lyman-$\alpha$ forests in the quasar spectra will be utilized to probe structure. Exposures of duration 1,800sec will guarantee a redshift determination efficiency of $90\%$ for ELGS and at least $50\%$ for LBGs. The survey will represent 100 nights per year for a 5-year MSE program. Finally, three ideas for additional projects of cosmological interest are proposed.

The aim of this work is to create a new catalog of reliable AGN candidates selected from the AKARI NEP-Deep field. Selection of the AGN candidates was done by applying a fuzzy SVM algorithm, which allows to incorporate measurement uncertainties into the classification process. The training dataset was based on the spectroscopic data available for selected objects in the NEP-Deep and NEP-Wide fields. The generalization sample was based on the AKARI NEP-Deep field data including objects without optical counterparts and making use of the infrared information only. A high quality catalog of previously unclassified 275 AGN candidates was prepared.

We present a machine learning analysis of five labelled galaxy catalogues from the Galaxy And Mass Assembly (GAMA): The SersicCatVIKING and SersicCatUKIDSS catalogues containing morphological features, the GaussFitSimple catalogue containing spectroscopic features, the MagPhys catalogue including physical parameters for galaxies, and the Lambdar catalogue, which contains photometric measurements. Extending work previously presented at the ESANN 2018 conference - in an analysis based on Generalized Relevance Matrix Learning Vector Quantization and Random Forests - we find that neither the data from the individual catalogues nor a combined dataset based on all 5 catalogues fully supports the visual-inspection-based galaxy classification scheme employed to categorise the galaxies. In particular, only one class, the Little Blue Spheroids, is consistently separable from the other classes. To aid further insight into the nature of the employed visual-based classification scheme with respect to physical and morphological features, we present the galaxy parameters that are discriminative for the achieved class distinctions.

We study projected underdensities in the cosmic galaxy density field known as 'troughs', and their overdense counterparts, which we call 'ridges'. We identify these regions using a bright sample of foreground galaxies from the photometric Kilo-Degree Survey (KiDS), specifically selected to mimic the spectroscopic Galaxy And Mass Assembly survey (GAMA). Using background galaxies from KiDS, we measure the weak gravitational lensing profiles of the troughs/ridges. We quantify the amplitude of their lensing strength $A$ as a function of galaxy density percentile rank $P$ and galaxy overdensity $\delta$, and find that the skewness in the galaxy density distribution is reflected in the total mass distribution measured by weak lensing. We interpret our results using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de l'Espai (MICE) simulation, and find a good agreement with our observations. Using signal-to-noise weights derived from the Scinet LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin, resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select troughs using a volume-limited sample of galaxies, split into two redshift bins between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$ = 1.9 Mpc/h, we find no significant difference in the comoving Excess Surface Density as a function of $P$ and $\delta$ between the low- and high-redshift sample. Using the MICE and SLICS mocks we predict that trough and ridge evolution could be detected with gravitational lensing using deeper and wider lensing surveys, such as those from the Large Synoptic Survey Telescope and Euclid.

We use the new release of the AKARI Far-Infrared all sky Survey matched with the NVSS radio database to investigate the local ($z<0.25$) far infrared-radio correlation (FIRC) of different types of extragalactic sources. To obtain the redshift information for the AKARI FIS sources we crossmatch the catalogue with the SDSS DR8. This also allows us to use emission line properties to divide sources into four categories: i) star-forming galaxies (SFGs), ii) composite galaxies (displaying both star-formation and active nucleus components), iii) Seyfert galaxies, and iv) low-ionization nuclear emission-line region (LINER) galaxies. We find that the Seyfert galaxies have the lowest FIR/radio flux ratios and display excess radio emission when compared to the SFGs. We conclude that FIRC can be used to separate SFGs and AGNs only for the most radio-loud objects.

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.

During its first observing run, in late 2015, the advanced LIGO facility announced 3 gravitational wave (GW) triggers to electromagnetic follow-up partners. Two of these have since been confirmed as being of astrophysical origin: both are binary black hole mergers at ~500 Mpc; the other trigger was later found not to be astrophysical. In this paper we report on the Swift follow up observations of the second and third triggers, including details of 21 X-ray sources detected; none of which can be associated with the GW event. We also consider the challenges that the next GW observing run will bring as the sensitivity and hence typical distance of GW events will increase. We discuss how to effectively use galaxy catalogues to prioritise areas for follow up, especially in the presence of distance estimates from the GW data. We also consider two galaxy catalogues and suggest that the high completeness at larger distances of the 2MASS Photometric Redshift Catalogue (2MPZ) makes it very well suited to optimise Swift follow-up observations.

We use the overlap between multiband photometry of the Kilo-Degree Survey (KiDS) and spectroscopic data based on the Sloan Digital Sky Survey (SDSS) and Galaxy And Mass Assembly (GAMA) to infer the colour-magnitude relation of red-sequence galaxies. We then use this inferred relation to select luminous red galaxies (LRGs) in the redshift range of $0.1<z<0.7$ over the entire KiDS Data Release 3 footprint. We construct two samples of galaxies with different constant comoving densities and different luminosity thresholds. The selected red galaxies have photometric redshifts with typical photo-z errors of $\sigma_z \sim 0.014 (1+z)$ that are nearly uniform with respect to observational systematics. This makes them an ideal set of galaxies for lensing and clustering studies. As an example, we use the KiDS-450 cosmic shear catalogue to measure the mean tangential shear signal around the selected LRGs. We detect a significant weak lensing signal for lenses out to $z \sim 0.7$.

We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $\Lambda$CDM Universe and assuming GR, we find $\Omega_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}.

Despite the high accuracy of photometric redshifts (zphot) derived using Machine Learning (ML) methods, the quantification of errors through reliable and accurate Probability Density Functions (PDFs) is still an open problem. First, because it is difficult to accurately assess the contribution from different sources of errors, namely internal to the method itself and from the photometric features defining the available parameter space. Second, because the problem of defining a robust statistical method, always able to quantify and qualify the PDF estimation validity, is still an open issue. We present a comparison among PDFs obtained using three different methods on the same data set: two ML techniques, METAPHOR (Machine-learning Estimation Tool for Accurate PHOtometric Redshifts) and ANNz2, plus the spectral energy distribution template fitting method, BPZ. The photometric data were extracted from the KiDS (Kilo Degree Survey) ESO Data Release 3, while the spectroscopy was obtained from the GAMA (Galaxy and Mass Assembly) Data Release 2. The statistical evaluation of both individual and stacked PDFs was done through quantitative and qualitative estimators, including a dummy PDF, useful to verify whether different statistical estimators can correctly assess PDF quality. We conclude that, in order to quantify the reliability and accuracy of any zphot PDF method, a combined set of statistical estimators is required.

We describe a tomographic dissection of the Planck CMB lensing data, cross-correlating this map with galaxies in different ranges of photometric redshift. We use the nearly all-sky 2MPZ and WISExSCOS catalogues for z<0.35, extending to z<0.6 using SDSS. We describe checks for consistency between the different datasets, and perform a test for possible leakage of thermal Sunyaev--Zel'dovich signal into our cross-correlation measurements. The amplitude of the cross-correlation allows us to estimate the evolution of density fluctuations as a function of redshift, thus providing a test of theories of modified gravity. Assuming the common parametrisation for the logarithmic growth rate, f_g=Omega_m(z)^gamma, we infer gamma=0.77 +/- 0.18 when Omega_m is fixed using external data. Thus CMB lensing tomography is currently consistent with Einstein gravity, where gamma=0.55 is expected. We discuss how such constraints may be expected to improve with future data.

New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-DR1, ii) direct modeling of KiDS cutouts and iii) positional offsets between different surveys (KiDS-vs-Gaia, Gaia-vs-2MASS), with purpose-built astrometric recalibrations. The first benchmark for the methods has been set by the recovery of known lenses. We are able to recover seven out of ten known lenses and pairs of quasars observed in the KiDS DR3 footprint, or eight out of ten with improved selection criteria and looser colour pre-selection. This success rate reflects the combination of all methods together, which, taken individually, performed significantly worse (four lenses each). One movelty of our analysis is that the comparison of the performances of the different methods has revealed the pros and cons of the approaches and, most of all, the complementarities. We finally provide a list of high-grade candidates found by one or more methods, awaiting spectroscopic follow-up for confirmation. Of these, KiDS 1042+0023 is to our knowledge the first confirmed lensed quasar from KiDS, exhibiting two quasar spectra at the same source redshift at either sides of a red galaxy, with uniform flux-ratio $f\approx1.25$ over the wavelength range $0.45\mu\mathrm{m}<\lambda<0.75\mu\mathrm{m}.$

The Deep Extragalactic VIsible Legacy Survey (DEVILS) is a large spectroscopic campaign at the Anglo-Australian Telescope (AAT) aimed at bridging the near and distant Universe by producing the highest completeness survey of galaxies and groups at intermediate redshifts ($0.3<z<1.0$). Our sample consists of $\sim$60,000 galaxies to Y$<$21.2mag, over $\sim$6deg$^{2}$ in three well-studied deep extragalactic fields (Cosmic Origins Survey field, COSMOS, Extended Chandra Deep Field South, ECDFS and the X-ray Multi-Mirror Mission Large-Scale Structure region, XMM-LSS - all Large Synoptic Survey Telescope deep-drill fields). This paper presents the broad experimental design of DEVILS. Our target sample has been selected from deep Visible and Infrared Survey Telescope for Astronomy (VISTA) Y-band imaging (VISTA Deep Extragalactic Observations, VIDEO and UltraVISTA), with photometry measured by ProFound. Photometric star/galaxy separation is done on the basis of NIR colours, and has been validated by visual inspection. To maximise our observing efficiency for faint targets we employ a redshift feedback strategy, which continually updates our target lists, feeding back the results from the previous night's observations. We also present an overview of the initial spectroscopic observations undertaken in late 2017 and early 2018.

We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $\delta z=-2e-4$ and scatter $\sigma_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $\mu$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $\delta z<4e-5$ and $\sigma_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation.

We derive updated constraints on the Integrated Sachs-Wolfe (ISW) effect through cross-correlation of the cosmic microwave background with galaxy surveys. We improve with respect to similar previous analyses in several ways. First, we use the most recent versions of extragalactic object catalogs: SDSS DR12 photometric redshift (photo-$z$) and 2MASS Photo-$z$ datasets, as well as employed earlier for ISW, SDSS QSO photo-$z$ and NVSS samples. Second, we use for the first time the WISE~$\times$~SuperCOSMOS catalog, which allows us to perform an all-sky analysis of the ISW up to $z\sim0.4$. Third, thanks to the use of photo-$z$s, we separate each dataset into different redshift bins, deriving the cross-correlation in each bin. This last step leads to a significant improvement in sensitivity. We remove cross-correlation between catalogs using masks which mutually exclude common regions of the sky. We use two methods to quantify the significance of the ISW effect. In the first one, we fix the cosmological model, derive linear galaxy biases of the catalogs, and then evaluate the significance of the ISW using a single parameter. In the second approach we perform a global fit of the ISW and of the galaxy biases varying the cosmological model. We find significances of the ISW in the range 4.7-5.0 $\sigma$ thus reaching, for the first time in such an analysis, the threshold of 5 $\sigma$. Without the redshift tomography we find a significance of $\sim$ 4.0 $\sigma$, which shows the importance of the binning method. Finally we use the ISW data to infer constraints on the Dark Energy redshift evolution and equation of state. We find that the redshift range covered by the catalogs is still not optimal to derive strong constraints, although this goal will be likely reached using future datasets such as from Euclid, LSST, and SKA.

We aim to investigate the galaxy environment in GAMA Galaxy Groups Catalogue (G3C) using a volume-limited galaxy sample from the Kilo Degree Survey Data Release 3. The k-Nearest Neighbour technique is adapted to take into account the probability density functions (PDFs) of photometric redshifts in our calculations. This algorithm was tested on simulated KiDS tiles, showing its capability of recovering the relation between galaxy colour, luminosity and local environment. The characterization of the galaxy environment in G3C groups shows systematically steeper density contrasts for more massive groups. The red galaxy fraction gradients in these groups is evident for most of group mass bins. The density contrast of red galaxies is systematically higher at group centers when compared to blue galaxy ones. In addition, distinct group center definitions are used to show that our results are insensitive to center definitions. These results confirm the galaxy evolution scenario which environmental mechanisms are responsible for a slow quenching process as galaxies fall into groups and clusters, resulting in a smooth observed colour gradients in galaxy systems.

We probe the isotropy of the Universe with the largest all-sky photometric redshift dataset currently available, namely WISE~$\times$~SuperCOSMOS. We search for dipole anisotropy of galaxy number counts in multiple redshift shells within the $0.10 < z < 0.35$ range, for two subsamples drawn from the same parent catalogue. Our results show that the dipole directions are in good agreement with most of the previous analyses in the literature, and in most redshift bins the dipole amplitudes are well consistent with $\Lambda$CDM-based mocks in the cleanest sample of this catalogue. In the $z<0.15$ range, however, we obtain a persistently large anisotropy in both subsamples of our dataset. Overall, we report no significant evidence against the isotropy assumption in this catalogue except for the lowest redshift ranges. The origin of the latter discrepancy is unclear, and improved data may be needed to explain it.

Taipan is a multi-object spectroscopic galaxy survey starting in 2017 that will cover 2pi steradians over the southern sky, and obtain optical spectra for about two million galaxies out to z<0.4. Taipan will use the newly-refurbished 1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN instrument, which includes an innovative 'Starbugs' positioning system capable of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up to 300 with a proposed upgrade) over the 6-deg diameter focal plane, and a purpose-built spectrograph operating from 370 to 870nm with resolving power R>2000. The main scientific goals of Taipan are: (i) to measure the distance scale of the Universe (primarily governed by the local expansion rate, H_0) to 1% precision, and the structure growth rate of structure to 5%; (ii) to make the most extensive map yet constructed of the mass distribution and motions in the local Universe, using peculiar velocities based on improved Fundamental Plane distances, which will enable sensitive tests of gravitational physics; and (iii) to deliver a legacy sample of low-redshift galaxies as a unique laboratory for studying galaxy evolution as a function of mass and environment. The final survey, which will be completed within 5 years, will consist of a complete magnitude-limited sample (i<17) of about 1.2x10^6 galaxies, supplemented by an extension to higher redshifts and fainter magnitudes (i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies. Observations and data processing will be carried out remotely and in a fully-automated way, using a purpose-built automated 'virtual observer' software and an automated data reduction pipeline. The Taipan survey is deliberately designed to maximise its legacy value, by complementing and enhancing current and planned surveys of the southern sky at wavelengths from the optical to the radio.

Using N-body simulations we study the impact of various systematic effects on the bulk flow (BF) and the Cosmic Mach Number (CMN). We consider two types of systematics: those related to survey properties and those induced by observer's location in the Universe. In the former category we model sparse sampling, velocity errors, and survey incompleteness. In the latter, we consider Local Group (LG) analogue observers, placed in a specific location within the Cosmic Web, satisfying various observational criteria. We differentiate such LG observers from Copernican ones, who are at random locations. We report strong systematic effects on the measured BF and CMN induced by sparse sampling, velocity errors and radial incompleteness. For BF most of these effects exceed 10\% for scales $R\leq100h^{-1}$Mpc. For CMN some of these systematics can be catastrophically large ($>50\%$) also on bigger scales. Moreover, we find that the position of the observer in the Cosmic Web significantly affects the locally measured BF (CMN), with effects as large as $\sim20\%$ ($30\%)$ at $R\leq50h^{-1}$Mpc for a LG-like observer as compared to a random one. This effect is comparable to the sample variance. To highlight the importance of these systematics, we additionally study a model of modified gravity (MG) with $\sim15\%$ enhanced growth rate. We found that the systematic effects can mimic the modified gravity signal. The worst-case scenario is realized for a case of a LG-like observer, when the effects induced by local structures are degenerate with the enhanced growth rate fostered by MG. Our results indicate that dedicated constrained simulations and realistic mock galaxy catalogs will be absolutely necessary to fully benefit from the statistical power of the forthcoming peculiar velocity data from surveys such as TAIPAN, WALLABY, Cosmic Flows-4 and SKA.

Using the almost all-sky 2MASS Photometric Redshift catalogue (2MPZ) we perform for the first time a tomographic analysis of galaxy angular clustering in the local Universe ($z<0.24$). We estimate the angular auto- and cross-power spectra of 2MPZ galaxies in three photometric redshift bins, and use dedicated mock catalogues to assess their errors. We measure a subset of cosmological parameters, having fixed the others at their Planck values, namely the baryon fraction $f_{b}=0.14^{+0.09}_{-0.06}$, the total matter density parameter $\Omega_{\rm m}=0.30\pm 0.06$, and the effective linear bias of 2MPZ galaxies $b_{\rm eff}$, which grows from $1.1^{+0.3}_{-0.4}$ at $\langle z \rangle=0.05$ up to $2.1^{+0.3}_{-0.5}$ at $\langle z\rangle=0.2$, largely because of the flux-limited nature of the dataset. The results obtained here for the local Universe agree with those derived with the same methodology at higher redshifts, and confirm the importance of the tomographic technique for next-generation photometric surveys such as Euclid or LSST.

Automatic source detection and classification tools based on machine learning (ML) algorithms are growing in popularity due to their efficiency when dealing with large amounts of data simultaneously and their ability to work in multidimensional parameter spaces. In this work, we present a new, automated method of outlier selection based on support vector machine (SVM) algorithm called one-class SVM (OCSVM), which uses the training data as one class to construct a model of 'normality' in order to recognize novel points. We test the performance of OCSVM algorithm on \textit{Wide-field Infrared Survey Explorer (WISE)} data trained on the Sloan Digital Sky Survey (SDSS) sources. Among others, we find $\sim 40,000$ sources with abnormal patterns which can be associated with obscured and unobscured active galactic nuclei (AGN) source candidates. We present the preliminary estimation of the clustering properties of these objects and find that the unobscured AGN candidates are preferentially found in less massive dark matter haloes ($M_{DMH}\sim10^{12.4}$) than the obscured candidates ($M_{DMH}\sim 10^{13.2}$). This result contradicts the unification theory of AGN sources and indicates that the obscured and unobscured phases of AGN activity take place in different evolutionary paths defined by different environments.

We investigate the nature of the extragalactic unresolved gamma-ray background (UGRB) by cross-correlating several galaxy catalogs with sky-maps of the UGRB built from 78 months of Pass 8 Fermi-Large Area Telescope data. This study updates and improves similar previous analyses in several aspects. Firstly, the use of a larger gamma-ray dataset allows us to investigate the energy dependence of the cross-correlation in more detail, using up to 8 energy bins over a wide energy range of [0.25-500] GeV. Secondly, we consider larger and deeper catalogs (2MASS Photometric-Redshift catalog, 2MPZ; WISE x SuperCOSMOS, WISC; and SDSS-DR12 photometric-redshift dataset) in addition to the ones employed in the previous studies (NVSS and SDSS-QSOs). Thirdly, we exploit the redshift information available for the above catalogs to divide them into redshift bins and perform the cross-correlation separately in each of them. Our results confirm, with higher statistical significance, the detection of cross-correlation signal between the UGRB maps and all the catalogs considered, on angular scales smaller than 1 degree. Significances range from 16.3 sigma for NVSS, 7 sigma for SDSS-DR12 and WISC, 5 sigma for 2MPZ and 4 sigma for SDSS-QSOs. Furthermore, including redshift tomography, the significance of the SDSS-DR12 signal strikingly rises up to 12 sigma and the one of WISC to 10.6 sigma. We offer a simple interpretation of the signal in the framework of the halo model. The precise redshift and energy information allows us to clearly detect a change over redshift in the spectral and clustering behavior of the gamma-ray sources contributing to the UGRB.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

Wide-angle photometric surveys of previously uncharted sky areas or wavelength regimes will always bring in unexpected sources whose existence and properties cannot be easily predicted from earlier observations: novelties or even anomalies. Such objects can be efficiently sought for with novelty detection algorithms. Here we present an application of such a method, called one-class support vector machines (OCSVM), to search for anomalous patterns among sources preselected from the mid-infrared AllWISE catalogue covering the whole sky. To create a model of expected data we train the algorithm on a set of objects with spectroscopic identifications from the SDSS DR13 database, present also in AllWISE. OCSVM detects as anomalous those sources whose patterns - WISE photometric measurements in this case - are inconsistent with the model. Among the detected anomalies we find artefacts, such as objects with spurious photometry due to blending, but most importantly also real sources of genuine astrophysical interest. Among the latter, OCSVM has identified a sample of heavily reddened AGN/quasar candidates distributed uniformly over the sky and in a large part absent from other WISE-based AGN catalogues. It also allowed us to find a specific group of sources of mixed types, mostly stars and compact galaxies. By combining the semi-supervised OCSVM algorithm with standard classification methods it will be possible to improve the latter by accounting for sources which are not present in the training sample but are otherwise well-represented in the target set. Anomaly detection adds flexibility to automated source separation procedures and helps verify the reliability and representativeness of the training samples. It should be thus considered as an essential step in supervised classification schemes to ensure completeness and purity of produced catalogues.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars. Here we present the third public data release (DR3) and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lensfit shear measurement packages are used. In DR3 stacked ugri images, weight maps, masks, and source lists for 292 new survey tiles (~300 sq.deg) are made available. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (447 sq.deg). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5 sigma in a 2 arcsec aperture) in ugri, respectively, and the typical r-band PSF size is less than 0.7 arcsec. The photometric homogenization scheme ensures accurate colors and an absolute calibration stable to ~2% for gri and ~3% in u. Separately released are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

We present the tomographic cross-correlation between galaxy lensing measured in the Kilo Degree Survey (KiDS-450) with overlapping lensing measurements of the cosmic microwave background (CMB), as detected by Planck 2015. We compare our joint probe measurement to the theoretical expectation for a flat $\Lambda$CDM cosmology, assuming the best-fitting cosmological parameters from the KiDS-450 cosmic shear and Planck CMB analyses. We find that our results are consistent within $1\sigma$ with the KiDS-450 cosmology, with an amplitude re-scaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck cosmology, we find our results are consistent within $2\sigma$, with $A_{\it Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases when the contamination to the signal by intrinsic galaxy alignments is accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic analysis of the galaxy lensing -- CMB lensing cross-correlation signal, and is based on five photometric redshift bins. We use this measurement as an independent validation of the multiplicative shear calibration and of the calibrated source redshift distribution at high redshifts. We find that constraints on these two quantities are strongly correlated when obtained from this technique, which should therefore not be considered as a stand-alone competitive calibration tool.

We measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.

We present The-wiZZ, an open source and user-friendly software for estimating the redshift distributions of photometric galaxies with unknown redshifts by spatially cross-correlating them against a reference sample with known redshifts. The main benefit of The-wiZZ is in separating the angular pair finding and correlation estimation from the computation of the output clustering redshifts allowing anyone to create a clustering redshift for their sample without the intervention of an "expert". It allows the end user of a given survey to select any sub-sample of photometric galaxies with unknown redshifts, match this sample's catalog indices into a value-added data file, and produce a clustering redshift estimation for this sample in a fraction of the time it would take to run all the angular correlations needed to produce a clustering redshift. We show results with this software using photometric data from the Kilo-Degree Survey (KiDS) and spectroscopic redshifts from the Galaxy and Mass Assembly (GAMA) survey and the Sloan Digital Sky Survey (SDSS). The results we present for KiDS are consistent with the redshift distributions used in a recent cosmic shear analysis from the survey. We also present results using a hybrid machine learning-clustering redshift analysis that enables the estimation of clustering redshifts for individual galaxies. The-wiZZ can be downloaded at http://github.com/morriscb/The-wiZZ/.

The Mg II emission line is visible in the optical band for intermediate redshift quasars (0.4 < z < 1.6) and it is thus an extremely important tool to measure the black hole mass and to understand the structure of the Broad Line Region. We aim to determine the substructure and the variability of the Mg II line with the aim to identify which part of the line comes from a medium in Keplerian motion. Using the Southern African Large Telescope (SALT) with the Robert Stobie Spectrograph (RSS) we performed ten spectroscopic observations of quasar HE 0435-4312 (z = 1.2231) over a period of three years (Dec 23/24, 2012 to Dec 7/8, 2015). Both the Mg II line and the Fe II pseudo-continuum increase with time. We clearly detect the systematic shift of the Mg II line with respect to the Fe II over the years, corresponding to the acceleration of 104 pm 14 km/s/year in the quasar rest frame. The Mg II line shape is clearly non-Gaussian but single-component, and the increase in line equivalent width and line shift is not accompanied with significant evolution of the line shape. We analyse the conditions in the Mg II and Fe II formation region and we note that the very large difference in the covering factor and the turbulent velocity also support the conclusion that the two regions are spatially separated. The measured acceleration of the line systematic shift is too large to connect it with the orbital motion at a distance of the Broad Line Region (BLR) in this source. It may imply a precessing inner disk illuminating the BLR. Further monitoring is still needed to better constrain the variability mechanism.

In this paper, we present the tools used to search for galaxy clusters in the Kilo Degree Survey (KiDS), and our first results. The cluster detection is based on an implementation of the optimal filtering technique that enables us to identify clusters as over-densities in the distribution of galaxies using their positions on the sky, magnitudes, and photometric redshifts. The contamination and completeness of the cluster catalog are derived using mock catalogs based on the data themselves. The optimal signal to noise threshold for the cluster detection is obtained by randomizing the galaxy positions and selecting the value that produces a contamination of less than 20%. Starting from a subset of clusters detected with high significance at low redshifts, we shift them to higher redshifts to estimate the completeness as a function of redshift: the average completeness is ~ 85%. An estimate of the mass of the clusters is derived using the richness as a proxy. We obtained 1858 candidate clusters with redshift 0 < z_c < 0.7 and mass 13.5 < log(M500/Msun) < 15 in an area of 114 sq. degrees (KiDS ESO-DR2). A comparison with publicly available Sloan Digital Sky Survey (SDSS)-based cluster catalogs shows that we match more than 50% of the clusters (77% in the case of the redMaPPer catalog). We also cross-matched our cluster catalog with the Abell clusters, and clusters found by XMM and in the Planck-SZ survey; however, only a small number of them lie inside the KiDS area currently available.

We present an analysis of the mid-infrared WISE sources seen within the equatorial GAMA G12 field, located in the North Galactic Cap. Our motivation is to study and characterize the behavior of WISE source populations in anticipation of the deep multi-wavelength surveys that will define the next decade, with the principal science goal of mapping the 3D large scale structures and determining the global physical attributes of the host galaxies. In combination with cosmological redshifts, we identify galaxies from their WISE W1 3.4um resolved emission, and by performing a star-galaxy separation using apparent magnitude, colors and statistical modeling of star-counts. The resultant galaxy catalog has ~590,000 sources in 60 deg^2, reaching a W1 5-sigma depth of 31 uJy. At the faint end, where redshifts are not available, we employ a luminosity function analysis to show that approximately 27% of all WISE extragalactic sources to a limit of 17.5 mag (31 uJy) are at high redshift, z > 1. The spatial distribution is investigated using two-point correlation functions and a 3D source density characterization at 5 Mpc and 20 Mpc scales. For angular distributions, we find brighter and more massive sources are strongly clustered relative to fainter and lower mass source; likewise, based on WISE colors, spheroidal galaxies have the strongest clustering, while late-type disk galaxies have the lowest clustering amplitudes. Along the radial direction, the strongest clustering is in the largest redshift shell, while the weakest is in the nearest redshift shell, consistent with the stellar mass and morphological type dependency results. In three dimensions, we find a number of distinct groupings, often bridged by filaments and super-structures. Using special visualization tools, we map these structures, exploring how clustering may play a role with stellar mass and galaxy type.

Verlinde (2016) proposed that the observed excess gravity in galaxies and clusters is the consequence of Emergent Gravity (EG). In this theory the standard gravitational laws are modified on galactic and larger scales due to the displacement of dark energy by baryonic matter. EG gives an estimate of the excess gravity (described as an apparent dark matter density) in terms of the baryonic mass distribution and the Hubble parameter. In this work we present the first test of EG using weak gravitational lensing, within the regime of validity of the current model. Although there is no direct description of lensing and cosmology in EG yet, we can make a reasonable estimate of the expected lensing signal of low redshift galaxies by assuming a background LambdaCDM cosmology. We measure the (apparent) average surface mass density profiles of 33,613 isolated central galaxies, and compare them to those predicted by EG based on the galaxies' baryonic masses. To this end we employ the ~180 square degrees overlap of the Kilo-Degree Survey (KiDS) with the spectroscopic Galaxy And Mass Assembly (GAMA) survey. We find that the prediction from EG, despite requiring no free parameters, is in good agreement with the observed galaxy-galaxy lensing profiles in four different stellar mass bins. Although this performance is remarkable, this study is only a first step. Further advancements on both the theoretical framework and observational tests of EG are needed before it can be considered a fully developed and solidly tested theory.

We report the discovery of a potentially major supercluster that extends across the Galactic Plane in the constellation of Vela, at a mean recessional velocity of ~18,000 km/s. Recent multi-object spectroscopic observations of this Vela Supercluster (VSCL), using AAOmega+2dF and the Southern African Large Telescope, confirm an extended galaxy overdensity in the Zone of Avoidance (ZOA) located where residual bulk flows predict a considerable mass excess. We present a preliminary analysis of ~4,500 new spectroscopic galaxy redshifts obtained in the ZOA centred on the Vela region (l=272.5+-20 deg, b=0+-10 deg). The presently sparsely-sampled dataset traces an overdensity that covers 25 deg in Galactic longitude on either side of the Plane, suggesting an extent of 25 deg x 20 deg, corresponding to ~115 x 90 $h_{70}$ Mpc at the supercluster redshift. In redshift space, the overdensity appears to consist of two merging wall-like structures, interspersed with clusters and groups. Both the velocity histogram and the morphology of the multi-branching wall structure are consistent with a supercluster classification. $K_s^o$ galaxy counts show an enhancement of ~1.2 over the survey area for galaxies brighter than $M_K^*$ at the VSCL distance, and a galaxy overdensity of $\delta=0.50\rm{-}0.77$ within a photometric redshift shell around the VSCL, when compared to various Two-Micron All-Sky Survey samples. Taking account of selection effects, the VSCL is estimated to contribute $v_\rm{LG} \gtrsim 50$ km/s to the motion of the Local Group.

We present a new training set for estimating empirical photometric redshifts of galaxies, which was created as part of the 2dFLenS project. This training set is located in a 700 sq deg area of the KiDS South field and is randomly selected and nearly complete at r<19.5. We investigate the photometric redshift performance obtained with ugriz photometry from VST-ATLAS and W1/W2 from WISE, based on several empirical and template methods. The best redshift errors are obtained with kernel-density estimation, as are the lowest biases, which are consistent with zero within statistical noise. The 68th percentiles of the redshift scatter for magnitude-limited samples at r<(15.5, 17.5, 19.5) are (0.014, 0.017, 0.028). In this magnitude range, there are no known ambiguities in the colour-redshift map, consistent with a small rate of redshift outliers. In the fainter regime, the KDE method produces p(z) estimates per galaxy that represent unbiased and accurate redshift frequency expectations. The p(z) sum over any subsample is consistent with the true redshift frequency plus Poisson noise. Further improvements in redshift precision at r<20 would mostly be expected from filter sets with narrower passbands to increase the sensitivity of colours to small changes in redshift.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

The two currently largest all-sky photometric datasets, WISE and SuperCOSMOS, were cross-matched by Bilicki et al. (2016) (B16) to construct a novel photometric redshift catalogue on 70% of the sky. Galaxies were therein separated from stars and quasars through colour cuts, which may leave imperfections because of mixing different source types which overlap in colour space. The aim of the present work is to identify galaxies in the WISExSuperCOSMOS catalogue through an alternative approach of machine learning. This allows us to define more complex separations in the multi-colour space than possible with simple colour cuts, and should provide more reliable source classification. For the automatised classification we use the support vector machines learning algorithm, employing SDSS spectroscopic sources cross-matched with WISExSuperCOSMOS as the training and verification set. We perform a number of tests to examine the behaviour of the classifier (completeness, purity and accuracy) as a function of source apparent magnitude and Galactic latitude. We then apply the classifier to the full-sky data and analyse the resulting catalogue of candidate galaxies. We also compare thus produced dataset with the one presented in B16. The tests indicate very high accuracy, completeness and purity (>95%) of the classifier at the bright end, deteriorating for the faintest sources, but still retaining acceptable levels of 85%. No significant variation of classification quality with Galactic latitude is observed. Application of the classifier to all-sky WISExSuperCOSMOS data gives 15 million galaxies after masking problematic areas. The resulting sample is purer than the one in B16, at a price of lower completeness over the sky. The automatic classification gives a successful alternative approach to defining a reliable galaxy sample as compared to colour cuts.

We assess the effect of the local large scale structure on the estimation of two-point statistics of the observed radial peculiar velocities of galaxies. A large N-body simulation is used to examine these statistics from the perspective of random observers as well as "Local Group (LG)-like" observers conditioned to reside in an environment resembling the observed universe within 20 Mpc. The local environment systematically distorts the shape and amplitude of velocity statistics with respect to ensemble-averaged measurements made by a Copernican (random) observer. The Virgo cluster has the most significant impact, introducing large systematic deviations in all the statistics. For a simple "top-hat" selection function, an idealized survey extending to $\sim 160h^{-1}\,{\rm Mpc}$ or deeper is needed to completely mitigate the effects of the local environment. Using shallower catalogues leads to systematic deviations of the order of $50$ to $200\%$ depending on the scale considered. For a flat redshift distribution similar to the one of the CosmicFlows-3 survey, the deviations are even more prominent in both the shape and amplitude at all separations considered $({\stackrel{<}{{}_\sim}} 100h^{-1}\,{\rm Mpc})$. Conclusions based on statistics calculated without taking into account the impact of the local environment should be revisited.

We cross-match the two currently largest all-sky photometric catalogs, mid-infrared WISE and SuperCOSMOS scans of UKST/POSS-II photographic plates, to obtain a new galaxy sample that covers 3pi steradians. In order to characterize and purify the extragalactic dataset, we use external GAMA and SDSS spectroscopic information to define quasar and star loci in multicolor space, aiding the removal of contamination from our extended-source catalog. After appropriate data cleaning we obtain a deep wide-angle galaxy sample that is approximately 95% pure and 90% complete at high Galactic latitudes. The catalog contains close to 20 million galaxies over almost 70% of the sky, outside the Zone of Avoidance and other confused regions, with a mean surface density of over 650 sources per square degree. Using multiwavelength information from two optical and two mid-IR photometric bands, we derive photometric redshifts for all the galaxies in the catalog, using the ANNz framework trained on the final GAMA-II spectroscopic data. Our sample has a median redshift of z_{med} = 0.2 but with a broad dN/dz reaching up to z>0.4. The photometric redshifts have a mean bias of |delta_z|~10^{-3}, normalized scatter of sigma_z = 0.033 and less than 3% outliers beyond 3sigma_z. Comparison with external datasets shows no significant variation of photo-z quality with sky position. Together with the overall statistics, we also provide a more detailed analysis of photometric redshift accuracy as a function of magnitudes and colors. The final catalog is appropriate for `all-sky' 3D cosmology to unprecedented depths, in particular through cross-correlations with other large-area surveys. It should also be useful for source pre-selection and identification in forthcoming surveys such as TAIPAN or WALLABY.

The WISE satellite has detected hundreds of millions sources over the entire sky. Classifying them reliably is however a challenging task due to degeneracies in WISE multicolour space and low levels of detection in its two longest-wavelength bandpasses. Here we aim at obtaining comprehensive and reliable star, galaxy and quasar catalogues based on automatic source classification in full-sky WISE data. This means that the final classification will employ only parameters available from WISE itself, in particular those reliably measured for a majority of sources. For the automatic classification we applied the support vector machines (SVM) algorithm, which requires a training sample with relevant classes already identified, and we chose to use the SDSS spectroscopic dataset for that purpose. By calibrating the classifier on the test data drawn from SDSS, we first established that a polynomial kernel is preferred over a radial one for this particular dataset. Next, using three classification parameters (W1 magnitude, W1-W2 colour, and a differential aperture magnitude) we obtained very good classification efficiency in all the tests. At the bright end, the completeness for stars and galaxies reaches ~95%, deteriorating to ~80% at W1=16 mag, while for quasars it stays at a level of ~95% independently of magnitude. Similar numbers are obtained for purity. Application of the classifier to full-sky WISE data, flux-limited to 16 mag (Vega) in the 3.4 {\mu}m channel, and appropriate a posteriori cleaning allowed us to obtain reliably-looking catalogues of star and galaxy candidates. However, the sources flagged by the classifier as `quasars' are in fact dominated by dusty galaxies but also exhibit contamination from sources located mainly at low ecliptic latitudes, consistent with Solar System objects. [abridged]

We describe the construction of an all-sky galaxy catalogue, using SuperCOSMOS scans of Schmidt photographic plates from the UKST and POSS2 surveys. The photographic photometry is calibrated using SDSS data, with results that are linear to 2% or better. All-sky photometric uniformity is achieved by matching plate overlaps and also by requiring homogeneity in optical-to-2MASS colours, yielding zero points that are uniform to 0.03 mag. or better. The typical AB depths achieved are B_J<21, R_F<19.5 and I_N<18.5, with little difference between hemispheres. In practice, the I_N plates are shallower than the B_J & R_F plates, so for most purposes we advocate the use of a catalogue selected in these two latter bands. At high Galactic latitudes, this catalogue is approximately 90% complete with 5% stellar contamination; we quantify how the quality degrades towards the Galactic plane. At low latitudes, there are many spurious galaxy candidates resulting from stellar blends: these approximately match the surface density of true galaxies at |b|=30 deg. Above this latitude, the catalogue limited in B_J & R_F contains in total about 20 million galaxy candidates, of which 75% are real. This contamination can be removed, and the sky coverage extended, by matching with additional datasets. This SuperCOSMOS catalogue has been matched with 2MASS and with WISE, yielding quasi-allsky samples of respectively 1.5 million and 18.5 million galaxies, to median redshifts of 0.08 and 0.20. This legacy dataset thus continues to offer a valuable resource for large-angle cosmological investigations.

The nature of dark energy, driving the accelerated expansion of the Universe, is one of the most important issues in modern astrophysics. In order to understand this phenomenon, we need precise astrophysical probes of the universal expansion spanning wide redshift ranges. Quasars have recently emerged as such a probe, thanks to their high intrinsic luminosities and, most importantly, our ability to measure their luminosity distances independently of redshifts. Here we report our ongoing work on observational reverberation mapping using the time delay of the Mg II line, performed with the South African Large Telescope (SALT).

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

Various aspects of cosmology require comprehensive all-sky mapping of the cosmic web to considerable depths. In order to probe the whole extragalactic sky beyond 100 Mpc, one must draw on multiwavelength datasets and state-of-the-art photometric redshift techniques. Here I summarize our dedicated program that employs the largest photometric all-sky surveys -- 2MASS, WISE and SuperCOSMOS -- to obtain accurate redshift estimates of millions of galaxies. The first outcome of these efforts -- the 2MASS Photometric Redshift catalog (2MPZ) -- was publicly released in 2013 and includes almost 1 million galaxies with a median redshift of z~0.1. I discuss how this catalog was constructed and how it is being used for various cosmological tests. I also present how combining the WISE mid-infrared survey with SuperCOSMOS optical data allowed us to push to depths over 1 Gpc on unprecedented angular scales. These photometric redshift samples, with about 20 million sources in total, provide access to volumes large enough to study observationally the Copernican Principle of universal homogeneity and isotropy, as well as to probe various aspects of dark energy and dark matter through cross-correlations with other data such as the cosmic microwave or gamma-ray backgrounds. Last but not least, they constitute a test-bed for forthcoming wide-angle multi-million galaxy samples expected from such instruments as the SKA, Euclid, or LSST.

We have undertaken a dedicated program of automatic source classification in the WISE database merged with SuperCOSMOS scans, comprehensively identifying galaxies, quasars and stars on most of the unconfused sky. We use the Support Vector Machines classifier for that purpose, trained on SDSS spectroscopic data. The classification has been applied to a photometric dataset based on all-sky WISE 3.4 and 4.6 $\mu$m information cross-matched with SuperCOSMOS B and R bands, which provides a reliable sample of $\sim170$ million sources, including galaxies at $z_{\rm med}\sim0.2$, as well as quasars and stars. The results of our classification method show very high purity and completeness (more than 96\%) of the separated sources, and the resultant catalogs can be used for sophisticated analyses, such as generating all-sky photometric redshifts.

We present the first results of our dedicated programme of automatised classification of galaxies, stars and quasars in the mid-infrared all-sky data from the WISE survey. We employ the Support Vector Machines (SVM) algorithm, which defines a hyperplane separating different classes of sources in a multidimensional space of arbitrarily chosen parameters. This approach consists of four general steps: 1) selection of the training sample, 2) selection of the optimal parameter space, 3) training of the classifier, 4) application to target data. Here, as the training set, we use sources from a cross-correlation of the WISE catalogue with the SDSS spectroscopic sample. The performance of the SVM classifier was tested as a function of size of the training set, dimension of the parameter space, WISE apparent magnitude and Galactic extinction. We find that our classifier provides promising results already for three classification parameters: magnitude, colour and differential aperture magnitude. Completeness and purity levels as high as 95% are obtained for quasars, while for galaxies and stars they vary between 80-95% depending on the magnitude, deteriorating for fainter sources.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

Using the 2MASS Photometric Redshift catalogue we perform a number of statistical tests aimed at detecting possible departures from statistical homogeneity and isotropy in the large-scale structure of the Universe. Making use of the angular homogeneity index, an observable proposed in a previous publication, as well as studying the scaling of the angular clustering and number counts with magnitude limit, we place constraints on the fractal nature of the galaxy distribution. We find that the statistical properties of our sample are in excellent agreement with the standard cosmological model, and that it reaches the homogeneous regime significantly faster than a class of fractal models with dimensions $D<2.75$. As part of our search for systematic effects, we also study the presence of hemispherical asymmetries in our data, finding no significant deviation beyond those allowed by the concordance model.

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

Quasars, as the most luminous persistent sources in the Universe, have broad applications for cosmological studies. In particular, they can be employed to directly measure the expansion history of the Universe, similarly to SNe Ia. The advantage of quasars is that they are numerous, cover a broad range of redshifts, up to $z = 7$, and do not show significant evolution of metallicity with redshift. The idea is based on the relation between the time delay of an emission line and the continuum, and the absolute monochromatic luminosity of a quasar. For intermediate redshift quasars, the suitable line is Mg II. Between December 2012 and March 2014, we performed five spectroscopic observations of the QSO CTS C30.10 ($z = 0.900$) using the South African Large Telesope (SALT), supplemented with photometric monitoring, with the aim of determining the variability of the line shape, changes in the total line intensity and in the continuum. We show that the method is very promising.

Quasars can be used as a complementary tool to SN Ia to probe the distribution of dark energy in the Universe by measuring the time delay of the emission line with respect to the continuum. The understanding of the Mg II emission line structure is important for cosmological application and for the black hole mass measurements of intermediate redshift quasars. The knowledge of the shape of Mg II line and its variability allows to determine which part of the line should be used to measure the time delay and the black hole mass. We thus aim at determination of the structure and the variability of the Mg II line as well as of the underlying Fe II pseudo-continuum. We performed five spectroscopic observations of a quasar CTS C30.10 (z = 0.9000) with the SALT telescope between December 2012 and March 2014, and we studied the variations of the spectral shape in the 2700 - 2900 A rest frame. We showed that the Mg II line in this source consists of two kinematic components, which makes the source representative of type B quasars. Both components were well modeled with a Lorentzian shape. Both components vary in a similar way. The Fe II contribution seems to be related only to the first (blue) Mg II component}. Broad band spectral fitting rather favor the use of the whole line profile. The contribution of the Narrow Line Region to Mg II is very low, below 2%. The Mg II variability is lower than the variability of the continuum, which is consistent with the simple reprocessing scenario. The variability level of CTS C30.10 and the measurement accuracy of the line and continuum is high enough to expect that further monitoring will allow to measure the time delay between the Mg II line and continuum.

Our view of the low-redshift Cosmic Web has been revolutionized by galaxy redshift surveys such as 6dFGS, SDSS and 2MRS. However, the trade-off between depth and angular coverage limits a systematic three-dimensional account of the entire sky beyond the Local Volume (z<0.05). In order to reliably map the Universe to cosmologically significant depths over the full celestial sphere, one must draw on multiwavelength datasets and state-of-the-art photometric redshift techniques. We have undertaken a dedicated program of cross-matching the largest photometric all-sky surveys -- 2MASS, WISE and SuperCOSMOS -- to obtain accurate redshift estimates of millions of galaxies. The first outcome of these efforts -- the 2MASS Photometric Redshift catalog (2MPZ, Bilicki et al. 2014a) -- has been publicly released and includes almost 1 million galaxies with a mean redshift of z=0.08. Here we summarize how this catalog was constructed and how using the WISE mid-infrared sample together with SuperCOSMOS optical data allows us to push to redshift shells of z~0.2--0.3 on unprecedented angular scales. Our catalogs, with ~20 million sources in total, provide access to cosmological volumes crucial for studies of local galaxy flows (clustering dipole, bulk flow) and cross-correlations with the cosmic microwave background such as the integrated Sachs-Wolfe effect or lensing studies.

While Norma (ACO3627) is the richest cluster in the Great Attractor (GA) region, its role in the local dynamics is poorly understood. The Norma cluster has a mean redshift (z_CMB) of 0.0165 and has been proposed as the "core" of the GA. We have used the Ks-band Fundamental Plane (FP) to measure Norma cluster's distance with respect to the Coma cluster. We report FP photometry parameters (effective radii and surface brightnesses), derived from ESO NTT SOFI images, and velocity dispersions, from AAT 2dF spectroscopy, for 31 early-type galaxies in the cluster. For the Coma cluster we use 2MASS images and SDSS velocity dispersion measurements for 121 early-type galaxies to generate the calibrating FP dataset. For the combined Norma-Coma sample we measure FP coefficients of a=1.465+/-0.059 and b=0.326+/-0.020. We find an rms scatter, in log(central velocity dispersion) of 0.08 dex which corresponds to a distance uncertainty of 28% per galaxy. The zero point offset between Norma's and Coma's FPs is 0.154+/-0.014 dex. Assuming that the Coma cluster is at rest with respect to the cosmic microwave background frame and z_CMB(Coma)=0.0240, we derive a distance to the Norma cluster of 5026+/-160 km/s, and the derived peculiar velocity is -72+/-170 km/s, i.e., consistent with zero. This is lower than previously reported positive peculiar velocities for clusters/groups/galaxies in the GA region and hence the Norma cluster may indeed represent the GA's "core".

Key cosmological applications require the three-dimensional galaxy distribution on the entire celestial sphere. These include measuring the gravitational pull on the Local Group, estimating the large-scale bulk flow and testing the Copernican principle. However, the largest all-sky redshift surveys -- the 2MRS and IRAS PSCz -- have median redshifts of only z=0.03 and sample the very local Universe. There exist all-sky galaxy catalogs reaching much deeper -- SuperCOSMOS in the optical, 2MASS in the near-IR and WISE in the mid-IR -- but these lack complete redshift information. At present, the only rapid way towards larger 3D catalogs covering the whole sky is through photometric redshift techniques. In this paper we present the 2MASS Photometric Redshift catalog (2MPZ) containing 1 million galaxies, constructed by cross-matching 2MASS XSC, WISE and SuperCOSMOS all-sky samples and employing the artificial neural network approach (the ANNz algorithm), trained on such redshift surveys as SDSS, 6dFGS and 2dFGRS. The derived photometric redshifts have errors nearly independent of distance, with an all-sky accuracy of sigma_z=0.015 and a very small percentage of outliers. In this way, we obtain redshift estimates with a typical precision of 12% for all the 2MASS XSC galaxies that lack spectroscopy. In addition, we have made an early effort towards probing the entire 3D sky beyond 2MASS, by pairing up WISE with SuperCOSMOS and training the ANNz on GAMA redshift data reaching currently to z_med~0.2. This has yielded photo-z accuracies comparable to those in the 2MPZ. These all-sky photo-z catalogs, with a median z~0.1 for the 2MPZ, and significantly deeper for future WISE-based samples, will be the largest and most complete of their kind for the foreseeable future.

Comparison of peculiar velocities of galaxies with their gravitational accelerations (induced by the density field) is one of the methods to constrain the redshift distortion parameter \beta=(\Omega_m^0.55)/b, where \Omega_m is the non-relativistic matter density parameter and b is the linear bias. In particular, one can use the motion of the Local Group (LG) for that purpose. Its peculiar velocity is known from the dipole component of the cosmic microwave background, whereas its acceleration can be estimated with the use of an all-sky galaxy catalog, from the so-called clustering dipole. At the moment, the biggest dataset of that kind is the Two Micron All Sky Survey Extended Source Catalog (2MASS XSC) containing almost 1 million galaxies and complete up to ~300 Mpc/h. We applied 2MASS data to measure LG acceleration and used two methods to estimate the beta parameter. Both of them yield \beta~0.4 with an error of several per cent, which is the most precise determination of this parameter from the clustering dipole to date.

High redshift quasars can be used to deduce the distribution of dark energy in the Universe, as a complementary tool to SN Ia. The method is based on determination of the size of the Broad Line Region from the emission line delay, determination of the absolute monochromatic luminosity either from the observed statistical relation or from a model of the formation of the Broad Line Region, and determination of the observed monochromatic flux from photometry. This allows to obtain the luminosity distance to a quasar independently from its redshift. The accuracy of the measurements is however, a key issue. We model the expected accuracy of the measurements by creating artificial quasar monochromatic lightcurves and responses from the Broad Line Region under various assumptions about the variability of a quasar, Broad Line Region extension, distribution of the measurements in time, accuracy of the measurements and the intrinsic line variability. We show that the five year monitoring based on Mg II line should give the accuracy of 0.06 - 0.32 magnitude in the distance modulus which allows to put interesting constraints on the cosmological models. Monitoring of higher redshift quasars based on CIV lines is problematic due to much higher level of the intrinsic variability of CIV in comparison with Mg II.

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

We analyse the possibility that our Universe could be described by the model recently proposed by Melia & Shevchuk (2012), where the Hubble scale R_h=c/H is at all times equal to the distance ct that light has travelled since the Big Bang. In such a model, the scale factor is proportional to cosmic time and there is neither acceleration nor deceleration of the expansion. We first point out problems with the very foundations of the model and its consequences for the evolution of the Universe. Next, we compare predictions of the model with observational data. As probes of the expansion we use distance data of supernovae type Ia, as well as Hubble rate data obtained from cosmic chronometers and radial baryon acoustic oscillations. We analyse the redshift evolution of the Hubble parameter and its redshift derivatives, together with the so-called O_m diagnostic and the deceleration parameter. To reliably estimate smooth functions and their derivatives from discrete data, we use the recently developed Gaussian Processes in Python package (GaPP). Our general conclusion is that the discussed model is strongly disfavoured by observations, especially at low redshifts (z<0.5). In particular, it predicts specific constant values for the deceleration parameter and for redshift derivatives of the Hubble parameter, which is ruled out by the data.

In this thesis, we use the motion of the Local Group of galaxies (LG) through the Universe to measure the cosmological parameter of non-relativistic matter density, Omega_m. For that purpose, we compare the peculiar velocity of the LG with its gravitational acceleration. The former is known from the dipole of the cosmic microwave background radiation and the latter is estimated here from the clustering dipole of galaxies in the Two Micron All Sky Survey (2MASS) Extended Source Catalog. We start by presenting the general framework of perturbation theory of gravitational instability in the expanding Universe and how it applies to the peculiar motion of the LG. Next, we study a particular effect for the dipole measurement, related to the fact that a nearby Local Void is partially hidden behind our Galaxy. We then describe in detail how we handled the 2MASS extragalactic data for the purpose of our analysis. Finally, we present two methods to estimate the density Omega_m, combined with the linear biasing b into the parameter beta = (Omega_m)^{0.55} / b, from the comparison of the LG velocity and acceleration. The first approach is to study the growth of the 2MASS clustering dipole with increased depth of the sample and compare it with theoretical expectations. The second is to apply the maximum-likelihood method in order to improve the precision of the measurement. With both these methods we find beta=0.4 and Omega_m=0.2, which is consistent with various independent estimates. We also briefly mention some future prospects in the field.

There is a long-standing controversy about the convergence of the dipole moment of the galaxy angular distribution (the so-called clustering dipole). We study the growth of the clustering dipole of galaxies as a function of the limiting flux of the sample from the Two Micron All Sky Survey (2MASS). Contrary to some earlier claims, we find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to 13.5 mag in the near-infrared K_s band (equivalent to an effective distance of 300 Mpc/h). We compare the observed growth of the dipole with the theoretically expected, conditional one (i.e., given the velocity of the Local Group relative to the CMB), for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of its theoretical counterpart once the proper observational window of the 2MASS flux-limited catalog is included. For a contrast, if the adopted window is a top-hat, then the predicted dipole grows significantly faster and converges to its final value for a distance of about 300 Mpc/h. By comparing the observational windows, we show that for a given flux limit and a corresponding distance limit, the 2MASS flux-weighted window passes less large-scale signal than the top-hat one. We conclude that the growth of the 2MASS dipole for effective distances greater than 200 Mpc/h is only apparent. On the other hand, for a distance of 80 Mpc/h (mean depth of the 2MASS Redshift Survey) and the LambdaCDM power spectrum, the true dipole is expected to reach only ~80% of its final value. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.04, which leads to an estimate of the density parameter Omega_m = 0.20+-0.08.

We study the growth of the clustering dipole of galaxies from the Two Micron All Sky Survey (2MASS). We find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to about 300 Mpc/h. We compare the observed growth of the dipole with the theoretically expected, conditional growth for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of the theoretical one, once the proper observational window of the 2MASS flux dipole is included. For a contrast, if the adopted window is a top hat, then the predicted dipole grows significantly faster and converges to its final value at a distance of about 200 Mpc/h. We study the difference between the top-hat window and the window for the flux-limited 2MASS survey and we conclude that the growth of the 2MASS dipole at effective distances greater than 200 Mpc/h is only apparent. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.05, which gives a rough estimate of Omega_m = 0.2+-0.1.

In measurements of the clustering dipole from all-sky surveys, an important problem is the lack of information about galaxy distribution in the so-called Zone of Avoidance (ZoA). The existence of the Local Void (LV) has a systematic effect on these measurements. If the ZoA is randomly filled with mock galaxies, then the calculated acceleration of the Local Group of galaxies (LG) has a spurious component, resulting from the lack of real galaxies in the intersection of the LV with the ZoA. This component affects both the misalignment angle between the clustering dipole and the CMB dipole, and the inferred value of mean matter density Omega_m. We calculate the amplitude of the spurious acceleration acting on the LG due to the LV. Its value depends on the geometry and size of the LV, as well as on its density contrast. However, under the simplest assumption of the LV being spherical and completely empty, within the linear theory this amplitude amounts only to about 45 km/s in units of velocity. The resulting change in the misalignment angle is smaller than 1 degree, and the fractional change in the deduced value of Omega_m is about 5%. Accounting for observationally indicated elongation of the LV and maintaining the maximising assumption of a complete lack of galaxies inside increases these numbers only moderately. Specifically, the amplitude of the spurious acceleration rises to about 60 km/s, the misalignment angle remains still smaller than 1 degree, and the fractional change in the deduced value of Omega_m is enhanced to about 7%. Thus, despite the overall importance of the Local Void for the motion of the Local Group, the influence of the intersection of the LV with the ZoA on measurements of the clustering dipole is found to be only a minor systematic effect.

We study the cosmic velocity-density relation using the spherical collapse model (SCM) as a proxy to non-linear dynamics. Although the dependence of this relation on cosmological parameters is known to be weak, we retain the density parameter Omega_m in SCM equations, in order to study the limit Omega_m -> 0. We show that in this regime the considered relation is strictly linear, for arbitrary values of the density contrast, on the contrary to some claims in the literature. On the other hand, we confirm that for realistic values of Omega_m the exact relation in the SCM is well approximated by the classic formula of Bernardeau (1992), both for voids (delta<0) and for overdensities up to delta ~ 3. Inspired by this fact, we find further analytic approximations to the relation for the whole range delta from -1 to infinity. Our formula for voids accounts for the weak Omega_m-dependence of their maximal rate of expansion, which for Omega_m < 1 is slightly smaller that 3/2. For positive density contrasts, we find a simple relation div v = 3 H_0 (Omega_m)^(0.6) [ (1+delta)^(1/6) - (1+delta)^(1/2) ], that works very well up to the turn-around (i.e. up to delta ~ 13.5 for Omega_m = 0.25 and neglected Omega_Lambda). Having the same second-order expansion as the formula of Bernardeau, it can be regarded as an extension of the latter for higher density contrasts. Moreover, it gives a better fit to results of cosmological numerical simulations.

A comparison of the 2MASS flux dipole to the CMB dipole can serve as a method to constrain a combination of the cosmological parameter Omega_m and the luminosity bias of the 2MASS survey. For this constraint to be as tight as possible, it is necessary to maximize the correlation between the two dipoles. This can be achieved by optimizing the survey window through which the flux dipole is measured. Here we explicitly construct such a window for the 2MASS survey. The optimization in essence reduces to excluding from the calculation of the flux dipole galaxies brighter than some limiting magnitude K_min of the near-infrared K_s band. This exclusion mitigates nonlinear effects and shot noise from small scales, which decorrelate the 2MASS dipole from the CMB dipole. Under the assumption of negligible shot noise we find that the optimal value of K_min is about five. Inclusion of shot noise shifts the optimal K_min to larger values. We present an analytical formula for shot noise for the 2MASS flux dipole, to be used in follow-up work with 2MASS data. The misalignment angle between the two dipoles is a sensitive measure of their correlation: the higher the correlation, the smaller the expectation value of the angle. A minimum of the misalignment is thus a sign of the optimal gravity window. We model analytically the distribution function for the misalignment angle and show that the misalignment estimated by Maller et al. is consistent with the assumed underlying model (though it is greater than the expectation value). We predict with about 90% confidence that the misalignment will decrease if 2MASS galaxies brighter than K_min = 5 mag are excluded from the calculation of the flux dipole. This prediction has been indirectly confirmed by the results of Erdogdu et al. (ABRIDGED)