A Leap-frog Algorithm for Stochastic Dynamics

Abstract

A third-order algorithm for stochastic dynamics (SD) simulations is proposed, identical to the powerful molecular dynamics leap-frog algorithm in the limit of infinitely small friction coefficient γ. It belongs to the class of SD algorithms, in which the integration time step Δt is not limited by the condition Δt ≤ γ^?1, but only by the properties of the systematic force. It is shown how constraints, such as bond length or bond angle constraints, can be incorporated in the computational scheme. It is argued that the third-order Verlet-type SD algorithm proposed earlier may be simplified without loosing its third-order accuracy. The leap-frog SD algorithm is proven to be equivalent to the verlet-type SD algorithm. Both these SD algorithms are slightly more economical on computer storage than the Beeman-type SD algorithm.     

Mixing of nanofluids: molecular dynamics simulations and modelling

In this paper we study a mixing scheme, which has recently been proposed for microfluids, on the nanoscale. We do this by performing a series of nonequilibrium molecular dynamics simulations. On the nanoscale the chaotic mixing regime is captured. We discover a new phenomenon where the two mixing fluids exchange positions after leaving the mixing intersection. The results from the molecular dynamics simulations also reveal complex spatio-temporal stream velocity profiles generated by the mixing device. We find that these profiles can be modelled through an approximate analytical solution to the Navier–Stokes equation.     

Monte Carlo Study of the Buckingham Exponential-six Fluid

Abstract

In this paper we report the results of extensive Monte Carlo simulations of a pure fluid of Buckingham modified exponential-six molecules. Data are presented for the configurational energy and pressure covering a wide range of temperatures and densities. These data are interpreted using the generalized van der Waals partition function with a novel separation into free volume and mean potential terms. We find, surprisingly, that the Buckingham fluid is described by a simple van der Waals-like equation of state provided that the b parameter is temperature dependent and chosen in a theoretically correct manner.     

Direct Evaluation of Vapour-Liquid Equilibria of Mixtures by Molecular Dynamics Using Gibbs-Duhem Integration

Abstract

We present an extension of the Gibbs-Duhem integration method that permits direct evaluation of vapour-liquid equilibria of mixtures by molecular dynamics. The Gibbs-Duhem integration combines the best elements of the Gibbs ensemble Monte Carlo technique and thermodynamic integration. Given conditions of coexistence of pure substances, simultaneous but independent molecular dynamics simulations of each phase at constant number of particles, constant pressure, constant temperature and constant fugacity fraction of species 2 are carried out in succession along coexistence lines. In each simulation, the coexistence pressure is adjusted to satisfy the Clapeyron-type equation. The Clapeyron-type equation is a first-order nonlinear differential equation that prescribes how the pressure must change with the fugacity fraction of species 2 to maintain coexistence at constant temperature. The Clapeyron-type equation is solved by the predictor-corrector method. Running averages of mole fraction and compressibility factor for the two phases are used to evaluate the right-hand side of the Clapeyron-type equation. The Gibbs-Duhem integration method is applied to three prototypes of binary mixtures of the two-centre Lennard-Jones fluid having various elongations. The starting points on the coexistence curve were taken from published data.     

避难所对七星瓢虫捕食桃蚜功能反应的影响

自然情况下的捕食是以植物为中介发生的。作者采用植物模型作为捕食环境,测试了单纯的物理因素(避难所)对七星瓢虫Coccinella septempunctata L.捕食桃蚜Myzus persicae(Sulzer)的功能反应的影响。结果表明,避难所的存在能造成七星瓢虫捕食量和捕食率极显著下降,对雄虫的影响显著比雌虫强,且在中等猎物密度时影响更为突出。Logistic回归分析发现,七星瓢虫雌虫和雄虫无论避难所是否存在,功能反应类型均为HollingⅡ型。避难所不存在的情况下,雌、雄成虫的圆盘方程分别为N a=1.0259 N0/(1+0.003078 N0)和N a=0.9246 N0/(1+0.002127 N0),避难所存在的情况下,其圆盘方程分别为N a=0.9500 N0/(1+0.00342 N0)和N a=0.6463N0/(1+0.002262 N0)。避难所的存在能导致两性成虫瞬时攻击率极显著下降,但仅雄虫处置时间极显著延长,说明雄虫在环境异质性较高时捕获猎物更困难。     

Life‐History Patterns of Cuban Poeciliid Fishes (Teleostei: Cyprinodontiformes)

