We study the properties of cosmological shock waves identified in high-resolution, N-body/hydrodynamic simulations of a Ã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 $ 100 or higher. In contrast, internal shocks have mostly low Mach numbers. For all shocks of M ! 1:5, the mean distance between shock surfaces over the entire computed volume is $4 h À1 Mpc at present, or $1 h À1 Mpc for internal shocks within nonlinear structures. Identified external shocks are more extensive, with their surface area $2 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 2dMd4 contribute about one-half 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 about one-half through the history of the universe. Our result supports scenarios in which the intracluster medium contains energetically significant populations of CRs. Subject headings: large-scale structure of universe -methods: numerical -shock waves

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 among 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 intracluster medium or the symmetry of the cluster. This result offers significant promise for precision cosmology using clusters of galaxies.

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 O VI 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 O VI in density and temperature exhibits two populations: one at T ∼ 10 5.5 K (collisionally ionized, 55% of total O VI) 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, O VI provides an important tool for a WHIM baryon census.

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 M pc/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 M < 5 typically trace mergers and complex flows, while 5 < M < 20 and M > 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.

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.Comment: 17 emulateapj pages, 17 figures, replaced with version accepted to Ap

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 > 5 × 10 7 K, ρ ∼ 10 −28 − 10 −27 g/cm 3 , with magnetic field strengths of 0.1 − 1.0µ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 use Enzo, a hybrid Eulerian AMR/N-body code including nongravitational 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.

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 3 h −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 ′ ) 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 > 5 ×10 13 M ⊙ 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.