The disruption of an icy moon in a collision with an interloping comet a few hundred million years ago is a simple way to create Saturn's rings. A ring parent moon with a mass comparable to Mimas could be trapped in mean motion resonance with Enceladus and Dione in an orbit near the current outer edge of the rings just beyond the Roche zone. I present collisional N-body simulations of cometary impacts that lead to the partial disruption of a differentiated moon with a rocky core and icy mantle. The core can survive largely intact while the debris from the mantle settles into a ring of predominantly ice particles straddling the orbital radius of the parent moon. The nascent ring spreads radially due to collisional viscosity while mass re-accretes onto the remnant rocky core to form a new moon that can be identified as Mimas. The icy debris that migrates into the Roche zone evolves into Saturn's ring system. Torques from tidal interaction with Saturn and resonant interactions with the rings push the recently formed Mimas outward to its current position on the same timescale of a few hundred million years. This scenario accounts for the high ice fraction observed in Saturn's rings and explains why the ring mass is comparable to the mass of Mimas. The prior existence of a ring parent moon in mean motion resonance results in a tidal heating rate for Enceladus in the recent past that is significantly larger than the current rate.

We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors.

Observations consistently show that elliptical galaxies follow a tight "fundamental plane" scaling relation between size, mean surface brightness and velocity dispersion, with the form $R \propto \sigma^{a}\mu^{b}$. This relation not only has very small (<0.05 dex) intrinsic scatter, but also has significantly different coefficients from the expected virial scaling (a "tilt"). We analyze hundreds of simulations of elliptical galaxies formed from mergers of spiral galaxies in groups to determine if the fundamental plane can emerge from multiple, mostly minor and hierarchical collisionless mergers. We find that these simulated ellipticals lie on a similar fundamental plane with $a \approx 1.7$ and $b \approx 0.3$. The scatter about this plane is not larger than observed, while the tilt is in the correct sense, although $a$ is larger than for typical observations. This supports the idea that collisionless mergers can contribute significantly to the tilt of the fundamental plane. The tilt is mainly driven by a mass-dependent dark matter fraction, such that more massive galaxies have larger dark matter fractions within $R_{e}$. We further discuss the origin of this mass-dependent dark matter fraction and its compatibility with strong lensing observations, as well as the links between the fundamental plane, dynamical masses and the virial theorem.

We present collisionless simulations of dry mergers in groups of three to twenty-five galaxies to test the hypothesis that elliptical galaxies form at the centers of such groups. Mock observations of the central remnants confirm their similarity to ellipticals, despite having no dissipational component. We vary the profile of the original spiral's bulge and find that ellipticals formed from spirals with exponential bulges have too low Sersic indices. Mergers of spirals with de Vaucouleurs (classical) bulges produce remnants with larger Sersic indices correlated with luminosity, as with SDSS ellipticals. Exponential bulge mergers are better fits to faint ellipticals, whereas classical bulge mergers better match luminous ellipticals. Similarly, luminous ellipticals are better reproduced by remnants undergoing many (>5) mergers, and fainter ellipticals by those with fewer mergers. The remnants follow tight size-luminosity and velocity dispersion-luminosity (Faber-Jackson) relations (<0.12 dex scatter), demonstrating that stochastic merging can produce tight scaling relations if the merging galaxies also follow tight scaling relations. The slopes of the size-luminosity and Faber-Jackson relations are close to observations but slightly shallower in the former case. Both relations' intercepts are offset - remnants are too large but have too low dispersions at fixed luminosity. Some remnants show substantial (v/\sigma > 0.1) rotational support, although most are slow rotators and few are very fast rotators (v/\sigma > 0.5). These findings contrast with previous studies concluding that dissipation necessary to produce ellipticals from binary mergers of spirals. Multiple, mostly minor and dry mergers can produce bright ellipticals, whereas significant dissipation could be required to produce faint, rapidly-rotating ellipticals.

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

