A simple model of human foveal ganglion cell responses to hyperacuity stimuli

We developed a physiologically plausible model of the first steps of spatial visual information processing in the fovea of the human retina. With the predictions of this model we could support the hypothesis that, for moderate contrasts ( 40%), hyperacuity is mediated by the magnocellular (MC-) pathway. Despite the lower sampling density in the MC pathway, as compared to the parvocellular (PC-) pathway, the information that is transferred by the MC ganglion cells is sufficient to achieve thresholds comparable to those of human subjects in psychophysical tasks. This is a result of the much higher signal-to-noise ratio of the MC pathway cell signals. The PC pathway cells do not transfer enough information for hyperacuity thresholds.     

蛋白质中原子与基团的可及性的一种新计算方法及其应用

实现了一种新的蛋白质中原子层次可及面积的计算方法－Ｍｏｎｔｅ－Ｃａｒｌｏ模拟方法。计算了溶菌酶各个原子与功能基团的可及性,并应用于蛋白质结构和功能的研究中。使得对蛋白质性能的预测更准确,提供了一种对蛋白质结构和功能研究的新思路。     

心肌电兴奋传导速度对心律失常影响的仿真研究

心电场是由心肌的电活动产生的。心肌细胞的电特性及心肌细胞间的传导关系决定了体表电位的分布及心电图的变化。心肌电兴奋传导速度则是影响心肌间兴奋传导关系的重要参数之一。由于很难通过实验方法来人为改变电兴奋传导速度,因而临床上有关该参数对心律影响的定量知识相当缺乏。本文采用真实三雏躯干模型及心脏模型,对心肌电兴奋传导速度与心律变化的关系进行定量仿真研究。结果表明,兴奋传导速度决定了整个心电图的变化,而局部普通心肌的传导速度在相当范围内变化似乎对心电图影响不明显,但传导速度超过一定范围后可能产生突变。     

The evolution of brood parasitism: the role of facultative parasitism

The hypothesis that facultative brood parasitism may serve asan intermediate step in the evolutionary transition from purelyparental reproduction to obligate parasitism was investigated.The population dynamics of a host-parasite complex were computer-simulatedin a model that incorporated different intensities of parasitismand host defense and considered a simplified semelparous birdspecies living in a homogeneous habitat The individuals usetwo different breeding strategies: provide parental care orparasitize the nest of those providing parental care. Underobligate parasitism, the parasites appeared unsuccessful, drovethe host population to extinction, or coexisted with the hostin stable or oscillating proportions. The behavior of the systemdepended on both the effectiveness of the parasite and the defenseof the host. Under facultative parasitism (making the best ofa bad job), the parasites reduced host numbers but did not reducethe population size below the number of breeding sites. Thus,facultative parasitism provides a better opportunity for thedevelopment of defense in the host. The population of a hostthat shows a certain level of defense can be more successfullyinvaded by obligate parasites so that stable coexistence ofhosts and parasites is possible.     

Mathematical simulation of the interactions among cyanobacteria, purple sulfur bacteria and chemotrophic sulfur bacteria in microbial mat communities

Abstract: A deterministic one-dimensional reaction diffusion model was constructed to simulate benthic stratification patterns and population dynamics of cyanobacteria, purple and colorless sulfur bacteria as found in marine microbial mats. The model involves the major biogeochemical processes of the sulfur cycle and includes growth metabolism and their kinetic parameters as described from laboratory experimentation. Hence, the metabolic production and consumption processes are coupled to population growth. The model is used to calculate benthic oxygen, sulfide and light profiles and to infer spatial relationships and interactions among the different populations. Furthermore, the model is used to explore the effect of different abiotic and biotic environmental parameters on the community structure. A strikingly clear pattern emerged of the interaction between purple and colorless sulfur bacteria: either colorless sulfur bacteria dominate or a coexistence is found of colorless and purple sulfur bacteria. The model predicts that purple sulfur bacteria only proliferate when the studied environmental parameters surpass well-defined threshold levels. However, once the appropriate conditions do occur, the purple sulfur bacteria are extremely successful as their biomass outweighs that of colorless sulfur bacteria by a factor of up to 17. The typical stratification pattern predicted closely resembles the often described bilayer communities which comprise a layer of purple sulfur bacteria below a cyanobacterial top-layer; colorless sulfur bacteria are predicted to sandwich in between both layers. The profiles of oxygen and sulfide shift on a diel basis similarly as observed in real systems.     

