We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

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

We present a detailed spectral and temporal study of the intermediate-type blazar ON 231 during the TeV outburst phase in 2008 June with observations performed by Swift and XMM-Newton. The X-ray flux of the source, which was significantly dominated by the soft photons (below $3-4$ keV), varies between 27$\%$ and 38$\%$ on day timescales, while mild variations were observed in the optical/UV emissions. We found a maximum soft lag of $\sim 1$ hr between the UV and soft X-ray bands, which can be understood if the magnetic field of the emitting region is $\sim 5.6~ \delta^{-1/3}$ G. The $0.6-10$ keV spectra can be well represented by a broken power-law model, which indicates the presence of both synchrotron and inverse Compton components in the studied X-ray regime. The synchrotron part of the SEDs constructed with simultaneous optical/UV and X-ray data follows a log-parabolic shape. A time-resolved spectral analysis shows that the break energy varies significantly between 2.4 and 7.3 keV with the changing flux state of the source, and the similar variations of the spectral slopes of the two components support the SSC scenario. The synchrotron tail, following a log-parabolic function, shows that the peak frequency ($\nu_{p}$) varies by two orders of magnitude ($\sim 10^{14}-10^{16}$ Hz) during the event. A significantly positive $E_{p}-\beta$ relation is observed from both SED and time-resolved spectral analyses. The most feasible scenario for the observed trend during the flaring event could be associated with a magnetic-field-driven stochastic process evolving toward an equilibrium energy level.

We present a detailed temporal and spectral study of the intermediate type blazar, ON 231 with observations taken by Swift and XMM-Newton satellites during an outburst phase in June 2008. The source shows multiple X-ray flares in that period when the flux amplitude varies in the range of 27-38%. The X-ray spectra are well fitted with a broken power-law model indicating the presence of both synchrotron and inverse Compton components. We find that the source shows strong and variable emission in the soft energy bands (below 3-4 keV), where the hard band emission is weak and stable. All the soft X-ray bands show correlated variability with zero lag, while a soft lag of -600 s between the 0.3-0.5 and 4-10 keV bands is observed with DCF analysis. A time-resolved spectral study of the flares gives a positive relation between the total fluxes and the break energies of the two emission components.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

Active galactic nuclei (AGNs) can be divided into two major classes, namely radio loud and radio quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff & Riley type I/II (FRI/II) radio galaxies. Following our previous work (Fan et al. 2011), we present a sample of 2400 sources with measured radio flux densities of the core and extended components. The sample contains 250 BL Lacs, 520 quasars, 175 Seyferts, 1178 galaxies, 153 Fanaroff & Riley type I and type II galaxies and 104 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also confirm that the correlation between core-dominance parameter and spectral index exists for a large sample presented in this work.

The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope which will be online in the early-2020s; and with the development of a new data archive facility for a range of telescopes at the National Astronomical Research Institute of Thailand. The Newton Fund enabled us to collaborate in designing a new versatile generic system that serves all purposes. In the end, we conclude that a single system would not meet the needs of all parties and only adopt similar front-ends while the back-ends are bespoke to our respective systems and data-flows.

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and latitude of 24.66^{\circ}, 5.28^{\circ}, 25 degrees away from the Galactic Center. The burst was found to be 43$\pm$5% linearly polarized with a rotation measure (RM) in the range -9 < RM < 12 rad m^{-2} (95% confidence level), consistent with zero. The burst was followed-up with 11 telescopes to search for radio, optical, X-ray, gamma-ray and neutrino emission. Neither transient nor variable emission was found to be associated with the burst and no repeat pulses have been observed in 17.25 hours of observing. The sightline to the burst is close to the Galactic plane and the observed physical properties of FRB 150215 demonstrate the existence of sight lines of anomalously low RM for a given electron column density. The Galactic RM foreground may approach a null value due to magnetic field reversals along the line of sight, a decreased total electron column density from the Milky Way, or some combination of these effects. A lower Galactic DM contribution might explain why this burst was detectable whereas previous searches at low latitude have had lower detection rates than those out of the plane.

