We describe a new 512-CPU Beowulf cluster with Teraflop performance dedicated to problems in computational astrophysics. The cluster incorporates a cubic network topology based on inexpensive commodity 24-port gigabit switches and point to point connections through the second gigabit port on each Linux server. This configuration has network performance competitive with more expensive cluster configurations and is scaleable to much larger systems using other network topologies. Networking represents only about 9% of our total system cost of USD$561K. The standard Top 500 HPL Linpack benchmark rating is 1.202 Teraflops on 512 CPUs so computing costs by this measure are $0.47/Megaflop. We also describe 4 different astrophysical applications using complex parallel algorithms for studying large-scale structure formation, galaxy dynamics, magnetohydrodynamic flows onto blackholes and planet formation currently running on the cluster and achieving high parallel performance. The MHD code achieved a sustained speed of 2.2 teraflops in single precision or 44% of the theoretical peak.

We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the accretion of discrete, stable subclusters. Any ``cooling flow'' present is overwhelmed by the velocity field within the cluster - the bulk flow of gas through the cluster typically has speeds up to about 2,000 km s^-1 and significant rotation is frequently present in the cluster core. The inclusion of consistent initial cosmological conditions for the cluster within its surrounding supercluster environment is crucial when simulating the evolution of cool cores in rich galaxy clusters. This new model for the hierarchical assembly of cool gas naturally explains the high frequency of cool cores in rich galaxy clusters despite the fact that a majority of these clusters show evidence of substructure which is believed to arise from recent merger activity. Furthermore, our simulations generate complex cluster cores in concordance with recent X-ray observations of cool fronts, cool ``bullets'', and filaments in a number of galaxy clusters. Our simulations were computed with a coupled N-body, Eulerian, adaptive mesh refinement, hydrodynamics cosmology code that properly treats the effects of shocks and radiative cooling by the gas. We employ up to seven levels of refinement to attain a peak resolution of 15.6 h^-1 kpc within a volume 256 h-1 Mpc on a side and assume a standard LambdaCDM cosmology.

We investigate the predicted present-day temperature profiles of the hot, X-ray emitting gas in galaxy clusters for two cosmological models - a current best-guess LCDM model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" of clusters are derived from high-resolution (15/h kpc) simulations which make use of a sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully captures the shocks which are essential for correctly modelling cluster temperatures. We show that the temperature structure on Mpc-scales is highly complex and non-isothermal. However, the temperature profiles of the simulated LCDM and SCDM clusters are remarkably similar and drop-off as $T +AFw-propto (1+-r/a_x)^{-+AFw-delta}$ where $a_x +AFw-sim r_{vir}/1.5$ and $+AFw-delta +AFw-sim 1.6$. This decrease is in good agreement with the observational results of Markevitch et al.(1998) but diverges, primarily in the innermost regions, from their fit which assumes a polytropic equation of state. Our result is also in good agreement with a recent sample of clusters observed by BeppoSAX though there is some indication of missing physics at small radii ($r<0.2 r_{vir}$). We discuss the interpretation of our results and make predictions for new x-ray observations that will extend to larger radii than previously possible. Finally, we show that, for $r>0.2 r_{vir}$, our universal temperature profile is consistent with our most recent simulations which include both radiative cooling and supernovae feedback.

In modern hierarchical theories of structure formation, rich clusters of galaxies form at the vertices of a weblike distribution of matter, with filaments emanating from them to large distances and with smaller objects forming and draining in along these filaments. The amount of mass contained in structure near the cluster can be comparable to the collapsed mass of the cluster itself. As the lensing kernel is quite broad along the line of sight around cluster lenses with typical redshifts near z=0.5, structures many Mpc away from the cluster are essentially at the same location as the cluster itself, when considering their effect on the cluster's weak lensing signal. We use large-scale numerical simulations of structure formation in a Lambda-dominated cold dark matter model to quantify the effect that large-scale structure near clusters has upon the cluster masses deduced from weak lensing analysis. A correction for the scatter in possible observed lensing masses should be included when interpreting mass functions from weak lensing surveys.