Evaluation of several methods for estimating phylogenetic trees when substitution rates differ over nucleotide sites

Several maximum likelihood and distance matrix methods for estimating phylogenetic trees from homologous DNA sequences were compared when substitution rates at sites were assumed to follow a gamma distribution. Computer simulations were performed to estimate the probabilities that various tree estimation methods recover the true tree topology. The case of four species was considered, and a few combinations of parameters were examined. Attention was applied to discriminating among different sources of error in tree reconstruction, i.e., the inconsistency of the tree estimation method, the sampling error in the estimated tree due to limited sequence length, and the sampling error in the estimated probability due to the number of simulations being limited. Compared to the least squares method based on pairwise distance estimates, the joint likelihood analysis is found to be more robust when rate variation over sites is present but ignored and an assumption is thus violated. With limited data, the likelihood method has a much higher probability of recovering the true tree and is therefore more efficient than the least squares method. The concept of statistical consistency of a tree estimation method and its implications were explored, and it is suggested that, while the efficiency (or sampling error) of a tree estimation method is a very important property, statistical consistency of the method over a wide range of, if not all, parameter values is prerequisite.     

Monte Carlo Study of Structural and Thermodynamic Properties of Liquid Chloroform Using a Five Site Model

A five site potential model combining Lennard–Jones plus Coulomb potential functions has been developed for chloroform molecule. The partial charges needed for Coulombic interactions were derived using the chelpg procedure implemented in the gaussian 92 program. These calculations were performed at the MP2 level with MC-311G* basis set for Cl and 6-311G** for C and H atoms. The parameters for the Lennard–Jones potentials were optimized to reproduce experimental values for the density and enthalpy of vaporization of the pure liquid at 298 K and 1 atm. The statistical mechanics calculations were performed with the Monte Carlo method in the isothermic and isobaric (NpT) ensemble. Besides the values obtained for density, ρ, and molar enthalpy of vaporization at constant pressure, Δ H_V, for liquid chloroform, results for molar volume, V_m, molar heat capacity, C_p, isobaric thermal expansivity, α_p, and isothermal compressibility, κ_T, for this pure liquid are also in very good agreement with experimental observations. Size effects on the values of thermodynamic properties were investigated. The potential model was also tested by computing the free energy for solvating one chloroform molecule into its own liquid at 298 K using a statistical perturbation approach. The result obtained compares well with the experimental value. Site–site pair correlation functions were calculated and are in good accordance with theoretical results available in the literature. Dipole–dipole correlation functions for the present five site model were also calculated at different carbon–carbon distances. These correlations were compared to those obtained using the four site model reported in the literature. An investigation of the solvent dependence of the relative free energy for cis/trans conversion of a hypothetical solute in TIP4P water and chloroform was accomplished. The results show strong interaction of water and chloroform molecules with the gauche conformer. The value obtained for the free energy barrier for cis/trans rotation in TIP4P water is higher than that for chloroform. This result is in agreement with the continuous theory for solvation as the conformer with higher dipole moment is more favoured by the solvent with higher dieletric constant. The results also show an increase in entropy as the solute goes from the cis to the trans geometry and this result is more appreciable in the aqueous solution. This revised version was published online in June 2006 with corrections to the Cover Date.     

中国若干森林水文要素地理分布规律的模拟

利用全国各森林植被地理分布区水文要素的实验观测数据,应用多元统计回归和主分量分析等数学方法建立了描述中国森林水文要素与地理要素（经度、纬度、海拔）之间函数关系Ｆ的三维地学模型（ＤＴＭ）,其中包括：全国年均降水分布规律模拟模型;全国林冠降水截留规律模拟模型;全国森林枯落物持水能力模拟模型;中国森林综合截留能力指数模型。并且利用地理信息系统的三维地学模型（ＤＴＭ）成图功能,实现全国范围内的森林水文要素     

一种模拟昆虫种群动态的改进的变维矩阵模型

提出了一种模拟昆虫种群动态的改进的变维矩阵模型,该模型以发有历期为维数,采用分解与合成的方法变维,并考虑了个体间的发育差异。经模拟检验,模型模拟结果略优于徐汝梅等（1981）变维矩阵模型的结果。