The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

Observations with the Chandra X-ray Observatory are used to examine the hot gas properties within a sample of 10 galaxy groups selected from the Galaxy And Mass Assembly survey's optical Friends-of-Friends group catalogue. Our groups have been screened to eliminate spurious and unrelaxed systems, and the effectiveness of this procedure is demonstrated by the detection of intergalactic hot gas in 80 per cent of our sample. However, we find that 9 of the 10 are X-ray underluminous by a mean factor of $\sim$4 compared to typical X-ray-selected samples. Consistent with this, the majority of our groups have gas fractions that are lower and gas entropies somewhat higher than those seen in typical X-ray-selected samples. Two groups, which have high 2{\sigma} lower limits on their gas entropy, are candidates for the population of high-entropy groups predicted by some active galactic nucleus feedback models.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Detailed measurements of the CMB lensing signal are an important scientific goal of ongoing ground-based CMB polarization experiments, which are mapping the CMB at high resolution over small patches of the sky. In this work we simulate CMB polarization lensing reconstruction for the $EE$ and $EB$ quadratic estimators with current-generation noise levels and resolution, and show that without boundary effects the known and expected zeroth and first order $N^{(0)}$ and $N^{(1)}$ biases provide an adequate model for non-signal contributions to the lensing power spectrum estimators. Small sky areas present a number of additional challenges for polarization lensing reconstruction, including leakage of $E$ modes into $B$ modes. We show how simple windowed estimators using filtered pure-$B$ modes can greatly reduce the mask-induced mean-field lensing signal and reduce variance in the estimators. This provides a simple method (used with recent observations) that gives an alternative to more optimal but expensive inverse-variance filtering.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

In recent years cross correlation of lensing of the Cosmic Microwave Background (CMB) with other large scale structure (LSS) tracers has been used as a method to detect CMB lensing. Current experiments are also becoming sensitive enough to measure CMB lensing without the help of auxiliary tracers. As data quality improves rapidly, it has been suggested that the CMB lensing-LSS cross correlation may provide new insights into parameters describing cosmological structure growth. In this work we perform forecasts that combine the lensing potential auto power spectrum from various future CMB experiments, the galaxy power spectrum from galaxy surveys, as well as the cross power spectrum between the two, marginalizing over a number galactic and non-galactic cosmological parameters. We find that the CMB lensing-LSS cross correlation contains significant information on parameters such as the redshift distribution and bias of LSS tracers. We also find that the cross correlation information will lead to independent probes of cosmological parameters such as neutrino mass and the reionization optical depth.

We have measured Johnson BV and Cousins RI magnitudes for comparison stars near five TeV blazars. We compare our values with published values, spanning 25 years in some cases, to identify those stars that are most likely proven stable. To avoid zero-point offsets mimicking long-term variability, we based our analysis on the standard deviation between measurements after a mean offset between data sets was removed. We found most stars to be stable at the 0.04 magnitude level. We confirm two stars as variable and identify two others as possibly being variable. In each of the five fields there are at least two stars, and typically many more, that show no evidence of variability.

Observatories constructed solely for photometric monitoring make it possible to understand the temporal nature of objects over time scales that historically have been difficult to achieve. We report on one such observatory, the Remote Observatory for Variable Object Research (ROVOR) which was constructed to enable both long-term and rapid cadence observations of brighter objects. ROVOR is an 0.4m optical telescope located in central Utah and commissioned for scientific observations in 2008. Principle research has been monitoring blazars, x-ray binaries, AGN, and an occasional gamma-ray burst afterglow. We describe the observatory, the control system, and its unique roof.

