Low-speed impacts between rubble piles modeled as collections of polyhedra

We present results of modeling rubble piles as collections of polyhedra. The use of polyhedra allows more realistic (irregular) shapes and interactions (e.g. collisions), particularly for objects of different sizes. Rotational degrees of freedom are included in the modeling, which may be important components of the motion. We solved the equations of rigid-body dynamics, including frictional/inelastic collisions, for collections of up to several hundred elements. As a demonstration of the methods and to compare with previous work by other researchers, we simulated low-speed collisions between km-scale bodies with the same general parameters as those simulated by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151]. High-speed collisions appropriate to present-day asteroid encounters require additional treatment of shock effects and fragmentation and are the subject of future work; here we study regimes appropriate to planetesimal accretion and re-accretion in the aftermath of catastrophic events. Collisions between equal-mass objects at low speeds () were simulated for both head-on and off-center collisions between rubble piles made of a power-law mass spectrum of sub-elements. Very low-speed head-on collisions produce single objects from the coalescence of the impactors. For slightly higher speeds, extensive disruption occurs, but re-accretion produces a single object with most of the total mass. For increasingly higher speeds, the re-accreted object has smaller mass, finally resulting in complete catastrophic disruption with all sub-elements on escape trajectories and only small amounts of mass in re-accreted bodies. Off-center collisions at moderately low speeds produce two re-accreted objects of approximately equal mass, separating at greater than escape speed. At high speed, complete disruption occurs as with the high-speed head-on collisions. Head-on collisions at low to moderate speeds result in objects of mostly oblate shape, while higher speed collisions produce mostly prolate objects, as do off-center collisions at moderate and high speeds. Collisions carried out with the same dissipative coefficients (coefficient of restitution ?n=0.8, zero friction) as used by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151] result in a value for specific energy for disruption , somewhat lower than the value of 2 J/kg found by them, while collisions with a lower coefficient of restitution and friction [?n=0.5, ?t=0, μ=0.5, similar to those used by Michel, et al. [Michel, P., Benz, W., Richardson, D.C., 2004. Planet. Space Sci. 52, 1109-1117] for SPH + N-body calculations] yield .     

Mass composition of the escaping plasma at Mars

Data from the Ion Mass Analyzer (IMA) sensor of the ASPERA-3 instrument suite on Mars Express have been analyzed to determine the mass composition of the escaping ion species at Mars. We have examined 77 different ion-beam events and we present the results in terms of flux ratios between the following ion species: CO⁺₂/O⁺ and O⁺₂/O⁺. The following ratios averaged over all events and energies were identified: CO⁺₂/O⁺ = 0.2 and O⁺₂/O⁺ = 0.9. The values measured are significantly higher, by a factor of 10 for O⁺₂/O⁺, than a contemporary modeled ratio for the maximum fluxes which the martian ionosphere can supply. The most abundant ion species was found to be O⁺, followed by O⁺₂ and CO⁺₂. We estimate the loss of CO⁺₂ to be by using the previous measurements of Phobos-2 in our calculations. The dependence of the ion ratios in relation to their energy ranges we studied, 0.3-3.0 keV, indicated that no clear correlation was found.     

Thermal interactions of the Huygens probe with the Titan environment: Constraint on near-surface wind

The Huygens probe lost heat to its cold environment during its descent through Titan's atmosphere and after landing. Here I report measurements of the probe's thermal behavior and comparison with ground tests (1) to provide a context for other scientific investigations, such as the release of volatiles from the landing site, and (2) to place constraints on Titan environmental parameters directly, such as the thermal conductivity of the surface material and the strength of winds at the surface. Near-surface winds are constrained to be less than 0.2 m s⁻¹, and probably much less.     

The peculiar case of the Agnia asteroid family

The Agnia asteroid family, a cluster of asteroids located near semimajor axis a=2.79 AU, has experienced significant dynamical evolution over its lifetime. The family, which was likely created by the breakup of a diameter D∼50 km parent body, is almost entirely contained within the high-order secular resonance z₁. This means that unlike other families, Agnia's full extent in proper eccentricity and inclination is a byproduct of the large-amplitude resonant oscillations produced by this resonance. Using numerical integration methods, we found that the spread in orbital angles observed among Agnia family members would have taken at least 40 Myr to create; this sets a lower limit on the family's age. To determine the upper bound on Agnia's age, we used a Monte Carlo model to track how the small members in the family evolve in semimajor axis by Yarkovsky thermal forces. Our results indicate the family is no more than 140 Myr old, with a best-fit age of 100⁺³⁰₋₂₀ Myr. Using two independent methods, we also determined that the D∼5 km fragments were ejected from the family-forming event at a velocity near 15 m/s. This velocity is consistent with results from numerical hydrocode simulations of asteroid impacts and observations of other similarly sized asteroid families. Finally, we found that 57% of known Agnia fragments were initially prograde rotators. The reason for this limited asymmetry is unknown, though we suspect it is a fluke produced by the stochastic nature of asteroid disruption events.     