We present results from a spectroscopic survey of the dwarf spheroidal And XXII and the two extended clusters EC1 and EC2. These three objects are candidate satellites of the Triangulum galaxy, M33, which itself is likely a satellite of M31. We use the DEep Imaging Multi-Object Spectrograph mounted on the Keck-II telescope to derive radial velocities for candidate member stars of these objects and thereby identify the stars that are most likely actual members. Eleven most probable stellar members (of 13 candidates) are found for AndXXII. We obtain an upper limit of sigma_v < 6.0 km s-1 for the velocity dispersion of AndXXII, [Fe/H] ~ -1.6 for its metallicity, and 255pc for the Plummer radius of its projected density profile. We construct a colour magnitude diagram for AndXXII and identify both the red giant branch and the horizontal branch. The position of the latter is used to derive a heliocentric distance to And XXII of 853 pm 26 kpc. The combination of the radial velocity, distance, and angular position of AndXXII indicates that it is a strong candidate for being the first known satellite of M33 and one of the very few examples of a galactic satellite of a satellite. N-body simulations imply that this conclusion is unchanged even if M31 and M33 had a strong encounter in the past few Gyr. We test the hypothesis that the extended clusters highlight tidally stripped galaxies by searching for an excess cloud of halo-like stars in their vicinity. We find such a cloud for the case of EC1 but not EC2. The three objects imply a dynamical mass for M33 that is consistent with previous estimates.

We present a series of hundreds of collisionless simulations of galaxy group mergers. These simulations are designed to test whether the properties of elliptical galaxies - including the key fundamental plane scaling relation, morphology and kinematics - can be simultaneously reproduced by dry multiple mergers in galaxy groups. Preliminary results indicate that galaxy group mergers can produce elliptical remnants lying on a tilted fundamental plane, even without a central dissipational component from a starburst. This suggests that multiple mergers in groups are an alternate avenue for the formation of elliptical galaxies which could well dominate for luminous ellipticals.

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $\mu_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an "S-shaped" appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

The dark matter halos in $\Lambda$CDM cosmological simulations are triaxial and highly flattened. In many cases, these triaxial equilibria are also tumbling slowly, typically about their short axes, with periods of order a Hubble time. Halos may therefore exert a slowly changing external torque on spiral galaxies that can affect their dynamical evolution in interesting ways. We examine the effect of the external torques exerted by a tumbling quadrupolar tidal field on the evolution of spiral galaxies using N-body simulations with realistic, disk galaxy models. We measure the amplitude of the external quadrupole moments of dark halos in cosmological simulations and use these to force disk galaxy models in a series of N-body experiments for a range of pattern speeds. We find that the torques are strong enough to induce long lived transient warps in disks similar to those observed in real spirals and also induce the bar instability at later times in some galaxy models that are otherwise stable for long periods of time in isolation. We also observe forced spiral structure near the edge of the disk where normally self-gravity is too weak to be responsible for such structure. This overlooked influence of dark halos may well be responsible for many of the peculiar aspects of disk galaxy dynamics.

In support of the new Sloan III survey, which will measure the baryon oscillation scale using the luminous red galaxies (LRGs), we have run the largest N-body simulation to date using $4120^3 = 69.9$ billion particles, and covering a volume of $(6.592 h^{-1} {\rm Gpc})^3$. This is over 2000 times the volume of the Millennium Run, and corner-to-corner stretches all the way to the horizon of the visible universe. LRG galaxies are selected by finding the most massive gravitationally bound, cold dark matter subhalos, not subject to tidal disruption, a technique that correctly reproduces the 3D topology of the LRG galaxies in the Sloan Survey.We have measured the covariance function, power spectrum, and the 3D topology of the LRG galaxy distribution in our simulation and made 32 mock surveys along the past lightcone to simulate the Sloan III survey. Our large N-body simulation is used to accurately measure the non-linear systematic effects such as gravitational evolution, redshift space distortion, past light cone space gradient, and galaxy biasing, and to calibrate the baryon oscillation scale and the genus topology. For example, we predict from our mock surveys that the baryon acoustic oscillation peak scale can be measured with the cosmic variance-dominated uncertainty of about 5% when the SDSS-III sample is divided into three equal volume shells, or about 2.6% when a thicker shell with $0.4<z<0.6$ is used. We find that one needs to correct the scale for the systematic effects amounting up to 5.2% to use it to constrain the linear theories. And the uncertainty in the amplitude of the genus curve is expected to be about 1% at 15 $h^{-1}$Mpc scale. We are making the simulation and mock surveys publicly available.

