We examine the integrated Sachs-Wolfe (ISW) imprint of voids and superclusters on the cosmic microwave background. We first study results from the Jubilee $N$-body simulation. From Jubilee, we obtain the full-sky ISW signal from structures out to redshift $z=1.4$ and a mock luminous red galaxy (LRG) catalogue. We confirm that the expected signal in the concordance \Lambda CDM model is very small and likely to always be much smaller than the anisotropies arising at the last scattering surface. Any current detections of such an imprint must, therefore, predominantly arise from something other than an ISW effect in a \Lambda CDM universe. Using the simulation as a guide, we then look for the signal using a catalogue of voids and superclusters from the Sloan Digital Sky Survey. We find a result that is consistent with the \Lambda CDM model, i.e. a signal consistent with zero.

We discuss the universality and self-similarity of void density profiles, for voids in realistic mock luminous red galaxy (LRG) catalogues from the Jubilee simulation, as well as in void catalogues constructed from the SDSS LRG and Main Galaxy samples. Voids are identified using a modified version of the ZOBOV watershed transform algorithm, with additional selection cuts. We find that voids in simulation are self-similar, meaning that their average rescaled profile does not depend on the void size, or -- within the range of the simulated catalogue -- on the redshift. Comparison of the profiles obtained from simulated and real voids shows an excellent match. The profiles of real voids also show a universal behaviour over a wide range of galaxy luminosities, number densities and redshifts. This points to a fundamental property of the voids found by the watershed algorithm, which can be exploited in future studies of voids.

The stacked density profile of cosmic voids in the galaxy distribution provides an important tool for the use of voids for precision cosmology. We study the density profiles of voids identified using the ZOBOV watershed transform algorithm in realistic mock luminous red galaxy (LRG) catalogues from the Jubilee simulation, as well as in void catalogues constructed from the SDSS LRG and Main Galaxy samples. We compare different methods for reconstructing density profiles scaled by the void radius and show that the most commonly used method based on counts in shells and simple averaging is statistically flawed as it underestimates the density in void interiors. We provide two alternative methods that do not suffer from this effect; one based on Voronoi tessellations is also easily able to account from artefacts due to finite survey boundaries and so is more suitable when comparing simulation data to observation. Using this method we show that voids in simulation are exactly self-similar, meaning that their average rescaled profile does not depend on the void size. Within the range of our simulation we also find no redshift dependence of the profile. Comparison of the profiles obtained from simulated and real voids shows an excellent match. The profiles of real voids also show a universal behaviour over a wide range of galaxy luminosities, number densities and redshifts. This points to a fundamental property of the voids found by the watershed algorithm, which can be exploited in future studies of voids.

The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case.

