Performance tuning of N-body codes on modern microprocessors: I. Direct integration with a hermite scheme on x86_64 architecture

The main performance bottleneck of gravitational N-body codes is the force calculation between two particles. We have succeeded in speeding up this pair-wise force calculation by factors between 2 and 10, depending on the code and the processor on which the code is run. These speed-ups were obtained by writing highly fine-tuned code for x86_64 microprocessors. Any existing N-body code, running on these chips, can easily incorporate our assembly code programs.In the current paper, we present an outline of our overall approach, which we illustrate with one specific example: the use of a Hermite scheme for a direct N² type integration on a single 2.0 GHz Athlon 64 processor, for which we obtain an effective performance of 4.05 Gflops, for double-precision accuracy. In subsequent papers, we will discuss other variations, including the combinations of N log N codes, single-precision implementations, and performance on other microprocessors.     

Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

Impact-induced N2 production from ammonium sulfate: Implications for the origin and evolution of N2 in Titan’s atmosphere

Chemical reactions and volatile supply through hypervelocity impacts may have played a key role for the origin and evolution of both planetary and satellite atmospheres. In this study, we evaluate the role of impact-induced N₂ production from reduced nitrogen-bearing solids proposed to be contained in Titan’s crust, ammonium sulfate ((NH₄)₂SO₄), for the replenishment of N₂ to the atmosphere in Titan’s history. To investigate the conversion of (NH₄)₂SO₄ into N₂ by hypervelocity impacts, we measured gases released from (NH₄)₂SO₄ that was exposed to hypervelocity impacts created by a laser gun. The sensitivity and accuracy of the measurements were enhanced by using an isotope labeling technique for the target. We obtained the efficiency of N₂ production from (NH₄)₂SO₄ as a function of peak shock pressure ranging from ∼8 to ∼45 GPa. Our results indicate that the initial and complete shock pressures for N₂ degassing from (NH₄)₂SO₄ are ∼10 and ∼25 GPa, respectively. These results suggest that cometary impacts on Titan (i.e., impact velocity v_i > ∼8 km/s) produce N₂ efficiently; whereas satellitesimal impacts during the accretion (i.e., v_i < 4 km/s) produce N₂ only inefficiently. Even when using the proposed small amount of (NH₄)₂SO₄ content in the crust (∼4 wt.%) (Fortes, A.D. et al., 2007. Icarus 188, 139-153), the total amount of N₂ provided through cometary impacts over 4.5 Ga reaches ∼2-6 times the present atmospheric N₂ (i.e., ∼7 × 10²⁰-2 × 10²¹ [mol]) based on the measured production efficiency and results of a hydrodynamic simulation of cometary impacts onto Titan. This implies that cometary impacts onto Titan’s crust have the potential to account for a large part of the present N₂ through the atmospheric replenishment after the accretion.     

N-Body simulations of growth from 1 km planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M_⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10⁵ particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.     

The Unusual Absorption Line Spectrum of Quasar SDSS J2220+0109

We report the Balmer broad absorption lines (BALs) in the quasar SDSS J2220 + 0109 discovered from the SDSS data, and present a detailed analysis of the peculiar absorption line spectrum, including the He I* multiplet at λλ3189, 3889 arising from the metastable 2³s-state helium and the Balmer Hα and Hβ lines from the excited hydrogen H I of n = 2 level, which are rarely seen in quasar spectra, as well as many absorption lines arising from the excited Fe II* of the levels 7 955 cm⁻¹, 13 474 cm⁻¹ and 13 673 cm⁻¹ in the wavelength range 3100∼3300 Å. Ca II H, K absorption line doublets also clearly appear in the SDSS spectrum. All absorption lines show a similar blueshifted velocity structure of Δv ≈ − 1500 ∼ 0 km·s⁻¹ relative to the quasar's systematic redshift determined from the emission lines. Detailed analysis suggests that the Balmer absorption lines should arise from the partially ionized region with a column density of N_HI ≈ 10²¹ cm⁻² for an electron density of ne ∼ 10⁶ cm⁻³; and that the hydrogen n = 2 level may be populated via collisional excitation with Lyα pumping.     

NN detectability in Pluto’s atmosphere