We report the results of recent numerical simulations of the head-on merger of a cooling flow cluster with an infalling subcluster of galaxies. These simulations examined the effects of different types of cluster mergers (with 16:1 and 4:1 mass ratios) on the evolution of cluster cooling flows (mass accretion rates of 100 and 400 M/year. The 2-dimensional simulations were performed with a Hydro/N-body code on a uniform grid with a resolution of 20 kpc. We find that the ram-pressure of the infalling gas is crucial in determining the fate of the cooling flow as disruption occurs when a substantial amount of subcluster gas reaches the primary's core. In such cases, the subcluster gas can increase the central cooling time by displacing the high-density cooling gas and by heating it via shocks and turbulent gas motions. Moreover, the fate of a merging cooling flow is also dependent on its initial cooling time. In cases where the initial cooling time is very small then, even if the flow is disrupted, the central cooling time will remain less than a Hubble time and the flow will likely re-establish itself. This has an important observational consequence because such clusters will be classified as cooling flows even though they have experienced a significant merger. In addition, there is a time delay (1-2 Gyrs) between core-crossing and the point at which the central cooling time of a disrupted flow becomes of order a Hubble time (abridged).

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback.

We present high resolution ROSAT x-ray and radio observations of three cooling flow clusters containing steep spectrum radio sources at their cores. All three systems exhibit strong signs of interaction between the radio plasma and the hot intracluster medium. Two clusters, A133 and A2626, show enhanced x-ray emission spatially coincident with the radio source whereas the third cluster, A2052, exhibits a large region of x-ray excess surrounding much of the radio source. Using 3-D numerical simulations, we show that a perturbed jet propagating through a cooling flow atmosphere can give rise to amorphous radio morphologies, particularly in the case where the jet was ``turned off'' and allowed to age passively. In addition, the simulated x-ray surface brightness produced both excesses and deficits as seen observationally.

New adaptive mesh refinement N-body + hydrodynamics numerical simulations are used to illustrate the complex and changing cluster environments in which many radio galaxies live and evolve. Groups and clusters of galaxies form at the intersections of filaments where they continue to accrete gas and dark matter to the present day. The accretion process produces shocks, turbulence, and transonic bulk flows forming a kind of stormy weather within the intracluster medium (ICM). Radio sources embedded within the stormy ICM form distorted, complex morphologies as observed in recent VLA cluster surveys. We show that the bending of wide-angle tailed radio sources can be understood as the result of recent cluster-subcluster mergers. We use new MHD simulations to illustrate how cluster radio halos can be formed by the shocks and turbulence produced during cluster mergers. Finally, we discuss new observations of distant Abell clusters that reveal a class of weak radio sources, probably starbursts, likely produced during the formation of the clusters as they accrete material from the supercluster environment.

A catalog of 732 optically selected, nearby poor clusters of galaxies covering the entire sky north of $\rm -3^{\circ}$ declination is presented. The poor clusters, called WBL clusters, were identified as concentrations of 3 or more galaxies with photographic magnitudes brighter than 15.7, possessing a galaxy surface overdensity of $10^{4/3}$. These criteria are consistent with those used in the identification of the original Yerkes poor clusters, and this new catalog substantially increases the sample size of such objects. These poor clusters cover the entire range of galaxy associations up to and including Abell clusters, systematically including poor and rich galaxy systems spanning over three orders of magnitude in the cluster mass function. As a result, this new catalog contains a greater diversity of richness and structures than other group catalogs, such as the Hickson or Yerkes catalogs. The information on individual galaxies includes redshifts and cross-references to other galaxy catalogs. The entries for the clusters include redshift (where available) and cross-references to other group and cluster catalogs.

With the availability of large-scale redshift survey data, it is now becoming possible to explore correlations between large-scale structure and the properties and morphologies of galaxy clusters. We investigate the spatial distributions of a 98% complete, volume-limited sample of nearby (z<0.1) Abell clusters with well-determined redshifts and find that cooling flow clusters with high mass-accretion rates have nearest neighbors which are much closer than those of other clusters in the sample (at the 99.8% confidence level), and reside in more crowded environments out to 30/h Mpc. Several possible explanations of this effect are discussed.

Hierarchical theories of structure formation predict that clusters of galaxies should be embedded in a web like structure, with filaments emanating from them to large distances. The amount of mass contained within such filaments near a cluster can be comparable to the collapsed mass of the cluster itself. Diffuse infalling material also contains a large amount of mass. Both these components can contribute to the cluster weak lensing signal. This ``projection bias'' is maximized if a filament lies close to the line-of-sight to a cluster. Using large--scale numerical simulations of structure formation in a cosmological constant dominated cold dark matter model, we show that the projected mass typically exceeds the actual mass by several tens of percent. This effect is significant for attempts to estimate cluster masses through weak lensing observations, and will affect weak lensing surveys aimed at constructing the cluster mass function.

In this letter, we present a new determination of the local (z<0.09) X-ray luminosity function (XLF) using a large, statistical sample of 294 Abell clusters and the ROSAT All-Sky-Survey. Given our large sample size, we have reduced errors by a factor of two for L(X)(0.5-2keV)>10^43 ergs/sec. We combine our data with previous work in order to explore possible constraints imposed by the shape of the XLF on cosmological models. A set of currently viable cosmologies is used to construct theoretical XLFs assuming Lx is proportional to M^p and a sigma_8-Omega_0 constraint (from Viana & Liddle 1996) based on the local X-ray temperature function. We fit these models to our observed XLF and verify that the simplest adiabatic, analytic scaling relation (e.g. Kaiser 1986) disagrees strongly with observations. If we assume that clusters can be described by the pre-heated, constant core-entropy models of Evrard & Henry (1991) then the observed XLF is consistent only with 0.1 < Omega_0 < 0.4 if the energy per unit mass in galaxies is roughly equal to the gas energy (ie if beta=1). (abridged)

The radio continuum properties of galaxy clusters with cooling flows are reviewed along with the relationship to the X-ray environment. We find that 60-70% of cD galaxies in cooling flows are radio-loud, a much higher fraction than the 14% found for typical cluster ellipticals. New ROSAT HRI observations reveal a variety of interesting correlations and anticorrelations between the X-ray structure within the inner cooling flows and the radio morphologies. It appears that the radio plasma can have a strong effect on the inner structure of cluster cooling cores as virtually all cooling flow clusters with central radio sources have non-symmetric X-ray structure. Numerical simulations of radio jets in cooling flow atmospheres are presented. We also discuss the prospects for destroying cooling flows via cluster-cluster mergers, using new hydro/N-body simulations which incorporate radiative cooling.

We present an observational analysis of numerical simulations of galaxy cluster mergers. We identify several observational signatures of recent merger activity, and quantitatively assess the uncertainty introduced into cluster mass estimates when invoking the commonly held assumptions of hydrostatic equilibrium, virial equilibrium, spherical symmetry and isothermality. We find that mergers result in multiple X-ray peaks, long-lived elongation of the X-ray emission as well as isophotal twisting and centroid shifting to a degree consistent with recent observations. We also find an enlargement of the X-ray core relative to the dark matter core. Mergers result in non-isothermal clusters exhibiting observable inhomogeneities in the emission-weighted X-ray temperature of several keV on linear scales of less than 0.5 Mpc. The resulting gas dynamics are extremely complex, and we present an example of what might be observed by a high resolution X-ray spectrograph. We further speculate that the gas dynamics, via shocks, bulk flows and turbulence, play an important role in the evolution of cluster galaxies and associated radio sources, particularly wide-angle tailed (WAT) sources and radio halos. We find that X-ray based by cluster mass estimates made under equilibrium assumptions can be uncertain 50\% or more in the first 2 Gyrs after a merger and by up to 25\% after 2 Gyrs depending on the details of the analysis and projection effects. Uncertainties can be considerably larger if the temperature is not well constrained.

From the poor cluster catalog of White et al. (1996), we define a sample of 71 optically-selected poor galaxy clusters. The surface-density enhance- ment we require for our clusters falls between that of the loose associations of Turner and Gott (1976) and the Hickson compact groups (Hickson, 1982). We review the selection biases and determine the statistical comleteness of the sample. For this sample, we report new velocity measurements made with the ARC 3.5-m Dual-Imaging spectrograph and the 2.3-m Steward Observatory MX fiber spectrograph. Combining our own measurements with those from the literature, we examine the velocity distributions, velocity dispersions, and 1-d velocity substructure for our poor cluster sample, and compare our results to other poor cluster samples. We find that approximately half of the sample may have significant 1-d velocity substructure. The optical morphology, large-scale environment, and velocity field of many of these clusters is indicative of young, dynamically evolving systems. In future papers, we will use this sample to derive the poor cluster X-ray luminosity function and gas mass function (see astro-ph/9606120), and will examine the optical/X-ray properties of the clusters in more detail.

We present the first X-ray Luminosity Function (XLF) for an optically-selected sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with z < 0.15. We have extended the measured cluster XLF by more than a factor of 10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog of groups with 0.01<z<0.03 composed of Zwicky galaxies. The X-ray emission was measured from the ROSAT all-sky survey. About 45% of the poor clusters were detected with 0.5-2.0 keV luminosities from (1.7-65)x 10^{41} h^{-2} ergs/sec. These are among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. For this sample, the poor cluster XLF was found to be smooth extrapolation of the rich cluster XLF. A new Hydro/N-body simulation of a Hot +Cold dark matter model with \Omega_{total}=1, \Omega_{\nu}=0.2, and a baryon fraction of 7.5% was used to model and understand our observational selection effects. We found that the observed cluster Gas Mass Function was consistent with our model.

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a Cold + Hot Dark Matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(hot) = 0.3, Omega(cold) = 0.6 and Omega(baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.