We examine the bar instability in galactic models with an exponential disk and a cuspy dark matter (DM) halo with a Navarro-Frenk-White (NFW) cosmological density profile. We construct equilibrium models from a 3-integral composite distribution function that are subject to the bar instability. We generate a sequence of models with a range of mass resolution from 1.8K to 18M particles in the disk and 10K to 100M particles in the halo along with a multi-mass model with an effective resolution of ~10^10 particles. We describe how mass resolution affects the bar instability, including its linear growth phase, the buckling instability, pattern speed decay through the resonant transfer of angular momentum to the DM halo, and the possible destruction of the halo cusp. Our higher resolution simulations show a converging spectrum of discrete resonance interactions between the bar and DM halo orbits. As the pattern speed decays, orbital resonances sweep through most of the DM halo phase space and widely distribute angular momentum among the halo particles. The halo does not develop a flat density core and preserves the cusp, except in the region dominated by gravitational softening. The formation of the bar increases the central stellar density and the DM is compressed adiabatically increasing the halo central density by 1.7X. Overall, the evolution of the bar displays a convergent behavior for halo particle numbers between 1M and 10M particles, when comparing bar growth, pattern speed evolution, the DM halo density profile and a nonlinear analysis of the orbital resonances. Higher resolution simulations clearly illustrate the importance of discrete resonances in transporting the angular momentum from the bar to the halo.

We examine the bar instability in models with an exponential disk and a cuspy NFW-like dark matter (DM) halo inspired by cosmological simulations. Bar evolution is studied as a function of numerical resolution in a sequence of models spanning 10K to 100M DM particles - including a multi-mass model with an effective resolution of 10G. The goal is to find convergence in dynamical behaviour. We characterize the bar growth, the buckling instability, pattern speed decay through resonant transfer of angular momentum, and possible destruction of the DM halo cusp. Overall, most characteristics converge in behaviour in detail for halos containing more than 10M particles. Notably, the formation of the bar does not destroy the density cusp in this case. These higher resolution simulations clearly illustrate the importance of discrete resonances in transporting angular momentum from the bar to the halo.

I begin with a brief history of N-body simulation and visualization and then go on to describe various methods for creating images and animations of modern simulations in cosmology and galactic dynamics. These techniques are incorporated into a specialized particle visualization software library called MYRIAD that is designed to render images within large parallel N-body simulations as they run. I present several case studies that explore the application of these methods to animations of star clusters, interacting galaxies and cosmological structure formation.

We explore the secular dynamical evolution of an N-body model of M31 in the presence of a population of 100 dark matter satellites over 10 Gyr. The satellite population has structural and kinematic characteristics modelled to follow the predictions of Lambda-CDM cosmological simulations. Vertical disk heating is a small effect despite many interactions with the satellite population with only a 20% increase in vertical velocity dispersion sigma_z and the disk scale height z_d at the equivalent solar radius R = 2.5R_d . However, the stellar disk is noticeably flared after 10 Gyr with z_d nearly doubling at the disk edge. Azimuthal disk heating is much larger with sigma_R and sigma_z both increasing by 1.7x. However, in a control experiment without satellites dispersion increases by 1.5x suggesting that most of the effect is due to heating through scattering off of spiral structure excited by swing-amplified noise. Surprisingly, direct impacts of satellites on the disk can excite spiral structure with a significant amplitude and in some cases impacts close to the disk center also induce the bar instability. The large number of dark matter satellite impacts expected over a galaxy's lifetime may be a significant source of external perturbations for driving disk secular evolution.

GRAVITAS is a self-published DVD that presents a visual and musical celebration of the beauty in a dynamic universe driven by gravity. Animations from supercomputer simulations of forming galaxies, star clusters, galaxy clusters, and galaxy interactions are presented as moving portraits of cosmic evolution. Billions of years of complex gravitational choreography are presented in 9 animations - each one interpreted with an original musical composition inspired by the exquisite movements of gravity. The result is an emotive and spiritually uplifting synthesis of science and art. The GRAVITAS DVD has been out for two years now but I am now making the DVD disk image freely available for personal and educational use through a bittorrent download. Download and burn at your leisure. The animations are also downloadable in various video formats.

We present an axisymmetric, equilibrium model for late-type galaxies which consists of an exponential disk, a Sersic bulge, and a cuspy dark halo. The model is specified by a phase space distribution function which, in turn, depends on the integrals of motion. Bayesian statistics and the Markov Chain Monte Carlo method are used to tailor the model to satisfy observational data and theoretical constraints. By way of example, we construct a chain of 10^5 models for the Milky Way designed to fit a wide range of photometric and kinematic observations. From this chain, we calculate the probability distribution function of important Galactic parameters such as the Sersic index of the bulge, the disk scale length, and the disk, bulge, and halo masses. We also calculate the probability distribution function of the local dark matter velocity dispersion and density, two quantities of paramount significance for terrestrial dark matter detection experiments. Though the Milky Way models in our chain all satisfy the prescribed observational constraints, they vary considerably in key structural parameters and therefore respond differently to non-axisymmetric perturbations. We simulate the evolution of twenty-five models which have different Toomre Q and Goldreich-Tremaine X parameters. Virtually all of these models form a bar, though some, more quickly than others. The bar pattern speeds are ~ 40 - 50 km/s/kpc at the time when they form and then decrease, presumably due to coupling of the bar with the halo. Since the Galactic bar has a pattern speed ~50 km/s/kpc we conclude that it must have formed recently.

