Future high spectral resolution telescopes will enable us to place tight constraints on turbulence in the intra-cluster medium through the line widths of strong emission lines. At the same time, these bright lines are the most prone to be optically thick. This requires us to separate the effects of resonant scattering from turbulence, both of which could broaden the lines. How this can be achieved has yet not been quantitatively addressed. In this paper, we propose a flexible new parametrization for the line profile, which allows these effects to be distinguished. The model has only 3 free parameters, which we calibrate with Monte-Carlo radiative transfer simulations. We provide fitting functions and tables that allow the results of these calculations to be easily incorporated into a fast spectral fitting package. In a mock spectral fit, we explicitly show that this parameterization allows us to correctly estimate the turbulent amplitude and metallicity of a cluster such as Perseus, which would otherwise give significantly biased results. We also show how the physical origin of the line shape can be understood analytically.

Upcoming high spectral resolution telescopes, particularly Astro-H, are expected to finally deliver firm quantitative constraints on turbulence in the intra-cluster medium (ICM). We develop a new spectral analysis technique which exploits not just the line width but the entire line shape, and show how the excellent spectral resolution of Astro-H can overcome its relatively poor spatial resolution in making detailed infer- ences about the velocity field. The spectrum is decomposed into distinct components, which can be quantitatively analyzed using Gaussian mixture models. For instance, bulk flows and sloshing produce components with offset means, while partial volume- filling turbulence from AGN or galaxy stirring leads to components with different widths. The offset between components allows us to measure gas bulk motions and separate them from small-scale turbulence, while component fractions and widths con- strain the emission weighted volume and turbulent energy density in each component. We apply mixture modeling to a series of analytic toy models as well as numerical simu- lations of clusters with cold fronts and AGN feedback respectively. From Markov Chain Monte Carlo and Fisher matrix estimates which include line blending and continuum contamination, we show that the mixture parameters can be accurately constrained with Astro-H spectra: at a \sim 10% level when components differ significantly in width, and a \sim 1% level when they differ significantly in mean value. We also study error scalings and use information criteria to determine when a mixture model is preferred. Mixture modeling of spectra is a powerful technique which is potentially applicable to other astrophysical scenarios.

The power spectrum of cosmic infrared background (CIB) anisotropies is sensitive to the connection between star formation and dark matter halos over the entire cosmic star formation history. Here we develop a model that associates star-forming galaxies with dark matter halos and their subhalos. The model is based on a parameterized relation between the dust-processed infrared luminosity and (sub)halo mass. By adjusting 3 free parameters, we attempt to simultaneously fit the 4 frequency bands of the Planck measurement of the CIB anisotropy power spectrum. To fit the data, we find that the star-formation efficiency must peak on a halo mass scale of ~ 5x10^12 solar mass and the infrared luminosity per unit mass must increase rapidly with redshift. By comparing our predictions with a well-calibrated phenomenological model for shot noise, and with a direct observation of source counts, we show that the mean duty cycle of the underlying infrared sources must be near unity, indicating that the CIB is dominated by long-lived quiescent star formation, rather than intermittent short "star bursts". Despite the improved flexibility of our model, the best simultaneous fit to all four Planck channels remains relatively poor. We discuss possible further extensions to alleviate the remaining tension with the data. Our model presents a theoretical framework for a future joint analysis of both background anisotropy and source count measurements.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

The gravitational waves (GWs) emitted by inspiraling binary black holes, expected to be detected by the Laser Interferometer Space Antenna (LISA), could be used to determine the luminosity distance to these sources with the unprecedented precision of <~ 1%. We study cosmological parameter constraints from such standard sirens, in the presence of gravitational lensing by large-scale structure. Lensing introduces magnification with a probability distribution function (PDF) whose shape is highly skewed and depends on cosmological parameters. We use Monte-Carlo simulations to generate mock samples of standard sirens, including a small intrinsic scatter, as well as the additional, larger scatter from lensing, in their inferred distances. We derive constraints on cosmological parameters, by simultaneously fitting the mean and the distribution of the residuals on the distance vs redshift (d_L - z) Hubble diagram. We find that for standard sirens at redshift z ~ 1, the sensitivity to a single cosmological parameter, such as the matter density Omega_m, or the dark energy equation of state w, is ~ 50%-80% tighter when the skewed lensing PDF is used, compared to the sensitivity derived from a Gaussian PDF with the same variance. When these two parameters are constrained simultaneously, the skewness yields a further enhanced improvement (by ~ 120%), owing to the correlation between the parameters. The sensitivity to the amplitude of the matter power spectrum, sigma_8 from the cosmological dependence of the PDF alone, however, is ~ 20% worse than that from the Gaussian PDF. At higher redshifts, the PDF resembles a Gaussian more closely, and the effects of the skewness become less prominent. These results highlight the importance of obtaining an accurate and reliable PDF of the lensing convergence, in order to realize the full potential of standard sirens as cosmological probes.