Squeezed primordial non-Gaussianity can strongly constrain early-universe physics, but it can only be observed on the CMB after it has been gravitationally lensed. We give a new simple non-perturbative prescription for accurately calculating the effect of lensing on any squeezed primordial bispectrum shape, and test it with simulations. We give the generalization to polarization bispectra, and discuss the effect of lensing on the trispectrum. We explain why neglecting the lensing smoothing effect does not significantly bias estimators of local primordial non-Gaussianity, even though the change in shape can be >~10%. We also show how tau_NL trispectrum estimators can be well approximated by much simpler CMB temperature modulation estimators, and hence that there is potentially a ~10-30% bias due to very large-scale lensing modes, depending on the range of modulation scales included. Including dipole sky modulations can halve the tau_NL error bar if kinematic effects can be subtracted using known properties of the CMB temperature dipole. Lensing effects on the g_NL trispectrum are small compared to the error bar. In appendices we give the general result for lensing of any primordial bispectrum, and show how any full-sky squeezed bispectrum can be decomposed into orthogonal modes of distinct angular dependence.

The recurrent nova RS Oph undergoes nova eruptions every ~ 10-20years as a result of thermonuclear runaway on the surface of a white dwarf close to the Chandrasekhar limit. Both the progress of the eruption, and its aftermath, depend on the (poorly known) composition of the red giant in the RS Oph system. Our aim is to understand better the effect of the giant secondary on the recurrent nova eruption. Synthetic spectra were computed for a grid of M-giant model atmospheres having a range of effective temperatures 3200 < Teff < 4400 K, gravities 0 < log g < 1 and abundances -4 <[Fe/H] < 0.5, and fit to infrared spectra of RS Oph as it returned to quiescence after its 2006 eruption. We have modelled the infrared spectrum in the range 1.4-2.5 micron to determine metallicity and effective temperature of the red giant. We find Teff = 4100 +/- 100 K, log g = 0.0 +/- 0.5, [Fe/H] = 0.0 +/- 0.5, [C/H] = -0.8 +/- 0.2, [N/H] = +0.6 +/- 0.3 in the atmosphere of the secondary, and demonstrate that that inclusion of some dust `veiling' in the spectra cannot improve our fits.

We study the geometry and topology of the large-scale structure traced by galaxy clusters in numerical simulations of a box of side 320 $h^{-1}$ Mpc, and compare them with available data on real clusters. The simulations we use are generated by the Zel'dovich approximation, using the same methods as we have used in the first three papers in this series. We consider the following models to see if there are measurable differences in the topology and geometry of the superclustering they produce: (i) the standard CDM model (SCDM); (ii) a CDM model with $\Omega_0=0.2$ (OCDM); (iii) a CDM model with a `tilted' power spectrum having $n=0.7$ (TCDM); (iv) a CDM model with a very low Hubble constant, $h=0.3$ (LOWH); (v) a model with mixed CDM and HDM (CHDM); (vi) a flat low-density CDM model with $\Omega_0=0.2$ and a non-zero cosmological $\Lambda$ term ($\Lambda$CDM). We analyse these models using a variety of statistical tests based on the analysis of: (i) the Euler-Poincar\'{e} characteristic; (ii) percolation properties; (iii) the Minimal Spanning Tree construction. Taking all these tests together we find that the best fitting model is $\Lambda$CDM and, indeed, the others do not appear to be consistent with the data. Our results demonstrate that despite their biased and extremely sparse sampling of the cosmological density field, it is possible to use clusters to probe subtle statistical diagnostics of models which go far beyond the low-order correlation functions usually applied to study superclustering.

We propose and investigate a new algorithm for quantifying the topological properties of cosmological density fluctuations. We first motivate this algorithm by drawing a formal distinction between two definitions of relevant topological characteristics, based on concepts, on the one hand, from differential topology and, on the other, from integral geometry. The former approach leads one to concentrate on properties of the contour surfaces which, in turn, leads to the algorithms CONTOUR2D and CONTOUR3D familiar to cosmologists. The other approach, which we adopt here, actually leads to much simpler algorithms in both two and three dimensions. (The 2D algorithm has already been introduced to the astronomical literature.) We discuss the 3D case in some detail and compare results obtained with it to analogous results using the CONTOUR3D algorithm.