Based on the vapor pressure behavior of Pluto’s surface ices, Pluto’s atmosphere is expected to be predominantly composed of N₂ gas. Measurement of the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere would provide important clues to the evolution of Pluto’s atmosphere from the time of formation to its present state. The most straightforward way of determining the N₂ isotopologue ¹⁵N/¹⁴N ratio in Pluto’s atmosphere is via spectroscopic observation of the ¹⁴N¹⁵N gas species. Recent calculations of the 80–100 nm absorption behavior of the ¹⁴N₂ and ¹⁴N¹⁵N isotopologues by Heays et al. (Heays, A.N. et al. [2011]. J. Chem. Phys. 135, 244301), Lewis et al. (Lewis, B.R., Heays, A.N., Gibson, S.T., Lefebvre-Brion, H., Lefebvre, R. [2008]. J. Chem. Phys. 129, 164306); Lewis et al. (Lewis, B.R., Gibson, S.T., Zhang, W., Lefebvre-Brion, H., Robbe, J.-M. [2005]. J. Chem. Phys. 122, 144302), and Haverd et al. (Haverd, V.E., Lewis, B.R., Gibson, S.T., Stark, G. [2005]. J. Chem. Phys. 123, 214304) show that the peak magnitudes of the ¹⁴N₂ and ¹⁴N¹⁵N absorption bandhead cross-sections are similar, but the locations of the bandhead peaks are offset in wavelength by ∼0.05–0.1 nm. These offsets make the segregation of the ¹⁴N₂ and ¹⁴N¹⁵N absorption signatures possible. We use the most recent N₂ isotopologue absorption cross-section calculations and the atmospheric density profiles resulting from photochemical models developed by Krasnopolsky and Cruickshank (Krasnopolsky, V.A., Cruickshank, D.P. [1999]. J. Geophys. Res. 104, 21979–21996) to predict the level of solar light that will be transmitted through Pluto’s atmosphere as a function of altitude during a Pluto solar occultation. We characterize the detectability of the isotopic absorption signature per altitude assuming ¹⁴N¹⁵N concentrations ranging from 0.1% to 2% of the ¹⁴N₂ density and instrumental spectral resolutions ranging from 0.01 to 0.3 nm. Our simulations indicate that optical depth of unity is attained in the key ¹⁴N¹⁵N absorption bands located between 85 and 90 nm at altitudes ∼1100–1600 km above Pluto’s surface. Additionally, an ¹⁴N¹⁵N isotope absorption depth ∼4–15% is predicted for observations obtained at these altitudes at a spectral resolution of ∼0.2–0.3 nm, if the N₂ isotopologue ¹⁵N/¹⁴N percent ratio is comparable to the 0.37–0.6% ratio observed at Earth, Titan and Mars. If we presume that the predicted absorption depth must be at least 25% greater than the expected observational uncertainty, then it follows that a statistically significant detection of these signatures and constraint of the N₂ isotopologue ¹⁴N/¹⁵N ratio within Pluto’s atmosphere will be possible if the attainable observational signal-to noise (S/N) ratio is ?9. The New Horizons (NH) Mission will be able to obtain high S/N, 0.27–0.35 nm full-width half-max 80–100 nm spectral observations of Pluto using the Alice spectrograph. Based on the NH/Alice specifications we have simulated 0.3 nm spectral resolution solar occultation spectra for the 1100–1600 km altitude range, assuming 30 s integration times. These simulations indicate that NH/Alice will obtain spectral observations within this altitude range with a S/N ratio ∼25–50, and should be able to reliably detect the ¹⁴N¹⁵N gas absorption signature between 85 and 90 nm if the ¹⁴N¹⁵N concentration is ∼0.3% or greater. This, additionally, implies that the non-detection of the ¹⁴N¹⁵N species in the 1100–1600 km range by NH/Alice may be used to reliably establish an upper limit to the N₂ isotopologue ¹⁵N/¹⁴N ratio within Pluto’s atmosphere. Similar results may be derived from 0.2 to 0.3 nm spectral resolution observations of any other N₂-rich Solar System or exoplanet atmosphere, provided the observations are attained with similar S/N levels.     

A momentum-conserving N-body scheme with individual time steps

