Planetary theories with the aid of the expansions of elliptic functions

For coplanar circular orbits, the mutual perturbations between two bodies can be expressed in term of the argument of Jacobian elliptic functions instead of the difference of the mean longitudes. For a given pair of planets, such a change of time variable improves the convergence of the developments. At the first order of planetary masses an integration of Lagrange's equations for the osculating elements is performed. When compared to classical developments the results are reduced by an important factor. The method is then extended to the mutual perturbations of Jupiter and Saturn, at any order of planetary masses, either with Fourier series with two arguments, or with one argument solely, taking advantage of the close commensurability of the mean motions.     

Numerical simulations of granular dynamics II: Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties.We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration.We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.     

The High-Resolution Coronal Imager (Hi-C)

The High-Resolution Coronal Imager (Hi-C) was flown on a NASA sounding rocket on 11 July 2012. The goal of the Hi-C mission was to obtain high-resolution (≈?0.3?–?0.4′′), high-cadence (≈?5 seconds) images of a solar active region to investigate the dynamics of solar coronal structures at small spatial scales. The instrument consists of a normal-incidence telescope with the optics coated with multilayers to reflect a narrow wavelength range around 19.3 nm (including the Fe xii 19.5-nm spectral line) and a 4096×4096 camera with a plate scale of 0.1′′?pixel^?1. The target of the Hi-C rocket flight was Active Region 11520. Hi-C obtained 37 full-frame images and 86 partial-frame images during the rocket flight. Analysis of the Hi-C data indicates the corona is structured on scales smaller than currently resolved by existing satellite missions.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

The effect of the dynamical state of clusters on gas expulsion and infant mortality

The star-formation efficiency (SFE) of a star cluster is thought to be the critical factor in determining if the cluster can survive for a significant (>50 Myr) time. There is an often quoted critical SFE of ～30% for a cluster to survive gas expulsion. I reiterate that the SFE is not the critical factor; rather, it is the dynamical state of the stars (as measured by their virial ratio) immediately before gas expulsion. If the stars in a star cluster are born in an (even slightly) cold dynamical state, then the survivability of a cluster can be greatly increased.     

Nanometer‐scale anatomy of entire Stardust tracks

Abstract– We have developed new sample preparation and analytical techniques tailored for entire aerogel tracks of Wild 2 sample analyses both on “carrot” and “bulbous” tracks. We have successfully ultramicrotomed an entire track along its axis while preserving its original shape. This innovation allowed us to examine the distribution of fragments along the entire track from the entrance hole all the way to the terminal particle. The crystalline silicates we measured have Mg‐rich compositions and O isotopic compositions in the range of meteoritic materials, implying that they originated in the inner solar system. The terminal particle of the carrot track is a ¹⁶O‐rich forsteritic grain that may have formed in a similar environment as Ca‐, Al‐rich inclusions and amoeboid olivine aggregates in primitive carbonaceous chondrites. The track also contains submicron‐sized diamond grains likely formed in the solar system. Complex aromatic hydrocarbons distributed along aerogel tracks and in terminal particles. These organics are likely cometary but affected by shock heating.     

Supergranule and Mesogranule Evolution

The MDI instrument on the SOHO satellite obtained a nearly continuous 45.5-hr run in high-resolution mode on 17–18 January 1997, collecting continuum, Dopplergram, and magnetogram images once per minute. This is one of the longest data sets yet obtained in this mode and shows significant evolution of the supergranulation pattern. After allowing for solar rotation within the fixed field of view, an area spanning 17° in latitude and 11° in longitude was extracted that covers the same area of the solar surface for the entire run. From the de-rotated continuum images, we computed flow maps of photospheric motions using local correlation techniques (LCT). Horizontal divergence maps constructed from the flow maps show local maxima of the size of mesogranules (5–10). We interpret these as mesogranules although the LCT flow map resolution (4.8 FWHM) may not completely resolve smaller mesogranules. Movies made from the divergence maps clearly show the outward convection (advection) of these mesogranules within each supergranule, and narrow boundaries of negative divergence outlining the supergranules. Several new supergranules are observed forming. These appear as areas of strong divergence that pop up between pre-existing supergranules and grow, pushing their neighbors apart. Others seem to perish between growing neighbors. We also computed the vertical component of vorticity from the flow maps. Movies of this vorticity do not show any obvious patterns.     

The Monitor project: rotation of low-mass stars in the open cluster NGC 2547

We report on the results of an I -band time-series photometric survey of NGC 2547 using the MPG/ESO 2.2-m telescope with Wide Field Imager, achieving better than 1 per cent photometric precision per data point over  14 ≲ I ≲ 18  . Candidate cluster members were selected from a V versus V − I colour–magnitude diagram over  12.5 < V < 24  (covering masses from  0.9 M_⊙  down to below the brown dwarf limit), finding 800 candidates, of which we expect ∼330 to be real cluster members, taking into account contamination from the field (which is most severe at the extremes of our mass range). Searching for periodic variations in these gave 176 detections over the mass range  0.1 ≲ M /M_⊙≲ 0.9  . The rotation period distributions were found to show a clear mass-dependent morphology, qualitatively intermediate between the distributions obtained from similar surveys in NGC 2362 and 2516, as would be expected from the age of this cluster. Models of the rotational evolution were investigated, finding that the evolution from NGC 2362 to 2547 was qualitatively reproduced (given the uncertainty in the age of NGC 2547) by solid body and core-envelope decoupled models from our earlier NGC 2516 study without need for significant modification.     

The dependence of galaxy formation on cosmological parameters: can we distinguish between the WMAP1 and WMAP3 parameter sets?

We combine N -body simulations of structure growth with physical modelling of galaxy evolution to investigate whether the shift in cosmological parameters between the first- and third-year results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) affects predictions for the galaxy population. Structure formation is significantly delayed in the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on the relevant scales. The decrease in dark matter clustering strength is, however, almost entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely altered. In both cosmologies, several combinations of physical parameters can reproduce observed, low-redshift galaxy properties; the star formation, supernova feedback and active galactic nucleus feedback efficiencies can be played off against each other to give similar results. Models which fit observed luminosity functions predict projected two-point correlation functions which scatter by about 10–20 per cent on large scale and by larger factors on small scale, depending both on cosmology and on details of galaxy formation. Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but careful treatment of the systematics is needed. Given present modelling uncertainties, it is not easy to distinguish between the WMAP1 and WMAP3 cosmologies on the basis of low-redshift galaxy properties. Model predictions diverge more dramatically at high redshift. Better observational data at   z > 2  will better constrain galaxy formation and perhaps also cosmological parameters.