Do binaries in clusters form in the same way as in the field?

We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities  (10⁴– 10⁵ M_⊙ pc⁻³)  almost all binaries with separations  >10³  au are destroyed after a few crossing times. At low densities [     ], many binaries with separations  >10³  au are destroyed, and no binaries with separations  >10⁴  au survive after a few crossing times. Therefore, the binary separations in clusters can be used as a tracer of the dynamical age and past density of a cluster.
We argue that the central region of the Orion nebula cluster was ∼100 times denser in the past with a half-mass radius of only 0.1–0.2 pc as (i) it is expanding, (ii) it has very few binaries with separations  >10³  au and (iii) it is well mixed and therefore dynamically old.
We also examine the origin of the field binary population. Binaries with separations  <10²  au are not significantly modified in any cluster, therefore at these separations the field reflects the sum of all star formation. Binaries with separations in the range  10²– 10⁴  au are progressively more and more heavily affected by dynamical disruption in increasingly dense clusters. If most star formation is clustered, these binaries must be overproduced relative to the field. Finally, no binary with a separation  >10⁴  au can survive in any cluster and so must be produced by isolated star formation, but only if all isolated star formation produces extremely wide binaries.     

High-inclination planets and asteroids in multistellar systems

The Kozai mechanism often destabilizes high-inclination orbits. It couples changes in the eccentricity and inclination, and drives high inclination, circular orbits to low inclination, eccentric orbits. In a recent study of the dynamics of planetesimals in the quadruple star system HD 98800, there were significant numbers of stable particles in circumbinary polar orbits about the inner binary pair which are apparently able to evade the Kozai instability.
Here, we isolate this feature and investigate the dynamics through numerical and analytical models. The results show that the Kozai mechanism of the outer star is disrupted by a nodal libration induced by the inner binary pair on a shorter time-scale. By empirically modelling the period of the libration, a criteria for determining the high-inclination stability limits in general triple systems is derived. The nodal libration feature is interesting and, although affecting inclination and node only, shows many parallels to the Kozai mechanism. This raises the possibility that high-inclination planets and asteroids may be able to survive in multistellar systems.     

Tailoring triaxial N-body models via a novel made-to-measure method

The made-to-measure N -body method slowly adapts the particle weights of an N -body model, whilst integrating the trajectories in an assumed static potential, until some constraints are satisfied, such as optimal fits to observational data. I propose a novel technique for this adaption procedure, which overcomes several limitations and shortcomings of the original method. The capability of the new technique is demonstrated by generating realistic N -body equilibrium models for dark matter haloes with prescribed density profile, triaxial shape and slowly outwardly growing radial velocity anisotropy.     

How well do starlab and nbody4 compare? I. Simple models

N -body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N -body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N -body codes starlab/kira and nbody4 are widely used for a range of applications, there is no thorough comparison study yet.
Here, we thoroughly compare basic quantities as derived from simulations performed either with starlab/kira or nbody4 .
We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields, etc.), evolve them in parallel with starlab/kira and nbody4 , analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we develop a bootstrap algorithm for functional dependencies.
We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies.
Our results reassure the comparability and the possibility to combine results from these two major N -body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. Further detailed comparison studies for more complex systems, e.g. including stellar/binary evolution, are required.     

Redistributing hot gas around galaxies: do cool clouds signal a solution to the overcooling problem?

We present a pair of high-resolution smoothed particle hydrodynamics simulations that explore the evolution and cooling behaviour of hot gas around Milky Way size galaxies. The simulations contain the same total baryonic mass and are identical other than their initial gas density distributions. The first is initialized with a low-entropy hot gas halo that traces the cuspy profile of the dark matter, and the second is initialized with a high-entropy hot halo with a cored density profile as might be expected in models with pre-heating feedback. Galaxy formation proceeds in dramatically different fashion depending on the initial setup. While the low-entropy halo cools rapidly, primarily from the central region, the high-entropy halo is quasi-stable for  ∼4 Gyr  and eventually cools via the fragmentation and infall of clouds from ∼100 kpc distances. The low-entropy halo's X-ray surface brightness is ∼100 times brighter than current limits and the resultant disc galaxy contains more than half of the system's baryons. The high-entropy halo has an X-ray brightness that is in line with observations, an extended distribution of pressure-confined clouds reminiscent of observed populations and a final disc galaxy that has half the mass and ∼50 per cent more specific angular momentum than the disc formed in the low-entropy simulation. The final high-entropy system retains the majority of its baryons in a low-density hot halo. The hot halo harbours a trace population of cool, mostly ionized, pressure-confined clouds that contain ∼10 per cent of the halo's baryons after 10 Gyr of cooling. The covering fraction for H  i and Mg  ii absorption clouds in the high-entropy halo is ∼0.4 and ∼0.6, respectively, although most of the mass that fuels disc growth is ionized, and hence would be under counted in H  i surveys.     

