Simulations of planet migration driven by planetesimal scattering

Evidence has mounted for some time that planet migration is an important part of the formation of planetary systems, both in the Solar System [Malhotra, R., 1993. Nature 365, 819-821] and in extrasolar systems [Mayor, M., Queloz, D., 1995. Nature 378, 355-359; Lin, D.N.C., Bodenheimer, P., Richardson, D.C., 1996. Nature 380, 606-607]. One mechanism that produces migration (the change in a planet's semi-major axis a over time) is the scattering of comet- and asteroid-size bodies called planetesimals [Fernandez, J.A., Ip, W.-H., 1984. Icarus 58, 109-120]. Significant angular momentum exchange can occur between the planets and the planetesimals during local scattering, enough to cause a rapid, self-sustained migration of the planet [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445]. This migration has been studied for the particular case of the four outer planets of the Solar System (as in Gomes et al. [Gomes, R.S., Morbidelli, A., Levison, H.F., 2004. Icarus 170, 492-507]), but is not well understood in general. We have used the Miranda [McNeil, D., Duncan, M., Levison, H.F., 2005. Astron. J. 130, 2884-2899] computer simulation code to perform a broad parameter-space survey of the physical variables that determine the migration of a single planet in a planetesimal disk. Migration is found to be predominantly inwards, and the migration rate is found to be independent of planet mass for low-mass planets in relatively high-mass disks. Indeed, a simple scaling relation from Ida et al. [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445] matches well with the dependencies of the migration rate:

(1)     

The clustering and host haloes of galaxy mergers at high redshift

High-resolution simulations of cosmological structure formation indicate that dark matter substructure in dense environments, such as groups and clusters, may survive for a long time. These dark matter subhaloes are the likely hosts of galaxies. We examine the small-scale spatial clustering of subhalo major mergers at high redshift using high-resolution N -body simulations of cosmological volumes. Recently merged, massive subhaloes exhibit enhanced clustering on scales  ∼100–300  h ⁻¹ kpc  , relative to all subhaloes of the same infall mass, for a short time after a major merger (<500 Myr). The small-scale clustering enhancement is smaller for lower mass subhaloes, which also show a deficit on scales just beyond the excess. Haloes hosting recent subhalo mergers tend to have more subhaloes; for massive subhaloes, the excess is stronger and it tends to increase for the most massive host haloes. The subhalo merger fraction is independent of halo mass for the scales we probe. In terms of satellite and central subhaloes, the merger increase in small-scale clustering for massive subhaloes arises from recently merged massive central subhaloes having an enhanced satellite population. Our mergers are defined via their parent infall mass ratios. Subhaloes experiencing major mass gains also exhibit a small-scale clustering enhancement, but these correspond to two-body interactions leading to two final subhaloes, rather than subhalo coalescence.     

Mass transfer in eccentric binaries: the new oil-on-water smoothed particle hydrodynamics technique

To measure the onset of mass transfer in eccentric binaries, we have developed a two-phase smoothed particle hydrodynamics (SPH) technique. Mass transfer is important in the evolution of close binaries, and a key issue is to determine the separation at which mass transfer begins. The circular case is well understood and can be treated through the use of the Roche formalism. To treat the eccentric case, we use a newly developed two-phase system. The body of the donor star is made up from high-mass water particles, whilst the atmosphere is modelled with low-mass oil particles. Both sets of particles take part fully in SPH interactions. To test the technique, we model circular mass-transfer binaries containing a  0.6 M_⊙  donor star and a  1 M_⊙  white dwarf; such binaries are thought to form cataclysmic variable (CV) systems. We find that we can reproduce a reasonable CV mass-transfer rate, and that our extended atmosphere gives a separation that is too large by approximately 16 per cent, although its pressure scale height is considerably exaggerated. We use the technique to measure the semimajor axis required for the onset of mass transfer in binaries with a mass ratio of   q = 0.6  and a range of eccentricities. Comparing to the value obtained by considering the instantaneous Roche lobe at pericentre, we find that the radius of the star required for mass transfer to begin decreases systematically with increasing eccentricity.     

Breaking halo occupation degeneracies with marked statistics

We show that a suitably defined marked correlation function can be used to break degeneracies in halo-occupation distribution modelling. The statistic can be computed on both three-dimensional and two-dimensional data sets and should be applicable to all upcoming galaxy surveys. A proof of principle, using mock catalogues created from N -body simulations, is given.     

Conceiving and marketing NASA’s Great Observatories

In early 1984, the astronomical community’s plans to launch a series of powerful telescopes able to study celestial sources at almost any wavelength were in trouble. The President of the United States had just declared his priority for a Space Station that was bound to be expensive, and Congress could not understand why yet another set of space observatories was needed when others were already being funded. To realize their aims, astronomers would have to advocate their needs much more effectively than in the past.     

B-Pol: detecting primordial gravitational waves generated during inflation

B-Pol is a medium-class space mission aimed at detecting the primordial gravitational waves generated during inflation through high accuracy measurements of the Cosmic Microwave Background polarization. We discuss the scientific background, feasibility of the experiment, and implementation developed in response to the ESA Cosmic Vision 2015-2025 Call for Proposals. See for the full list of collaboration members and a full copy of the B-Pol proposal.     

Predictions for the LCROSS mission

Abstract— We describe the results of a variety of model calculations for predictions of observable results of the LCROSS mission to be launched in 2009. Several models covering different aspects of the event are described along with their results. Our aim is to bracket the range of expected results and produce a useful guide for mission planning. In this paper, we focus on several different questions, which are modeled by different methods. The questions include the size of impact crater, the mass, velocity, and visibility of impact ejecta, and the mass and temperature of the initial vapor plume. The mass and velocity profiles of the ejecta are of primary interest, as the ejecta will be the main target of the S‐S/C observations. In particular, we focus on such quantities as the amount of mass that reaches various heights. A height of 2 km above the target is of special interest, as we expect that the EDUS impact will take place on the floor of a moderate‐sized crater ?30 km in diameter, with a rim height of 1–2 km. The impact ejecta must rise above the crater rim at the target site in order to scatter sunlight and become visible to the detectors aboard the S‐S/C. We start with a brief discussion of crater scaling relationships as applied to the impact of the EDUS second stage and resulting estimated crater diameter and ejecta mass. Next we describe results from the RAGE hydrocode as applied to modeling the short time scale (t 0.1 s) thermal plume that is expected to occur immediately after the impact. We present results from several large‐scale smooth‐particle hydrodynamics (SPH) calculations, along with results from a ZEUS‐MP hydrocode model of the crater formation and ejecta mass‐velocity distribution. We finish with two semi‐analytic models, the first being a Monte Carlo model of the distribution of expected ejecta, based on scaling models using a plausible range of crater and ejecta parameters, and the second being a simple model of observational predictions for the shepherding spacecraft (S‐S/C) that will follow the impact for several minutes until its own impact into the lunar surface. For the initial thermal plume, we predict an initial expansion velocity of ?7 km s^?1, and a maximum temperature of ?1200 K. Scaling relations for crater formation and the SPH calculation predict a crater with a diameter of ?15 m, a total ejecta mass of ?10⁶kg, with ?10⁴kg reaching an altitude of 2 km above the target. Both the SPH and ZEUS‐MP calculations predict a maximum ejecta velocity of ?1 km s^?1. The semi‐analytic Monte Carlo calculations produce more conservative estimates (by a factor of ?5) for ejecta at 2 km, but with a large dispersion in possible results.     

Stability at high pressure, elastic behavior and pressure-induced structural evolution of “Al5BO9”, a mullite-type ceramic material

Elastic behavior and pressure-induced structural evolution of synthetic boron-mullite “Al₅BO₉” (a = 5.678(2) Å, b = 15.015(4) Å and c = 7.700(3) Å, space group Cmc2₁, Z = 4) were investigated up to 7.4 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase transition or anomalous compressional behavior occurred within the investigated P range. Fitting the P–V data with a truncated second-order (in energy) Birch-Murnaghan Equation-of-State (BM-EoS), using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 656.4(3) Å³ and K _T0 = 165(7) GPa (β _V0 = 0.0061(3) GPa^?1). The evolution of the Eulerian finite strain versus normalized stress (f _E–F _E plot) leads to an almost horizontal trend, showing that a truncated second-order BM-EoS is appropriate to describe the elastic behavior of “Al₅BO₉” within the investigated P range. The weighted linear regression through the data points gives: F _E(0) = 159(11) GPa. Axial compressibility coefficients yielded: β _a  = 1.4(2) × 10^?3 GPa^?1, β _b  = 3.4(4) × 10^?3 GPa^?1, and β _c  = 1.7(3) × 10^?3 GPa^?1 (β _a :β _b :β _c  = 1:2.43:1.21). The highest compressibilities observed in this study within (100) can be ascribed to the presence of voids represented by five-membered rings of polyhedra: Al1–Al3–Al4–Al1–Al3, which allow accommodating the effect of pressure by polyhedral tilting. Polyhedral tilting around the voids also explains the higher compressibility along [010] than along [001]. The stiffer crystallographic direction observed here might be controlled by the infinite chains of edge-sharing octahedra running along [100], which act as “pillars”, making the structure less compressible along the a-axis than along the b- and c-axis. Along [100], compression can only be accommodated by deformation of the edge-sharing octahedra (and/or by compression of the Al–O bond lengths), as no polyhedral tilting can occur. In addition, a comparative elastic analysis among the mullite-type materials is carried out.