We investigate the evolution of a population of 100 dark matter satellites orbiting in the gravitational potential of a realistic model of M31. We find that after 10 Gyr, seven subhalos are completely disrupted by the tidal field of the host galaxy. The remaining satellites suffer heavy mass loss and overall, 75% of the mass initially in the subhalo system is tidally stripped. Not surprisingly, satellites with pericentric radius less than 30 kpc suffer the greatest stripping and leave a complex structure of tails and streams of debris around the host galaxy. Assuming that the most bound particles in each subhalo are kinematic tracers of stars, we find that the halo stellar population resulting from the tidal debris follows an r^{-3.5} density profile at large radii. We construct B-band photometric maps of stars coming from disrupted satellites and find conspicuous features similar both in morphology and brightness to the observed Giant Stream around Andromeda. An assumed star formation efficiency of 5-10% in the simulated satellite galaxies results in good agreement with the number of M31 satellites, the V-band surface brightness distribution, and the brightness of the Giant Stream. During the first 5 Gyr, the bombardment of the satellites heats and thickens the disk by a small amount. At about 5 Gyr, satellite interations induce the formation of a strong bar which, in turn, leads to a significant increase in the velocity dispersion of the disk.

We describe a new set of self-consistent, equilibrium disk galaxy models that incorporate an exponential disk, a Hernquist model bulge, an NFW halo and a central supermassive black hole. The models are derived from explicit distribution functions for each component and the large number of parameters permit detailed modeling of actual galaxies. We present techniques that use structural and kinematic data such as radial surface brightness profiles, rotation curves and bulge velocity dispersion profiles to find the best-fit models for the Milky Way and M31. Through N-body realizations of these models we explore their stability against the formation of bars. The models permit the study of a wide range of dynamical phenomenon with a high degree of realism.

We report the search for intracluster light in four Abell Type II/III (non-cD) galaxy clusters: Abell 801, 1234, 1553, & 1914. We find on average that these clusters contain $\sim$ 10% of their detected stellar luminosity in a diffuse component. We show that for two of the clusters the intracluster light closely follows the galaxy distribution, but in the other two cases, there are noticeable differences between the spatial distribution of the galaxies and the intracluster light. We report the results of a search for intracluster tidal debris in each cluster, and note that Abell 1914 in particular has a number of strong tidal features likely due to its status as a recent cluster merger. One of the Abell 1914 features appears to be spatially coincident with an extension seen in weak lensing maps, implying the feature traces a large amount of mass. We compare these results to numerical simulations of hierarchically-formed galaxy clusters, and find good general agreement between the observed and simulated images, although we also find that our observations sample only the brightest features of the intracluster light. Together, these results suggest that intracluster light can be a valuable tool in determining the evolutionary state of galaxy clusters.

Weak lensing measurements are entering a precision era to statistically map the distribution of matter in the universe. The most common measurement has been of the variance of the projected surface density of matter, which corresponds to the induced correlation in alignments of background galaxies. This measurement of the fluctuations is insensitive to the total mass content, like using waves on the ocean to measure its depths. But when the depth is shallow as happens near a beach, waves become skewed. Similarly, a measurement of skewness in the projected matter distribution directly measures the total matter content of the universe. While skewness has already been convincingly detected, its constraint on cosmology is still weak. We address optimal analyses for the CFHT Legacy Survey in the presence of noise. We show that a compensated Gaussian filter with a width of 2.5 arc minutes optimizes the cosmological constraint, yielding $\Delta \Omega_m/\Omega_m\sim 10%$. This is significantly better than other filters which have been considered in the literature. This can be further improved with tomography and other sophisticated analyses.