We use the lensing potential map from Planck CMB lensing reconstruction analysis and the "Public Cosmic Void Catalog" to measure the stacked void lensing potential. In this profile, four parameters are needed to describe the shape of voids with different characteristic radii $R_V$. However, we have found that after reducing the background noise by subtracting the average background, there is a residue lensing power left in the data. The inclusion of the environment shifting parameter, $\gamma_V$, is necessary to get a better fit to the data with the residue lensing power. We divide the voids into two redshift bins: cmass1 ($0.45 < z < 0.5$) and cmass2 ($0.5 < z < 0.6$). Our best-fit parameters are $\alpha = 1.989\pm0.149$, $\beta = 12.61\pm0.56$, $\delta_c=-0.697\pm0.025$, $R_S/R_V=1.039\pm0.030$, $\gamma_v=(-7.034\pm0.150) \times 10^{-2}$ for the cmass1 sample with 123 voids and $\alpha = 1.956\pm0.165$, $\beta = 12.91\pm0.60$, $\delta_c=-0.673\pm0.027$, $R_S/R_V=1.115\pm0.032$, $\gamma_v=(-4.512\pm0.114) \times 10^{-2}$ for the cmass2 sample with 393 voids at 68% C.L. The addition of the environment parameter is consistent with the conjecture that the Sloan Digital Sky Survey voids reside in an underdense region.