The dynamical architecture and habitable zones of the quintuplet planetary system 55 Cancri

We perform numerical simulations to study the secular orbital evolution and dynamical structure of the quintuplet planetary system 55 Cancri with the self-consistent orbital solutions by Fischer and coworkers. In the simulations, we show that this system can be stable for at least 108 yr. In addition, we extensively investigate the planetary configuration of four outer companions with one terrestrial planet in the wide region of 0.790AU 〈 a 〈 5.900AU to examine the existence of potential asteroid structure and Habitable Zones （HZs）. We show that there are unstable regions for orbits about 4：1, 3：1 and 5：2 mean motion resonances （MMRs） of the outermost planet in the system, and several stable orbits can remain at 3：2 and 1：1 MMRs, which resembles the asteroid belt in the solar system. From a dynamical viewpoint, proper HZ candidates for the existence of more potential terrestrial planets reside in the wide area between 1.0 AU and 2.3 AU with relatively low eccentricities.     

Dependence of the accuracy of precipitation and cloud simulation on temporal and spatial scales

Precipitation and associated cloud hydrometeors have large temporal and spatial variability, which makes accurate quantitative precipitation forecasting difficult. Thus, dependence of accurate precipitation and associated cloud simulation on temporal and spatial scales becomes an important issue. We report a cloud-resolving modeling analysis on this issue by comparing the control experiment with experiments perturbed by initial temperature, water vapor, and cloud conditions. The simulation is considered to be accurate only if the root-mean-squared difference between the perturbation experiments and the control experiment is smaller than the standard deviation. The analysis may suggest that accurate precipitation and cloud simulations cannot be obtained on both fine temporal and spatial scales simultaneously, which limits quantitative precipitation forecasting. The accurate simulation of water vapor convergence could lead to accurate precipitation and cloud simulations on daily time scales, but it may not be beneficial to precipitation and cloud simulations on hourly time scales due to the dominance of cloud processes.     

Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years

We describe the evolutionary response of northern and southern hemisphere summer monsoons to orbital forcing over the past 280,000 years using a fully coupled general circulation ocean-atmosphere model in which the orbital forcing is accelerated by a factor of 100. We find a strong and positive response of northern (southern) summer monsoon precipitation to northern (southern) summer insolation forcing. On average, July (January) precipitation maxima and JJA (DJF) precipitation maxima have high coherence and are approximately in phase with June (December) insolation maxima, implying an average lag between forcing and response of about 30° of phase at the precession period. The average lag increases to over 40° for 4-month precipitation averages, JJAS (DJFM). The phase varies from region to region. The average JJA (DJF) land temperature maxima also lag the June orbital forcing maxima by about 30° of phase, whereas ocean temperature maxima exhibit a lag of about 60° of phase at the precession period. Using generalized measures of the thermal and hydrologic processes that produce monsoons, we find that the summer monsoon precipitation indices for the six regions all fall within the phase limits of the process indices for the respective hemispheres. Selected observational studies from four of the six monsoon regions report approximate in-phase relations of summer monsoon proxies to summer insolation. However other observational studies report substantial phase lags of monsoon proxies and a strong component of forcing associated with glacial-age boundary conditions or other factors. An important next step will be to include glacial-age boundary condition forcing in long, transient paleoclimate simulations, along with orbital forcing.     

Numerical modelling of a point-absorbing wave energy converter in irregular and extreme waves

Based on the Navier-Stokes (RANS) equations, a three-dimensional (3-D) mathematical model for the hydrodynamics and structural dynamics of a floating point-absorbing wave energy converter (WEC) with a stroke control system in irregular and extreme waves is presented. The model is validated by a comparison of the numerical results with the wave tank experiment results of other researchers. The validated model is then utilized to examine the effect of wave height on structure displacements and connection rope tension. In the examined cases, the differences in WEC’s performance exhibited by an inviscid fluid and a viscous fluid can be neglected. Our results also reveal that the differences in behavior predicted by boundary element method (BEM) and the RANS-based method can be significant and vary considerably, depending on wave height.