Dynamic lattice model

The Dynamic Lattice Model of Surfaces. The Experimental Facts. Let us now review what is known about the structure of clean surfaces without adsorbates. There are three phenomena that have been identified relaxation, reconstruction and the presence of steps and kinks. 

There is, of course, a mass of rather direct evidence on orientation at the liquid-vapor interface, much of which is at least implicit in this chapter and in Chapter IV. The methods of statistical mechanics are applicable to the calculation of surface orientation of assymmetric molecules, usually by introducing an angular dependence to the inter-molecular potential function (see Refs. 67, 68, 77 as examples). Widom has applied a mean-held approximation to a lattice model to predict the tendency of AB molecules to adsorb and orient perpendicular to the interface between phases of AA and BB [78]. In the case of water, a molecular dynamics calculation concluded that the surface dipole density corresponded to a tendency for surface-OH groups to point toward the vapor phase [79]. 

The parameter /r tunes the stiffness of the potential. It is chosen such that the repulsive part of the Leimard-Jones potential makes a crossing of bonds highly improbable (e.g., k= 30). This off-lattice model has a rather realistic equation of state and reproduces many experimental features of polymer solutions. Due to the attractive interactions the model exhibits a liquid-vapour coexistence, and an isolated chain undergoes a transition from a self-avoiding walk at high temperatures to a collapsed globule at low temperatures. Since all interactions are continuous, the model is tractable by Monte Carlo simulations as well as by molecular dynamics. Generalizations of the Leimard-Jones potential to anisotropic pair interactions are available e.g., the Gay-Beme potential [29]. This latter potential has been employed to study non-spherical particles that possibly fomi liquid crystalline phases. 

Despite their simplicity, certainly compared to the all-atom potentials used in molecular dynamics studies, these contact energy functions enable the exploration of different interaction scenarios. This diversity is achieved by changing the heterogeneity of the sequence, by altering the number N of different types of residues that are being used. The most elementary lattice model involves only two types of monomers hydrophobic... 

The function / incorporates the screening effect of the surfactant, and is the surfactant density. The exponent x can be derived from the observation that the total interface area at late times should be proportional to p. In two dimensions, this implies R t) oc 1/ps and hence x = /n. The scaling form (20) was found to describe consistently data from Langevin simulations of systems with conserved order parameter (with n = 1/3) [217], systems which evolve according to hydrodynamic equations (with n = 1/2) [218], and also data from molecular dynamics of a microscopic off-lattice model (with n= 1/2) [155]. The data collapse has not been quite as good in Langevin simulations which include thermal noise [218]. 

These apparent restrictions in size and length of simulation time of the fully quantum-mechanical methods or molecular-dynamics methods with continuous degrees of freedom in real space are the basic reason why the direct simulation of lattice models of the Ising type or of solid-on-solid type is still the most popular technique to simulate crystal growth processes. Consequently, a substantial part of this article will deal with scientific problems on those time and length scales which are simultaneously accessible by the experimental STM methods on one hand and by Monte Carlo lattice simulations on the other hand. Even these methods, however, are too microscopic to incorporate the boundary conditions from the laboratory set-up into the models in a reahstic way. Therefore one uses phenomenological models of the phase-field or sharp-interface type, and finally even finite-element methods, to treat the diffusion transport and hydrodynamic convections which control a reahstic crystal growth process from the melt on an industrial scale. 

The Monte Carlo method as described so far is useful to evaluate equilibrium properties but says nothing about the time evolution of the system. However, it is in some cases possible to construct a Monte Carlo algorithm that allows the simulated system to evolve like a physical system. This is the case when the dynamics can be described as thermally activated processes, such as adsorption, desorption, and diffusion. Since these processes are particularly well defined in the case of lattice models, these are particularly well suited for this approach. The foundations of dynamical Monte Carlo (DMC) or kinetic Monte Carlo (KMC) simulations have been discussed by Eichthom and Weinberg (1991) in terms of the theory of Poisson processes. The main idea is that the rate of each process that may eventually occur on the surface can be described by an equation of the Arrhenius type ... 

KovalyovEV, Resnyanskii ED, Elokhin VI, Bal zhinimaev BS, Myshlyavtsev AV. 2003. Novel statistical lattice model for the supported nanoparticle. Features of the reaction performance influenced by the dynamically changed shape and surface morphology of the supported active particle. Phys Chem Chem Phys 5 784-790. 

Coarse-grained polymer models neglect the chemical detail of a specific polymer chain and include only excluded volume and topology (chain connectivity) as the properties determining universal behavior of polymers. They can be formulated for the continuum (off-lattice) as well as for a lattice. For all coarse-grained models, the repeat unit or monomer unit represents a section of a chemically realistic chain. MD techniques are employed to study dynamics with off-lattice models, whereas MC techniques are used for the lattice models and for efficient equilibration of the continuum models.36 2 A tutorial on coarse-grained modeling can be found in this book series.43... 

Since none of the lattice models is now clearly superior, the choice for interpretation of spin relaxation in polymers is arbitrary. Familiarity leads us to select the Jones and Stockmayer model so we will now consider application of this model to several well studied polymer systems in order to compare dynamics from polymer to polymer. Also the equations required to consider anisotropic Internal rotation of substituent groups and overall molecular tumbling as independent motions in addition to backbone rearrangements caused by the three-bond jump are available for the Jones and Stockmayer model (13). 

A. Jones Yes, I didn t show all the data. Most of the data is from the literature. Polystyrene data and polyisobutylene data at two field strengths can be accounted for by these lattice models based upon the three bond jump. Some recent data by Dr. Bovey on the polybutene was more difficult to fit with regard to the frequency dependence. (F. C. Schilling, R. E. Cais and F. A. Bovey, Macromolecules ]A, 325 (19T8).) Frequency dependent data I think is important information to acquire when trying to understand the dynamics. 

These two methods are different and are usually employed to calculate different properties. Molecular dynamics has a time-dependent component, and is better at calculating transport properties, such as viscosity, heat conductivity, and difftisivity. Monte Carlo methods do not contain information on kinetic energy. It is used more in the lattice model of polymers, protein stmcture conformation, and in the Gibbs ensemble for phase equilibrium. 

C.J. Perez, A. Corral, A. Diaz Guilera, K. Christensen, and A. Arenas. On self-organized criticality and synchronization in lattice models of coupled dynamical systems, http //xxx.lanl.gov, paperno. cond-mat/9601102, 1996. 

There has been a long history in theoretical efforts to understand H + H/Cu(lll) and its isotopic analogs because it represents the best studied prototype of an ER/HA reaction. These have evolved from simple 2D collinear quantum dynamics on model PES [386] to 6D quasi-classical dynamics on PES fit to DFT calculations [380,387,388], and even attempts to include lattice motion on ER/HA reactions [389]. These studies show that there is little reflection of incident H because of the deep well and energy scrambling upon impact, i.e., a % 1. Although some of the... 

---