Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated. We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard-Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Phase transitions from the solid state of monatomic elements

The principle aims of this thesis include the development of models of sublimation and melting from first principles and the application of these models to the rare gases.

A simple physical model is constructed to represent the sublimation of monatomic elements. According to this model, the solid and gas phases are two states of a single physical system. The nature of the phase transition is clearly revealed, and the relations between the vapor pressure, the latent heat, and the transition temperature are derived. The resulting theory is applied to argon, krypton, and xenon, and good agreement with experiment is found.

For the melting transition, the solid is represented by an anharmonic model and the liquid is described by the Percus-Yevick approximation. The behavior of the liquid at high densities is studied on the isotherms kT/∈ = 1.3, 1.8, and 2.0, where k is Boltzmann's constant, T is the temperature, and e is the well depth of the Lennard-Jones 12-6 pair potential. No solutions of the PercusYevick equation were found for ρσ³ above 1.3, where ρ is the particle density and σ is the radial parameter of the Lennard-Jones potential. The liquid structure is found to be very different from the solid structure near the melting line. The liquid pressures are about 50 percent low for experimental melting densities of argon. This discrepancy gives rise to melting pressures up to twice the experimental values.

伦纳德琼斯薄膜材料导热系数的分子动力学研究

提出从微观的角度,借助计算机工具,将薄膜破坏发展的细节展现出来的分子动力学研究的思想。使得实验上难以观察的现象变得形象而便于理解。应用分子动力学理论,使用伦纳德琼斯势函数,采用预校正积分法和虚拟外力约束标定方法,模拟薄膜体系的传热系数受体系的密度、温度的影响,同时结合体系粒子的径向分布函数和长程分布函数分析了相应的系统结构特性。另外,采用不同的模拟尺寸获得了低维材料所特有的“高温尺寸效应”。结果显示,导热系数随密度的增加变大,随温度的上升而变大。这些数据现有测量手段是难以得到的,这类模拟可以为研究提供一些

Efficient sampling of atomic configurational spaces.

We describe a method to explore the configurational phase space of chemical systems. It is based on the nested sampling algorithm recently proposed by Skilling (AIP Conf. Proc. 2004, 395; J. Bayesian Anal. 2006, 1, 833) and allows us to explore the entire potential energy surface (PES) efficiently in an unbiased way. The algorithm has two parameters which directly control the trade-off between the resolution with which the space is explored and the computational cost. We demonstrate the use of nested sampling on Lennard-Jones (LJ) clusters. Nested sampling provides a straightforward approximation for the partition function; thus, evaluating expectation values of arbitrary smooth operators at arbitrary temperatures becomes a simple postprocessing step. Access to absolute free energies allows us to determine the temperature-density phase diagram for LJ cluster stability. Even for relatively small clusters, the efficiency gain over parallel tempering in calculating the heat capacity is an order of magnitude or more. Furthermore, by analyzing the topology of the resulting samples, we are able to visualize the PES in a new and illuminating way. We identify a discretely valued order parameter with basins and suprabasins of the PES, allowing a straightforward and unambiguous definition of macroscopic states of an atomistic system and the evaluation of the associated free energies.

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.

高温高压气体状态方程研究及钱学森方程改进

钱学森根据分子动理论的普遍原则,结合Lennard-Jones Devanshire(LJD)液体理论的结果,给出了一个适用于高温高压气体的普适状态方程.该方程有物理基础,且具有计算精度高、形式简单及显含温度的优点.本文用分子动力学方法检验了LJD理论的适用性,并给出精确的高温高压气体状态方程分子动力学数值解,在此基础上修订了钱学森状态方程.改进后的钱学森状态方程在高密度范围内与MD结果平均绝对偏差率小于10%.

原子力显微镜Tapping模式探针动力学与操作方法研究

自从1986年G. Binnig 等发明原子力显微镜(AFM)以来，纳米科技研究得到了快速发展，纳米科技研究的最终目标是从纳米甚至分子与原子尺度上制造功能器件或系统，而实现此目标的手段之一是研制能够在纳米尺度上进行精确可控操作的装置、方法与技术。对此，学者们对基于 AFM 的纳米操作进行了大量研究。基于AFM-Tapping模式的纳米观测与操作具有成本低、易实现性等特点，因而近年来成为纳米技术研究的热点。但基于AFM探针的纳米观测和操作的基本理论仍然很不完善，许多情况下缺乏基于探针观测与操作的纳观作用力的合理解释与描述，因此限制了基于探针的纳米技术发展。为了解决这一问题，本论文在国家自然科学基金项目“基于机器人化的纳米操作方法研究(60575060)”的支持下，广泛收集和查阅国内外相关资料，重点开展了AFM-Tapping模式下探针动力学解释方法与建模研究，并进行了基于AFM-Tapping的纳米操作方法研究。主要研究工作如下： 1)针对探针动力学建模问题，开展了基于Hamilton原理和变分原理的有边界条件的探针一维悬臂梁建模方法研究。利用模态叠加法和模态函数正交性对探针各阶模态进行了解耦分析，得到了解耦的探针动力学方程。 2)开展了纳观尺度下探针-样品作用力研究。分析与阐述了毛细力、接触力和摩擦力等纳观尺度作用力；开展了基于Lennard-Jones势模型的不同形状探针与样品间的相互作用势描述方法研究；并根据势能-力之间的关系，研究了可描述半球形探针与样品间相互作用力的数学表达式。 3)开展了基于弹簧振子模型的纳观非线性力作用下探针振动问题分析研究。首次提出了基于探针弹性力与非线性纳观力平衡关系的探针振动状态区域划分方法。在探针振动能量平衡分析的基础上，研究了AFM-Tapping模式中的双稳态现象，进而得到了描述双稳态现象的解析表达式。最后采用Lyapunov指数分析的方法，证明了探针在非线性振动中确实存在着混沌现象。 4)研究了基于AFM-Tapping模式的纳米操作方法。通过研究探针振幅调节控制和预编程作业轨迹规划策略，提出了一种具有增强探针概念的纳米作业方法。该方法可以在参数设定和预编程条件下自动完成纳米尺度的平移和旋转操作，从而可以大大提高纳米操作的效率和精度。实验结果验证了该方法的有效性。本文的工作为纳米操作方法的研究提供了可以借鉴的理论与实验经验，有助于推动基于AFM的纳米操作与制造技术的发展。