Semi-analytic simulations of galactic winds: volume filling factor, ejection of metals and parameter study

We present a semi-analytic treatment of galactic winds within high-resolution, large-scale cosmological N -body simulations of a Λ cold dark matter (ΛCDM) universe. The evolution of winds is investigated by following the expansion of supernova-driven superbubbles around the several hundred thousand galaxies that form in an approximately spherical region of space with diameter 52  h ⁻¹ Mpc and mean density close to the mean density of the universe. We focus our attention on the impact of winds on the diffuse intergalactic medium. Initial conditions for mass loss at the base of winds are taken from Shu, Mo & Mao. Results are presented for the volume filling factor and the mass fraction of the intergalactic medium (IGM) affected by winds, and their dependence on the model parameters is carefully investigated. The mass-loading efficiency of bubbles is a key factor to determine the evolution of winds and their global impact on the IGM: the higher the mass loading, the later the IGM is enriched with metals. Galaxies with 10⁹ < M _★ < 10¹⁰ M_⊙ are responsible for most of the metals ejected into the IGM at   z = 3  , while galaxies with   M _★ < 10⁹ M_⊙   give a non-negligible contribution only at higher redshifts, when larger galaxies have not yet assembled. We find a higher mean IGM metallicity than Lyα forest observations suggest, and we argue that the discrepancy may be explained by the high temperatures of a large fraction of the metals in winds, which may not leave detectable imprints in absorption in the Lyα forest.     

X-ray synchrotron emission from the oblique shock in the jet of the powerful radio galaxy 3C 346

We report the first detection, with Chandra , of X-ray emission from the jet of the powerful narrow-line radio galaxy 3C 346. X-rays are detected from the bright radio and optical knot at which the jet apparently bends by approximately 70°. The Chandra observation also reveals a bright galaxy-scale atmosphere within the previously known cluster and provides a good X-ray spectrum for the bright core of 3C 346. The X-ray emission from the knot is synchrotron radiation, as seen in lower-power sources. In common with these sources, there is evidence of morphological differences between the radio/optical and X-ray structures, and the spectrum is inconsistent with a one-component continuous-injection model. We suggest that the X-ray-bright knot is associated with a strong oblique shock in a moderately relativistic, light jet, at ∼ 20° to the line of sight, and that this shock is caused by the jet interacting with the wake in the cluster medium behind the companion galaxy of 3C 346. The general jet curvature can result from pressure gradients in the cluster atmosphere.     

A statistical study of 233 pulsar proper motions

We present and analyse a catalogue of 233 pulsars with proper motion measurements. The sample contains a wide variety of pulsars including recycled objects and those associated with globular clusters or supernova remnants. After taking the most precise proper motions for those pulsars for which multiple measurements are available, the majority of the proper motions (58 per cent) are derived from pulsar timing methods, 41 per cent using interferometers and the remaining 1 per cent using optical telescopes. Many of the one-dimensional (1D) and two-dimensional (2D) speeds (referring to speeds measured in one coordinate only and the magnitudes of the transverse velocities, respectively) derived from these measurements are somewhat lower than earlier estimates because of the use of the most recent electron density model in determining pulsar distances. The mean 1D speeds for the normal and recycled pulsars are 152(10) and 54(6) km s⁻¹, respectively. The corresponding mean 2D speeds are 246(22) and 87(13) km s⁻¹. PSRs B2011+38 and B2224+64 have the highest inferred 2D speeds of  ∼1600 km s⁻¹  . We study the mean speeds for different subsamples and find that, in general, they agree with previous results. Applying a novel deconvolution technique to the sample of 73 pulsars with characteristic ages less than 3 Myr, we find the mean three-dimensional (3D) pulsar birth velocity to be 400(40) km s⁻¹. The distribution of velocities is well described by a Maxwellian distribution with  1D rms σ= 265 km s⁻¹  . There is no evidence for a bimodal velocity distribution. The proper motions for PSRs B1830−08 and B2334+61 are consistent with their proposed associations with the supernova remnants W41 and G114.3+0.3, respectively.