grommet: an N-body code for high-resolution simulations of individual galaxies

This paper presents a fast, economical particle-multiple-mesh N -body code optimized for large- N modelling of collisionless dynamical processes, such as black hole wandering or bar–halo interactions, occurring within isolated galaxies. The code has been specially designed to conserve linear momentum. Despite this, it also has variable softening and an efficient block-time-step scheme: the force between any pair of particles is calculated using the finest mesh that encloses them both (respecting Newton's third law) and is updated only on the longest time-step of the two (which conserves momentum). For realistic galaxy models with   N ≳ 10⁶  , it is faster than the fastest comparable tree code by factors ranging from ∼2 (using single time-steps) to ∼10 (multiple time-steps in a concentrated galaxy).     

General relativistic MHD simulations of black hole accretion disks and jets

Accretion disks orbiting black holes power high-energy systems such as X-ray binaries and Active Galactic Nuclei. Observations are providing increasingly detailed quantitative information about such systems. This data has been interpreted using standard toy-models that rely on simplifying assumptions such as regular flow geometry and a parameterized stress. Global numerical simulations offer a way to investigate the basic physical dynamics of accretion flows without these assumptions and, in principle, lead to a genuinely predictive theory. In recent years we have developed a fully three-dimensional general relativistic magnetohydrodynamic simulation code that evolves time-dependent inflows into Kerr black holes. Although the resulting global simulations of black hole accretion are still somewhat simplified, they have brought to light a number of interesting results. These include the formation of electro-magnetically dominated jets powered by the black hole’s rotation, and the presence of strong stresses in the plunging region of the accretion flow. The observational consequences of these features are gradually being examined. Increasing computer power and increasingly sophisticated algorithms promise a bright future for the computational approach to black hole accretion.     

Graded-Density Reservoirs for Accessing High Stress Low Temperature Material States

In recently developed laser-driven shockless compression experiments an ablatively driven shock in a primary target is transformed into a ramp compression wave in a secondary target via unloading followed by stagnation across an intermediate vacuum gap. Current limitations on the achievable peak longitudinal stresses are limited by the ability of shaping the temporal profile of the ramp compression pulse. We report on new techniques using graded density reservoirs for shaping the loading profile and extending these techniques to high peak pressures.     

Enceladus: Present internal structure and differentiation by early and long-term radiogenic heating

Pre-Cassini images of Saturn's small icy moon Enceladus provided the first indication that this satellite has undergone extensive resurfacing and tectonism. Data returned by the Cassini spacecraft have proven Enceladus to be one of the most geologically dynamic bodies in the Solar System. Given that the diameter of Enceladus is only about 500 km, this is a surprising discovery and has made Enceladus an object of much interest. Determining Enceladus' interior structure is key to understanding its current activity. Here we use the mean density of Enceladus (as determined by the Cassini mission to Saturn), Cassini observations of endogenic activity on Enceladus, and numerical simulations of Enceladus' thermal evolution to infer that this satellite is most likely a differentiated body with a large rock-metal core of radius about 150 to 170 km surrounded by a liquid water-ice shell. With a silicate mass fraction of 50% or more, long-term radiogenic heating alone might melt most of the ice in a homogeneous Enceladus after about 500 Myr assuming an initial accretion temperature of about 200 K, no subsolidus convection of the ice, and either a surface temperature higher than at present or a porous, insulating surface. Short-lived radioactivity, e.g., the decay of ²⁶Al, would melt all of the ice and differentiate Enceladus within a few million years of accretion assuming formation of Enceladus at a propitious time prior to the decay of ²⁶Al. Long-lived radioactivity facilitates tidal heating as a source of energy for differentiation by warming the ice in Enceladus so that tidal deformation can become effective. This could explain the difference between Enceladus and Mimas. Mimas, with only a small rock fraction, has experienced relatively little long-term radiogenic heating; it has remained cold and stiff and less susceptible to tidal heating despite its proximity to Saturn and larger eccentricity than Enceladus. It is shown that the shape of Enceladus is not that of a body in hydrostatic equilibrium at its present orbital location and rotation rate. The present shape could be an equilibrium shape corresponding to a time when Enceladus was closer to Saturn and spinning more rapidly, or more likely, to a time when Enceladus was spinning more rapidly at its present orbital location. A liquid water layer on Enceladus is a possible source for the plume in the south polar region assuming the survivability of such a layer to the present. These results could place Enceladus in a category similar to the large satellites of Jupiter, with the core having a rock-metal composition similar to Io, and with a deep overlying ice shell similar to Europa and Ganymede. Indeed, the moment of inertia factor of a differentiated Enceladus, C/MR², could be as small as that of Ganymede, about 0.31.     

An improved shadow measurement technique for constraining the morphometry of simple impact craters

Abstract— The lengths of the shadows cast within simple, bowl‐shaped impact craters have been used to constrain their depths on a variety of planetary bodies. This technique, however, only yields the “true” crater depth if the shadow transects the crater center where the floor is deepest. In the past, attempts have been made to circumvent this limitation by choosing only craters where the shadow tip lies very near the crater center; but this approach may introduce serious artifacts that adversely affect the slope of the regressed depth vs. diameter data and its variance. Here we introduce an improved method for deriving depth information from shadow measurements that considers three basic shape variations of simple craters: paraboloidal, conical, and flat‐floored. We show that the shape of the cast shadow can be used to constrain crater shape and we derive improved equations for finding the depths of these simple craters.     

Oort Cloud Comet Perihelion Asymmetries: Galactic Tide, Shower or Observational Bias?

We investigate the distribution of Oort cloud comet perihelia. The data considered includes comets having orbital elements of the two highest quality classes with original energies designated as new or young. Perihelion directions are determined in galactic, ecliptic and geocentric equatorial coordinates. Asymmetries are detected in the scatter and are studied statistically for evidence of adiabatic galactic tidal dynamics, an impulse-induced shower and observational bias. The only bias detected is the well-known deficiency of observations with perihelion distances q > 2.5 AU. There is no significant evidence of a seasonal dependence. Nor is there a substantive hemispherical bias in either ecliptic or equatorial coordinates. There is evidence for a weak stellar shower previously detected by Biermann which accounts for ≈ 10% of the total observations. Both the q bias and the Biermann star track serve to weaken the evidence for a galactic tidal imprint. Nevertheless, statistically significant asymmetries in galactic latitude and longitude of perihelia remain. A latitude asymmetry is produced by a dominant tidal component perpendicular to the galactic disk. The longitude signal implies that ≈ 20% of new comets need an additional dynamical mechanism. Known disk non-uniformities and an hypothetical bound perturber are discussed as potential explanations. We conclude that the detected dynamical signature of the galactic tide is real and is not an artifact of observational bias, impulsive showers or poor data. This revised version was published online in July 2006 with corrections to the Cover Date.     

The mean magnetic field of the Sun: Method of observation and relation to the interplanetary magnetic field

The mean solar magnetic field as measured in integrated light has been observed since 1968. Since 1970 it has been observed both at Hale Observatories and at the Crimean Astrophysical Observatory. The observing procedures at both observatories and their implications for mean field measurements are discussed. A comparison of the two sets of daily observations shows that similar results are obtained at both observatories. A comparison of the mean field with the interplanetary magnetic polarity shows that the IMF sector structure has the same pattern as the mean field polarity.     

FOUR NEW IRON METEORITE FINDS

Four new irons are described; Buenaventura (IIIB) from Chihuahua, Mexico: mass 114 kg; Denver City (anomalous) from Texas, USA: mass 26.1 kg; Kinsella (IIIB) from Alberta, Canada: mass 3.7 kg; and Tacoma (IA) from Washington, USA: mass 17 g. Denver City is unique, i.e., not related to any other known iron. Tacoma is the smallest iron meteorite recorded. All were purchased for the UCLA collection following a publicity drive for new meteorites     

DETERMINANTS OF STRUCTURE IN AQUATIC NON-PELAGIC PROTOZOAN COMMUNITIES

Although microbes have traditionally been used as models for testing ecological theory, research on the organization of these communities has largely been isolated from other areas of community ecology. Early studies indicated that microbial populations can form dynamic, interacting assemblages. Observations of distributional patterns suggest that certain deterministic forces regulate community structure. It appears that protozoans are tolerant of a wide range of environmental conditions although evidence for this is largely based on gross environmental analyses. Laboratory studies have suggested the importance of negative and positive biotic interactions in determining community membership, but little field evidence for this exists. The importance of transport processes in controlling community composition is being increasingly recognized. All three types of processes likely act to regulate the colonization and successional dynamics of these communities. A simple model is presented to promote the devel