Equilibrium polymerization of cyclic carbonate oligomers. II. Role of multiple active sites

Ring opening polymerization of bisphenol A polycarbonate is studied by Monte Carlo simulations of a model comprising a fixed number of Lennard-Jones particles and harmonic bonds [J. Chem. Phys. 115, 3895 (2001)]. Bond interchanges produced by a low concentration (0.10%less than or equal toc(a)less than or equal to0.36%) of chemically active particles lead to equilibrium polymerization. There is a continuous transition in both 2D and 3D from unpolymerized cyclic oligomers at low density to a system of linear chains at high density, and the polymeric phase is much more stable in three dimensions than in two. The steepness of the polymerization transition increases rapidly as c(a) decreases, suggesting that it is discontinuous in the limit c(a)-->0. The transition is entropy driven, since the average potential energy increases systematically upon polymerization, and there is a steady decline in the degree of polymerization as the temperature is lowered. The mass distribution functions for open chains and for rings are unimodal, with exponentially decaying tails that can be fitted by Zimm-Schulz functions and simpler exponential forms. (C) 2002 American Institute of Physics.

Simulation of interfaces between room temperature ionic liquids and other liquids

The structure and properties of the interfaces between the room temperature ionic liquid dimethylimidazolium chloride ([dmim]Cl) and different Lennard-Jones fluids and between ionic liquid and water have been studied by molecular dynamics simulations, and compared to the ionic liquid-vapour interface. Two contrasting types of interface were investigated, thermodynamically stable interfaces between ionic liquid and vapour and between ionic liquid and Lennard-Jones fluids, and diffusing interfaces between miscible phases of different compositions involving water. The density profiles of different species through the interface are presented. The cations and water molecules near the former type of interface are aligned relative to the surface, but no orientational preference was found near or in the broad diffusing interface. The ionic liquid has a negative electrostatic potential relative to vapour or Lennard-Jones fluid, but is more positive than pure water. This contrast is explained in terms of the relative importance of orientation and concentration differences in the two types of interface.

Demonstrating the Transferability and the Descriptive Power of Sketch-Map

Increasingly, it is recognized that new automated forms of analysis are required to understand the high-dimensional output obtained from atomistic simulations. Recently, we introduced a new dimensionality reduction algorithm, sketch-map, that was designed specifically to work with data from molecular dynamics trajectories. In what follows, we provide more details on how this algorithm works and on how to set its parameters. We also test it on two well-studied Lennard-Jones clusters and show that the coordinates we extract using this algorithm are extremely robust. In particular, we demonstrate that the coordinates constructed for one particular Lennard-Jones cluster can be used to describe the configurations adopted by a second, different cluster and even to tell apart different phases of bulk Lennard-Jonesium.

Exploring the free energy surfaces of clusters using reconnaissance metadynamics

A new approach is proposed for exploring the low-energy structures of small to medium-sized aggregates of atoms and molecules. This approach uses the recently proposed reconnaissance metadynamics method [G. A. Tribello, M. Ceriotti, and M. Parrinello. Proc. Natl. Acad. Sci. U.S.A. 107(41), 17509 (2010)] in tandem with collective variables that describe the average structure of the coordination sphere around the atoms/molecules. We demonstrate this method on both Lennard-Jones and water clusters and show how it is able to quickly find the global minimum in the potential energy surface, while exploring the finite temperature free energy surface. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3628676]

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

An application of cell-cluster theory to a rare gas crystal /|n[by] K. Westera. - 260 St. Catharines [Ont.] : Dept. of Physics, Brock University,

The Lennard-Jones Devonshire 1 (LJD) single particle theory for liquids is extended and applied to the anharmonic solid in a high temperature limit. The exact free energy for the crystal is expressed as a convergent series of terms involving larger and larger sets of contiguous particles called cell-clusters. The motions of all the particles within cell-clusters are correlated to each other and lead to non-trivial integrals of orders 3, 6, 9, ... 3N. For the first time the six dimensional integral has been calculated to high accuracy using a Lennard-Jones (6-12) pair interaction between nearest neighbours only for the f.c.c. lattice. The thermodynamic properties predicted by this model agree well with experimental results for solid Xenon.