Computing the effects of YORP on the spin rate distribution of the NEO population

The overall change of NEO spin rate due to planetary encounters and YORP is evaluated by using a Monte Carlo model. A large sample of test objects mimicking a source population is evolved over a timescale comparable with the Solar System age until they reach a steady state spin distribution that should reproduce the current NEO distribution. The spin change due to YORP is computed for each body according to a simplified model based on Scheeres [Scheeres, D.J., 2007a. Icarus 188, 430-450].The steady state cumulative distribution of NEO spin rates obtained from our simulation nicely reproduces the observed one, once our results are biased to match the diameter distribution of the sample of objects included in the observational database. The excellent agreement strongly suggests that YORP is responsible for the concentration of spin at low rotation rates. In fact, in the absence of YORP the steady state population significantly deviates from the observed one. The spin evolution due to YORP is also so rapid for NEOs that the initial rotation rate distribution of any source population is quickly relaxed to that of the observed population. This has profound consequences for the study of NEO origin since we cannot trace the sources of NEOs from their rotation rate only.     

Considerations on the magnitude distributions of the Kuiper belt and of the Jupiter Trojans

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the ‘Nice model’. This leads us to predict that the HD of the hot population should have the same slope of the HD of the cold population for 6.5<H<9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H∼9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. [Bernstein, G.M., and 5 colleagues, 2004. Astron. J. 128, 1364-1390] and Fuentes and Holman [Fuentes, C.I., Holman, M.J., 2008. Astron. J. 136, 83-97]. Finally, we show that the constraint on the total mass of the primordial trans-neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.     

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Trace constituents of Europa's atmosphere

Europa is bombarded by intense radiation that erodes the surface, launching molecules into a thin “atmosphere” representative of surface composition. In addition to atoms and molecules created in the mostly water ice surface such as H₂O, O₂, H₂, the atmosphere is known to have species representative of trace surface materials. These trace species are carried off with the 10-10⁴ H₂O molecules ejected by each energetic heavy ion, a process we have simulated using molecular dynamics. Using the results of those simulations, we found that a neutral mass spectrometer orbiting ∼100 km above the surface could detect species with surface concentrations above ∼0.03%. We have also modeled the atmospheric spatial structure of the volatile species CO₂ and SO₂ under a variety of assumptions. Detections of these species with moderate time and space resolution would allow us to constrain surface composition, chemistry and to study space weathering processes.     

The interior structure of Mercury and its core sulfur content

Depth-dependent interior structure models of Mercury are calculated for several plausible chemical compositions of the core and of the mantle. For those models, we compute the associated libration amplitude, obliquity, tidal deformation, and tidal changes in the external potential. In particular we study the relation between the interior structure parameters for five different mantle mineralogies and two different temperature profiles together with two extreme crust density values. We investigate the influence of the core light element concentration, temperature, and melting law on core state and inner and outer core size. We show that a sulfur concentration above 10 wt% is unlikely if the temperature at the core-mantle boundary is above 1850 K and the silicate shell at least 240 km thick. The interior models can only have an inner core if the sulfur weight fraction is below 5 wt% for core-mantle boundary temperature in the 1850-2200 K range. Within our modeling hypotheses, we show that with the expected precision on the moment of inertia the core size can be estimated to a precision of about 50 km and the core sulfur concentration with an error of about 2 wt%. This uncertainty can only be reduced when more information on the mantle mineralogy of Mercury becomes available. However, we show that the uncertainty on the core size estimation can be greatly reduced, to about 25 km, if tidal surface displacements and tidal variations in the external potential are considered.     

3D triangulation of a Sun-grazing comet

The bright Kreutz Comet C/2007 L3 (SOHO) entered the fields of view of the twin Solar Terrestrial Relations Observatory (STEREO) COR1 telescopes on 7–8 June 2007. The 12° separation between the two spacecraft at the time afforded the opportunity to derive the position of the comet's tail in three-dimensional space using direct triangulation. The track of the comet's orbit is compared against more traditional orbital calculations using observations from the STEREO COR2 telescopes, and from the Large Angle and Spectrometric Coronagraph (LASCO) aboard the Solar and Heliospheric Observatory (SOHO). The shape of the comet's tail shows that it is composed of dust particles released when the comet was between 18 and 22 solar radii, with no significant dust production after that. The comet did not survive perihelion passage, but a rare faint remnant of the comet tail persisted for several hours after the break-up, and was seen by both the SOHO and STEREO coronagraphs to drift slowly away from the Sun. This tail remnant was found to be composed of particles far back from the head of the comet. The motion of the tail remnant shows a loss of angular momentum during the passage through the solar corona. Atmospheric drag is estimated to account for a significant fraction of this change in angular momentum, but indications are that other mechanisms may be required to completely account for the total amount of change.     

iVINE – Ionization in the parallel tree/sph code VINE: first results on the observed age-spread around O-stars

We present a three-dimensional, fully parallelized, efficient implementation of ionizing ultraviolet (UV) radiation for smoothed particle hydrodynamics ( sph ) including self-gravity. Our method is based on the sph / tree code vine . We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from, for example, young massive stars on their turbulent parental molecular clouds. In this paper, we describe the concept and the numerical implementation of the radiative transfer for a plane-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several 10⁵ yr earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of protostellar envelopes by ionizing radiation.     

Crash-testing the cauldron code for joint lensing and dynamics analysis of early-type galaxies

We apply the joint lensing and dynamics code for the analysis of early-type galaxies, 'Combined Algorithm for Unified Lensing and Dynamics ReconstructiON ( cauldron )', to a rotating N -body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial symmetry supported by a two-integral distribution function. The goal is to study how cauldron performs in an extreme case, and to determine which galaxy properties can still be robustly recovered. Three data sets, corresponding to orthogonal lines of sight, are generated from the N -body system and analysed with the identical procedure followed in the study of real lens galaxies, adopting an axisymmetric power-law total density distribution. We find that several global properties of the N -body system are recovered with remarkable accuracy, despite the fact that the adopted power-law model is too simple to account for the lack of symmetry of the true density distribution. In particular, the logarithmic slope of the total density distribution is robustly recovered to within less than 10 per cent (with the exception of the ill-constrained very inner regions), the inferred angle-averaged radial profile of the total mass closely follows the true distribution, and the dark matter fraction of the system (inside the effective radius) is correctly determined within ∼10 per cent of the total mass. Unless the line-of-sight direction is almost parallel to the total angular momentum vector of the system, reliably recovered quantities also include the angular momentum, the   V /σ  ratio and the anisotropy parameter δ. We conclude that the cauldron code can be safely and effectively applied to real early-type lens galaxies, also providing reliable information for the systems that depart significantly from the method's assumptions.     

High-order 3D Voronoi tessellation for identifying isolated galaxies, pairs and triplets

A geometric method based on the high-order 3D Voronoi tessellation is proposed for identifying single galaxies, pairs and triplets. This approach allows us to select small galaxy groups and isolated galaxies in different environments and to find the isolated systems. The volume-limited sample of galaxies from the Sloan Digital Sky Survey Data Release 5 spectroscopic survey was used. We conclude that in such small groups as pairs and triplets, segregation by luminosity is clearly observed: galaxies in isolated pairs and triplets are on average two times more luminous than isolated galaxies. We consider the dark matter content in different systems. The median values of mass-to-luminosity ratio are  12 M_⊙/L_⊙  for the isolated pairs and  44 M_⊙/L_⊙  for the isolated triplets, and 7 (8)  M_⊙/L_⊙  for the most compact pairs (triplets). We also found that systems in denser environments have greater rms velocity and mass-to-luminosity ratio.