Scaling relations among galaxy cluster observables, which will become available in large future samples of galaxy clusters, could be used to constrain not only cluster structure, but also cosmology. We study the utility of this approach, employing a physically motivated parametric model to describe cluster structure, and applying it to the expected relation between the Sunyaev-Zel'dovich decrement (S_\nu) and the emission-weighted X-ray temperature (T_ew). The slope and normalization of the entropy profile, the concentration of the dark matter potential, the pressure at the virial radius, and the level of non-thermal pressure support, as well as the mass and redshift-dependence of these quantities are described by free parameters. With a suitable choice of fiducial parameter values, the cluster model satisfies several existing observational constraints. We employ a Fisher matrix approach to estimate the joint errors on cosmological and cluster structure parameters from a measurement of S_\nu vs. T_ew in a future survey. We find that different cosmological parameters affect the scaling relation differently: predominantly through the baryon fraction (\Omega_m and \Omega_b), the virial overdensity (w_0 and w_a for low-z clusters) and the angular diameter distance (w_0, w_a for high-z clusters; \Omega_DE and h). We find that the cosmology constraints from the scaling relation are comparable to those expected from the number counts (dN/dz) of the same clusters. The scaling relation approach is relatively insensitive to selection effects and it offers a valuable consistency check; combining the information from the scaling relation and dN/dz is also useful to break parameter degeneracies and help disentangle cluster physics from cosmology.

In the absence of H_2 molecules, the primordial gas in early dark matter halos with virial temperatures just above T_vir >~ 10^4 K cools by collisional excitation of atomic H. Although it cools efficiently, this gas remains relatively hot, at a temperature near T ~ 8000 K, and consequently might be able to avoid fragmentation and collapse directly into a supermassive black hole (SMBH). In order for H_2--formation and cooling to be strongly suppressed, the gas must be irradiated by a sufficiently intense ultraviolet (UV) flux. We performed a suite of three--dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of gas collapse in three different protogalactic halos with T_vir >~ 10^4 K, irradiated by a UV flux with various intensities and spectra. We determined the critical specific intensity, Jcrit, required to suppress H_2 cooling in each of the three halos. For a hard spectrum representative of metal--free stars, we find (in units of 10^{-21} erg s^{-1} Hz^{-1} sr^{-1} cm^{-2}) 10^4<Jcrit<10^5, while for a softer spectrum, which is characteristic of a normal stellar population, and for which H^{-} --dissociation is important, we find 30<Jcrit<300. These values are a factor of 3--10 lower than previous estimates. We attribute the difference to the higher, more accurate H_2 collisional dissociation rate we adopted. The reduction in Jcrit exponentially increases the number of rare halos exposed to super--critical radiation. When H_2 cooling is suppressed, gas collapse starts with a delay, but it ultimately proceeds more rapidly. The infall velocity is near the increased sound speed, and an object as massive as M ~ 10^5 solar mass may form at the center of these halos, compared to the M ~ 10^2 solar mass stars forming when H_2--cooling is efficient.

We present the first results of a survey of 14 low redshift galaxy clusters using Suzaku. Although luminous (L_x> 1\times 10^{43} erg s^{-1} (0.1-2.4 keV)), these clusters have no prior pointed X-ray data. Together with 47 other systems they form a flux limited sample (f_x >1.0 \times 10^{-11} erg s^{-1}cm^{-2} (0.1-2.4 keV)) with z<0.1 in the northern celestial hemisphere. Using this total sample we evaluate the local $L-T$ relationship and the local cluster temperature function. Suzaku temperature measurements appear to be in accord with those of other missions. General agreement is found with other published estimates of the low redshift cluster temperature function, however the sample used here exhibits slightly lower space densities at gas temperatures below 4-5 keV. We find a corresponding deficit in the number of clusters with temperatures between approximately 4.5 and 5.5 keV. Although at low significance, a similar feature exists in previous low redshift cluster datasets. We suggest that low-redshift cluster samples, while crucial for calibrating precision cosmology measurements, must be used with caution due to their limited size and susceptibility to the effects of cosmic variance.

Energetic feedback processes during the formation of galaxy clusters may have heated and ionized a large fraction of the intergalactic gas in proto-cluster regions. When such a highly ionized hot ``super-bubble'' falls along the sightline to a background quasar, it would be seen as a large void, with little or no absorption, in the Lyman alpha forest. We examine the spectra of 137 quasars in the Sloan Digital Sky Survey, to search for such voids, and find no clear evidence of their existence. The size distribution of voids in the range 5-70 Angstrom (corresponding to physical sizes of approximately 3-35 comoving Mpc/h) is consistent with the standard model for the Lyman alpha forest without additional hot bubbles. We adapt a physical model for HII bubble growth during cosmological reionization (Furlanetto, Zaldarriaga and Hernquist 2004), to describe the expected size-distribution of hot super-bubbles at redshift around z = 3. This model incorporates the conjoining of bubbles around individual neighboring galaxies. Using the non-detection of voids, we find that models in which the volume filling factor of hot bubbles exceeds approximately 20 percent at z=3 can be ruled out, primarily because they overproduce the number of large (40-50 Angstrom) voids. We conclude that any pre-heating mechanism that explains galaxy cluster observations must avoid heating the low-density gas in the proto-cluster regions, either by operating relatively recently (z<3) or by depositing entropy in the high-density regions.