Surface dynamics, simulation

Small metal clusters are also of interest because of their importance in catalysis. Despite the fact that small clusters should consist of mostly surface atoms, measurement of the photon ionization threshold for Hg clusters suggest that a transition from van der Waals to metallic properties occurs in the range of 20-70 atoms per cluster [88] and near-bulk magnetic properties are expected for Ni, Pd, and Pt clusters of only 13 atoms [89] Theoretical calculations on Sin and other semiconductors predict that the stmcture reflects the bulk lattice for 1000 atoms but the bulk electronic wave functions are not obtained [90]. Bartell and co-workers [91] study beams of molecular clusters with electron dirfraction and molecular dynamics simulations and find new phases not observed in the bulk. Bulk models appear to be valid for their clusters of several thousand atoms (see Section IX-3). 

Alejandre J, Tildesley D J and Chapela G A 1995 Molecular dynamics simulation of the orthobaric densities and surface tension of water J. Chem. Phys. 102 4574-83... 

Weakliem P C and Carter E A 1993 Surface chemical reactions studied via ab /n/f/o-derived molecular dynamics simulations fluorine etching of Si(IOO) J. Chem Phys. 98 737-45... 

Helmer B A and Graves D B 1998 Molecular dynamics simulations of Ar" and Cl" Impacts onto silicon surfaces distributions of reflected energies and angles J. Vac. Sc/. Technol. A 16 3503-14... 

Juffer, A.H., Berendsen, H.J.C. Dynamic surface boundary conditions A simple boundary model for molecular dynamics simulations. Mol. Phys. 79 (1993) 623-644. 

The correct treatment of boundaries and boundary effects is crucial to simulation methods because it enables macroscopic properties to be calculated from simulations using relatively small numbers of particles. The importance of boundary effects can be illustrated by considering the following simple example. Suppose we have a cube of volume 1 litre which is filled with water at room temperature. The cube contains approximately 3.3 X 10 molecules. Interactions with the walls can extend up to 10 molecular diameters into the fluid. The diameter of the water molecule is approximately 2.8 A and so the number of water molecules that are interacting with the boundary is about 2 x 10. So only about one in 1.5 million water molecules is influenced by interactions with the walls of the container. The number of particles in a Monte Carlo or molecular dynamics simulation is far fewer than 10 -10 and is frequently less than 1000. In a system of 1000 water molecules most, if not all of them, would be within the influence of the walls of the boundary. Clecirly, a simulation of 1000 water molecules in a vessel would not be an appropriate way to derive bulk properties. The alternative is to dispense with the container altogether. Now, approximately three-quarters of the molecules would be at the surface of the sample rather than being in the bulk. Such a situation would be relevcUit to studies of liquid drops, but not to studies of bulk phenomena. 

Fig. 9.19 Schematic illustration of an energy surface. A high-temperature molecular dynamics simulation may be ah to ooercome very high energy barriers and so explore conformational space. On minimisation, the appropriate minimum energy conformation is obtained (arrcrws).

To calculate the properties of a molecule, you need to generate a well-defined structure. A calculation often requires a structure that represents a minimum on a potential energy surface. HyperChem contains several geometry optimizers to do this. You can then calculate single point properties of a molecule or use the optimized structure as a starting point for subsequent calculations, such as molecular dynamics simulations. 

The relative intensities of the bands in the transmission and RAIR spectra were used to determine the orientation of the long axis of the 4-MPP molecules with respect to the normal to the gold surface. It was found that this tilt angle was about 21°, a value that was similar to that obtained from molecular dynamics simulations [11]. 

In thermal building-dynamics simulation codes, outdoor conditions are mostly input by the so-called weather data file, containing (usually hourly) data for air temperature, wind speed and direction, air humidity, and global and diffuse solar radiation on horizontal surfaces. 

In order to run a dynamic simulation, weather data files are needed, providing (usually hourly) data on the following quantities air temperature, air humidity, ground temperature, global and diffuse irradiation on horizontal surfaces, and wind speed and direction. 

T. Matsuda, G. D. Smith, R. G. Winkler, D. Y. Yoon. Stochastic dynamics simulations of n-alkane melts confined between solid surfaces Influence of surface properties and comparison with Schetjens-Fleer theory. Macromolecules 28 65- 13, 1995. 

Importantly, all biological procedures are operating at a temperature of 310 Kelvin, not at 0 Kelvin as the potential energy is calculated by the force fields. The kinetic energy must also be considered. Molecules and proteins at room temperature change the conformation at least at the surface and in loop region. Molecular dynamics simulation (MD) is an approach to tackle these kinetic and stability problems. 

---