The following work provides basic information about the life history of 10 Cuban species of the family Poeciliidae. Adult fish stocks were captured in their natural habitat, and litters obtained from them were raised and maintained in captivity for 19 weeks. For each species, we present the mean value of newborn length (TL_o), age at sexual maturity (AM), total length at sexual maturity (TLM), as well as the patterns of postnatal growth in aquarium conditions, which were described using size–age curves and nonlinear regression equations (Richards model). There are differences in growth dynamics among species. In general, growth rates differ for both sexes in all poeciliids studied, males maturing earlier than females, who reach higher values of total length at the 19th week (TL_f). Sexual size dimorphism could be explained by the specific roles of each sex (fecundity in females and early maturity in males) while differences in growth among species could be related to their distribution patterns in the wild. The data summarized in this contribution can be useful for the conservation of these fish species. Zoo Biol 32:251–256, 2013. © 2012 Wiley Periodicals, Inc.     

THE EFFECTS OF EXTENDED WORKING HOURS ON HEALTH AND SOCIAL WELL-BEING—A COMPARATIVE ANALYSIS OF FOUR INDEPENDENT SAMPLES

Using structural equation modeling, it can be shown that long weekly working hours and work on weekends, nights, and in shifts have detrimental effects on psychovegetative health. Employees' reported subjective work-life balance also decreases with increasing number of hours worked/week, days worked on weekends, or at nights, and with working shifts. A decrease in work-life balance in turn increases the risk of psychovegetative impairments (PVIs). Thus, long and unusual working hours increase the risk of psychovegetative health impairments both directly and indirectly, moderated by the subjective work-life balance. In fact, the indirect effects of working time on PVIs via the work-life balance seem to be stronger than the direct effects. Results of a cross-validation study of four independent and representative samples from Germany and the European Union (N?>?50,000) indicate high structural stability of these results and thus an increased validity and range for generalization. (Author correspondence: anna.wirtz@gawo-ev.de)     

Relative roles of dispersal dynamics and competition in determining the isotopic niche breadth of a wetland fish

1. The niche variation hypothesis predicts that among‐individual variation in niche use will increase in the presence of intraspecific competition and decrease in the presence of interspecific competition. We sought to determine whether the local isotopic niche breadth of fish inhabiting a wetland was best explained by competition for resources and the niche variation hypothesis, by dispersal of individuals from locations with different prey resources or by a combination of the two. We analysed stable isotopes of carbon and nitrogen as indices of feeding niche and compared metrics of within‐site spread to characterise site‐level isotopic niche breadth. We then evaluated the explanatory power of competing models of the direct and indirect effects of several environmental variables spanning gradients of disturbance, competition strength and food availability on among‐individual variation of the eastern mosquitofish (Gambusia holbrooki). 2. The Dispersal model posits that only the direct effect of disturbance (i.e. changes in water level known to induce fish movement) influences among‐individual variation in isotopic niche. The Partitioning model allows for only direct effects of local food availability on among‐individual variation. The Combined model allows for both hypotheses by including the direct effects of disturbance and food availability. 3. A linear regression of the Combined model described more variance than models limited to the variables of either the Dispersal or Partitioning models. Of the independent variables considered, the food availability variable (per cent edible periphyton) explained the most variation in isotopic niche breadth, followed closely by the disturbance variable (days since last drying event). 4. Structural equation modelling provided further evidence that the Combined model was best supported by the data, with the Partitioning and the Dispersal models only modestly less informative. Again, the per cent edible periphyton was the variable with the largest direct effect on niche variability, with other food availability variables and the disturbance variable only slightly less important. Indirect effects of heterospecific and conspecific competitor densities were also important, through their effects on prey density. 5. Our results support the Combined hypotheses, although partitioning mechanisms appear to explain the most diet variation among individuals in the eastern mosquitofish. The results also support some predictions of the niche variation hypothesis, although both conspecific and interspecific competition appeared to increase isotopic niche breadth in contrast to predictions that interspecific competition would decrease it. We think this resulted from high diet overlap of co‐occurring species, most of which consume similar macroinvertebrates.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

Allometric disparity in rodent evolution

In this study, allometric trajectories for 51 rodent species, comprising equal representatives from each of the major clades (Ctenohystrica, Muroidea, Sciuridae), are compared in a multivariate morphospace (=allometric space) to quantify magnitudes of disparity in cranial growth. Variability in allometric trajectory patterns was compared to measures of adult disparity in each clade, and dietary habit among the examined species, which together encapsulated an ecomorphological breadth. Results indicate that the evolution of allometric trajectories in rodents is characterized by different features in sciurids compared with muroids and Ctenohystrica. Sciuridae was found to have a reduced magnitude of inter‐trajectory change and growth patterns with less variation in allometric coefficient values among members. In contrast, a greater magnitude of difference between trajectories and an increased variation in allometric coefficient values was evident for both Ctenohystrica and muroids. Ctenohystrica and muroids achieved considerably higher adult disparities than sciurids, suggesting that conservatism in allometric trajectory modification may constrain morphological diversity in rodents. The results provide support for a role of ecology (dietary habit) in the evolution of allometric trajectories in rodents.