High-resolution N-body Simulations of Galactic Cannibalism: The Magellanic Stream

Hierarchical clustering represents the favoured paradigm for galaxy formation throughout the Universe; due to its proximity, the Magellanic system offers one of the few opportunities for astrophysicists to decompose the full six-dimensional phase-space history of a satellite in the midst of being cannibalised by its host galaxy. The availability of improved observational data for the Magellanic Stream and parallel advances in computational power has led us to revisit the canonical tidal model describing the disruption of the Small Magellanic Cloud and the consequent formation of the Stream. We suggest improvements to the tidal model in light of these recent advances.Comment: 6 pages, 4 figures, LaTeX (gcdv.sty). Refereed contribution to the 5th Galactic Chemodynamics conference held in Swinburne, July 2003. Accepted for publication in PASA. Version with high resolution figures available at http://astronomy.swin.edu.au/staff/tconnors/publications.htm

Galactic Cannibalism: the Origin of the Magellanic Stream

We are in a privileged location in the Universe which allows us to observe galactic interactions from close range -- the merger of our two nearest dwarf satellite galaxies, the LMC and SMC. It is important to understand the local merger process before we can have confidence in understanding mergers at high redshift. We present high resolution Nbody+SPH simulations of the disruption of the LMC and SMC and the formation of the Magellanic Stream, and discuss the implications for galaxy formation and evolution.Comment: 2 pages, 1 figure, to appear in "The Evolution of Galaxies II: Basic Building Blocks", (2002) ed. M. Sauvage et al. (Kluwer

Simulating a White Dwarf-dominated Galactic Halo

Observational evidence has suggested the possibility of a Galactic halo which is dominated by white dwarfs (WDs). While debate continues concerning the interpretation of this evidence, it is clear that an initial mass function (IMF) biased heavily toward WD precursors (1 < m/Msol < 8), at least in the early Universe, would be necessary in generating such a halo. Within the framework of homogeneous, closed-box models of Galaxy formation, such biased IMFs lead to an unavoidable overproduction of carbon and nitrogen relative to oxygen (as measured against the abundance patterns in the oldest stars of the Milky Way). Using a three-dimensional Tree N-body smoothed particle hydrodynamics code, we study the dynamics and chemical evolution of a galaxy with different IMFs. Both invariant and metallicity-dependent IMFs are considered. Our variable IMF model invokes a WD-precursor-dominated IMF for metallicities less than 5% solar (primarily the Galactic halo), and the canonical Salpeter IMF otherwise (primarily the disk). Halo WD density distributions and C,N/O abundance patterns are presented. While Galactic haloes comprised of ~5% (by mass) of WDs are not supported by our simulations, mass fractions of ~1-2% cannot be ruled out. This conclusion is consistent with the present-day observational constraints.Comment: accepted for publication in MNRA

Interactions of Satellite Galaxies in Cosmological Dark Matter Halos

We present a statistical analysis of the interactions between satellite galaxies in cosmological dark matter halos taken from fully self-consistent high-resolution simulations of galaxy clusters. We show that the number distribution of satellite encounters has a tail that extends to as many as 3-4 encounters per orbit. On average 30% of the substructure population had at least one encounter (per orbit) with another satellite galaxy. However, this result depends on the age of the dark matter host halo with a clear trend for more interactions in younger systems. We also report a correlation between the number of encounters and the distance of the satellites to the centre of the cluster: satellite galaxies closer to the centre experience more interactions. However, this can be simply explained by the radial distribution of the substructure population and merely reflects the fact that the density of satellites is higher in those regions. In order to find substructure galaxies we applied (and present) a new technique based upon the N-body code MLAPM. This new halo finder MHF (MLAPM's-Halo-Finder) acts with exactly the same accuracy as the N-body code itself and is therefore free of any bias and spurious mismatch between simulation data and halo finding precision related to numerical effects.Comment: 6 pages, 4 figures, accepted by PASA (refereed contribution to the 5th Galactic Chemodynamics workshop, July 2003

Mapping Substructures in Dark Matter Halos

We present a detailed study of the real and integrals-of-motion space distributions of a satellite obtained from a self-consistent high-resolution simulation of a galaxy cluster and re-simulated using various analytical halo potentials. We found that the disrupted satellite appears as a coherent structure in integrals-of-motion space in all models (``live'' and analytical potential) although the distribution is significantly smeared for the live host halo. Further the primary mechanism for this smearing is the mass growth of the host, which changes both the energy and angular momentum of the satellite debris. Hence, this must be considered when searching for (stellar) streams with future observational experiments such as RAVE and GAIA.Comment: 5 pages, 6 figures, MNRAS accepted - minor editing without changing the conclusions, a high-resolution version of the paper is available from http://astronomy.swin.edu.au/~sgill/downloads/downloads.htm

ROAM: a Radial-basis-function Optimization Approximation Method for diagnosing the three-dimensional coronal magnetic field

The Coronal Multichannel Polarimeter (CoMP) routinely performs coronal polarimetric measurements using the Fe XIII 10747 $\AA$ and 10798 $\AA$ lines, which are sensitive to the coronal magnetic field. However, inverting such polarimetric measurements into magnetic field data is a difficult task because the corona is optically thin at these wavelengths and the observed signal is therefore the integrated emission of all the plasma along the line of sight. To overcome this difficulty, we take on a new approach that combines a parameterized 3D magnetic field model with forward modeling of the polarization signal. For that purpose, we develop a new, fast and efficient, optimization method for model-data fitting: the Radial-basis-functions Optimization Approximation Method (ROAM). Model-data fitting is achieved by optimizing a user-specified log-likelihood function that quantifies the differences between the observed polarization signal and its synthetic/predicted analogue. Speed and efficiency are obtained by combining sparse evaluation of the magnetic model with radial-basis-function (RBF) decomposition of the log-likelihood function. The RBF decomposition provides an analytical expression for the log-likelihood function that is used to inexpensively estimate the set of parameter values optimizing it. We test and validate ROAM on a synthetic test bed of a coronal magnetic flux rope and show that it performs well with a significantly sparse sample of the parameter space. We conclude that our optimization method is well-suited for fast and efficient model-data fitting and can be exploited for converting coronal polarimetric measurements, such as the ones provided by CoMP, into coronal magnetic field data.Comment: 23 pages, 12 figures, accepted in Frontiers in Astronomy and Space Science