We present an N-body code called Taichi for galactic dynamics and controlled numerical experiments. The code includes two high-order hierarchical multipole expansion methods: the Barnes-Hut (BH) tree and the fast multipole method (FMM). For the time integration, the code can use either a conventional adaptive KDK or a Hamiltonian splitting integrator. The combination of FMM and the Hamiltonian splitting integrator leads to a momentum-conserving N-body scheme with individual time steps. We find Taichi performs well in the typical applications in galactic dynamics. In the isolated and interacting galaxies tests, the momentum conserving scheme produces the same result as a conventional BH tree code. But for similar force accuracies, FMM significantly speeds up the simulations compared to the monopole BH tree. In the cold collapse test, we find the inner structure after relaxation can be sensitive to the force accuracies. Taichi is ready to incorporate special treatment of close encounters thanks to the Hamiltonian splitting integrator, suitable for studying dynamics around central massive bodies.     

Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

Revised vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ Northern Spring Equinox from 2006 to 2008: Application of new methane absorption data and comparison with Neptune

Long-slit spectroscopy observations of Uranus by the United Kingdom InfraRed Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s Northern Spring Equinox in December 2007.These spectra were analysed and presented by Irwin et al. (Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2009]. Icarus 203, 287-302), but since publication, a new set of methane absorption data has become available (Karkoschka, E., Tomasko, M. [2010]. Methane absorption coefficients for the jovian planets from laboratory, Huygens, and HST data. Icarus 205, 674-694.), which appears to be more reliable at the cold temperatures and high pressures of Uranus’ deep atmosphere. We have fitted k-coefficients to these new methane absorption data and we find that although the latitudinal variation and inter-annual changes reported by Irwin et al. (2009) stand, the new k-data place the main cloud deck at lower pressures (2-3 bars) than derived previously in the H-band of ∼3-4 bars and ∼3 bars compared with ∼6 bars in the J-band. Indeed, we find that using the new k-data it is possible to reproduce satisfactorily the entire observed centre-of-disc Uranus spectrum from 1 to 1.75 μm with a single cloud at 2-3 bars provided that we make the particles more back-scattering at wavelengths less than 1.2 μm by, for example, increasing the assumed single-scattering albedo from 0.75 (assumed in the J and H-bands) to near 1.0. In addition, we find that using a deep methane mole fraction of 4% in combination with the associated warm ‘F’ temperature profile of Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92, 14987-15001), the retrieved cloud deck using the new (Karkoschka and Tomasko, 2010) methane absorption data moves to between 1 and 2 bars.The same methane absorption data and retrieval algorithm were applied to observations of Neptune made during the same programme and we find that we can again fit the entire 1-1.75 μm centre-of-disc spectrum with a single cloud model, providing that we make the stratospheric haze particles (of much greater opacity than for Uranus) conservatively scattering (i.e. ω = 1) and we also make the deeper cloud particles, again at around the 2 bar level more reflective for wavelengths less than 1.2 μm. Hence, apart from the increased opacity of stratospheric hazes in Neptune’s atmosphere, the deeper cloud structure and cloud composition of Uranus and Neptune would appear to be very similar.     

Carbon monoxide fluorescence from Titan's atmosphere

We report on the discovery of emissions due to carbon monoxide from Titan's atmosphere, from mid-infrared observations with the ISAAC spectrometer at the Very Large Telescope and covering the 4.50-4.85 μm range. We detected about 45 emission lines coinciding with CO ro-vibrational lines, including CO(1-0) (P18 to R11) and CO(2-1) (P11 to R11). We show that these emissions cannot be generated thermally but occur in non-LTE conditions, due to radiative de-excitation from the v=1 and v=2 CO levels after excitation at 4.7 and 2.3 μm by solar radiation. A complete fluorescence model is then developed, allowing to compute the state populations of the two most abundant CO isotopes and N₂(1). It includes absorption by CO and CH₄, and vibrational-thermal and vibrational-vibrational collisional exchanges with CO, N₂, CH₄, and H₂. Emerging radiances at the top of the atmosphere are evaluated with a line-by-line code and compared to observations. Contribution functions show that the CO emissions sound Titan's stratosphere: while the (1-0) lines generally probe two layers, located respectively at 100-250 km and 300-550 km, the (2-1) lines are sensitive to the intermediate layer at 150-300 km. A sensitivity study is performed to establish the effect of the main model parameters (temperature profile, collisional scenario, and CO stratospheric abundance) on the results. Models reproduce the essential structure of the observed emissions. The (1-0) fundamental band is generally well fit with a nominal CO mixing ratio of 32 ppm—as inferred in the troposphere from observations at 4.80-5.10 μm (Lellouch et al., 2003, Icarus 162, 126-143). However, this band is only weakly dependent on the CO abundance, and given temperature and collisional scenario uncertainties, it constrains the CO stratospheric mixing ratio only to within a factor of ∼3. In addition, the nominal model with 32 ppm CO underestimates the first hot (2-1) transition by approximately a factor of 2. This discrepancy can be resolved by a combined adjustment of collisional rates and an increased CO stratospheric ratio of 60 ppm, consistent with the determination of Gurwell and Muhleman (2000, Icarus 145, 653-656). In contrast, the CO vertical profile suggested by Hidayat et al. (1998, Icarus 133, 109-133), strongly depleted above 200 km, cannot match the data for any realistic collisional scenario, and is therefore not supported by our results.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

