Liquid state molecular dynamics

By employing a very strong external field, a gedankexperiment may be set up whereby the natural thermal motion of the molecules is put in competition with the aligning effect of the field. This method reveals some properties of the molecular liquid state which are otherwise hidden. In order to explain the observable effects of the applied fields, it is necessary to use equations of motion more generally valid than those of Benoit. These equations may be incorporated within the general structure of reduced model theory " (RMT) and illustrate the use of RMT in the context of liquid-state molecular dynamics. (Elsewhere in this volume RMT is applied to problems in other fields of physics where consideration of stochastic processes is necessary.) In this chapter modifications to the standard methods are described which enable the detailed study of field-on molecular dynamics. 

Both vibrational and rotovibrational relaxation can be described analyti-caDy as multiplicative stochastic processes. For these processes, RMT is equivalent to the stochastic Liouville equation of Kubo, with the added feature that RMT takes into account the back-reaction from the molecule imder consideration on the thermal bath. The stochastic Liouville equation has been used successfully to describe decoupling in the transient field-on condition and the effect of preparation on decay. When dealing with liquid-state molecular dynamics, RMT provides a rigorous justification for itinerant oscillator theory, widely applied to experimental data by Evans and coworkers. This implies analytically that decoupling effects should be exhibited in molecular liquids treated with strong fields. In the absence of experimental data, the computer runs described earlier amount to an independent means of verifying Grigolini s predictions. In this context note that the simulation of Oxtoby and coworkers are semistochastic and serve a similar purpose. 

The coherent motion initiated by an excitation pulse can be monitored by variably delayed, ultrashort probe pulses. Since these pulses may also be shorter in duration than the vibrational period, individual cycles of vibrational oscillation can be time resolved and spectroscopy of vibrationally distorted species (and other unstable species) can be carried out. In the first part of this section, the mechanisms through which femtosecond pulses may initiate and probe coherent lattice and molecular vibrational motion are discussed and illustrated with selected experimental results. Next, experiments in the areas of liquid state molecular dynamics and chemical reaction dynamics are reviewed. These important areas can be addressed incisively by coherent spectroscopy on the time scale of individual molecular collisions or half-collisions. 

To examine the solid as it approaches equilibrium (atom energies of 0.025 eV) requires molecular dynamic simulations. Molecular dynamic (MD) simulations follow the spatial and temporal evolution of atoms in a cascade as the atoms regain thermal equilibrium in about 10 ps. By use of MD, one can follow the physical and chemical effects that influence the final cascade state. Molecular dynamics have been used to study a variety of cascade phenomena. These include defect evolution, recombination dynamics, liquid-like core effects, and final defect states. MD programs have also been used to model sputtering processes. 

I. Benjamin,/. Chem. Phys., 103, 2459 (1995). Photodissociation of ICN in Liquid Chloroform Molecular Dynamics of Ground and Excited State Recombination, Cage Escape and Hydrogen Abstraction Reaction. 

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... 

Maxwell J C 1874 Van der Waals on the continuity of the gaseous and liquid states Nature 10 477 Maxwell J C 1875 On the dynamical evidence of the molecular constitution of bodies Nature 11 357... 

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. 

In Chapter 2, a brief discussion of statistical mechanics was presented. Statistical mechanics provides, in theory, a means for determining physical properties that are associated with not one molecule at one geometry, but rather, a macroscopic sample of the bulk liquid, solid, and so on. This is the net result of the properties of many molecules in many conformations, energy states, and the like. In practice, the difficult part of this process is not the statistical mechanics, but obtaining all the information about possible energy levels, conformations, and so on. Molecular dynamics (MD) and Monte Carlo (MC) simulations are two methods for obtaining this information... 

Molecular mechanics methods have been used particularly for simulating surface-liquid interactions. Molecular mechanics calculations are called effective potential function calculations in the solid-state literature. Monte Carlo methods are useful for determining what orientation the solvent will take near a surface. Molecular dynamics can be used to model surface reactions and adsorption if the force held is parameterized correctly. 

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... 

Models for description of liquids should provide us with an understanding of the dynamic behavior of the molecules, and thus of the routes of chemical reactions in the liquids. While it is often relatively easy to describe the molecular structure and dynamics of the gaseous or the solid state, this is not true for the liquid state. Molecules in liquids can perform vibrations, rotations, and translations. A successful model often used for the description of molecular rotational processes in liquids is the rotational diffusion model, in which it is assumed that the molecules rotate by small angular steps about the molecular rotation axes. One quantity to describe the rotational speed of molecules is the reorientational correlation time T, which is a measure for the average time elapsed when a molecule has rotated through an angle of the order of 1 radian, or approximately 60°. It is indirectly proportional to the velocity of rotational motion. 

In this article, methods for extracting ground state structural, dynamical and molecular electronic properties for liquid crystals have been outlined. It is clear that these methods have been applied only to a small number of systems... 

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]. 

---