In this paper we investigate how the halo mass function evolves with redshift, based on a suite of very large (with N_p = 3072^3 - 6000^3 particles) cosmological N-body simulations. Our halo catalogue data spans a redshift range of z = 0-30, allowing us to probe the mass function from the dark ages to the present. We utilise both the Friends-of-Friends (FOF) and Spherical Overdensity (SO) halofinding methods to directly compare the mass function derived using these commonly used halo definitions. The mass function from SO haloes exhibits a clear evolution with redshift, especially during the recent era of dark energy dominance (z < 1). We provide a redshift-parameterised fit for the SO mass function valid for the entire redshift range to within ~20% as well as a scheme to calculate the mass function for haloes with arbitrary overdensities. The FOF mass function displays a weaker evolution with redshift. We provide a `universal' fit for the FOF mass function, fitted to data across the entire redshift range simultaneously, and observe redshift evolution in our data versus this fit. The relative evolution of the mass functions derived via the two methods is compared and we find that the mass functions most closely match at z=0. The disparity at z=0 between the FOF and SO mass functions resides in their high mass tails where the collapsed fraction of mass in SO haloes is ~80% of that in FOF haloes. This difference grows with redshift so that, by z>20, the SO algorithm finds a ~50-80% lower collapsed fraction in high mass haloes than does the FOF algorithm, due in part to the significant over-linking effects known to affect the FOF method.

We present the first results from the JUropa huBbLE volumE (Jubilee) project, based a large N-body, dark matter-only cosmological simulation with a volume of $V=(6 h^{-1}\mathrm{Gpc})^3$, containing 6000$^3$ particles, performed within the concordance $\Lambda$CDM cosmological model. The simulation volume is sufficient to probe extremely large length scales in the universe, whilst at the same time the particle count is high enough so that dark matter haloes down to $1.5\times10^{12} h^{-1}\mathrm{M}_\odot$ can be resolved. At $z = 0$ we identify over 400 million haloes. The cluster mass function is derived using three different halofinders and compared to fitting functions in the literature. The distribution of clusters of maximal mass across redshifts agrees well with predicted masses of extreme objects, and we explicitly confirm that the Poisson distribution is very good at describing the distribution of rare clusters. The Poisson distribution also matches well the level to which cosmic variance can be expected to affect number counts of high mass clusters. We find that objects like the Bullet cluster exist in the far-tail of the distribution of mergers in terms of relative collisional speed. We also derive the number counts of voids in the simulation box for $z = 0$, $0.5$ and $1$.

We present initial results from the Jubilee ISW project, which models the expected \LambdaCDM Integrated Sachs-Wolfe (ISW) effect in the Jubilee simulation. The simulation volume is (6 Gpc/h)^3, allowing power on very large-scales to be incorporated into the calculation. Haloes are resolved down to a mass of 1.5x10^12 M_sun/h, which allows us to derive a catalogue of mock Luminous Red Galaxies (LRGs) for cross-correlation analysis with the ISW signal. We find the ISW effect observed on a projected sky to grow stronger at late times with the evolution of the ISW power spectrum matching expectations from linear theory. Maps of the gravitational lensing effect, including the convergence and deflection fields, are calculated using the same potential as for the ISW. We calculate the redshift dependence of the ISW-LRG cross-correlation signal for a full sky survey with no noise considerations. For l < 30, the signal is strongest for lower redshift bins (z ~ 0.2 to 0.5), whereas for l > 30 the signal is best observed with surveys covering z ~ 0.6-1.0.

The XMM Cluster Survey (XCS) is a serendipitous search for galaxy clusters using all publicly available data in the XMM-Newton Science Archive. Its main aims are to measure cosmological parameters and trace the evolution of X-ray scaling relations. In this paper we present the first data release from the XMM Cluster Survey (XCS-DR1). This consists of 503 optically confirmed, serendipitously detected, X-ray clusters. Of these clusters, 255 are new to the literature and 356 are new X-ray discoveries. We present 464 clusters with a redshift estimate (0.06 < z < 1.46), including 261 clusters with spectroscopic redshifts. In addition, we have measured X-ray temperatures (Tx) for 402 clusters (0.4 < Tx < 14.7 keV). We highlight seven interesting subsamples of XCS-DR1 clusters: (i) 10 clusters at high redshift (z > 1.0, including a new spectroscopically-confirmed cluster at z = 1.01); (ii) 67 clusters with high Tx (> 5 keV); (iii) 131 clusters/groups with low Tx (< 2 keV); (iv) 27 clusters with measured Tx values in the SDSS `Stripe 82' co-add region; (v) 78 clusters with measured Tx values in the Dark Energy Survey region; (vi) 40 clusters detected with sufficient counts to permit mass measurements (under the assumption of hydrostatic equilibrium); (vii) 105 clusters that can be used for applications such as the derivation of cosmological parameters and the measurement of cluster scaling relations. The X-ray analysis methodology used to construct and analyse the XCS-DR1 cluster sample has been presented in a companion paper, Lloyd-Davies et al. (2010).

Basic aspects of the relationship between the magnetic field and polarized maser radiation are described with the emphasis on interpreting the observed spectra. Special attention is given to three issues--the limitations on the applicability of the classic solutions of Goldreich, Keeley & Kwan (1973), inferring the strength of the magnetic field from the circular polarization when the Zeeman splitting is much less than the spectral linebreadth (especially for SiO masers), and the significance of the absence of components of the Zeeman triplet in the spectra of OH masers in regions of star formation.

The anisotropy of velocities in MHD turbulence is demonstrated explicitly by calculating the velocity gradients as a function of direction in representative simulations of decaying turbulence. It follows that the optical depths of spectral lines are anisotropic when there is MHD turbulence, and that this anisotropy influences the polarization characteristics of the emergent radiation. We calculate the linear polarization that results for the microwave lines of the CO molecule in star-forming gas and show that it is comparable to the polarization that is observed. This and our earlier result--that the anisotropy of MHD turbulence may be the cause for the absence of the Zeeman $\pi $-components in the spectra of OH mainline masers--are the first demonstrations of the occurrence of anisotropy in the optical depths caused by MHD turbulence. A non-local approximation is developed for the radiative transfer and the results are compared with those from a local (LVG) approximation.