Turbulent viscosity and lifetime of Saturn's rings

The viscosity (the angular momentum flux) in the disk of mutually gravitating particles of Saturn's rings is investigated. The hydrodynamic theory of the gravitational Jeans-type instability of small gravity perturbations (e.g., those produced by spontaneous disturbances) of the disk is developed. It is suggested that in such a system the hydrodynamic turbulence may arise as a result of the instability. The turbulence is related to stochastic motions of “fluid” elements. The objective of this paper is to show that in the Jeans-unstable Saturnian ring disk the turbulent viscosity may exceed the ordinary microscopic viscosity substantially. The main result of local N-body simulations of planetary rings by Daisaka et al. (2001. Viscosity in a dense planetary ring with self-gravitating particles. Icarus 154, 296-312) is explained: in the presence of the gravitationally unstable density waves, the effective turbulent viscosity ν_eff is given as ν_eff=CG²Σ²/Ω³, where G, Σ, and Ω are the gravitational constant, the surface mass density of a ring, and the angular velocity, respectively, and the nondimensional correction factor C≈10. We argue that both Saturn's main rings and their irregular of the order of 100 m or even less fine-scale structure (being recurrently created and destroyed on the time scale of an order of Keplerian period ) are not likely much younger than the solar system.     

Accretion among preplanetary bodies: The many faces of runaway growth

When preplanetary bodies reach proportions of ∼1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. The numerical model is novel as it attempts to preserve the individual particle nature of the bodies (like N-body codes); yet it is statistical in nature since it must incorporate the very large number of planetesimals. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R₀. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We associate the initial runaway growth phase with increasing GF-factors for the largest body. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high-mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two-component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. We identify the peak of GF with the transition between the runaway growth and oligarchy accretion stages. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius, especially at the outer disk. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion, which could lead to a very fast growth.     

Analysis of the Hungaria asteroid population

There are approximately 5000 known asteroids in the Hungaria orbital space, a region defined by orbits with high inclination (16° < i < 34°), low eccentricities (e < 0.18), and semi-major axes 1.78 < a < 2.0 AU. We argue that this region is populated by a large number of asteroids formed after a catastrophic collision involving (434) Hungaria, the presumptive largest fragment of the Hungaria collisional family. The remaining objects form a background population that share orbital characteristics with the family members. Due to the general dynamic stability of the region, it is likely that most asteroids in Hungaria space (the Hungaria “group”) have been in this region since the formation of the Solar System or at least since the planets assumed their current orbital configuration. Our examination of the Hungaria group included comparing rotation rates, taxonomic classification, and orbital dynamics to determine the characteristics of the family and background populations. We first found there is an excess of slow rotators among the group but, otherwise, the distribution of spin frequencies is essentially uniform, i.e., that a plot of the cumulative number of objects over the range of 1 d⁻¹ < f < 9 d⁻¹ is nearly a straight line or, put another way, if the distribution over the range is binned by equal intervals of f (1-2 d⁻¹, 2-3 d⁻¹, etc.), the number of objects in each bin is statistically the same.There is a distinct family within the Hungaria group, centered at a semi-major axis of 1.940 AU, with a dispersion range that increases with decreasing size of members, as expected of an evolved collisional family. The larger members with well-determined taxonomic class, including (434) Hungaria itself, have flat spectra, mostly likely type E or similar. The degree of spreading versus size of family members is consistent with that expected from Yarkovsky thermal drift in roughly 0.5 Gyr, suggesting that age for the family. The Asteroid (434) Hungaria is displaced in semi-major axis by 0.004 AU from the center of the Hungaria family. The collision event that produced the family should not have left the largest body displaced by more than 0.001 AU from the original orbit, thus we infer that the displacement of (434) Hungaria is mainly due to Yarkovsky drift, and is consistent with the expected drift for that size body in ∼0.5 Gyr. Below ∼1.93 AU heliocentric distance the Hungaria family is perturbed by at least two secular resonances, 2g − g₅ − g₆ and one of the family of 4th or 6th order secular resonances near s ∼ −22.25 ″/year. Their combined effect results in larger inclination dispersion of the family members.     

