Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new 400 Mpc/h DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.

We present the nonlinear 2D galaxy power spectrum, $P(k,\mu)$, in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual $\mu$ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the $\mu<0.2$ simulation data, which we show is not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of $\sim 5$% at $k<0.6 h$Mpc$^{-1}$. This use of individual $\mu $ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low $\mu$ simulation data to constrain the nonlinear bias, and $\mu\ge0.2$ to constrain the growth rate and show that $f$ can be constrained to $\sim 26\, (22)$% to a $k_{\rm max}< 0.4\, (0.6) h$Mpc$^{-1}$ from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual $\mu $ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each $\mu$ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy-dark matter connection.

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 Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

X-ray observations were used to examine the complex temperature structure of Abell 85, a cool-core galaxy cluster. Temperature features can provide evidence of merging events which shock heat the intracluster gas. Temperature maps were made from both \textit{Chandra} and \textit{XMM-Newton} obervations. The combination of a new, long-exposure \textit{XMM} observation and an improved temperature map binning technique produced the highest fidelity temperature maps of A85 to date. Hot regions were detected near the subclusters to the South and Southwest in both the \textit{Chandra} and \textit{XMM} temperature maps. The presence of these structures implies A85 is not relaxed. The hot regions may indicate the presence of shocks. The Mach numbers were estimated to be $\sim$1.9 at the locations of the hot spots. Observational effects will tend to systematically reduce temperature jumps, so the measured Mach numbers are likely underestimated. Neither temperature map showed evidence for a shock in the vicinity of the presumed radio relic near the Southwest subcluster. However, the presence of a weak shock cannot be ruled out. There was tension between the temperatures measured by the two instruments.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Simulations of galaxy clusters have a difficult time reproducing the radial gas-property gradients and red central galaxies observed to exist in the cores of galaxy clusters. Thermal conduction has been suggested as a mechanism that can help bring simulations of cluster cores into better alignment with observations by stabilizing the feedback processes that regulate gas cooling, but this idea has not yet been well tested with cosmological numerical simulations. Here we present cosmological simulations of ten galaxy clusters performed with five different levels of isotropic Spitzer conduction, which alters both the cores and outskirts of clusters, but not dramatically. In the cores, conduction flattens central temperature gradients, making them nearly isothermal and slightly lowering the central density but failing to prevent a cooling catastrophe there. Conduction has little effect on temperature gradients outside of cluster cores because outward conductive heat flow tends to inflate the outer parts of the intracluster medium (ICM) instead of raising its temperature. In general, conduction tends reduce temperature inhomogeneity in the ICM, but our simulations indicate that those homogenizing effects would be extremely difficult to observe in ~5 keV clusters. Outside the virial radius, our conduction implementation lowers the gas densities and temperatures because it reduces the Mach numbers of accretion shocks. We conclude that despite the numerous small ways in which conduction alters the structure of galaxy clusters, none of these effects are significant enough to make the efficiency of conduction easily measurable unless its effects are more pronounced in clusters hotter than those we have simulated.

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.

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.

The acceleration of electrons at shock fronts is thought to be responsible for radio relics, extended radio features in the vicinity of merging galaxy clusters. By combining high resolution Adaptive Mesh Refinement Hydro/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. From these simulations, we produce radio, SZE and X-ray images for a large sample of galaxy clusters along with radio luminosity functions and scaling relationships. We find that with upcoming radio arrays, 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 our understanding of electron shock acceleration and kinematic structure of galaxy clusters.

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.

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.

Recently, Suzaku has produced temperature and entropy profiles, along with profiles of gas density, gas fraction, and mass, for multiple galaxy clusters out to ~r_200 (~= virial radius). In this paper, we compare these novel X-ray observations with results from N-body + hydrodynamic adaptive mesh refinement cosmological simulations using the Enzo code. There is excellent agreement in the temperature, density, and entropy profiles between a sample of 27 mostly substructure-free massive clusters in the simulated volume and the observed clusters. This supports our previous contention that clusters have "universal" outer temperature profiles. Furthermore, it appears that the simplest adiabatic gas physics used in these Enzo simulations is adequate to model the outer regions of these clusters without other mechanisms (e.g., non-gravitational heating, cooling, magnetic fields, or cosmic rays). However, the outskirts of these clusters are not in hydrostatic equilibrium. There is significant bulk flow and turbulence in the outer intracluster medium created by accretion from filaments. Thus, the gas is not fully supported by thermal pressure. The implications for mass estimation from X-ray data are discussed.

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.

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.