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 .     

Carbon cycling in a high-arctic marine ecosystem – Young Sound, NE Greenland

Young Sound is a deep-sill fjord in NE Greenland (74°N). Sea ice usually begins to form in late September and gains a thickness of 1.5 m topped with 0–40 cm of snow before breaking up in mid-July the following year. Primary production starts in spring when sea ice algae begin to flourish at the ice–water interface. Most biomass accumulation occurs in the lower parts of the sea ice, but sea ice algae are observed throughout the sea ice matrix. However, sea ice algal primary production in the fjord is low and often contributes only a few percent of the annual phytoplankton production. Following the break-up of ice, the immediate increase in light penetration to the water column causes a steep increase in pelagic primary production. Usually, the bloom lasts until August–September when nutrients begin to limit production in surface waters and sea ice starts to form. The grazer community, dominated by copepods, soon takes advantage of the increased phytoplankton production, and on an annual basis their carbon demand (7–11 g C m⁻²) is similar to phytoplankton production (6–10 g C m⁻²). Furthermore, the carbon demand of pelagic bacteria amounts to 7–12 g C m⁻² yr⁻¹. Thus, the carbon demand of the heterotrophic plankton is approximately twice the estimated pelagic primary production, illustrating the importance of advected carbon from the Greenland Sea and from land in fuelling the ecosystem.In the shallow parts of the fjord (<40 m) benthic primary producers dominate primary production. As a minimum estimate, a total of 41 g C m⁻² yr⁻¹ is fixed by primary production, of which phytoplankton contributes 15%, sea ice algae <1%, benthic macrophytes 62% and benthic microphytes 22%. A high and diverse benthic infauna dominated by polychaetes and bivalves exists in these shallow-water sediments (<40 m), which are colonized by benthic primary producers and in direct contact with the pelagic phytoplankton bloom. The annual benthic mineralization is 32 g C m⁻² yr⁻¹ of which megafauna accounts for 17%. In deeper waters benthic mineralization is 40% lower than in shallow waters and megafauna, primarily brittle stars, accounts for 27% of the benthic mineralization. The carbon that escapes degradation is permanently accumulated in the sediment, and for the locality investigated a rate of 7 g C m⁻² yr⁻¹ was determined.A group of walruses (up to 50 adult males) feed in the area in shallow waters (<40 m) during the short, productive, ice-free period, and they have been shown to be able to consume <3% of the standing stock of bivalves (Hiatella arctica, Mya truncata and Serripes Groenlandicus), or half of the annual bivalve somatic production. Feeding at greater depths is negligible in comparison with their feeding in the bivalve-rich shallow waters.     

Measurement of Impact Ejecta from Regolith Targets in Oblique Impacts

This paper reports the results of experiments on projectile impact into regolith targets at various impact angles. Copper projectiles of 240 mg are accelerated to 197 to 272 m s⁻¹ using an electromagnetic gun. The ejecta are detected by thin Al foil targets as secondary targets, and the resulting holes on the foil are measured to derive the spatial distribution of the ejecta. The ejecta that penetrated the foil are concentrated toward the downrange azimuths of impacting projectiles in oblique impacts. In order to investigate the ejecta velocity distribution, the nondimensional volume of ejecta with velocities higher than a given value is calculated from the spatial distribution. In the case of the vertical impact of the projectile, most ejecta have velocities lower than 24% of the projectile speed (∼50 m s⁻¹), and there are only several ejecta with velocities higher than 72 m s⁻¹. This result confirms the existence of an upper limit to the ejection velocity in the ejecta velocity distribution (Hartmann cutoff velocity) (W. K. Hartmann, 1985, Icarus63, 69-98). On the other hand, it is found that, in the oblique impacts, there are a large number of ejecta with velocities higher than the Hartmann cutoff velocity. The relative quantity of ejecta above the Hartmann cutoff velocity increases as the projectile impact angle decreases. Taking these results with the results of S. Yamamoto and A. M. Nakamura (1997, Icarus128, 160-170) from impact experiments using an impact angle of 30°, it can be concluded that the ejecta from these regolith targets exhibit a bimodal velocity distribution. Below a few tens of m s⁻¹, we see the expected velocity distribution of ejecta, but above this velocity we see a separate group of high-velocity ejecta.     

Electron-impact rotational excitation of linear molecular ions

Molecular R -matrix calculations are performed to give rotational excitation rates for electron collisions with linear molecular ions. Results are presented for CO⁺, HCO⁺, NO⁺ and H₂⁺ up to electron temperatures of 10 000 K. De-excitation rates and critical electron densities are also given. It is shown that the widely used Coulomb–Born approximation is valid for Δ j =1 transitions when the molecular ion has a dipole greater than about 2D, but otherwise is not reliable for studying electron-impact rotational excitation. In particular, transitions with Δ j >1 are found to have appreciable rates and are found to be entirely dominated by short-range effects.     

Quantum peculiarities of the primordial hydrogen recombination

Some problems of the two-photon decay of highly excited atomic hydrogen and helium levels are considered. The necessary accuracy of calculating the transition probabilities is shown to require abandoning the purely dipole approximation in the theory of atomic hydrogen radiation. The range of applicability of the perturbation theory to this process is discussed. We suggest a new mechanism of the Lyman photon redistribution due to very close locations of the 20P and 19P levels in the hydrogen and deuterium atoms, respectively (hydrogen-deuterium resonance).