Molecular dynamics atomic liquids

Electrons residing in molecular clusters can be viewed as microscopic probes of both the local liquid structure and the molecular dynamics of liquids, and as such their transitory existence becomes a theoretical and experimental metaphor for one of the major fundamental and contemporary problems in chemical and molecular physics, that is, how to describe the transition between the microscopic and macroscopic realms of physical laws in the condensed phase. Since this chapter was completed in the Spring of 1979, several new and important observations have been made on the dynamics and structure of e, which, as a fundamental particle interacting with atoms and molecules in a fundamental way, serves to assist that transformation for electronic states in disordered systems. In a sense, disorder has become order on the subpicosecond time-scale, as we study events whose time duration is shorter than, or comparable to, the period during which the atoms or molecules retain some memory of the initial quantum state, or of the velocity or phase space correlations of the microscopic system. This approach anticipated the new wave of theoretical and experimental interest in developing microscopic theories of... 

In Fig. III-7 we show a molecular dynamics computation for the density profile and pressure difference P - p across the interface of an argonlike system [66] (see also Refs. 67, 68 and citations therein). Similar calculations have been made of 5 in Eq. III-20 [69, 70]. Monte Carlo calculations of the density profile of the vapor-liquid interface of magnesium how stratification penetrating about three atomic diameters into the liquid [71]. Experimental measurement of the transverse structure of the vapor-liquid interface of mercury and gallium showed structures that were indistinguishable from that of the bulk fluids [72, 73]. 

It was noted in connection with Eq. III-56 that molecular dynamics calculations can be made for a liquid mixture of rare gas-like atoms to obtain surface tension versus composition. The same calculation also gives the variation of density for each species across the interface [88], as illustrated in Fig. Ill-13b. The density profiles allow a calculation, of course, of the surface excess quantities. 

It is possible to use the quantum states to predict the electronic properties of the melt. A typical procedure is to implement molecular dynamics simulations for the liquid, which pemiit the wavefiinctions to be detemiined at each time step of the simulation. As an example, one can use the eigenpairs for a given atomic configuration to calculate the optical conductivity. The real part of tire conductivity can be expressed as... 

Figure B3.3.12. Sulphur atoms in liquid iron at the Earth s core conditions, simnlated by first-principle Car-Parrinello molecular dynamics, (a) Initial conditions, showing a mannally-prepared initial cluster of snlphur atoms, (b) A short tune later, indicating spontaneous dispersal of the snlphur atoms, which mingle with the surroundmg iron atoms. Thanks are dne to D Alfe and M J Gillan for this figure. For fiirtlier details see [210. 211].

Further support for this approach is provided by modern computer studies of molecular dynamics, which show that much smaller translations than the average inter-nuclear distance play an important role in liquid state atom movement. These observations have conhrmed Swalin s approach to liquid state diffusion as being very similar to the calculation of the Brownian motion of suspended particles in a liquid. The classical analysis for this phenomenon was based on the assumption that the resistance to movement of suspended particles in a liquid could be calculated by using the viscosity as the frictional force in the Stokes equation... 

The only feasible procedure at the moment is molecular dynamics computer simulation, which can be used since most systems are currently essentially controlled by classical dynamics even though the intermolecular potentials are often quantum mechanical in origin. There are indeed many intermolecular potentials available which are remarkably reliable for most liquids, and even for liquid mixtures, of scientific and technical importance. However potentials for the design of membranes and of the interaction of fluid molecules with membranes on the atomic scale are less well developed. 

To make evaluations more definite, we use optical and microwave experimental data, as well as calculations of molecular dynamics of certain simple liquids which usually fit the experiment. Rotation is everywhere considered as classical, and the objects are two-atomic and spherical molecules, as well as hard ellipsoids. 

Force fields split naturally into two main classes all-atom force fields and united atom force fields. In the former, each atom in the system is represented explicitly by potential functions. In the latter, hydrogens attached to heavy atoms (such as carbon) are removed. In their place single united (or extended) atom potentials are used. In this type of force field a CH2 group would appear as a single spherical atom. United atom sites have the advantage of greatly reducing the number of interaction sites in the molecule, but in certain cases can seriously limit the accuracy of the force field. United atom force fields are most usually required for the most computationally expensive tasks, such as the simulation of bulk liquid crystal phases via molecular dynamics or Monte Carlo methods (see Sect. 5.1). 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

Theories and simulation of the operation of AFM in liquid have been attempted [102-104], In principle, molecular dynamics or NEMD may be a suitable method to mimic the operation of a scanning tip. The time scale, however, precludes simulating a long-enough scan to see a complete atom. Most studies, therefore, were made with equilibrium conditions and a fixed position of the AFM tip. Explicit consideration of electrolytes and electrostatic effects has not been modeled. 

In tfiis chapter we address first the electrochemical application of the more familiar method of molecular (or atom) dynamics, and later turn to consider Monte Carlo methods, in each case giving a short introduction that should motivate the reader to pursue reading more specific works. Although the present research field is relatively new, the investigations are already too extensive to review in detail in a single chapter. For this reason, we discuss here the more extended research branches in the field and present a few representative examples. The application of simulations applied to nanostructuring problems is discussed in Chapter 36 liquid-liquid interfaces have been addressed by I. Benjamin (1997). 

The studies of molecular dynamics is now not only limited to the systems of atom and small molecules in dilute gas phase, the intense current interest is occurring in the condensed phase, i.e. in liquid solutions and in the solid state. 

To date, the only applications of these methods to the solution/metal interface have been reported by Price and Halley, who presented a simplified treatment of the water/metal interface. Briefly, their model involves the calculation of the metal s valence electrons wave function, assuming that the water molecules electronic density and the metal core electrons are fixed. The calculation is based on a one-electron effective potential, which is determined from the electronic density in the metal and the atomic distribution of the liquid. After solving the Schrddinger equation for the wave function and the electronic density for one configuration of the liquid atoms, the force on each atom is ciculated and the new positions are determined using standard molecular dynamics techniques. For more details about the specific implementation of these general ideas, the reader is referred to the original article. ... 

---