We examine the observable properties of simulated barred galaxies including radial mass profiles, edge-on structure and kinematics, bar lengths and pattern speed evolution for detailed comparison to real systems. We have run several simulations in which bars are created through inherent instabilities in self-consistent simulations of a realistic disc+halo galaxy model with a disc-dominated, flat rotation curve. These simulations were run at high (N=20M particles) and low (N=500K) resolution to test numerical convergence. We determine the pattern speeds in simulations directly from the phase angle of the bar versus time and the Tremaine-Weinberg method. Fundamental dynamics do not change between the high and low resolution, suggesting that convergence has been reached in this case. We find the higher resolution is needed to simulate structural and kinematic properties accurately. The edge-on view of the higher-resolution system clearly shows the bending instability and formation of a peanut-shaped bulge. We determined bar lengths by different means to determine the simulated bar is fast, with a corotation to bar length ratio under 1.5. Simulated bars in these models form with pattern speeds slower than those observed and slow down during their evolution. Dynamical friction between the bar and dark halo is responsible for this deceleration, as revealed by the transfer of angular momentum between the disc and the halo. However, even though the pattern speed is reduced at later times, the instantaneous scale length of the disc has grown sufficiently for the bar motion to agree with many observations. By using a different model and simulation technique than other authors, we are able to compare the robustness of these methods.

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 describe a parallel, cosmological N-body code based on a hybrid scheme using the particle-mesh (PM) and Barnes-Hut (BH) oct-tree algorithm. We call the algorithm GOTPM for Grid-of-Oct-Trees-Particle-Mesh. The code is parallelized using the Message Passing Interface (MPI) library and is optimized to run on Beowulf clusters as well as symmetric multi-processors. The gravitational potential is determined on a mesh using a standard PM method with particle forces determined through interpolation. The softened PM force is corrected for short range interactions using a grid of localized BH trees throughout the entire simulation volume in a completely analogous way to P$^3$M methods. This method makes no assumptions about the local density for short range force corrections and so is consistent with the results of the P$^3$M method in the limit that the treecode opening angle parameter, $\theta \to 0$. (abridged)

Weak gravitational lensing provides a direct statistical measure of the dark matter distribution. The variance is easiest to measure, which constrains the degenerate product \sigma_8\Omega^0.6. The degeneracy is broken by measuring the skewness arising from the fact that densities must remain positive, which is not possible when the initially symmetric perturbations become non-linear. Skewness measures the non-linear mass scale, which in combination with the variance measures \Omega directly. We present the first detection of dark matter skewness from the Virmos-Decart survey. We have measured the full three point function, and its projections onto windowed skewness. We separate the lensing mode and the B mode. The lensing skewness is detected for a compensated Gaussian on scales of 5.37 arc minutes to be \kappa^3=1.06+/-0.06x10^-6. The B-modes are consistent with zero at this scale. The variance for the same window function is \kappa^2= 5.32+/-0.62+/-0.98x10^-5, resulting in S_3=375^{+342}_{-124}. Comparing to N-body simulations, we find \Omega_0<0.5 at 90% confidence. The Canada-France-Hawaii-Telescope legacy survey and newer simulations should be able to improve significantly on the constraint.

We conduct an analysis of the shapes of molecular cloud cores using recently compiled catalogs of observed axis ratios of individual cores mapped in ammonia or through optical selection. We apply both analytical and statistical techniques to deproject the observed axis ratios in order to determine the true distribution of cloud core shapes. We find that neither pure oblate nor pure prolate cores can account for the observed distribution of core shapes. Intrinsically triaxial cores produce distributions which agree with observations. The best-fit triaxial distribution contains cores which are more nearly oblate than prolate.

After studying the orientation of a warp generated by a companion satellite, we show that the Galactic Warp would be oriented in a different way if the Magellanic Clouds were its cause. We have treated the problem analytically, and complemented it with numerical N-Body simulations.

We discuss how to use tidal streams from globular clusters to measure the mass distribution of the Milky Way. Recent proper motion determinations for globular clusters from plate measurements and Hipparcos astrometry provide several good candidates for Galactic mass determinations in the intermediate halo, far above the Galactic disk, including Pal 5, NGC 4147, NGC 5024 (M53) and NGC 5466; the remaining Hipparcos clusters provide candidates for measurements several kpc above and below the disk. These clusters will help determine the profile and shape of the inner halo. To aid this effort, we present two methods of mass determination: one, a generalization of rotation-curve mass measurements, which gives the mass and potential from complete position-velocity observations for stream stars; and another using a simple chi^2 estimator, which can be used when only projected positions and radial velocities are known for stream stars. We illustrate the use of the latter method using simulated tidal streams from Pal 5 and find that fairly accurate mass determinations are possible even for relatively poor data sets. Follow-up observations of clusters with proper motion determinations may reveal tidal streams; obtaining radial velocity measurements would enable accurate measurements of the mass distribution in the inner Galaxy.

