Galaxy clusters, the most massive collapsed structures, have been routinely used to determine cosmological parameters. When using clusters for cosmology, the crucial assumption is that they are relaxed. However, subarcminute resolution Sunyaev-Zel'dovich (SZ) effect images compared with high resolution X-ray images of some clusters show significant offsets between the two peaks. We have carried out self-consistent N-body/hydrodynamical simulations of merging galaxy clusters using FLASH to study these offsets quantitatively. We have found that significant displacements result between the SZ and X-ray peaks for large relative velocities for all masses used in our simulations as long as the impact parameters were about 100-250 kpc. Our results suggest that the SZ peak coincides with the peak in the pressure times the line-of-sight characteristic length and not the pressure maximum (as it would for clusters in equilibrium). The peak in the X-ray emission, as expected, coincides with the density maximum of the main cluster. As a consequence, the morphology of the SZ signal and therefore the offset between the SZ and X-ray peaks change with viewing angle. As an application, we compare the morphologies of our simulated images to observed SZ and X-ray images and mass surface densities derived from weak lensing observations of the merging galaxy cluster CL0152-1357. We find that a large relative velocity of 4800 km/s is necessary to explain these observations. We conclude that an analysis of the morphologies of multi-frequency observations of merging clusters can be used to put meaningful constraints on the initial parameters of the progenitors.

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

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

Asteroids can be eclipsed by other bodies in the Solar System, but no direct observation of an asteroid eclipse has been reported to date. We describe a statistical method to predict an eclipse for an asteroid based on the analysis of the orbital elements covariance matrix. By propagating a set of Virtual Asteroids to an epoch correspondent to a close approach with a Solar System planet or natural satellite, it is possible to estimate the probability of a partial or total eclipse. The direct observation of an eclipse can provide data useful to improve the asteroid orbit, especially for dim asteroids typically observed only for a few days. We propose two different methods: the first, based on the inclusion of the apparent magnitude residuals into the orbit's least squares minimization process, capable of improving the asteroid's nominal orbit and the related covariance matrix; the second, based on weighting different Virtual Asteroids in relation to their apparent magnitude during the eclipse, useful for recovery purposes. As an application, we have numerically investigated the possibility of a Near Earth Asteroid eclipsed by the Moon or the Earth in the 1990-2050 period. A total of 74 distinct eclipses have been found, involving 59 asteroids. In particular, the asteroid (99942) Apophis has a probability of about 74% to enter the Moon's penumbra cone and a probability of about 6% to enter the umbra cone on April 14, 2029, less than six hours after a very close approach to Earth.

We summarize the results of numerical simulations of colliding gas-rich disk galaxies in which the impact velocity is set parallel to the spin axes of the two galaxies. The effects of varying the impact speed are studied with particular attention to the resulting gaseous structures and shockwave patterns, and the time needed to produce these structures. The simulations employ an N-body treatment of the stars and dark matter, together with an SPH treatment of the gas, in which all components of the models are gravitationally active. The results indicate that for such impact geometries, collisions can lead to the very rapid formation of a central, rapidly rotating, dense gas disk, and that in all cases extensive star formation is predicted by the very high gas densities and prevalence of shocks, both in the nucleus and out in the galactic disks. As the dense nucleus is forming, gas and stars are dispersed over very large volumes, and only fall back towards the nucleus over long times. In the case of low impact velocities, this takes an order of magnitude more time than that needed for the formation of a dense nucleus.

We resolve the extended X-ray emission from the prototypical ultraluminous infrared galaxy Arp 220. Extended, faint edge-brightened, soft X-ray lobes outside the optical galaxy are observed to a distance of 10 to 15 kpc on each side of the nuclear region. Bright plumes inside the optical isophotes coincide with the optical line emission and extend 11 kpc from end to end across the nucleus. The data for the plumes cannot be fit by a single temperature plasma, and display a range of temperatures from 0.2 to 1 keV. The plumes emerge from bright, diffuse circumnuclear emission in the inner 3 kpc centered on the Halpha peak, which is displaced from the radio nuclei. There is a close morphological correspondence between the Halpha and soft X-ray emission on all spatial scales. We interpret the plumes as a starburst-driven superwind, and discuss two interpretations of the emission from the lobes in the context of simulations of the merger dynamics of Arp 220.

Colliding galaxies are excellent laboratories for studying galactic evolution and global star formation. Computer simulations of galaxy collisions, in which at least one galaxy has a significant gaseous component, show the production of density enhancements and shock waves in the interstellar medium. These high-density regions coincide with the locations of recent, large-scale star formation in observations of some real colliding galaxies. Thus, combined n-body/hydrodynamic computer simulations can be used to explore the history and conditions of star forming regions in colliding galaxies. The galaxy system Arp 119 (CPG 29) contains a southern member, Arp 119S (Mrk 984), that has an extremely disturbed appearance, and a northern member, Arp 119N, a gas-poor elliptical. The morphology of both members can be fit well by a simulation in which a gas-rich disk galaxy has been impacted by an equal mass elliptical that had a trajectory approximately perpendicular to the plane of the disk and passed through the disk slightly off-center. We compare multi-wavelength observations of the Arp 119 system with a combined n-body/SPH simulation of colliding galaxies. We find that the progression of recent, large-scale star formation in this galaxy can be accounted for by a single, outwardly propagating collision-induced density wave in the gas.

We present a brief overview of some recent observations of colliding galaxies and relevant numerical simulations. These are compared, and details of the locations and history of collision induced star formation are explored, with possible application to star formation at earlier epochs.

In this paper we demonstrate for the first time the connection between the spatial and temporal progression of star formation and the changing locations of the very dense regions in the gas of a massive disk galaxy (NGC 1144) in the aftermath of its collision with a massive elliptical (NGC 1143). These two galaxies form the combined object Arp 118, a collisional ring galaxy system. The results of 3D, time-dependent, numerical simulations of the behavior of the gas, stars, and dark matter of a disk galaxy and the stars and dark matter in an elliptical during a collision are compared with multiwavelength observations of Arp 118. The collision that took place approximately 22 Myr ago generated a strong, non-linear density wave in the stars and gas in the disk of NGC 1144, causing the gas to became clumped on a large scale. This wave produced a series of superstarclusters along arcs and rings that emanate from the central point of impact in the disk. The locations of these star forming regions match those of the regions of increased gas density predicted the time sequence of models. The models also predict the large velocity gradients observed across the disk of NGC 1144. These are due to the rapid radial outflow of gas coupled to large azimuthal velocities in the expanding ring, caused by the impact of the massive intruder.