Using the youngest asteroid clusters to constrain the space weathering and gardening rate on S-complex asteroids

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10⁶ years and to newly identified asteroid pairs with ages <5 × 10⁵ years. We have identified three S-complex asteroid clusters amongst the set of clusters with ages in the range 10^5-6 years—(1270) Datura, (21509) Lucascavin and (16598) 1992 YC2. The average color of the objects in these clusters agrees with the color predicted by the space weathering model of Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). SDSS five-filter photometry of the members of the very young asteroid pairs with ages <10⁵ years was used to determine their taxonomic classification. Their types are consistent with the background population near each object. The average color of the S-complex pairs is PC₁ = 0.49 ± 0.03, over 5σ redder than predicted by Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). This may indicate that the most likely pair formation mechanism is a gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. If this is the case then our color measurement allows us to set an upper limit of ∼64% on the amount of surface disturbed in the separation process. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic time scales of τ_w = 960 ± 160 Myr and τ_g = 2000 ± 290 Myr respectively. The first principal component color for fresh S-complex material is PC₁ = 0.37 ± 0.01 while the maximum amount of local reddening is ΔPC₁ = 0.33 ± 0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of τ_g ∼ 270 Myr based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a ‘honeycomb’ mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.     

Clouds, haze, and CH4, CH3D, HCN, and C2H2 in the atmosphere of Titan probed via 3 μm spectroscopy

Using synthetic spectra derived from an updated model atmosphere together with a continuum model that includes contributions from haze, cloud and ground, we have re-analyzed the recently published (Geballe et al., 2003, Astrophys. J. 583, L39-L42) high-resolution 3 μm spectrum of Titan which contains newly-detected bands of HCN (in emission) and C₂H₂ and CH₃D (in absorption), in addition to previously detected bands of CH₄. In the 3.10-3.54 μm interval the analysis yields strong evidence for the existence of a cloud deck or optically thick haze layer at about the 10 mbar (∼ 100 km) level. The haze must extend well above this altitude in order to mask the strong CH₄ lines at 3.20-3.50 μm. These cloud and haze components must be transparent at 2.87-2.92 μm, where analysis of the CH₃D spectrum demonstrates that Titan's surface is glimpsed through a second cloud deck at about the 100 mbar (∼ 50 km) level. Through a combination of areal distribution and optical depth this cloud deck has an effective transmittance of ∼ 20%. The spectral shape of Titan's continuum indicates that the higher altitude cloud and haze particles responsible for suppressing the CH₄ absorptions have a largely organic make-up. The rotational temperature of the HCN ranges from 140 to 180 K, indicating that the HCN emission occurs over a wide range of altitudes. This emission, remodeled using an improved collisional deactivation rate, implies mesospheric mixing ratio curves that are consistent with previously predictions. The stratospheric and mesospheric C₂H₂ mixing ratios are ∼¹⁰⁻⁵, considerably less than previous model predictions (Yung et al., 1984), but approximately consistent with recent observational results. Upper limits to mixing ratios of HC₃N and C₄H₂ are derived from non-detections of those species near 3.0 μm.     

Interplanetary survival probability of Aspergillus terreus spores under simulated solar vacuum ultraviolet irradiation

