A Generalized Expression Optimization Hook for C++ on High-Performance Architectures

C++ has gained broad acceptance as an object-oriented evolutionary extension to the C language, but it severely constrains methods for operating on class objects by forcing all data manipulation through an interface which assumes that all basic operations can be implemented as they are written: as unary or binary operators. C++ allows great flexibility in the creation of complex data structures which can perform the same functionality as built-in types of many other languages, but unfortunately it does not allow an equivalent level of flexibility so that operators acting on those data types can achieve the same level of efficiency as their counterparts in language-level implementations. This limitation becomes even more pronounced on high-performance computers whose advanced features require considerable cooperation between the algorithm, the compiler, and the architecture for maximum performance. This paper describes a language enhancement proposal: a special member function which interacts with overload operators as complete expressions

Hierarchical Tree-Structures as Adaptive Meshes

Introduction: Two basic types of simulations exist for modeling systems of many particles: grid-based (point particles indirectly interacting with one another through the potential calculated from equivalent particle densities on a mesh) and particle-based (point particles directly interacting with a one another through potentials at their positions calculated from the other particles in the system). Grid-based solvers traditionally model continuum problems, such as fluid and gas systems, and mixed particle-continuum systems. Particle-based solvers find more use modeling discrete systems such as stars within galaxies or other rarefied gases. Many different physical systems, including electromagnetic interactions, gravitational interactions, and fluid vortex interactions, all are governed by Poisson\u27s Equation..

A 2MASS All-Sky View of the Sagittarius Dwarf Galaxy: IV. Modeling the Sagittarius Tidal Tails

M giants recovered from the Two Micron All-Sky Survey (2MASS) have recently been used to map the position and velocity distributions of tidal debris from the Sagittarius (Sgr) dwarf spheroidal galaxy entirely around the Galaxy. We compare this data set to both test particle orbits and N-body simulations of satellite destruction run within a variety of rigid Milky Way potentials and find that the mass of the Milky Way within 50 kpc of its center should be 3.8-5.6 x 10^11 Msun in order for any Sgr orbit to simultaneously fit the velocity gradient in the Sgr trailing debris and the apocenter of the Sgr leading debris. Orbital pole precession of young debris and leading debris velocities in regions corresponding to older debris provide contradictory evidence in favor of oblate/prolate Galactic halo potentials respectively, leading us to conclude that the orbit of Sgr has evolved over the past few Gyr. Based upon the velocity dispersion and width along the trailing tidal stream we estimate the current bound mass of Sgr to be M_Sgr = 2 - 5 x 10^8 Msun independant of the form of the Galactic potential; this corresponds to a range of mass to light ratios (M/L)_Sgr = 14 - 36 (M/L)_Sun for the Sgr core. Models with masses in this range best fit the apocenter of leading Sgr tidal debris when they orbit with a radial period of roughly 0.85 Gyr and have periGalactica and apoGalactica of about 15 kpc and 60 kpc respectively. These distances will scale with the assumed distance to the Sgr dwarf and the assumed depth of the Galactic potential. The density distribution of debris along the orbit in these models is consistent with the M giant observations, and debris at all orbital phases where M giants are obviously present is younger (i.e. was lost more recently from the satellite) than the typical age of a Sgr M giant star.Comment: 42 pages, 13 figures; Accepted for publication by ApJ (October 08, 2004; originally submitted May 10, 2004). Fixed typos and added references. PDF file with high resolution figures may be downloaded from http://www.astro.caltech.edu/~drlaw/Papers/Sgr_paper4.pd

Evidence for a Triaxial Milky Way Dark Matter Halo from the Sagittarius Stellar Tidal Stream

Observations of the lengthy tidal streams produced by the destruction of the Sagittarius dwarf spheroidal (Sgr dSph) are capable of providing strong constraints on the shape of the Galactic gravitational potential. However, previous work, based on modeling different stream properties in axisymmetric Galactic models has yielded conflicting results: while the angular precession of the Sgr leading arm is most consistent with a spherical or slightly oblate halo, the radial velocities of stars in this arm are only reproduced by prolate halo models. We demonstrate that this apparent paradox can be resolved by instead adopting a triaxial potential. Our new Galactic halo model, which simultaneously fits all well-established phase space constraints from the Sgr stream, provides the first conclusive evidence for, and tentative measurement of, triaxiality in an individual dark matter halo. The Milky Way halo within ~ 60 kpc is best characterized by a minor/major axis ratio of the isovelocity contours c/a ~ 0.67, intermediate/major axis ratio b/a ~ 0.83, and triaxiality parameter T ~ 0.56. In this model, the minor axis of the dark halo is coincident with the Galactic X axis connecting the Sun and the Galactic Center to within ~ 15 degrees, while the major axis also lies in the Galactic plane, approximately along the Galactic Y axis.Comment: Accepted for publication in ApJ Letters. 3 figure

