The Populations of Comet-Like Bodies in the Solar system

A new classification scheme is introduced for comet-like bodies in the Solar system. It covers the traditional comets as well as the Centaurs and Edgeworth-Kuiper belt objects. At low inclinations, close encounters with planets often result in near-constant perihelion or aphelion distances, or in perihelion-aphelion interchanges, so the minor bodies can be labelled according to the planets predominantly controlling them at perihelion and aphelion. For example, a JN object has a perihelion under the control of Jupiter and aphelion under the control of Neptune, and so on. This provides 20 dynamically distinct categories of outer Solar system objects in the Jovian and trans-Jovian regions. The Tisserand parameter with respect to the planet controlling perihelion is also often roughly constant under orbital evolution. So, each category can be further sub-divided according to the Tisserand parameter. The dynamical evolution of comets, however, is dominated not by the planets nearest at perihelion or aphelion, but by the more massive Jupiter. The comets are separated into four categories -- Encke-type, short-period, intermediate and long-period -- according to aphelion distance. The Tisserand parameter categories now roughly correspond to the well-known Jupiter-family comets, transition-types and Halley-types. In this way, the nomenclature for the Centaurs and Edgeworth-Kuiper belt objects is based on, and consistent with, that for comets.Comment: MNRAS, in press, 11 pages, 6 figures (1 available as postscript, 5 as gif). Higher resolution figures available at http://www-thphys.physics.ox.ac.uk/users/WynEvans/preprints.pd

Ruthenium/Iridium Ratios in the Cretaceous-tertiary Boundary Clay: Implications for Global Dispersal and Fractionation Within the Ejecta Cloud