We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainter magnitudes, optical UVX selection is still required to reach our $g{\approx}22.5$ limit. Using both these techniques we observed quasar redshifts at sky densities up to $90$ deg$^{-2}$. By comparing 2QDESp with other surveys (SDSS, 2QZ and 2SLAQ) we find that quasar clustering is approximately luminosity independent, with results for all four surveys consistent with a correlation scale of $r_{0}{=}6.1{\pm}0.1 \: h^{-1}$Mpc, despite their decade range in luminosity. We find a significant redshift dependence of clustering, particularly when BOSS data with $r_{0}{=}7.3{\pm}0.1 \: h^{-1}$Mpc are included at $z{\approx}2.4$. All quasars remain consistent with having a single host halo mass of ${\approx}2{\pm}1{\times}10^{12} \: h^{-1}M_\odot$. This result implies that either quasars do not radiate at a fixed fraction of the Eddington luminosity or AGN black hole and dark matter halo masses are weakly correlated. No significant evidence is found to support fainter, X-ray selected quasars at low redshift having larger halo masses as predicted by the `hot halo' mode AGN model of Fanidakis et al. 2013. Finally, although the combined quasar sample reaches an effective volume as large as that of the original SDSS LRG sample, we do not detect the BAO feature in these data.

We investigate the potential of using cosmic voids as a probe to constrain cosmological parameters through the gravitational lensing effect of the cosmic microwave background (CMB) and make predictions for the next generation surveys. By assuming the detection of a series of $\approx 5 - 10$ voids along a line of sight within a square-degree patch of the sky, we found that they can be used to break the degeneracy direction of some of the cosmological parameter constraints (for example $\omega_b$ and $\Omega_\Lambda$) in comparison with the constraints from random CMB skies with the same size area for a survey with extensive integration time. This analysis is based on our current knowledge of the average void profile and analytical estimates of the void number function. We also provide combined cosmological parameter constraints between a sky patch where series of voids are detected and a patch without voids (a randomly selected patch). The full potential of this technique relies on an accurate determination of the void profile to $\approx 10$% level. For a small-area CMB observation with extensive integration time and a high signal-to-noise ratio, CMB lensing with such series of voids will provide a complementary route to cosmological parameter constraints to the CMB observations. Example of parameter constraints with a series of five voids on a $1.0^{\circ} \times 1.0^{\circ}$ patch of the sky are $100\omega_b = 2.20 \pm 0.27$, $\omega_c = 0.120 \pm 0.022$, $\Omega_\Lambda = 0.682 \pm 0.078$, $\Delta_{\mathcal{R}}^2 = \left(2.22 \pm 7.79\right) \times 10^{-9}$, $n_s = 0.962 \pm 0.097$ and $\tau = 0.925 \pm 1.747$ at 68% C.L.

The prime evidence underpinning the standard LCDM cosmological model is the CMB power spectrum as observed by WMAP and other microwave experiments. But Sawangwit & Shanks (2010) have recently shown that the WMAP CMB power spectrum is highly sensitive to the beam profile of the WMAP telescope. Here, we use the source catalogue from the Planck Early Data Release to test further the WMAP beam profiles. We confirm that stacked beam profiles at Q, V and particularly at W appear wider than expected when compared to the Jupiter beam, normalised either directly to the radio source profiles or using Planck fluxes. The same result is also found based on WMAP-CMBfree source catalogues and NVSS sources. The accuracy of our beam profile measurements is supported by analysis of CMB sky simulations. However the beam profiles from WMAP7 at the W band are narrower than previously found in WMAP5 data and the rejection of the WMAP beam is now only at the ~3sigma level. We also find that the WMAP source fluxes demonstrate possible non-linearity with Planck fluxes. But including ground-based and Planck data for the bright Weiland et al (2011) sources may suggest that the discrepancy is a linear offset rather than a non-linearity. Additionally, we find that the stacked Sunyaev-Zel'dovich (SZ) decrements of ~151 galaxy clusters observed by Planck are in agreement with the WMAP data. We find that there is no evidence for a WMAP SZ deficit as has previously been reported. In the particular case of Coma we find evidence for the presence of an ~0.1mK downwards CMB fluctuation. We conclude that beam profile systematics can have significant effects on both the amplitude and position of the acoustic peaks, with potentially important implications for cosmology parameter fitting.

We present a clustering analysis of Luminous Red Galaxies in SDSS Stripe 82. We study the angular 2-point correlation function, of 130,000 LRG candidates via colour-cut selections in izK with the K band coverage coming from UKIDSS LAS. We have used the cross-correlation technique of Newman (2008) to establish the LRG redshift distribution. Cross-correlating with SDSS QSOs, MegaZ-LRGs and DEEP2 galaxies implies an average LRG redshift of z~1 with space density, n_g~3.2 +/-0.16 x10^-4 h^3 Mpc^-3. For theta<10', w(theta) significantly deviates from a single power-law. A double power-law with a break at r_b~2.4 h^-1 Mpc fits the data better, with best-fit scale length, r_0,1=7.63+/-0.27 h^-1Mpc and slope gamma_1=2.01 +/-0.02 at small scales and r_0,2=9.92 +/-0.40 h^-1 Mpc and gamma_2=1.64 +/-0.04 at large scales. Due to the flat slope at large scales, we find that a standard LambdaCDM linear model is accepted only at 2-3sigma, with the best-fit bias factor, b=2.74+/-0.07. We also fitted HOD models and estimate an effective halo mass of M_eff=3.3 +/-0.6x10^13 h^-1 M_sun. But at large scales, the current HOD models did not help explain the power excess in the clustering signal. We then compare the w(theta) results to those of Sawangwit et al. (2011) from 3 samples of LRGs at lower redshifts to measure clustering evolution. We find that a long-lived model may be a poorer fit than at lower redshifts, although this assumes that the Stripe 82 LRGs are luminosity-matched to the AAOmega LRGs. We find stronger evidence for evolution in the form of the z~1 LRG correlation function, with the above flat 2-halo slope maintaining to r>50 h^-1 Mpc. Applying the cross-correlation test of Ross et al. (2011), we find little evidence that the result is due to systematics. Otherwise it may provide evidence for primordial non-Gaussianity in the matter distribution, with f^local_NL=90+/-30.[Abridged]

Our goals are (i) to search for BAO and large-scale structure in current QSO survey data and (ii) to use these and simulation/forecast results to assess the science case for a new, >10x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ surveys to form a survey of ~60000 QSOs. We find a hint of a peak in the QSO 2-point correlation function, xi(s), at the same scale (~105h^-1 Mpc) as detected by Eisenstein et al (2005) in their sample of DR5 LRGs but only at low statistical significance. We then compare these data with QSO mock catalogues from the Hubble Volume simulation used by Hoyle et al (2002) and find that both routes give statistical error estimates that are consistent at ~100h^-1 Mpc scales. Mock catalogues are then used to estimate the nominal survey size needed for a 3-4 sigma detection of the BAO peak. We find that a redshift survey of ~250000 z<2.2 QSOs is required over ~3000 deg^2. This is further confirmed by static log-normal simulations where the BAO are clearly detectable in the QSO power spectrum and correlation function. The nominal survey would on its own produce the first detection of, for example, discontinuous dark energy evolution in the so far uncharted 1<z<2.2 redshift range. A survey with ~50% higher QSO sky densities and 50% bigger area will give an ~6sigma BAO detection, leading to an error ~60% of the size of the BOSS error on the dark energy evolution parameter, w_a. Another important aim for a QSO survey is to place new limits on primordial non-Gaussianity at large scales, testing tentative evidence we have found for the evolution of the linear form of the combined QSO xi(s) at z~1.6. Such a QSO survey will also determine the gravitational growth rate at z~1.6 via z-space distortions, allow lensing tomography via QSO magnification bias while also measuring the exact luminosity dependence of small-scale QSO clustering.

We cross match the NVSS and FIRST surveys with three large photometric catalogues of luminous red galaxies (LRGs) to define radio-loud samples. These have median redshifts 0.35, 0.55 and 0.68 and, by matching rest-frame optical and radio properties, we construct uniform samples across the three surveys. This paper is concerned with the clustering properties of these samples derived from the angular correlation function. The primary aim is to characterise any evolution in the clustering amplitude of radio galaxies bellow z~0.68. We find no evidence for evolution in the large-scale (~1-50h^{-1}Mpc) clustering amplitude. Our radio galaxy autocorrelations are consistent with previous findings indicating little-to-no evolution in the redshift range 0.68 to 0 (~6Gyr of time). We also cross correlate radio galaxies with the parent LRG samples to increase the precision of our results and again find no evidence for evolution. Our results are inconsistent with a long-lived model for the clustering evolution that assumes radio sources randomly sample the LRG population. A model where the halo mass is constant with redshift is consistent with the data. This is similar to QSOs that have clustering amplitudes consistent with a single halo mass at all redshifts. Given that the brightest radio sources show stronger evolution in space density compared to fainter radio sources we restrict our samples to include only objects with L>10^{26}W/Hz and repeat the analysis. Again we find no evidence for evolution in the comoving correlation amplitude. These radio sources appear to inhabit the same mass halos as fainter radio galaxies (~9x10^{13}h^{-1}M_\odot). These halos are ~twice as massive as those of the general LRG population and ~30 times as massive as optical AGN/QSOs.

We present the angular correlation function measured from photometric samples comprising 1,562,800 luminous red galaxies (LRGs). Three LRG samples were extracted from the Sloan Digital Sky Survey (SDSS) imaging data, based on colour-cut selections at redshifts, z ~ 0.35, 0.55 and 0.7 as calibrated by the spectroscopic surveys, SDSS-LRG, 2dF-SDSS LRG and QSO (2SLAQ), and the AAOmega LRG survey. The galaxy samples cover ~7600 s.q.deg of the sky, probing a total cosmic volume of ~5.5 h^-3} Gpc^3. The small and intermediate scale correlation functions generally show significant deviations from a single power-law fit with a well-detected break at ~1 h^-1 Mpc, consistent with the transition scale between the 1- and 2-halo terms in halo occupation models. For galaxy separations 1-20 h^-1 Mpc and at fixed luminosity, we see virtually no evolution of the clustering with redshift and the data is consistent with a simple high peaks biasing model where the comoving LRG space density is constant with z. At fixed z, the LRG clustering amplitude increases with luminosity in accordance with the simple high peaks model, with a typical LRG dark matter halo mass 10^13-10^14 h^-1 M_sol. For r < 1 h^-1 Mpc, the evolution is slightly faster and the clustering decreases towards high redshift consistent with a virialised clustering model. However, assuming the HOD and LCDM halo merger frameworks, ~2-3 per cent per Gyr of the LRGs is required to merge in order to explain the small scales clustering evolution, consistent with previous results [abridge].

QSOs from SDSS, 2QZ and 2SLAQ covering an order of magnitude in luminosity at fixed redshift exhibit similar amplitudes of clustering. In addition, QSO clustering evolution at z>0.5 is well fitted by a model that assumes a fixed host halo mass, implying that QSOs may occur in a relatively narrow range of halo and BH mass. We argue that the slow evolution of early-type galaxies out to z~1-2 may also provide support for a slow evolution of QSO host BH masses. The result would mean that if high-z QSOs radiate at Eddington rates then low-z SyI must radiate at ~100x less than Eddington. We conclude that models where QSOs radiate at L_Edd require M_BH and M_halo to be decoupled to circumvent the clustering results. While single BH mass and flickering models fit the z>0.5 clustering results, they appear to be rejected by the z~0, M_BH-L relation from reverberation mapping. We find that the inclusion of z<0.5 QSO clustering data improves the fit of a long-lived QSO model and suggest that the predictions of a PLE model for QSO BH masses agree reasonably with UV-bump and reverberation estimates (abridged).

Using the published WMAP 5-year data, we first show how sensitive the WMAP power spectra are to the form of the WMAP beam. It is well known that the beam profile derived from observations of Jupiter is non-Gaussian and indeed extends, in the W band for example, well beyond its 12.'6 FWHM core out to more than 1 degree in radius. This means that even though the core width corresponds to wavenumber l\approx1800, the form of the beam still significantly affects the WMAP results even at l\approx200 which is the scale of the first acoustic peak. The difference between the beam convolved C_l and the final C_l is \approx70% at the scale of the first peak, rising to \approx400% at the scale of the second. New estimates of the Q, V and W-band beam profiles are then presented, based on a stacking analysis of the WMAP5 radio source catalogue and temperature maps. The radio sources show a significantly (3-4\sigma) broader beam profile on scales of 10'-30' than that found by the WMAP team whose beam analysis is based on measurements of Jupiter. Beyond these scales the beam profiles from the radio sources are too noisy to give useful information. Furthermore, we find tentative evidence for a non-linear relation between WMAP and ATCA/IRAM 95 GHz source fluxes. We discuss whether the wide beam profiles could be caused either by radio source extension or clustering and find that neither explanation is likely. We also argue against the possibility that Eddington bias is affecting our results. The reasons for the difference between the radio source and the Jupiter beam profiles are therefore still unclear. If the radio source profiles were then used to define the WMAP beam, there could be a significant change in the amplitude and position of even the first acoustic peak. It is therefore important to identify the reasons for the differences between these two beam profile estimates.

We first discuss the sensitivity of the WMAP CMB power spectrum to systematic errors by calculating the raw CMB power spectrum from WMAP data. We find that the power spectrum is surprisingly sensitive to the WMAP radiometer beam profile even at the position of the first acoustic peak on ~1 degree scales. Although the WMAP beam profile core is only 12.6arcmin FWHM at W, there is a long power-law tail to the beam due to side-lobes and this causes significant effects even at the first peak position. We then test the form of the beam-profile used by the WMAP team which is based on observations of Jupiter. We stacked radio source beam profiles as observed in each WMAP band and found that they showed a wider profile in Q, V, W than the Jupiter profile. We have now checked that this is not due to any Eddington or other bias in our sample by showing that the same results are obtained when radio sources are selected at 1.4GHz and that our methods retrieve the Jupiter beam when it is employed in simulations. Finally, we show that the uncertainty in the WMAP beam profile allows the position as well as the amplitude of the first peak to be changed and how this could allow simpler cosmologies than standard Lambda-CDM to fit the CMB data.

We present the cross-correlation of the density map of LRGs and the temperature fluctuation in the CMB as measured by the WMAP5 observations. The LRG samples were extracted from imaging data of the SDSS based on two previous spectroscopic redshift surveys, the SDSS-LRG and the 2SLAQ surveys at average redshifts z~0.35 and z~0.55. In addition we have added a higher-redshift photometric LRG sample based on the selection of the AAOmega LRG redshift survey at z~0.7. The total LRG sample thus comprises 1.5 million galaxies, sampling a redshift range of 0.2 < z < 0.9 over ~7600 square degrees of the sky, probing a total cosmic volume of ~5.5 h^{-3} Gpc^3. We find that the new LRG sample at z~0.7 shows very little positive evidence for the ISW effect. Indeed, the cross-correlation is negative out to ~1 deg. The standard LCDM model is rejected at ~2-3% significance by the new LRG data. We then performed a new test on the robustness of the LRG ISW detections at z~0.35 and 0.55. We made 8 rotations through 360deg of the CMB maps with respect to the LRG samples around the galactic pole. We find that in both cases there are stronger effects at angles other than zero. This implies that the z~0.35 and 0.55 ISW detections may still be subject to systematic errors which combined with the known sizeable statistical errors may leave these ISW detections looking unreliable. We have further made the rotation test on several other samples where ISW detections have been claimed and find that they also show peaks when rotated. We conclude that in the samples we have tested the ISW effect may be absent and we argue that this result may not be in contradiction with previous results.

We investigate the use of simple colour cuts applied to the SDSS optical imaging to perform photometric selections of emission line galaxies out to z<1. From colour-cuts using the SDSS g, r and i bands, we obtain mean photometric redshifts of z=0.32+-0.08, z=0.44+-0.12 and z=0.65+-0.21. We further calibrate our high redshift selection using spectroscopic observations with the AAOmega spectrograph on the 4m Anglo-Australian Telescope (AAT), observing ~50-200 galaxy candidates in 4 separate fields. With just 1-hour of integration time and with seeing of ~1.6", we successfully determined redshifts for ~65% of the targeted candidates. We calculate the angular correlation functions of the samples and find correlation lengths of r0=2.64 h-1 Mpc, r0=3.62 h-1 Mpc and r0=5.88 h-1 Mpc for the low, mid and high redshift samples respectively. Comparing these results with predicted dark matter clustering, we estimate the bias parameter for each sample to be b=0.70, b=0.92 and b=1.46. We calculate the 2-point redshift-space correlation function at z~0.6 and find a clustering amplitude of s0=6.4 h-1 Mpc. Finally, we use our photometric sample to search for the Integrated Sachs-Wolfe signal in the WMAP 5yr data. We cross-correlate our three redshift samples with the WMAP W, V, Q and K bands and find an overall trend for a positive signal similar to that expected from models. However, the signal in each is relatively weak. Combining all three galaxy samples we find a signal of wTg(<100')=0.20+-0.12 microK in the WMAP W-band, a significance of 1.7sigma.

We have measured the bias of QSOs as a function of QSO luminosity at fixed redshift (z<1) by cross-correlating them with LRGs in the same spatial volume, hence breaking the degeneracy between QSO luminosity and redshift. We use three QSO samples from 2SLAQ, 2QZ and SDSS covering a QSO absolute magnitude range, -24.5<M_{b_J}<-21.5, and cross-correlate them with 2SLAQ (z~0.5) and AAOmega (z~0.7) photometric and spectroscopic LRGs in the same redshift ranges. The 2-D and 3-D cross-clustering measurements are generally in good agreement. Our (2SLAQ) QSO-LRG clustering amplitude (r_0=6.8_{-0.3}^{+0.1}h^{-1}Mpc) as measured from the semi-projected cross-correlation function appears similar to the (2SLAQ) LRG-LRG auto-correlation amplitude (r_0=7.45\pm0.35h^{-1}Mpc) and both are higher than the (2QZ+2SLAQ) QSO-QSO amplitude (r_0\simeq5.0h^{-1}Mpc). Our measurements show remarkably little QSO-LRG cross-clustering dependence on QSO luminosity. If anything, the results imply that brighter QSOs may be less highly biased than faint QSOs, the opposite direction expected from simple high peaks biasing models. Assuming a standard LCDM model and values for b_{LRG} measured from LRG autocorrelation analyses, we find b_Q=1.45\pm0.11 at M_{b_J}\approx-24 and b_Q=1.90\pm0.16 at M_{b_J}~-22. We also find consistent results for the QSO bias from a z-space distortion analysis of the QSO-LRG cross-clustering at z~0.55. The dynamical infall results give \beta _Q=0.55\pm0.10, implying b_Q=1.4\pm0.2. Thus both the z-space distortion and the amplitude analyses yield b_Q~1.5 at M_{b_J}~-23. The implied DM halo mass inhabited by QSOs at z~0.55 is \sim10^{13}h^{-1}M_{\sun}, again approximately independent of QSO luminosity.