Sixth- and eighth-order Hermite integrator for N-body simulations

We present sixth- and eighth-order Hermite integrators for astrophysical N-body simulations, which use the derivatives of accelerations up to second-order (snap) and third-order (crackle). These schemes do not require previous values for the corrector, and require only one previous value to construct the predictor. Thus, they are fairly easy to implement. The additional cost of the calculation of the higher-order derivatives is not very high. Even for the eighth-order scheme, the number of floating-point operations for force calculation is only about two times larger than that for traditional fourth-order Hermite scheme. The sixth-order scheme is better than the traditional fourth-order scheme for most cases. When the required accuracy is very high, the eighth-order one is the best. These high-order schemes have several practical advantages. For example, they allow a larger number of particles to be integrated in parallel than the fourth-order scheme does, resulting in higher execution efficiency in both general-purpose parallel computers and GRAPE systems.     

Mass transfer in double white dwarf binaries

I review the implications of our recent numerical simulations of the hydrodynamics of Roche lobe overflow in double white dwarf binaries, with special emphasis on the stability and survival of such systems after the onset of dynamically unstable mass transfer. We find that some of these systems may survive and enter a secular phase of orbital expansion, driven by gravitational radiation, undergoing slow and decreasing mass transfer, and be observable as AM Canum Venaticorum binaries.     

Dark matter halo profiles in scale-free cosmologies

We explore the dependence of the central logarithmic slope of dark matter halo density profiles α on the spectral index n of the linear matter power spectrum P ( k ) using cosmological N -body simulations of scale-free models [i.e. P ( k ) ∝ k ⁿ ]. These simulations are based on a set of clear, reproducible and physically motivated criteria that fix the appropriate starting and stopping times for runs, and allow one to compare haloes across models with different spectral indices and mass resolutions. For each of our simulations we identify samples of well-resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parametrizing the mass profile using a 'generalized' Navarro, Frenk & White profile in which the central logarithmic slope α is allowed to vary while preserving the r ⁻³ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between α and n , such that α becomes shallower as n becomes steeper. However, if we normalize our mass profiles by r ₋₂, the radius at which the logarithmic slope of the density profile is −2, we find that these differences are no longer present. This is apparent if we plot the maximum slope     as a function of r / r ₋₂– we find that the profiles are similar for haloes forming in different n models. This reflects the importance of concentration, and reveals that the concentrations of haloes forming in steep- n cosmologies tend to be smaller than those of haloes forming in shallow- n cosmologies. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.     

The influence of gas expulsion and initial mass segregation on the stellar mass function of globular star clusters

Recently, De Marchi, Paresce & Pulone studied a sample of 20 globular clusters and found that all clusters with high concentrations have steep stellar mass functions while clusters with low concentration have comparatively shallow mass functions. No globular clusters were found with a flat mass function and high concentration. This seems curious since more concentrated star clusters are believed to be dynamically more evolved and should have lost more low-mass stars via evaporation, which would result in a shallower mass function in the low-mass part.
We show that this effect can be explained by residual-gas expulsion from initially mass segregated star clusters, and is enhanced further through unresolved binaries. If gas expulsion is the correct mechanism to produce the observed trend in the   c –α  -plane, then observation of these parameters would allow to constrain cluster starting conditions such as star formation efficiency and the time-scale of gas expulsion.     

Multimass schemes for collisionless N-body simulations

We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f ₀. Our method uses importance sampling to find the sampling DF f _s that minimizes the mean-square formal errors in a given set of projections of the DF f ₀. The result is a multimass N -body realization of the galaxy model in which 'interesting' regions of phase space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher mass particles.
As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N -body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor of ∼100 reduction in the variance at the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N -body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N -body simulations closer to the collisionless ideal.     

On the structure of tidal tails

We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster.     

The evolution of the galaxy red sequence in simulated clusters and groups

N -body/hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a Λ cold dark matter (ΛCDM) cosmology are used in order to follow the building-up of the colour–magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the red sequence (RS) as they get aged over times and eventually set upon a 'dead sequence' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However, galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z ∼ 1, after which major events only involve the brightest cluster galaxies.
The emerging scenario is that the evolution of the colour–magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter haloes assembly. Galaxies progressively quenching their star formation settle to a very sharp 'red and dead' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z ∼ 1.5) and environment.
Differently from the DS, the operatively defined RS evolves more evidently with z , the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star-forming one: this goes from z ≃ 1 in clusters down to 0.4 in normal groups.     

CFD 在建筑物对气象台站测风影响分析中的应用

受城市化、探测环境变化等人为因素影响,部分气象站风速序列存在不均一性,在使用这类气象站的风速资料时必须进行订正。本文以东山气象站为例,应用测站受建筑物影响前后各5 a的风速资料及探测环境历史沿革资料,利用计算流体力学方法（CFD）研究建筑物对气象站测风数据的影响。结果表明：建筑物与测风杆的相对位置不同,建筑物对测风的影响程度亦不同,测站北面建筑物是影响该站偏北风风速观测的主要因素;研究不同风速下建筑的影响,可得模拟结果与输入风速间呈线性关系,并由此建立风速订正关系式,订正后基本消除了建筑物对N方位风速的影响;订正受建筑物影响较小的NNE方位风速需考虑其他障碍物的影响。本文通过试验,验证了CFD方法可用于定量评估建筑物对测风的影响,从而重建受建筑物影响台站的风速序列。