This work is a part of ESA/EU SURE project aiming to quantify the survival probability of fungal spores in space under solar irradiation in the vacuum ultraviolet (VUV) (110-180 nm) spectral region. The contribution and impact of VUV photons, vacuum, low temperature and their synergies on the survival probability of Aspergillus terreus spores is measured at simulated space conditions on Earth. To simulate the solar VUV irradiation, the spores are irradiated with a continuous discharge VUV hydrogen photon source and a molecular fluorine laser, at low and high photon intensities at 10¹⁵ photon m⁻² s⁻¹ and 3.9×10²⁷ photons pulse⁻¹ m⁻² s⁻¹, respectively. The survival probability of spores is independent from the intensity and the fluence of photons, within certain limits, in agreement with previous studies. The spores are shielded from a thin carbon layer, which is formed quickly on the external surface of the proteinaceous membrane at higher photon intensities at the start of the VUV irradiation. Extrapolating the results in space conditions, for an interplanetary direct transfer orbit from Mars to Earth, the spores will be irradiated with 3.3×10²¹ solar VUV photons m⁻². This photon fluence is equivalent to the irradiation of spores on Earth with 54 laser pulses with an experimental ∼92% survival probability, disregarding the contribution of space vacuum and low temperature, or to continuous solar VUV irradiation for 38 days in space near the Earth with an extrapolated ∼61% survival probability. The experimental results indicate that the damage of spores is mainly from the dehydration stress in vacuum. The high survival probability after 4 days in vacuum (∼34%) is due to the exudation of proteins on the external membrane, thus preventing further dehydration of spores. In addition, the survival probability is increasing to ∼54% at 10 K with 0.12 K/s cooling and heating rates.     

Radar and optical observations and physical modeling of triple near-Earth Asteroid (136617) 1994 CC

We report radar, photometric, and spectroscopic observations of near-Earth Asteroid (136617) 1994 CC. The radar measurements were obtained at Goldstone (8560 MHz, 3.5 cm) and Arecibo (2380 MHz, 12.6 cm) on 9 days following the asteroid’s approach within 0.0168 AU on June 10, 2009. 1994 CC was also observed with the Panchromatic Robotic Optical Monitoring and Polarimetry Telescopes (PROMPT) on May 21 and June 1-3. Visible-wavelength spectroscopy was obtained with the 5-m Hale telescope at Palomar on August 25. Delay-Doppler radar images reveal that 1994 CC is a triple system; along with (153591) 2001 SN263, this is only the second confirmed triple in the near-Earth population. Photometry obtained with PROMPT yields a rotation period for the primary P = 2.38860 ± 0.00009 h and a lightcurve amplitude of ∼0.1 mag suggesting a shape with low elongation. Hale telescope spectroscopy indicates that 1994 CC is an Sq-class object. Delay-Doppler radar images and shape modeling reveal that the primary has an effective diameter of 0.62 ± 0.06 km, low pole-on elongation, few obvious surface features, and a prominent equatorial ridge and sloped hemispheres that closely resemble those seen on the primary of binary near-Earth Asteroid (66391) 1999 KW4. Detailed orbit fitting reported separately by Fang et al. (Fang, J., Margot, J.-L., Brozovic, M., Nolan, M.C., Benner, L.A.M., Taylor, P.A. [2011]. Astron. J. 141, 154-168) gives a mass of the primary of 2.6 × 10¹¹ kg that, coupled with the effective diameter, yields a bulk density of 2.1 ± 0.6 g cm⁻³. The images constrain the diameters of the inner and outer satellites to be 113 ± 30 m and 80 ± 30 m, respectively. The inner satellite has a semimajor axis of ∼1.7 km (∼5.5 primary radii), an orbital period of ∼30 h, and its Doppler dispersion suggests relatively slow rotation, 26 ± 12 h, consistent with spin-orbit lock. The outer satellite has an orbital period of ∼9 days and a rotation period of 14 ± 7 h, establishing that the rotation is not spin-orbit locked. Among all binary and triple systems observed by radar, at least 25% (7/28) have a satellite that rotates more rapidly than its orbital period. This suggests that asynchronous configurations with P_rotation < P_orbital are relatively common among multiple systems in the near-Earth population. 1994 CC’s outer satellite has an observed maximum separation from the primary of ∼5.7 km (∼18.4 primary radii) that is the largest separation relative to primary radius seen to date among all 36 known binary and triple NEA systems. 1994 CC, (153591) 2001 SN263, and 1998 ST27 are the only triple and binary systems known with satellite separations >10 primary radii, suggesting either a detection bias, or that such widely-separated satellites are relatively uncommon in NEA multiple systems.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.