N-body simulation for self-gravitating collisional systems with a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of ∼20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved ∼90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N ∼ 10⁵ on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.     

Impact of the blockage ratio on the transport of sediment in the presence of a hydrokinetic turbine: Numerical modeling of the interaction sediment and turbine

Interactions between hydrokinetic turbines and near scale and far scale bed sediment particles are considered a critical area of assessment; however, a limited number of research studies have been published to address this issue. The current paper explores the influence of the blockage ratio on the morphology evolution downstream from a hydrokinetic turbine. A modeling framework is derived to predict the significant transport induced by a turbine installed on the erodible fluvial sandy bed surface, such as the Eulerian multi-phase model for the sediment and the Blade Element Momentum Theory (BEMT) for the turbine, using the open source platform OpenFOAM. Two configurations of different blockage ratios are considered. The configurations show differences in the bed morphology and the sediment transport. It has been shown for both configurations that the scour capabilities are enhanced below the turbine due to the acceleration of the flow and increasing local shear stress on the sediment.     

The magnetic field near Mars: A comparison between a hybrid model, Mars Global Surveyor and Mars Express observations

Mars Express (MEX) does not carry its own magnetometer which complicates interpretation of ASPERA-3/MEX ion measurements. The direction of the interplanetary magnetic field (IMF) is especially important because it, among other things, determines the direction of the convective electric field and orientation of the cross tail current sheet and tail lobes. In this paper we present a case study to show the properties of the magnetic field near Mars in a quasi-neutral hybrid (QNH) model at the orbits where the Mars Global Surveyor (MGS) has made measurements, present a method to derive the IMF clock angle by comparing fields in a hybrid model and the direction of the magnetic field measured by MGS by deriving the IMF clock angle. We also use H⁺ ring velocity distribution observations upstream of the bow shock measured by the IMA/ASPERA-3 instrument on board MEX spacecraft. These observations are used to indirectly provide the orientation of the IMF. We use a QNH model (HYB-Mars) where ions are modeled as particles while electrons form a mass-less charge neutralizing fluid. We found that the direct MGS and non-direct IMA observations of the orientation magnetic field vectors in non-crustal magnetic field regions are consistent with the global magnetic field draping pattern predicted by the global model.     

Intermediate-mass black holes and ultraluminous X-ray sources in the Cartwheel ring galaxy

Chandra and XMM–Newton observations of the Cartwheel galaxy show ∼17 bright X-ray sources  (≳5 × 10³⁸ erg s⁻¹)  , all within the gas-rich outer ring. We explore the hypothesis that these X-ray sources are powered by intermediate-mass black holes (IMBHs) accreting gas or undergoing mass transfer from a stellar companion. To this purpose, we run N -body/smoothed particle hydrodynamics simulations of the galaxy interaction which might have led to the formation of Cartwheel, tracking the dynamical evolution of two different IMBH populations: halo and disc IMBHs. Halo IMBHs cannot account for the observed X-ray sources, as only a few of them cross the outer ring. Instead, more than half of the disc IMBHs are pulled in the outer ring as a consequence of the galaxy collision. However, also in the case of disc IMBHs, accretion from surrounding gas clouds cannot account for the high luminosities of the observed sources. Finally, more than 500 disc IMBHs are required to produce ≲15 X-ray sources via mass transfer from very young stellar companions. Such number of IMBHs is very large and implies extreme assumptions. Thus, the hypothesis that all the observed X-ray sources in Cartwheel are associated with IMBHs is hardly consistent with our simulations, even if it is still possible that IMBHs account for the few (≲1–5) brightest ultraluminous X-ray sources.     

Are ring galaxies the ancestors of giant low surface brightness galaxies?

We simulate the collisional formation of a ring galaxy and we integrate its evolution up to 1.5 Gyr after the interaction. About 100–200 Myr after the collision, the simulated galaxy is very similar to observed ring galaxies (e.g. Cartwheel). After this stage, the ring keeps expanding and fades. Approximately 0.5–1 Gyr after the interaction, the disc becomes very large (∼100 kpc) and flat. Such extended discs have been observed only in giant low surface brightness galaxies (GLSBs). We compare various properties of our simulated galaxies (surface brightness profile, morphology, H  i spectrum and rotation curve) with the observations of four well-known GLSBs (UGC 6614, Malin 1, Malin 2 and NGC 7589). The simulations match quite well the observations, suggesting that ring galaxies could be the progenitors of GLSBs. This result is crucial for the cold dark matter (CDM) model, as it was very difficult, so far, to explain the formation of GLSBs within the CDM scenario.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

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.