Ruthenium (Ru) and iridium (Ir) are the least mobile platinum group elements (PGE's) within the Cretaceous-Tertiary (K-T) boundary clay (BC). The Ru/Ir ratio is, therefore, the most useful PGE interelement ratio for distinguishing terrestrial and extraterrestrial contributions to the BC. The Ru/Ir ratio of marine K-T sections (1.77 +/- 0.53) is statistically different from that of the continental sections (0.93 +/- 0.28). The marine Ru/Ir ratios are chondritic (C1 = 1.48 +/- 0.09), but the continental ratios are not. We discovered an inverse correlation of shocked quartz size (or distance from the impact site) and Ru/Ir ratio. This correlation may arise from the difference in Ru and Ir vaporization temperature and/or fractionation during condensation from the ejecta cloud. Postsedimentary alteration, remobilization, or terrestrial PGE input may be responsible for the Ru/Ir ratio variations within the groups of marine and continental sites studied. The marine ratios could also be attained if approximately 15 percent of the boundary metals were contributed by Deccan Trap emissions. However, volcanic emissions could not have been the principal source of the PGE's in the BC because mantle PGE ratios and abundances are inconsistent with those measured in the clay. The Ru/Ir values for pristine Tertiary mantle xenoliths (2.6 +/- 0.48), picrites (4.1 +/- 1.8), and Deccan Trap basalt (3.42 +/- 1.96) are all statistically distinct from those measured in the K-T BC

New Velocity Distribution in the Context of the Eddington Theory

Exotic dark matter together with the vacuum energy (associated with the cosmological constant) seem to dominate the Universe. Thus its direct detection is central to particle physics and cosmology. Supersymmetry provides a natural dark matter candidate, the lightest supersymmetric particle (LSP). One essential ingredient in obtaining the direct detection rates is the density and velocity distribution of the LSP. The detection rate is proportional to this density in our vicinity. Furthermore, since this rate is expected to be very low, one should explore the two characteristic signatures of the process, namely the modulation effect, i.e. the dependence of the event rate on the Earth's motion and the correlation of the directional rate with the motion of the sun. Both of these crucially depend on the LSP velocity distribution. In the present paper we study simultaneously density profiles and velocity distributions based on the Eddington theory.Comment: 40 LaTex pages, 19 figures and one table. The previous version was expanded to include new numerical solutions to Poisson's equation. Sheduled to appear in vol. 588, ApJ, May 1, 300

Numerical Evidence for Divergent Burnett Coefficients

In previous papers [Phys. Rev. A {\bf 41}, 4501 (1990), Phys. Rev. E {\bf 18}, 3178 (1993)], simple equilibrium expressions were obtained for nonlinear Burnett coefficients. A preliminary calculation of a 32 particle Lennard-Jones fluid was presented in the previous paper. Now, sufficient resources have become available to address the question of whether nonlinear Burnett coefficients are finite for soft spheres. The hard sphere case is known to have infinite nonlinear Burnett coefficients (ie a nonanalytic constitutive relation) from mode coupling theory. This paper reports a molecular dynamics caclulation of the third order nonlinear Burnett coefficient of a Lennard-Jones fluid undergoing colour flow, which indicates that this term is diverges in the thermodynamic limit.Comment: 12 pages, 9 figure

2D and 3D Dense-Fluid Shear Flows via Nonequilibrium Molecular Dynamics. Comparison of Time-and-Space-Averaged Tensor Temperature and Normal Stresses from Doll's, Sllod, and Boundary-Driven Shear Algorithms

Homogeneous shear flows (with constant strainrate du/dy) are generated with the Doll's and Sllod algorithms and compared to corresponding inhomogeneous boundary-driven flows. We use one-, two-, and three-dimensional smooth-particle weight functions for computing instantaneous spatial averages. The nonlinear stress differences are small, but significant, in both two and three space dimensions. In homogeneous systems the sign and magnitude of the shearplane stress difference, P(xx) - P(yy), depend on both the thermostat type and the chosen shearflow algorithm. The Doll's and Sllod algorithms predict opposite signs for this stress difference, with the Sllod approach definitely wrong, but somewhat closer to the (boundary-driven) truth. Neither of the homogeneous shear algorithms predicts the correct ordering of the kinetic temperatures, T(xx) > T(zz) > T(yy).Comment: 34 pages with 12 figures, under consideration by Physical Review

The beginnings of geography teaching and research in the University of Glasgow: the impact of J.W. Gregory

J.W. Gregory arrived in Glasgow from Melbourne in 1904 to take up the post of foundation Professor of Geology in the University of Glasgow. Soon after his arrival in Glasgow he began to push for the setting up of teaching in Geography in Glasgow, which came to pass in 1909 with the appointment of a Lecturer in Geography. This lecturer was based in the Department of Geology in the University's East Quad. Gregory's active promotion of Geography in the University was matched by his extensive writing in the area, in textbooks, journal articles and popular books. His prodigious output across a wide range of subject areas is variably accepted today, with much of his geomorphological work being judged as misguided to varying degrees. His 'social science' publications - in the areas of race, migration, colonisation and economic development of Africa and Australia - espouse a viewpoint that is unacceptable in the twenty-first century. Nonetheless, that viewpoint sits squarely within the social and economic traditions of Gregory's era, and he was clearly a key 'Establishment' figure in natural and social sciences research in the first half of the twentieth century. The establishment of Geography in the University of Glasgow remains enduring testimony of J.W. Gregory's energy, dedication and foresight

Axisymmetric Three-Integral Models for Galaxies

We describe an improved, practical method for constructing galaxy models that match an arbitrary set of observational constraints, without prior assumptions about the phase-space distribution function (DF). Our method is an extension of Schwarzschild's orbit superposition technique. As in Schwarzschild's original implementation, we compute a representative library of orbits in a given potential. We then project each orbit onto the space of observables, consisting of position on the sky and line-of-sight velocity, while properly taking into account seeing convolution and pixel binning. We find the combination of orbits that produces a dynamical model that best fits the observed photometry and kinematics of the galaxy. A key new element of this work is the ability to predict and match to the data the full line-of-sight velocity profile shapes. A dark component (such as a black hole and/or a dark halo) can easily be included in the models. We have tested our method, by using it to reconstruct the properties of a two-integral model built with independent software. The test model is reproduced satisfactorily, either with the regular orbits, or with the two-integral components. This paper mainly deals with the technical aspects of the method, while applications to the galaxies M32 and NGC 4342 are described elsewhere (van der Marel et al., Cretton & van den Bosch). (abridged)Comment: minor changes, accepted for publication in the Astrophysical Journal Supplement