The survival of the Sgr dwarf galaxy and the flatness of the rotation curve of the Galaxy

How has the ``fluffy'' core of the Sgr dwarf galaxy survived multiple strong shocks from the tidal force of the Galactic halo and disc since the formation of the core a Hubble time ago? A scenario that Sgr was deflected to its current orbit by the Magellanic Clouds after a rendezvous on the north Galactic pole 2-3 Gyrs ago is examined. It is shown that the conditions of the collision fix both the sense of circulation of Sgr and the LMC around the Galaxy and the slope of the Galactic rotation curve. The model argues that the two orthogonal polar circles traced by a dozen or so Galactic halo dwarf galaxies and globular clusters (LMC-SMC-Magellanic Stream-Draco-Ursa Minor along $l \approx 270^o$ and M54-Ter 7-Ter 8-Arp 2-NGC 2419-Pal 15 along $l \approx 0^o$) are streams of tidal relics from two ancient galaxies which was captured on two intersecting polar rosette orbits by the Galaxy. Our results favor the interpretation of microlensing towards the LMC being due to source or lens stars in tidal features of the Magellanic Clouds. We discuss direct and indirect observations to test the collision scenario.Comment: 11 pages including 3 ps figures and an added appendix for deriving the strength and probability of the encounte

Detection of a Corrugated Velocity Pattern in the Spiral Galaxy NGC 5427

Here we report the detection, in Halpha emission, of a radial corrugation in the velocity field of the spiral galaxy NGC 5427. The central velocity of the Halpha line displays coherent, wavy-like variations in the vicinity of the spiral arms. The spectra along three different arm segments show that the maximum amplitude of the sinusoidal line variations are displaced some 500 pc from the central part of the spiral arms. The peak blueshifted velocities appear some 500 pc upstream the arm, whereas the peak redshifted velocities are located some 500 pc downstream the arm. This kinematical behavior is similar to the one expected in a galactic bore generated by the interaction of a spiral density wave with a thick gaseous disk, as recently modeled by Martos & Cox (1998).Comment: Accepted for publication in Ap

The Northern wraps of the Sagittarius Stream as traced by Red Clump stars: distances, intrinsic widths and stellar densities

We trace the tidal Stream of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using Red Clump stars from the catalog of the Sloan Digital Sky Survey - Data Release 6, in the range 150{\deg} < RA < 220{\deg}, corresponding to the range of orbital azimuth 220{\deg} < Lambda < 290{\deg}. Substructures along the line of sight are identified as significant peaks in the differential star count profiles (SCP) of candidate Red Clump stars. A proper modeling of the SCPs allows us to obtain: (a) <10% accurate, purely differential distances with respect to the main body of Sgr, (b) estimates of the FWHM along the line of sight, and (c) estimates of the local density, for each detected substructure. In the range 255{\deg} < Lambda < 290{\deg} we cleanly and continuously trace various coherent structures that can be ascribed to the Stream, in particular: the well known northern portion of the leading arm, running from d~43 kpc at Lambda ~ 290{\deg} to d ~ 30 kpc at Lambda ~ 255{\deg}, and a more nearby coherent series of detections lying at constant distance d ~ 25 kpc, that can be identified with a wrap of the trailing arm. The latter structure, predicted by several models of the disruption of Sgr dSph, was never traced before; comparison with existing models indicates that the difference in distance between these portions of the leading and trailing arms may provide a powerful tool to discriminate between theoretical models assuming different shapes of the Galactic potential. A further, more distant wrap in the same portion of the sky is detected only along a couple of lines of sight.[abridged]Comment: 31 pages, 28 figures, accepted for publication in the Astrophysical Journal, a version with figures at full resolution can be downloaded at the following URL: http://www.bo.astro.it/SGR

The Vertical Structure of the Outer Milky Way HI Disk

We examine the outer Galactic HI disk for deviations from the b=0 plane by constructing maps of disk surface density, mean height, and thickness. We find that the Galactic warp is well described by a vertical offset plus two Fourier modes of frequency 1 and 2, all of which grow with Galactocentric radius. Adding the m=2 mode accounts for the large asymmetry between the northern and southern warps. We use a Morlet wavelet transform to investigate the spatial and frequency localization of higher frequency modes; these modes are often referred to as "scalloping." We find that the m=10 and 15 scalloping modes are well above the noise, but localized; this suggests that the scalloping does not pervade the whole disk, but only local regions.Comment: Accepted by ApJ. 17 pages, 18 figures. Color maps are available at http://astron.berkeley.edu/~elevin