The advent of sensitive low frequency radio observations has revealed a number of diffuse radio objects with peculiar properties that are challenging our understanding about the physics of the intracluster medium. Here, we report the discovery of a steep spectrum radio halo surrounding the central Brightest Cluster Galaxy (BCG) in the galaxy cluster SPT-CL J2031-4037. This cluster is morphologically disturbed yet has a weak cool core, an example of cool core/non-cool core transition system, which harbours a radio halo of $\sim 0.7$ Mpc in size. The halo emission detected at 1.7 GHz is less extended compared to that in the 325 MHz observation, and the spectral index of the part of the halo visible at 325 MHz to 1.7 GHz frequencies was found to be $-1.35 \pm 0.07$. Also, $P_{1.4\ \mathrm{GHz}}$ was found to be $0.77 \times 10^{24}$ W Hz$^{-1}$ which falls in the region where radio mini-halos, halo upper limits and ultra-steep spectrum (USS) halos are found in the $P_{1.4\ \mathrm{GHz}} - L_\mathrm{X}$ plane. Additionally, simulations presented in the paper provide support to the scenario of the steep spectrum. The diffuse radio emission found in this cluster may be a steep spectrum "intermediate" or "hybrid" radio halo which is transitioning into a mini-halo.

Cool core galaxy clusters are considered to be dynamically relaxed clusters with regular morphology and highly X-ray luminous central region. However, cool core clusters can also be sites for merging events that exhibit cold fronts in X-ray and mini-halos in radio. We present recent radio/X-ray observations of the Phoenix Cluster or SPT-CL J2344-4243 at the redshift of $z=0.596$. Using archival {\it Chandra} X-ray observations, we detect spiraling cool gas around the cluster core as well as discover two cold fronts near the core. It is perhaps the most distant galaxy cluster to date known to host cold fronts. Also, we present JVLA\footnote{Jansky Very Large Array\\ \url{https://science.nrao.edu/facilities/vla}} 1.52 GHz observations of the minihalo, previously discovered at 610 MHz with GMRT\footnote{Giant Metrewave Radio Telescope\\ \url{http://www.gmrt.ncra.tifr.res.in}} observations in the center of the Phoenix galaxy cluster. The minihalo flux density at 1.52 GHz is $9.65 \pm 0.97$ mJy with the spectral index between 610 MHz and 1.52 GHz being $-0.98 \pm 0.16$\footnote{$S_{\nu} = \nu^{\alpha}$ where $S_{\nu}$}. A possible origin of these radio sources is turbulence induced by sloshing of the gas in the cluster core.

ClusterPyXT is a new software pipeline to generate spectral temperature, X-ray surface brightness, pressure, and density maps from X-ray observations of galaxy clusters. These data products help elucidate the physics of processes occurring within clusters of galaxies, including turbulence, shock fronts, nonthermal phenomena, and the overall dynamics of cluster mergers. ClusterPyXT automates the creation of these data products with minimal user interaction, and allows for rapid analyses of archival data with user defined parameters and the ability to straightforwardly incorporate additional observations. In this paper, we describe in detail the use of this code and release it as an open source Python project on GitHub.

We present results from an X-ray and radio study of the merging galaxy cluster Abell 115. We use the full set of 5 Chandra observations taken of A115 to date (360 ks total integration) to construct high-fidelity temperature and surface brightness maps. We also examine radio data from the Very Large Array at 1.5 GHz and the Giant Metrewave Radio Telescope at 0.6 GHz. We propose that the high X-ray spectral temperature between the subclusters results from the interaction of the bow shocks driven into the intracluster medium by the motion of the subclusters relative to one another. We have identified morphologically similar scenarios in Enzo numerical N-body/hydrodynamic simulations of galaxy clusters in a cosmological context. In addition, the giant radio relic feature in A115, with an arc-like structure and a relatively flat spectral index, is likely consistent with other shock-associated giant radio relics seen in other massive galaxy clusters. We suggest a dynamical scenario that is consistent with the structure of the X-ray gas, the hot region between the clusters, and the radio relic feature.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

X-ray astronomy is an important tool in the astrophysicist's toolkit to investigate high-energy astrophysical phenomena. Theoretical numerical simulations of astrophysical sources are fully three-dimensional representations of physical quantities such as density, temperature, and pressure, whereas astronomical observations are two-dimensional projections of the emission generated via mechanisms dependent on these quantities. To bridge the gap between simulations and observations, algorithms for generating synthetic observations of simulated data have been developed. We present an implementation of such an algorithm in the yt analysis software package. We describe the underlying model for generating the X-ray photons, the important role that yt and other Python packages play in its implementation, and present a detailed workable example of the creation of simulated X-ray observations.

The galaxy cluster Abell 3667 is an ideal laboratory to study the plasma processes in the intracluster medium (ICM). High resolution Chandra X-ray observations show a cold front in Abell 3667. At radio wavelengths, Abell 3667 reveals a double radio-relic feature in the outskirts of the cluster. These suggest multiple merger events in this cluster. In this paper, we analyze the substantial archival X-ray observations of Abell 3667 from ChandraX-ray Observatory and compare these with existing radio observations as well as state-of-the-art AMR (Adaptive Mesh Refinement) MHD cosmological simulations using Enzo. We have used two temperature map making techniques, Weighted Voronoi Tessellation and Adaptive Circular Binning, to produce the high resolution and largest field-of-view temperature maps of Abell 3667. These high fidelity temperature maps allow us to study the X-ray shocks in the cluster using a new 2-dimensional shock-finding algorithm. We have also estimated the Mach numbers from the shocks inferred from previous ATCA radio observations. The combined shock statistics from the X-ray and radio data are in agreement with the shock statistics in a simulated MHD cluster. We have also studied the profiles of the thermodynamic properties across the cold front using 447 ksec from the combined Chandra observations on Abell 3667. Our results show that the stability of the cold front in Abell 3667 can be attributed to the suppression of the thermal conduction across the cold front by a factor of 100 - 700 compared to the classical Spitzer value.

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

Flux limited X-ray surveys of galaxy clusters show that clusters come in two roughly equally proportioned varieties: "cool core" clusters (CCs) and non-"cool core" clusters (NCCs). In previous work, we have demonstrated using cosmological $N$-body + Eulerian hydrodynamic simulations that NCCs are often consistent with early major mergers events that destroy embryonic CCs. In this paper we extend those results and conduct a series of simulationsusing different methods of gas cooling, and of energy and metal feedback from supernovae, where we attempt to produce a population of clusters with realistic central cooling times, entropies, and temperatures. We find that the use of metallicity-dependent gas cooling is essential to prevent early overcooling,and that adjusting the amount of energy and metal feedback can have a significant impact on observable X-ray quantities of the gas. We are able to produce clusters with more realistic central observable quantities than have previously been attained. However, there are still significant discrepancies between the simulated clusters and observations, which indicates that a different approach to simulating galaxies in clusters is needed. We conclude by looking towards a promising subgrid method of modeling galaxy feedback in clusters which may help to ameliorate the discrepancies between simulations and observations.

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu G and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

We track the histories of massive clusters of galaxies formed within a cosmological hydrodynamic simulation. Specifically, we track the time evolution of the energy in random bulk motions of the intracluster medium and X-ray measures of cluster structure and their relationship to cluster mergers. We aim to assess the viability of the turbulent re-acceleration model for the generation of giant radio halos by comparing the level of turbulent kinetic energy in simulated clusters with the observed properties of radio halo clusters, giving particular attention to the association of radio halos to clusters with disturbedX-ray structures. The evolution of X-ray cluster structure and turbulence kinetic energy, k, in simulations can then inform us about the expected lifetime of radio halos and the fraction of clusters as a function of redshift expected to host them. We find strong statistical correlation of disturbed structure measures and the presence of enhancements in k. Specifically, quantitatively "disturbed", radio halo-like X-ray morphology in our sample indicates a 92% chance of the cluster in question having k elevated to more than twice its minimum value over the cluster's life. The typical lifetime of episodes of elevated turbulence is on the order of 1 Gyr, though these periods can last 5 Gyrs or more. This variation reflects the wide range of cluster histories; while some clusters undergo complex and repeated mergers spending a majority of their time in elevated k states, other clusters are relaxed over nearly their entire history. We do not find a bimodal relationship between cluster X-ray luminosity and the total energy in turbulence that might account directly for a bimodal L_X-P_{1.4 GHz} relation. However, our result may be consistent with the observed bimodality, as here we are not including a full treatment of cosmic rays sources and magnetic fields.

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees.

Cosmological shocks are a critical part of large-scale structure formation, and are responsible for heating the intracluster medium in galaxy clusters. In addition, they are also capable of accelerating non-thermal electrons and protons. In this work, we focus on the acceleration of electrons at shock fronts, which is thought to be responsible for radio relics - extended radio features in the vicinity of merging galaxy clusters. By combining high resolution AMR/N-body cosmological simulations with an accurate shock finding algorithm and a model for electron acceleration, we calculate the expected synchrotron emission resulting from cosmological structure formation. We produce synthetic radio maps of a large sample of galaxy clusters and present luminosity functions and scaling relationships. With upcoming long wavelength radio telescopes, we expect to see an abundance of radio emission associated with merger shocks in the intracluster medium. By producing observationally motivated statistics, we provide predictions that can be compared with observations to further improve our understanding of magnetic fields and electron shock acceleration.

We perform a series of cosmological simulations using Enzo, an Eulerian adaptive-mesh refinement, N-body + hydrodynamical code, applied to study the warm/hot intergalactic medium. The WHIM may be an important component of the baryons missing observationally at low redshift. We investigate the dependence of the global star formation rate and mass fraction in various baryonic phases on spatial resolution and methods of incorporating stellar feedback. Although both resolution and feedback significantly affect the total mass in the WHIM, all of our simulations find that the WHIM fraction peaks at z ~ 0.5, declining to 35-40% at z = 0. We construct samples of synthetic OVI absorption lines from our highest-resolution simulations, using several models of oxygen ionization balance. Models that include both collisional ionization and photoionization provide excellent fits to the observed number density of absorbers per unit redshift over the full range of column densities (10^13 cm-2 <= N_OVI <= 10^15 cm^-2). Models that include only collisional ionization provide better fits for high column density absorbers (N_OVI > 10^14 cm^-2). The distribution of OVI in density and temperature exhibits two populations: one at T ~ 10^5.5 K (collisionally ionized, 55% of total OVI) and one at T ~ 10^4.5 K (photoionized, 37%) with the remainder located in dense gas near galaxies. While not a perfect tracer of hot gas, OVI provides an important tool for a WHIM baryon census.

We examine the incidence of cold fronts in a large sample of galaxy clusters extracted from a (512h^-1 Mpc) hydrodynamic/N-body cosmological simulation with adiabatic gas physics computed with the Enzo adaptive mesh refinement code. This simulation contains a sample of roughly 4000 galaxy clusters with M > 10^14 M_sun at z=0. For each simulated galaxy cluster, we have created mock 0.3-8.0 keV X-ray observations and spectroscopic-like temperature maps. We have searched these maps with a new automated algorithm to identify the presence of cold fronts in projection. Using a threshold of a minimum of 10 cold front pixels in our images, corresponding to a total comoving length L_cf > 156h^-1 kpc, we find that roughly 10-12% of all projections in a mass-limited sample would be classified as cold front clusters. Interestingly, the fraction of clusters with extended cold front features in our synthetic maps of a mass-limited sample trends only weakly with redshift out to z=1.0. However, when using different selection functions, including a simulated flux limit, the trending with redshift changes significantly. The likelihood of finding cold fronts in the simulated clusters in our sample is a strong function of cluster mass. In clusters with M>7.5x10^14 M_sun the cold front fraction is 40-50%. We also show that the presence of cold fronts is strongly correlated with disturbed morphology as measured by quantitative structure measures. Finally, we find that the incidence of cold fronts in the simulated cluster images is strongly dependent on baryonic physics.

Cosmological simulations of galaxy formation often rely on prescriptions for star formation and feedback that depend on halo properties such as halo mass, central over-density, and virial temperature. In this paper we address the convergence of individual halo properties, based on their number of particles N, focusing in particular on the mass of halos near the resolution limit of a simulation. While it has been established that the halo mass function is sampled on average down to N~30 particles, we show that individual halo properties exhibit significant scatter, and some systematic biases, as one approaches the resolution limit. We carry out a series of cosmological simulations using the Gadget2 and Enzo codes with N_p=64^3 to N_p=1024^3 total particles, keeping the same large-scale structure in the simulation box. We consider boxes from l_{box} = 8 Mpc/h to l_{box} = 512 Mpc/h to probe different halo masses and formation redshifts. We cross-identify dark matter halos in boxes at different resolutions and measure the scatter in their properties. The uncertainty in the mass of single halos depends on the number of particles (scaling approximately as N^{-1/3}), but the rarer the density peak, the more robust its identification. The virial radius of halos is very stable and can be measured without bias for halos with N>30. In contrast, the average density within a sphere containing 25% of the total halo mass is severely underestimated (by more than a factor 2) and the halo spin is moderately overestimated for N<100. If sub-grid physics is implemented upon a cosmological simulation, we recommend that rare halos (~3sigma peaks) be resolved with N>100 particles and common halos (~1sigma peaks) with N>400 particles to avoid excessive numerical noise and possible systematic biases in the results.

In our ongoing multi-wavelength study of cluster AGN, we find ~75% of the spectroscopically identified cluster X-ray point sources (XPS) with L(0.3-8.0keV)>10^{42} erg s^{-1} and cluster radio galaxies with P(1.4 GHz) > 3x10^{23} W Hz^{-1} in 11 moderate redshift clusters (0.2<z<0.4) are located within 500 kpc from the cluster center. In addition, these sources are much more centrally concentrated than luminous cluster red sequence (CRS) galaxies. With the exception of one luminous X-ray source, we find that cluster XPSs are hosted by passive red sequence galaxies, have X-ray colors consistent with an AGN power-law spectrum, and have little intrinsic obscuring columns in the X-ray (in agreement with previous studies). Our cluster radio sources have properties similar to FR1s, but are not detected in X-ray probably because their predicted X-ray emission falls below our sensitivity limits. Based on the observational properties of our XPS population, we suggest that the cluster XPSs are low-luminosity BL Lac objects, and thus are beamed low-power FR 1s. Extrapolating the X-ray luminosity function of BL Lacs and the Radio luminosity function of FR 1s down to fainter radio and X-ray limits, we estimate that a large fraction, perhaps all CRSs with L>L* possess relativistic jets which can inject energy into the ICM, potentially solving the uniform heating problem in the central region of clusters.

Using Chandra imaging spectroscopy and VLA L-band maps, we have identified radio galaxies at P(1.4 GHz) >= 3x10^{23} W Hz^{-1} and X-ray point sources (XPSs) at L(0.3-8 keV) >= 10^{42} ergs s^{-1} in 11 moderate redshift (0.2<z<0.4) clusters of galaxies. Each cluster is uniquely chosen to have a total mass similar to predicted progenitors of the present-day Coma Cluster. Within a projected radius of 1 Mpc we detect 20 radio galaxies and 8 XPSs confirmed to be cluster members above these limits. 75% of these are detected within 500 kpc of the cluster center. This result is inconsistent with a random selection from bright, red sequence ellipticals at the > 99.999% level. All but one of the XPSs are hosted by luminous ellipticals which otherwise show no other evidence for AGN activity. These objects are unlikely to be highly obscured AGN since there is no evidence for large amounts of X-ray or optical absorption. The most viable model for these sources are low luminosity BL Lac Objects. The expected numbers of lower luminosity FR 1 radio galaxies and BL Lacs in our sample converge to suggest that very deep radio and X-ray images of rich clusters will detect AGN in a large fraction of bright elliptical galaxies in the inner 500 kpc. Because both the radio galaxies and the XPSs possess relativistic jets, they can inject heat into the ICM. Using the most recent scalings of P_jet ~ L_radio^{0.5} from Birzan et al. (2008), radio sources weaker than our luminosity limit probably contribute the majority of the heat to the ICM. If a majority of ICM heating is due to large numbers of low power radio sources, triggered into activity by the increasing ICM density as they move inward, this may be the feedback mechanism necessary to stabilize cooling in cluster cores.

To better constrain models of cool core galaxy cluster formation, we have used X-ray observations taken from the Chandra and ROSAT archives to examine the properties of cool core and non-cool core clusters, especially beyond the cluster cores. We produced X-ray images, surface brightness profiles, and hardness ratio maps of 30 nearby rich Abell clusters (17 cool cores and 13 non-cool cores). We show that the use of double beta-models with cool core surface brightness profiles and single beta-models for non-cool core profiles yield statistically significant differences in the slopes (i.e., beta values) of the outer surface brightness profiles, but similar cluster core radii, for the two types of clusters. Hardness ratio profiles as well as spectroscopically-fit temperatures suggest that non-cool core clusters are warmer than cool core clusters of comparable mass beyond the cluster cores. We compared the properties of these clusters with the results from analogously reduced simulations of 88 numerical clusters created by the AMR Enzo code. The simulated surface brightness profiles have steeper beta-model fits in the outer cluster regions for both cool cores and non-cool cores, suggesting additional ICM heating is required compared to observed cluster ICMs. Temperature and surface brightness profiles reveal that the simulated clusters are over-cooled in their cores. As in the observations, however, simulated hardness ratio and temperature profiles indicate that non-cool core clusters are warmer than cool core clusters of comparable mass far beyond the cluster cores. The general similarities between observations and simulations support a model described in Burns et al. 2008 suggesting that non-cool core clusters suffered early major mergers destroying nascent cool cores.

Detection of the Warm-Hot Intergalactic Medium (WHIM) using Sunyaev-Zeldovich effect (SZE) surveys is an intriguing possibility, and one that may allow observers to quantify the amount of "missing baryons" in the WHIM phase. We estimate the necessary sensitivity for detecting low density WHIM gas with the South Pole Telescope (SPT) and Planck Surveyor for a synthetic 100 square degree sky survey. This survey is generated from a very large, high dynamic range adaptive mesh refinement cosmological simulation performed with the Enzo code. We find that for a modest increase in the SPT survey sensitivity (a factor of 2-4), the WHIM gas makes a detectable contribution to the integrated sky signal. For a Planck-like satellite, similar detections are possible with a more significant increase in sensitivity (a factor of 8-10). We point out that for the WHIM gas, the kinematic SZE signal can sometimes dominate the thermal SZE where the thermal SZE decrement is maximal (150 GHz), and that using the combination of the two increases the chance of WHIM detection using SZE surveys. However, we find no evidence of unique features in the thermal SZE angular power spectrum that may aid in its detection. Interestingly, there are differences in the power spectrum of the kinematic SZE, which may not allow us to detect the WHIM directly, but could be an important contaminant in cosmological analyses of the kSZE-derived velocity field. Corrections derived from numerical simulations may be necessary to account for this contamination.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 Mpc/h)^3, cosmological volume and study the shock Mach number (M) distribution as a function of pre-shock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that results from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with Mach<5 typically trace mergers and complex flows, while 5<Mach<20 and Mach>20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic ray protons. The acceleration of cosmic ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.

We use Enzo, a hybrid Eulerian AMR/N-body code including non-gravitational heating and cooling, to explore the morphology of the X-ray gas in clusters of galaxies and its evolution in current generation cosmological simulations. We employ and compare two observationally motivated structure measures: power ratios and centroid shift. Overall, the structure of our simulated clusters compares remarkably well to low-redshift observations, although some differences remain that may point to incomplete gas physics. We find no dependence on cluster structure in the mass-observable scaling relations, T_X-M and Y_X-M, when using the true cluster masses. However, estimates of the total mass based on the assumption of hydrostatic equilibrium, as assumed in observational studies, are systematically low. We show that the hydrostatic mass bias strongly correlates with cluster structure and, more weakly, with cluster mass. When the hydrostatic masses are used, the mass-observable scaling relations and gas mass fractions depend significantly on cluster morphology, and the true relations are not recovered even if the most relaxed clusters are used. We show that cluster structure, via the power ratios, can be used to effectively correct the hydrostatic mass estimates and mass-scaling relations, suggesting that we can calibrate for this systematic effect in cosmological studies. Similar to observational studies, we find that cluster structure, particularly centroid shift, evolves with redshift. This evolution is mild but will lead to additional errors at high redshift. Projection along the line of sight leads to significant uncertainty in the structure of individual clusters: less than 50% of clusters which appear relaxed in projection based on our structure measures are truly relaxed.

Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian hydrodynamic simulations, including radiative cooling and heating, are used to address this question as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. These adaptive mesh refinement simulations produce both CC and NCC clusters in the same volume. They have a peak resolution of 15.6 h^{-1} kpc within a (256 h^{-1} Mpc)^3 box. Our simulations suggest that there are important evolutionary differences between CC clusters and their NCC counterparts. Many of the numerical CC clusters accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CC's survived the collisions. By contrast, NCC clusters experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevented CC re-formation. As a result, our simulations predict observationally testable distinctions in the properties of CC and NCC beyond the core regions in clusters. In particular, we find differences between CC versus NCC clusters in the shapes of X-ray surface brightness profiles, between the temperatures and hardness ratios beyond the cores, between the distribution of masses, and between their supercluster environs. It also appears that CC clusters are no closer to hydrostatic equilibrium than NCC clusters, an issue important for precision cosmology measurements.

We present the first results from a new generation of simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have simulated a very large (512^3h^{-3}Mpc^3) volume with unprecedented dynamic range. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish SZE surveys. Smaller beam surveys (~1 arcmin) have more than one massive cluster within a beam diameter 5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at the maximum decrement. We find that there is a contribution from gas outside clusters of ~16% per object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. Finally, we find that in images where objects with M > 5x10^{13} M_{\odot} have had their SZE signatures removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming generation of SZE surveys.

We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse effects resulting from use of X-ray fitted beta-model parameters with Sunyaev-Zeldovich effect (SZE) data. There is a fundamental incompatibility between model fits to X-ray surface brightness profiles and those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X-ray parameters are often used in SZE analysis. We show that this leads to biased estimates of the integrated Compton y-parameter inside r_{500} calculated from clusters. We suggest a simple correction of the method, using a non-isothermal beta-model modified by a universal temperature profile, which brings these calculated quantities into closer agreement with the true values.

We present preliminary results from simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. These simulations are the most advanced calculations of their kind. The simulated sky surveys allow a direct comparison of large N-body/hydro cosmological simulations to current and pending sky surveys. Our synthetic surveys provide an indispensable guide for observers in the interpretation of large area sky surveys, and will develop the tools necessary to discriminate between models for cluster baryonic physics, and to accurately determine cosmological parameters.

Why do some clusters have cool cores while others do not? In this paper, cosmological simulations, including radiative cooling and heating, are used to examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. Numerical CC clusters at z=0 accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CCs survived the collisions. By contrast, NCC clusters of similar mass experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevent CC re-formation. We discuss observational consequences.

Withdrawn at the current time due to errors.

We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numerically simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect (SZE), we reduce the observations to obtain measurements of the cluster gas mass. We are thus able to quantify the possible sources of uncertainty and systematic bias associated with the common simplifying assumptions used in reducing real cluster observations including isothermality and hydrostatic equilibrium. We find that intrinsic variations in clusters limit the precision of observational gas mass estimation to ~10% to 1 sigma confidence excluding instrumental effects. Gas mass estimates performed via all methods surprisingly show little or no trending in the scatter as a function of cluster redshift. For the full cluster sample, methods that use SZE profiles out to roughly the virial radius are the simplest, most accurate, and unbiased way to estimate cluster mass. X-ray methods are systematically more precise mass estimators than are SZE methods if merger and cool core systems are removed, but X-ray methods slightly overestimate (5-10%) the cluster gas mass on average. We find that cool core clusters in our samples are particularly poor candidates for observational mass estimation, even when excluding emission from the core region. The effects of cooling in the cluster gas alter the radial profile of the X-ray and SZE surface brightness outside the cool core region, leading to poor gas mass estimates in cool core clusters. Finally, we find that methods using a universal temperature profile estimate cluster masses to higher precision than those assuming isothermality.

Recent X-ray observations reveal growing evidence for heating by active galactic nuclei (AGN) in clusters and groups of galaxies. AGN outflows play a crucial role in explaining the riddle of cooling flows and the entropy problem in clusters. Here we study the effect of AGN on the intra-cluster medium in a cosmological simulation using the adaptive mesh refinement FLASH code. We pay particular attention to the effects of conductivity and viscosity on the dissipation of weak shocks generated by the AGN activity in a realistic galaxy cluster. Our 3D simulations demonstrate that both viscous and conductive dissipation play an important role in distributing the mechanical energy injected by the AGN, offsetting radiative cooling and injecting entropy to the gas. These processes are important even when the transport coefficients are at a level of 10% of the Spitzer value. Provided that both conductivity and viscosity are suppressed by a comparable amount, conductive dissipation is likely to dominate over viscous dissipation. Nevertheless, viscous effects may still affect the dynamics of the gas and contribute a significant amount of dissipation compared to radiative cooling. We also present synthetic Chandra observations. We show that the simulated buoyant bubbles inflated by the AGN, and weak shocks associated with them, are detectable with the Chandra observatory.

We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four different samples that differ only in their assumed baryonic physics. We show that the strength of the thermal SZE integrated over a significant fraction of the virialized region of the clusters is relatively insensitive to the detailed heating and cooling processes in the cores of clusters by demonstrating that the derived scaling relations are nearly identical between the four cluster samples considered. For our synthetic images, the central Comptonization parameter shows significant boosting during transient merging events, but the integrated SZE appears to be relatively insensitive to these events. Most importantly, the integrated SZE closely tracks the underlying cluster mass. Observations through the thermal SZE allow a strikingly accurate mass estimation from relatively simple measurements that do not require either parametric modeling or geometric deprojection and thus avoid assumptions regarding the physics of the ICM or the symmetry of the cluster. This result offers significant promise for precision cosmology using clusters of galaxies.

We present Chandra observations of the cool cluster A168, for which previous X-ray imaging and optical studies indicated a merger of two subclusters nearly in the plane of the sky. We derive a temperature map for A168, which shows that the merger has proceeded beyond the core passage and is near subcluster turnaround. It also reveals an unusual feature -- the gas core of one of the subclusters forms a tongue-like structure extending ahead (in the direction of motion) of the subcluster center. The coolest cluster gas is found in a crescent-shaped region at the tip of this tongue, and forms a cold front in pressure equilibrium with the external gas. In contrast with this feature's forward location, previously observed merger cold fronts (e.g., A3667, 1E0657--56) lagged behind their host subclusters, as expected in the presense of ram pressure. We propose that A168 illustrates a much later stage in the evolution of a cold front, when its host subcluster approaches the apocenter of the merger orbit where the ram pressure on its gas drops sharply. As a result, a large chunk of the subcluster gas ``slingshots'' past the dark matter center, becomes unbound from the subcluster and expands adiabatically, as seen in some recent hydrodynamic simulations.

We study the properties of cosmological shock waves identified in high-resolution, N-body/hydrodynamic simulations of a $\Lambda$CDM universe and their role on thermalization of gas and acceleration of nonthermal, cosmic ray (CR) particles. External shocks form around sheets, filaments and knots of mass distribution when the gas in void regions accretes onto them. Within those nonlinear structures, internal shocks are produced by infall of previously shocked gas to filaments and knots, and during subclump mergers, as well as by chaotic flow motions. Due to the low temperature of the accreting gas, the Mach number of external shocks is high, extending up to $M\sim 100$ or higher. In contrast, internal shocks have mostly low Mach numbers. For all shocks of $M\ge1.5$ the mean distance between shock surfaces over the entire computed volume is $\sim4 h^{-1}$ Mpc at present, or $\sim 1 h^{-1}$ Mpc for internal shocks within nonlinear structures. Identified external shocks are more extensive, with their surface area $\sim2$ times larger than that of identified internal shocks at present. However, especially because of higher preshock densities, but also due to higher shock speeds, internal shocks dissipate more energy. Hence, the internal shocks are mainly responsible for gas thermalization as well as CR acceleration. In fact, internal shocks with $2 \la M \la 4$ contribute $\sim 1/2$ of the total dissipation. Using a nonlinear diffusive shock acceleration model for CR protons, we estimate the ratio of CR energy to gas thermal energy dissipated at cosmological shock waves to be $\sim1/2$ through the history of the universe. Our result supports scenarios in which the intracluster medium contains energetically significant populations of CRs.

Shocks are a ubiquitous consequence of cosmic structure formation, and they play an essential role in heating galaxy cluster media. Virtually all of the gas in clusters has been processed by one or more shocks of at least moderate strength. These are collisionless shocks, so likely sites for diffusive shock acceleration of high energy particles. We have carried out numerical simulations of cosmic structure formation that directly include acceleration and transport of nonthermal protons, as well as primary and secondary electrons. Nonthermal emissions have also been computed from the resulting particle spatial and energy distributions. Here we outline some of our current findings, showing that nonthermal protons may contribute a significant pressure in cluster media, and that expected radio, X-ray and $\gamma$-ray emissions from these populations should be important cluster diagnostics.