The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled baryonic and dark matter components in massive halos, and its knowledge will enhance cosmological parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present a computational study of intrinsic covariance in cluster properties using halo populations derived from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 Mpc/h Millennium Simulation performed with gas dynamics under two different physical treatments: shock heating driven by gravity only (GO) and a second treatment with cooling and preheating (PH). We examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) signals for samples of thousands of halos with M_200 > 5 \times 10^{13} Msun/h and z < 2. While the X-ray scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas observable--mass relations and ~10% departures from self-similar redshift evolution for 10^14 Msun/h halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs display mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction (r=0.88/0.69 for PH/GO at z=0) or X-ray temperature (r=0.62/0.83). We discuss halo mass selection by signal pairs, and find a minimum mass scatter of 4% in the \PH model by combining thermal SZ and gas fraction measurements.

Cosmological tests based on cluster counts require accurate calibration of the space density of massive halos, but most calibrations to date have ignored complex gas physics associated with halo baryons. We explore the sensitivity of the halo mass function to baryon physics using two pairs of gas-dynamic simulations that are likely to bracket the true behavior. Each pair consists of a baseline model involving only gravity and shock heating, and a refined physics model aimed at reproducing the observed scaling of the hot, intracluster gas phase. One pair consists of billion-particle re-simulations of the original 500 Mpc/h Millennium Simulation of Springel et al. (2005), run with the SPH code Gadget-2 and using a refined physics treatment approximated by preheating (PH) at high redshift. The other pair are high-resolution simulations from the adaptive-mesh refinement code ART, for which the refined treatment includes cooling, star formation, and supernova feedback (CSF). We find that, although the mass functions of the gravity-only (GO) treatments are consistent with the recent calibration of Tinker et al. (2008), both pairs of simulations with refined baryon physics show significant deviations. Relative to the GO case, the masses of ~10^{14} Msun/h halos in the PH and CSF treatments are shifted by averages of -15 \pm 1 percent and +12 \pm 5 percent, respectively. These mass shifts cause ~ 30% deviations in number density relative to the Tinker function, significantly larger than the 5% statistical uncertainty of that calibration.

We present a new measurement of the scaling relation between X-ray luminosity and total mass for 17,000 galaxy clusters in the maxBCG cluster sample. Stacking sub-samples within fixed ranges of optical richness, N_200, we measure the mean 0.1-2.4 keV X-ray luminosity, <L_X>, from the ROSAT All-Sky Survey. The mean mass, <M_200>, is measured from weak gravitational lensing of SDSS background galaxies (Johnston et al. 2007). For 9 <= N_200 < 200, the data are well fit by a power-law, <L_X>/10^42 h^-2 erg/s = (12.6+1.4-1.3 (stat) +/- 1.6 (sys)) (<M_200>/10^14 h^-1 M_sun)^1.65+/-0.13. The slope agrees to within 10% with previous estimates based on X-ray selected catalogs, implying that the covariance in L_X and N_200 at fixed halo mass is not large. The luminosity intercent is 30%, or 2\sigma, lower than determined from the X-ray flux-limited sample of Reiprich & Bohringer (2002), assuming hydrostatic equilibrium. This difference could arise from a combination of Malmquist bias and/or systematic error in hydrostatic mass estimates, both of which are expected. The intercept agrees with that derived by Stanek et al. (2006) using a model for the statistical correspondence between clusters and halos in a WMAP3 cosmology with power spectrum normalization sigma_8 = 0.85. Similar exercises applied to future data sets will allow constraints on the covariance among optical and hot gas properties of clusters at fixed mass.

We explore how the behavior of galaxy cluster scaling relations are affected by flux-limited selection biases and intrinsic covariance among observable properties. Our models presume log-normal covariance between luminosity (L) and temperature (T) at fixed mass (M), centered on evolving, power-law mean relations as a function of host halo mass. Selection can mimic evolution; the \lm and \lt relations from shallow X-ray flux-limited samples will deviate from mass-limited expectations at nearly all scales while the relations from deep surveys ($10^{-14} \cgsflux$) become complete, and therefore unbiased, at masses above $\sims 2 \times 10^{14} \hinv \msol$. We derive expressions for low-order moments of the luminosity distribution at fixed temperature, and show that the slope and scatter of the \lt relation observed in flux-limited samples is sensitive to the assumed \lt correlation coefficient. In addition, \lt covariance affects the redshift behavior of halo counts and mean luminosity in a manner that is nearly degenerate with intrinsic population evolution.

We investigate the relationship between soft \xray luminosity and mass for low redshift clusters of galaxies by comparing observed number counts to expectations of $\Lambda$CDM cosmologies. We use a three-parameter model for the conditional probability of luminosity given mass and epoch, described as a log-normal distribution of fixed width centered on a power-law scaling relation, $L \spropto M^p\rhoc^s(z)$. We use an ensemble of simulated clusters to argue that the observed, intrinsic variance in the temperature--luminosity relation is directly indicative of mass--luminosity variance, and derive $\sigm \se 0.43 \pm 0.06$ from HIFLUGCS data. Adding this to the likelihood analysis results in best-fit estimates $p \se 1.59 \pm 0.05$, $\lnlf \se 1.34 \pm 0.09$, and $\sigm \se 0.37 \pm 0.05$ for self-similar redshift evolution in a concordance ($\Omega_m \se 0.3$, $\Omega_\Lambda \se 0.7$, $\sigma_8 \se0.9$) universe. We show that the present-epoch intercept is very sensitive to power spectrum normalization, $\lnlf \spropto \sigate^{-4}$, and the slope is weakly sensitive to the matter density, $p \spropto \Omega_m^{1/2}$. The intercept derived here is dimmer by a factor 2, and slope slightly steeper, than the L-M relation published using hydrostatic mass estimates of the HIFLUGCS sample. We show that this discrepancy is largely due to Malmquist bias of the \xray\ flux-limited sample. In light of new WMAP constraints, we discuss the interplay between parameters and sources of systematic error, and offer a compromise model with $\Omega_m \se 0.24$, $\sigma_8 \se 0.85$, and somewhat lower scatter $\sigm \se 0.25$, in which hydrostatic mass estimates remain accurate to $\ssim 15%$. We stress the need for independent calibration of the L-M relation via weak gravitational lensing.

We investigate the stellar disk properties of a sample of 19 nearby spiral galaxies with low inclination and late Hubble type (Scd or later). We combine our high-resolution HST I-band observations with existing ground-based optical images to obtain surface brightness profiles that cover a high dynamic range of galactic radius. Most of these galaxies contain a nuclear star cluster, as discussed in a separate paper. The main goal of the present work is to constrain the properties of stellar bulges at these extremely late Hubble types. We find that the surface brightness profiles of the latest-type spirals are complex, with a wide range in shapes. We have sorted our sample in a sequence, starting with ``pure'' disk galaxies (approximately 30% of the sample). These galaxies have exponential stellar disks that extend inwards to within a few tens of pc from the nucleus, where the light from the nuclear cluster starts to dominate. They appear to be truly bulge-less systems. Progressing along the sequence, the galaxies show increasingly prominent deviations from a simple exponential disk model on kpc scales. Traditionally, such deviations have prompted ``bulge-disk'' decompositions. Indeed, the surface brightness profiles of these galaxies are generally well fit by adding a second (exponential) bulge component. However, we find that most surface brightness profiles can be fit equally well (or better) with a single Sersic-type R^{1/n} profile over the entire radial range of the galaxy, without requiring a separate ``bulge'' component. We warn in a general sense against identification of bulges solely on the basis of single-band surface brightness profiles. (abridged)