We present polarization observations in DR21(OH) from thermal dust emission at 3 mm and from CO J=1-0 line emission. The observations were obtained using the Berkeley-Illinois-Maryland Association array. Previous observations of this region at 1.3 mm for the polarized continuum emission, and also the CO J=2-1 polarized emission have been made before. Our continuum polarization results are consistent with those results. However, the direction of the linear polarization for the J=1-0 is perpendicular to that of the CO J=2-1 polarization. This unexpected result was explored by obtaining numerical solutions to the multilevel, radiative transfer equations for a gas with anisotropic optical depths. We find that in addition to the anisotropic optical depths, anisotropic excitation due to a source of radiation that is external to the CO is needed to understand the orthogonality in the directions of polarization. The continuum emission by dust grains at the core of DR21(OH) is sufficient to provide this external radiation. The CO polarization must arise in relatively low density (nH2 ~ 100 [cm^-3]) envelope gas. We infer B ~ 10 [micro Gauss] in this gas, which implies that the envelope is subcritical.

(abridged) The polarized thermal radiation at far infrared and submillimeter wavelengths from dust grains in interstellar clouds with irregular magnetic fields is simulated. We show that the polarization properties of quiescent clouds and of star-forming regions are consistent with Kolmogorov-like turbulent magnetic fields that are comparable in magnitude to the uniform component of the magnetic fields. If the beam size is much smaller than the correlation length $L_{\rm corr}$ of the fields, the calculated percentage polarization $p$ decreases to an asymptotic value when the number of correlation lengths $N_{\rm corr}$ through the cloud exceeds a few tens. For these values of $N_{\rm corr}$, the dispersion in the position angles $\sigma_\alpha$ is still appreciable--decreasing to only about $20^{\circ}$. However, when the finite size of a telescope beam is taken into account, the asymptotic value of $p$ is reached for fewer correlation lengths (smaller $N_{\rm corr}$) due to averaging over the beam; $\sigma_\alpha$ becames much smaller and consistent with the observational data. The smoothing of the polarization properties due to the combined effect of the thickness of the cloud and the finite size of the beam can be described by a single variable which we designate as the generalized number of correlation lengths. In addition, we study various factors that may contribute to the decrease in the linear polarization percentage with increasing intensity that is observed at submillimeter and far infrared wavelengths in many, though not in all, dark clouds (the "polarization hole" effect).

We consider astrophysical maser radiation that is created in the presence of mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on the OH masers for which the magnetic field is strong enough that the separations of the Zeeman components are greater than the spectral linebreadths. A longstanding puzzle has been the absence of the Zeeman $\pi $ components and the high circular polarization in the observed spectra of these masers. We first argue that the elongation of eddies along the field that has recently been recognized in MHD turbulence will enhance the optical depth parallel to the magnetic field in comparison with that perpendicular to the magnetic field. We then simulate maser emission with a numerical model of MHD turbulence to demonstrate quantitatively how the intensities of the linearly polarized $\pi $ components are suppressed and the intensities of the nearly circularly polarized $\sigma $ components are enhanced. This effect is also generic in the sense that most spectral lines in MHD turbulence with Mach number $M\sim 1$ should have larger optical depths parallel to the magnetic field than perpendicular. The effect is reduced considerably when $M<1$. The simulations also demonstrate that the velocity and magnetic field variations due to the turbulence can (but do not necessarily) cause one of the $\sigma $ components to be much more intense than the other, as is often observed for mainline OH masers.

Calculations are performed for the spectral line profiles and images of astrophysical maser radiation that emerges from isolated spheres and thin disks viewed edge-on. In contrast to previous investigations where various approximations are made, the full equations are solved here for the frequency-dependent radiative transport that includes the thermal motion of the molecules. The spectral line profiles for spheres and disks are found to rebroaden to the full thermal Doppler breadth with increasing saturation in essentially the same way as is well known to occur for a linear maser. The variation with frequency in the apparent angular sizes of masing spheres and thin disks is found to be negligible at frequencies within the spectral line where the flux is significant. Calculations also are performed for spherical and disk masers that are not isolated, but for which the seed radiation for the masers is incident from one side as would occur when a strong continuum source is on the far side of the masers. Again, the spectral line profiles are found to rebroaden to the full thermal breadths with increasing saturation and there are no significant variations in the apparent angular sizes with frequency. However, the full rebroadening does occur at somewhat higher saturation and the variation of the apparent angular sizes as a function of the degree of saturation is quite different from that of the isolated masers. Spheres and disks have served as idealized geometries with which to examine possible deviations from the linear approximation for astrophysical masers.

Calculations are performed to demonstrate the deviations from Gaussian that occur in the spectral line profiles of a linear maser as a result of the amplification process. Near-Gaussian profiles are presented for bright, interstellar 22 GHz water masers obtained from high resolution Very Long Baseline Array (VLBA) observations of W3 IRS 5. For the profiles to be so close to Gaussian, the calculations indicate that these masers must originate in quite hot gas with temperatures greater than 1200 K -- a conclusion that is supportive of C-type shocks as the origin of these masers. In addition, the degree of saturation of these masers must be less than approximately one-third, from which it follows that the beaming angles are less than about 10^{-4} ster and the actual luminosities are modest. If spectral profiles that are as close to Gaussian as the profiles presented in this initial investigation are found to occur widely, they can be valuable diagnostics for the environments of astrophysical masers.