The evolution of galaxies is driven strongly by dynamical processes including internal instabilities, tidal interactions and mergers. The cluster environment is a useful laboratory for studying these effects. I present recent results on simulations of interacting populations of spiral and elliptical galaxies in the cosmological collapse of a cluster, showing the formation of the central, brightest cluster galaxies through merging, the effect of tidal interactions and merging on galaxy morphology over cosmic history, and the distribution and kinematics of the tidal debris field. (See http://www.cita.utoronto.ca/~dubinski/rutgers98 for extensive figures and animations.)

We present an extensive parameter survey to study the influence of halo mass profiles on the development of tidal tails in interacting disk galaxies. We model the galaxies using a fixed exponential disk with a central bulge and vary the halo potential over a range of parameters using both the Hernquist and NFW mass distributions, probing the effect of the halo mass, extent and concentration. We examine the consistency of the results against both observational and theoretical constraints on halo profiles and comment on the failures and weaknesses of different models. Galaxies with rising or flat rotation curves dominated by the halo are inhibited from forming tidal tails unless the halo is abruptly cut off just beyond the disk edge. Conversely, models with declining rotation curves -- resulting either from compact, low mass halos, or from massive disk components in low concentration dark halos -- produce tidal tails very similar to those observed in well-studied interacting systems. As argued by Springel & White (1998), a unifying, quantitative relation for all cases is that the ratio of the escape velocity to the circular velocity at around the solar radius must be v_e/v_c<2.5 for tidal tails to be produced. The galaxy models which appear to fit most of the observational constraints are those which have disk-dominated rotation curves and low concentration halos. We discuss our results in a cosmological context using recent studies which link halo properties to cosmological models.

Most clusters and groups of galaxies contain a giant elliptical galaxy in their centres which far outshines and outweighs normal ellipticals. The origin of these brightest cluster galaxies is intimately related to the collapse and formation of the cluster. Using an N-body simulation of a cluster of galaxies in a hierarchical cosmological model, we show that galaxy merging naturally produces a massive, central galaxy with surface brightness and velocity dispersion profiles similar to observed BCG's. To enhance the resolution of the simulation, 100 dark halos at $z=2$ are replaced with self-consistent disk+bulge+halo galaxy models following a Tully-Fisher relation using 100000 particles for the 20 largest galaxies and 10000 particles for the remaining ones. This technique allows us to analyze the stellar and dark matter components independently. The central galaxy forms through the merger of several massive galaxies along a filament early in the cluster's history. Galactic cannibalism of smaller galaxies through dynamical friction over a Hubble time only accounts for a small fraction of the accreted mass. The galaxy is a flattened, triaxial object whose long axis aligns with the primordial filament and the long axis of the cluster galaxy distribution agreeing with observed trends for galaxy-cluster alignment.

Massive and extended dark halos can inhibit the formation of long tidal tails in galaxy collisions. We examine this effect using an extensive survey of simulations with different dark halo potentials to constrain halo properties of interacting galaxies. These constraints are compared to other observational limits and theoretical predictions of halo structure. The dark halos predicted by $\Omega=1$ cosmological models like CDM are too massive and extended to produce the long tidal tails seen in nearby galaxy collisions. There is also a conflict with the halo potentials inferred from satellite kinematics; such halos would likewise inhibit tail formation in galaxy collisions.

Secondary infall in galaxies could cause the angular momentum of the outer halo to change its orientation. The generation of warps due to this effect is studied using N-body simulations.

We conduct gravitational microlensing experiments in a galaxy taken from a cosmological N-body simulation. Hypothetical observers measure the optical depth and event rate toward hypothetical LMCs and compare their results with model predictions. Since we control the accuracy and sophistication of the model, we can determine how good it has to be for statistical errors to dominate over systematic ones. Several thousand independent microlensing experiments are performed. When the ``best-fit'' triaxial model for the mass distribution of the halo is used, the agreement between the measured and predicted optical depths is quite good: by and large the discrepancies are consistent with statistical fluctuations. If, on the other hand, a spherical model is used, systematic errors dominate. Even with our ``best-fit'' model, there are a few rare experiments where the deviation between the measured and predicted optical depths cannot be understood in terms of statistical fluctuations. In these experiments there is typically a clump of particles crossing the line of sight to the hypothetical LMC. These clumps can be either gravitationally bound systems or transient phenomena in a galaxy that is still undergoing phase mixing. Substructure of this type, if present in the Galactic distribution of Machos, can lead to large systematic errors in the analysis of microlensing experiments. We also describe how hypothetical WIMP and axion detection experiments might be conducted in a simulated N-body galaxy.

We examine the effect of different dark matter halo potentials on the morphology and kinematics of tidals tails in a merger model of NGC 7252. We find that models of merging galaxies with low halo masses of Mh = 4-8 M{disk+bulge} (Mdb) can fit the observed morphology and kinematics of the NGC 7252 tails while galaxies with high mass halos (Mh = 16-32 Mdb) fail in this respect. In high mass models, the deep potential only allows weakly bound disk material (stars or gas) at R > 5 disk scale lengths to be ejected in tidal tails which tend to fall back onto the parent galaxies before the final merger. Galaxies with massive, low density halos are somewhat more successful at ejecting tidal debris during mergers, but still have difficulties recreating the thin, gas-rich tails observed in NGC 7252. Our models suggest upper limits for the dark halo masses in the NGC 7252 progenitor galaxies of roughly Mh <= 10 Mdb. We note, however, that our calculations have focused on the rather idealized case of the isolated merging of galaxies with distinct dark matter halos; calculations which employ more realistic (cosmological) initial conditions are needed to fully explore the use of tidal tails in constraining dark matter in galaxies.

We describe PTreeSPH, a gravity treecode combined with an SPH hydrodynamics code designed for massively parallel supercomputers having distributed memory. Our computational algorithm is based on the popular TreeSPH code of Hernquist & Katz (1989). PTreeSPH utilizes a domain decomposition procedure and a synchronous hypercube communication paradigm to build self-contained subvolumes of the simulation on each processor at every timestep. Computations then proceed in a manner analogous to a serial code. We use the Message Passing Interface (MPI) communications package, making our code easily portable to a variety of parallel systems. PTreeSPH uses individual smoothing lengths and timesteps, with a communication algorithm designed to minimize exchange of information while still providing all information required to accurately perform SPH computations. We have additionally incorporated cosmology, periodic boundary conditions with forces calculated using a quadrupole Ewald summation method, and radiative cooling and heating from a parameterized ionizing background following Katz, Weinberg & Hernquist (1996). The addition of other physical processes, such as star formation, is straightforward. A cosmological simulation from z=49 to z=2 with 64^3 gas particles and 64^3 dark matter particles requires ~6000 node-hours on a Cray T3D, with a communications overhead of ~10% and is load balanced to a ~90% level. When used on the new Cray T3E, this code will be capable of performing cosmological hydrodynamical simulations down to z=0 with ~2x10^6 particles, or to z=2 with ~10^7 particles, in a reasonable amount of time. Even larger simulations will be practical in situations where the matter is not highly clustered or when periodic boundaries are not required.

We describe a new implementation of a parallel N-body tree code. The code is load-balanced using the method of orthogonal recursive bisection to subdivide the N-body system into independent rectangular volumes each of which is mapped to a processor on a parallel computer. On the Cray T3D, the load balance in the range of 70-90\% depending on the problem size and number of processors. The code can handle simulations with $>$ 10 million particles roughly a factor of 10 greater than allowed in vectorized tree codes.

We use simulations of merging galaxies to explore the sensitivity of the morphology of tidal tails to variations of the halo mass distributions in the parent galaxies. Our goal is to constrain the mass of dark halos in well-known merging pairs. We concentrate on prograde encounters between equal mass galaxies which represent the best cases for creating tidal tails, but also look at systems with different relative orientations, orbital energies and mass ratios. As the mass and extent of the dark halo increase in the model galaxies, the resulting tidal tails become shorter and less massive, even under the most favorable conditions for producing these features. Our simulations imply that the observed merging galaxies with long tidal tails ($\sim 50-100$ kpc) such as NGC 4038/39 (the Antennae) and NGC 7252 probably have halo:disk+bulge mass ratios less than 10:1. These results conflict with the favored values of the dark halo mass of the Milky Way derived from satellite kinematics and the timing argument which give a halo:disk+bulge mass ratio of $\sim 30:1$. However, the lower bound of the estimated dark halo mass in the Milky Way (mass ratio $\sim 10:1$) is still consistent with the inferred tidal tail galaxy masses. Our results also conflict with the expectations of $\Omega=1$ cosmologies such as CDM which predict much more massive and extended dark halos.

We describe methods for setting up self-consistent disk-bulge-halo galaxy models. The bulge and halo distribution functions (\df) are functions of $E$ and $L_z$ only. The halo's flattening and rotation can be specified. The disk \df\ is a function of $E$ and $L_z$ and a third ``integral'', $E_z$, the vertical energy, which is approximately conserved in a warm disk with vertical extent. The models also have finite extent making them suitable for N-body simulation. A simulation of a sample model shows that in practice the models are very close to equilibrium making them ideal for experiments on instabilities in galactic disks. We also present a sequence of models closely resembling the Milky Way mass distribution with 5 exponential scale radii and varying halo mass and radial extent.

We generate random Gaussian turbulent velocity fields with a Kolmogorov spectrum and use these to obtain synthetic line-of-sight velocity profiles. The profiles are found to be similar to line profiles observed in molecular clouds. We suggest methods for analysing measured line profiles to test whether they might arise from Gaussian Kolmogorov turbulence.

When a galaxy forms, the disk may initially be tilted with respect to a flattened dark halo. The misalignment between the disk and the halo is a common explanation for galactic disk warps, since in this state, disks have precessing bending modes which resemble real warps. A dubious assumption in this theory is that the gravitational response of the halo is negligible. We therefore calculate the response of an oblate halo to a precessing inclined exponential disk using a variety of techniques. We construct models with a rigid exponential disk precessing in a particle halo, a particle disk precessing inside a static bulge/halo potential, and a self-consistent model with a particle disk, bulge and halo. When the disk:halo mass ratio is small ($\sim$10\%) within 5 exponential scale radii, the disk settles to the equatorial plane of the halo within 5 orbital times. When the disk and halo mass are comparable, the halo rapidly aligns with the disk within a few orbital times, while the disk inclination drops. The rapid response of the halo to a inclined disk suggests that the warps seen in galactic disks are not likely the bending modes of a precessing disk inside of triaxial halos. If a galaxy forms inclined to the principal plane of a dark halo, either the disk will settle to a principal plane or the inner halo will twist to align with the disk. The outer halo will remain misaligned for a much longer time and therefore may still exert a torque. Warped bending modes may still exist if the misalignment of the outer halo persists for a Hubble time.

We construct models of highly-inclined (``polar'') rings in an external potential including both self--gravity and dissipation due to a drag force. We adopt the oblate spheroidal scale--free logarithmic potential with axis ratio $q=0.85$ and an initial inclination of 80$^\circ$ for the self--gravitating rings. Stellar (dissipationless) rings suffer from mass loss during their evolution which drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than $\sim$ 0.02 of the enclosed galaxy mass, rings remain bound and and are not destroyed by differential precession. The rings precess at a constant rate for more than a precession period $\tau_p$ finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self--gravity if the characteristic viscous inflow time scale $\tau_{vi}$ is shorter than $\sim 25 \tau_p$. Rings with $\tau_{vi}/\tau_p \ \sles \ 25$ collapse toward the nucleus of the potential within one precession period independent of the amount of self--gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.

We describe an analytic distribution function of a finite, oblate stellar system that is useful for the practical modelling of dark halos. The function is determined by lowering Evans's (1993) distribution function of a flattened, cored isothermal system in analogy to the lowering of the singular, isothermal sphere in the definition of the King (1966) models. We derive analytic expressions of the density, maximal streaming velocity and velocity dispersion profiles as a function of the potential. As for the King models, the potential must be calculated numerically. We also present a recipe for generating N-body realizations of this distribution function and examine the stability in three models with dimensionless spins $\lambda=$0.0, 0.05 and 0.18 using N-body simulations with 50,000 particles. The $\lambda=0.18$ model is unstable to the formation of a triaxial bar within $\sim 5$ King radii while the other models appear stable. We conclude that the slowly rotating systems are useful for modelling flattened dark halos.

The dissipative infall of gas during the formation of a galaxy modifies the density profile and shape of the dark halo. Gas dissipates energy radiatively and sinks to the center of the dark halo forming the luminous part of a galaxy. The resulting central density enhancement can alter the halo's orbital distribution. We simulate dissipative infall inside of an initially triaxial N-body dark halo by slowly growing a potential in the center of the particle distribution. The dark halo transforms from a prolate-triaxial halo ($T \sim 0.8$) to an oblate-triaxial halo ($T \sim 0.5$) while approximately preserving the flattening ($c/a \sim 0.5$). The main implication is that dark halos are rounder and more oblate than previous predictions of purely collisionless simulations with the new constraint that $b/a \gapp 0.7$. At the same time, the distribution of intrinsic flattenings of dark halos ($\langle c/a \rangle = 0.5$, $\langle (c/a)^2 \rangle ^{1/2} = 0.15$) is preserved during the period of baryonic infall. The oval distortions of disk galaxies should therefore be slightly less than original predictions from collisionless dark halos. The predicted distribution of shapes for dark halos from cosmological N-body simulations is in better agreement with distribution of ellipticities of elliptical galaxies if we assume all halos and galaxies become oblate-triaxial in response to baryonic dissipation. The observed distribution of kinematic misalignment angles is also consistent with this shape distribution if