Metastability and exponential approach to equilibrium for low-temperature stochastic Ising models

We present a proof of the exponential convergence to equilibrium of single-spin-flip stochastic dynamics for the two-dimensional Ising ferromagnet in the low-temperature case with not too small external magnetic fieldh uniformly in the volume and in the boundary conditions.     

Stochastic modeling of a billiard in a gravitational field: Power law behavior of Lyapunov exponents

We consider the motion of a point particle (billiard) in a uniform gravitational field constrained to move in a symmetric wedge-shaped region. The billiard is reflected at the wedge boundary. The phase space of the system naturally divides itself into two regions in which the tangent maps are respectively parabolic and hyperbolic. It is known that the system is integrable for two values of the wedge half-angle ₁ and ₂ and chaotic for ₁<< ₂. We study the system at three levels of approximation: first, where the deterministic dynamics is replaced by a random evolution; second, where, in addition, the tangent map in each region is, replaced by its average; and third, where the tangent map is replaced by a single global average. We show that at all three levels the Lyapunov exponent exhibits power law behavior near ₁ and ₂ with exponents 1/2 and 1, respectively. We indicate the origin of the exponent 1, which has not been observed in unaccelerated billiards.     

Computer simulation of some dynamical properties of the Lorentz gas

We carried out molecular dynamics simulations of a Lorentz gas, consisting of a lone hydrogen molecule moving in a sea of stationary argon atoms. A Lennard-Jones form was assumed for the H₂-Ar potential. The calculations were performed at a reduced temperatureK ^* =kT/_{H 2–Ar} = 4.64 and at reduced densities ^*= _{Ar Ar} ³ in the range 0.074–0.414. The placement of Ar atoms was assumed to be random rather than dictated by equilibrium considerations. We followed the trajectories of many H₂ molecules, each of which is assigned in turn a velocity given by the Maxwell-Boltzmann distribution at the temperature of the simulation. Solving the equations of motion classically, we obtained the translational part of the incoherent dynamic structure factor for the H₂ molecule,S _tr(q, ). This was convoluted with the rotational structure factorS _rot(q, ) calculated assuming unhindered rotation to obtain the total structure factorS(q, ). Our results agree well with experimental data on this function obtained by Egelstaffet al. At the highest density ( ^*=0.414) we studied the dependence ofS(q, ) on system size (number of Ar atoms), number of H₂ molecules for which trajectories are generated, and the length of time over which these trajectories are followed.     

Dynamics of the current-driven Josephson junction

The irreversible macroscopic dynamics of the Josephson junction coupled to external wires acting as a current source is derived rigorously from the underlying microscopic Hamiltonian quantum mechanics. The external systems are treated in the singular coupling limit. The use of this limit is explicitly justified via an interpretation of the singular coupling limit in terms of the relative magnitudes of system, reservoir, and coupling energies. The qualitative behavior of the macroscopic dynamical equations is shown to depend sensitively and crucially on the interaction between the wires and the superconductors and on the size of the wires: the dc Josephson effect only happens when one lets Cooper pairs be driven into the junction by collective (i.e., small) reservoirs.     

Ar—Kr溶液扩散系数的分子动力学模拟及其与温度的关系

用分子动力学模拟方法研究确定Ａｒ－Ｋｒ溶液的自扩散系数Ｄ１、Ｄ２和互扩散系数Ｄ１２以及它们随温度变化的规律。结果表明，分别用Ｇｒｅｅｎ－Ｋｕｂｏ法和Ｅｉｎｓｔｅｉｎ法得到的扩散系数在数值上一致；该溶液的３种扩散系数均满足Ｄ＝Ｄ０ｅ＾Ｅ／ＲＴ关系。     

液体和非晶态Ni64B36合金结构的从头算分子动力学模拟

运用从头算分子动力学模拟了液体以及猝冷后形成的非晶态Ni₆₄B₃₆合金体系, 得到了它们的对相关函数、结构因子、键对分析方面的结构信息, 与实验结果相当一致; 结果表明, 猝冷得到的合金性质与液体合金性质相似, 为非晶态结构. B原子多数以B—B双原子成键形式分散于Ni原子构成的骨架中. 电子态密度分析表明, Ni 3d电子最活泼, 因此在合金中Ni为活性位. 轨道电荷分析从电子结构角度揭示了在NiB 催化剂中B作为修饰剂的机理.     

Local elevation: A method for improving the searching properties of molecular dynamics simulation

Summary The concept of memory has been introduced into a molecular dynamics algorithm. This was done so as to persuade a molecular system to visit new areas of conformational space rather than be confined to a small number of low-energy regions. The method is demonstrated on a simple model system and the 11-residue cyclic peptide cyclosporin A. For comparison, calculations were also performed using simulated temperature annealing and a potential energy annealing scheme. Although the method can only be applied to systems with a small number of degrees of freedom, it offers the chance to generate a multitude of different low-energy structures, where other methods only give a single one or few. This is clearly important in problems such as drug design, where one is interested in the conformational spread of a system.     

电化学斑图中空间耦合的时空动力学

斑图是在空间或时间上具有某种规律性的非均匀宏观结构,普遍存在于自然界,事实上大多数电化学体系在适当参数范围内也可以形成各种电化学斑图。本文介绍了电化学斑图中空间耦合的三种方式,重点总结了在不同机理作用下可能出现的斑图种类,并简述了电化学斑图模拟方面的发展。     

Structure-activity relationships of pyrethroid insecticides. Part 2. The use of molecular dynamics for conformation searching and average parameter calculation

Summary Molecular dynamics simulations have been performed on a number of conformationally flexible pyrethroid insecticides. The results indicate that molecular dynamics is a suitable tool for conformational searching of small molecules given suitable simulation parameters. The structures derived from the simulations are compared with the static conformation used in a previous study. Various physicochemical parameters have been calculated for a set of conformations selected from the simulations using multivariate analysis. The averaged values of the parameters over the selected set (and the factors derived from them) are compared with the single conformation values used in the previous study.     

Pu（^7Fg）＋H2（X^1∑g^＋,0,0）的分子反应动力学

基于PuH_2分子基态(X~7A_1)的分析势能函数，用准经典的Monte-Carlo轨线法 对Pu(~7Fg)+H_2(X~1∑_g~+，0，0)的分子反应动力学过程进行了计算。结果表明 ：Pu(~7F_g)与H_2(X~1∑_g~+，0，0)碰撞是弹性碰撞。