The relationship between the magnetic field and the circular polarization of astrophysical maser radiation due to the Zeeman effect under idealized conditions is investigated when the Zeeman splitting is much smaller than the spectral linebreadth and when radiative saturation is significant. The description of the circular polarization, as well as inferences about the magnetic field from the observations, are clearest when the rate for stimulated emission is much less than the Zeeman splitting. The calculations here are performed in this regime, which is relevant for some (if not most) observations of astrophysical masers. We demonstrate that Stokes-$V$ is proportional to the Zeeman splitting and the fractional linear polarization is independent of the Zeeman splitting when the ratio of the Zeeman splitting to the spectral line breadth is small---less than about 0.1. In contrast to its behavior for ordinary spectral lines, the circular polarization for masers that are at least partially saturated does not decrease with increasing angle between the magnetic field and the line-of-sight until they are nearly perpendicular.

Theoretical research on extragalactic masers (and hence this review) is almost entirely directed toward 22 GHz water masers in circumnuclear disks with the focus being the near ideal masing disk at the nucleus of the galaxy NGC4258. The discussion here is organized around (1) the excitation and conditions for the masers, (2) the spatial and spectral appearance of the maser emission, and (3) the mass accretion rate and structure of the disk. In addition, a summary is given of the basic physics (spectral linebeadths and maser polarization) that underlies certain interpretations of these and other astrophysical masers.

Calculations are performed for the circular polarization of maser radiation from a turbulent, Keplerian disk that is intended to represent the sub-parsec disk at the nucleus of the galaxy NGC4258. The polarization in the calculations is a result of the Zeeman effect in the regime in which the Zeeman splitting is much less than the spectral linebreadth. Plausible configurations for turbulent magnetic and velocity fields in the disk are created by statistical methods. This turbulence, along with the Keplerian velocity gradients and the blending of the three hyperfine components to form the $6_{16} - 5_{23}$ masing transition of water, are key ingredients in determining the appearance of the polarized spectra that are calculated. These spectra are quite different from the polarized spectra that would be expected for a two-level transition where there is no hyperfine structure. The effect of the hyperfine structure on the polarization is most striking in the calculations for the maser emission that represents the central (or systemic) features of NGC4258. Information about magnetic fields is inferred from observations for polarized maser radiation and bears on the structure of accretion disks.

Calculations are performed for the linear polarization of starlight due to extinction by aligned dust grains when the starlight traverses a medium with irregular magnetic fields. This medium is intended to represent the optically thick components of interstellar clouds which are observed to make little, if any, contribution to the polarization of starlight. In agreement with the observations, we find that the polarization properties of the starlight---the average fractional polarization and the dispersion in position angles---can be essentially unchanged. For this, the rms of the irregular component must be greater than the average magnetic field, which in turn tends to imply that the turbulent velocities in these interstellar clouds are super-Alfvenic.

The strengths of magnetic fields in interstellar gas clouds are obtained through observations of the circular polarization of spectral line radiation. Irregularities in this magnetic field may be present due to turbulence, waves or perhaps other causes, and may play an essential role in the structure and evolution of the gas clouds. To infer information about these irregularities from the observational data, we develop statistical relationships between the rms values of the irregular component of the magnetic field and spatial variations in the circular polarization of the spectral line radiation. The irregularities are characterized in analogy with descriptions of turbulence---by a sum of Fourier waves having a power spectrum with a slope similar to that of Kolmogorov turbulence. For comparison, we also perform computations in which turbulent magnetic and velocity fields from representative MHD simulations by others are utilized. Although the effects of the variations about the mean value of the magnetic field along the path of a ray tend to cancel, a significant residual effect in the polarization of the emergent radiation remains for typical values of the relevant parameters. A map of observed spectra of the 21 cm line toward Orion A is analyzed and the results are compared with our calculations in order to infer the strength of the irregular component of the magnetic field. The rms of the irregular component is found to be comparable in magnitude to the mean magnetic field within the cloud. Hence, the turbulent and Alfven velocities should also be comparable.

The linear polarization that is observed, together with likely changes in the orientation of the magnetic field along the line of sight and hence of the optical axes of the medium, can lead to the circular polarization that is observed in the radiation of the circumstellar SiO masers. A magnetic field greater than only about 30 mG is required, in contrast to 10-100 G that would be implied by the Zeeman interpretation. To assess quantitatively the likely changes in orientation of the magnetic field, calculations are performed with representative field configurations that are created by statistical sampling using a Kolmogorov-like power spectrum.