基于应变率场的全球地震活动特征模拟探索

提出利用全球应变率资料模拟全球地震活动特征的基本思想,并作了初步探索。利用GSRM的全球应变率场结果,初步设计了模拟全球地震活动时空分布特征的细胞自动机模型。该模型把地球考虑为一个自组织的整体系统,制定了细胞自动机的演化规则,获得了模拟的人工地震目录。初步的模拟结果基本反映了全球地震活动的主要分布特征,体现了全球构造活动强弱的主要格局,初步达到了利用GPS等实测资料计算的应变率作为细胞自动机网格状态及其改变量来模拟复杂的全球地震活动特征的实验目的。     

Temporal stability of soil water content as affected by climate and soil hydraulic properties: a simulation study

Temporal stability of soil water content (TS SWC) is an often‐observed phenomenon, which characterization finds multiple applications. Climate and variability in soil properties are usually mentioned as factors of TS SWC, but their effects are far from clear. The objective of this work was to use SWC modeling to evaluate the effects of climate and soil hydraulic properties on the TS of soil water at different measurement schedules. We selected four representative climates found in USA and simulated the multiyear SWC dynamics for sandy loam, loam, and silty clay loam soils, all having the lognormal spatial distribution of the saturated hydraulic conductivity. The CLIMGEN and the HYDRUS6 codes were used to generate weather patterns and to simulate SWC, respectively. Four different methods were applied to select the representative location (RL). The low probability of having the same variability of mean relative differences of soil water under different climates was found in most of the cases. The probability that the variance of mean relative differences depended on sampling frequency was generally higher than 91% for the three soils. The interannual difference in mean relative differences variation from short and intensive summer campaigns was highly probable for all climates and soils. The RLs changed as climate and measurement scheduling changed, and they were less pronounced for coarse‐textured soils. The RL selection methods based solely on bias provided more consistency as compared with other methods. The TS appears to be the result of the interplay between climate, soil properties, and survey protocols. One implication of this factor interaction effect on TS SWC is that a simulation study can be useful to decide on the feasibility of including a search for TS‐based RLs for a specific site. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

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.     

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.