Speed of adaptation in structured populations

Adaptation of populations takes place with the occurrence and subsequent fixation of mutations that confer some selective advantage to the individuals which acquire it. For this reason, the study of the process of fixation of advantageous mutations has a long history in the population genetics literature. Particularly, the previous investigations aimed to find out the main evolutionary forces affecting the strength of natural selection in the populations. In the current work, we investigate the dynamics of fixation of beneficial mutations in a subdivided population. The subpopulations (demes) can exchange migrants among their neighbors, in a migration network which is assumed to have either a random graph or a scale-free topology. We have observed that the migration rate drastically affects the dynamics of mutation fixation, despite of the fact that the probability of fixation is invariant on the migration rate, accordingly to Maruyama's conjecture. In addition, we have noticed a topological dependence of the adaptive evolution of the population when clonal interference becomes effective.     

热动力学问题的数值计算

 对含热传导的流体动力学方程组，用有限元方法进行数值求解。采用傅里叶热传导计算热流、用热流连续条件计算单元间接触面的温度、用三角形传热法计算体单元表面的热流，考虑各向同性弹塑性流体材料模型、三项式物态方程和导热系数与状态的相关性，给出了傅里叶热传导、接触传热、热应力应变效应、以及混合物冲击压缩特性等算例。对混合物冲击温度的数值模拟表明，小颗粒混合物在冲击压缩过程中，颗粒间的温度有差别、稍有波动，并随时间趋向于一致，以至热平衡。     

自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。     

Thermal helium clusters at 3.2 Kelvin in classical and semiclassical simulations

The thermodynamic stability of⁴He_4–13 at 3.2 K is investigated with the classical Monte Carlo method, with the semiclassical path-integral Monte Carlo (PIMC) method, and with the semiclassical all-order many-body method. In the all-order many-body simulation the dipole-dipole approximation including short-range correction is used. The resulting stability plots are discussed and related to recent TOF experiments by Stephens and King. It is found that with classical Monte Carlo of course the characteristics of the measured mass spectrum cannot be resolved. With PIMC, switching on more and more quantum mechanics. by raising the number of virtual time steps results in more structure in the stability plot, but this did not lead to sufficient agreement with the TOF experiment. Only the all-order many-body method resolved the characteristic structures of the measured mass spectrum, including magic numbers. The result shows the influence of quantum statistics and quantum mechanics on the stability of small neutral helium clusters.     

Electron-stimulated desorption of D (H) from condensed D2O (H2O) films

The electron-stimulated desorption (ESD) of D⁻ and H⁻ ions from condensed D₂O and H₂O films is investigated. Three low-energy peaks are observed in the ESD anion yield, which are identified as arising from excitation of ²B₁, ²A₁ and ²B₂ dissociative electron attachment (DEA) resonances. Additional structure is observed between 18 and 32 eV, which may be due to ion pair formation or to DEA resonances involving the 2a₁ orbital. The ion yield resulting from excitation of the ²B₁ resonance increases as the film is heated. We attribute the increase in the ion yield to thermally induced hydrogen bond breaking near the surface, which enhances the lifetimes of the excited states that lead to desorption.     

A RAPID APPROXIMATION FOR PREDICTING THE HYDROCARBON DISCOVERY RATE: PART 1. ASSESSING THE ACCURACY

The hydrocarbon discovery prediction problem is important to firms having to make decisions about the deployment of scarce exploration resources. Traditional methods for estimating the discovery rate rely on the completion of time consuming simulation experiments. A rapid approximation that does not require the completion of simulation exists and has been shown to have some promise as a prediction tool. This paper investigates the accuracy of the approximation method under a wide variety of distributional and drilling efficiency assumptions. The results indicate that the approximation produces predictions close to those of simulation under most of the tested conditions. This suggests that resource exploration firms could conveniently use the method for a wide variety of planning purposes without incurring the same costs in time and personnel required for simulation.     

微裂纹扩展中的分形现象及其计算机仿真

本文探讨了采用仿真技术研究微裂纹连接过程和生成分形断口表面的可行性，发展了用于仿真微裂纹随机生长行为的ＧＯ模型的基本构相，并进行了计算机仿真实验，由此得出了一此有的结论。     

曙光一号自由电子激光器的理论计算

系统总结曙光一号自由电子激光器理论计算的主要结果:包括曙光一号装置主要参数的选取和理解;磁场失谐曲线的计算;常参数摇摆器和变参数摇摆器的主要结果;高阶波导模的贡献;电子束参数扰动对激光性能的影响;空间电荷效应等。计算结果表明,常参数摇摆器激光输出功率可达80MW,效率约50％;变参数摇摆器激光输出功率可达250MW左右,效率约16％。     

高功率振荡器锁相的粒子模拟

用粒子(PIC,Particle-in-Cell)模拟方法研究了当前颇感兴趣的高功率振荡器的锁相问题。给出了模拟计算结果。讨论了锁相(Phase Locking)、Priming和InjectionSeeding等概念的区别。     

三维射频加速腔程序包MAFIA的开发和应用

介绍了三维射频加速器计算机辅助设计(CAD)程序包MAFIA的特点、功能及其应用。用CAD的方法实现了射频加速器特别是非轴对称带耦合腔的加速器的三维腔形设计。同时介绍了节省内存、节约机时的模型方法,使得在中小型机器上运行大程序(八万条)成为可能。文章主要对L波段的加速腔进行了数值模拟,并和实验进行了比较,结果符合较好;通过三维数值模拟及图形模拟的方法比较了各种形状谐振腔的性能;计算分析了输入耦合腔及高阶模耦合器对加速腔电磁场分布等参量的影响;提出了加耦合器后谐振频率改变的补偿方法。