Molecular dynamics microcanonical ensembles

In the above discussion it was assumed that the barriers are low for transitions between the different confonnations of the fluxional molecule, as depicted in figure A3.12.5 and therefore the transitions occur on a timescale much shorter than the RRKM lifetime. This is the rapid IVR assumption of RRKM theory discussed in section A3.12.2. Accordingly, an initial microcanonical ensemble over all the confonnations decays exponentially. However, for some fluxional molecules, transitions between the different confonnations may be slower than the RRKM rate, giving rise to bottlenecks in the unimolecular dissociation [4, ]. The ensuing lifetime distribution, equation (A3.12.7), will be non-exponential, as is the case for intrinsic non-RRKM dynamics, for an mitial microcanonical ensemble of molecular states. 

The first classical trajectory study of iinimoleciilar decomposition and intramolecular motion for realistic anhannonic molecular Hamiltonians was perfonned by Bunker [12,13], Both intrinsic RRKM and non-RRKM dynamics was observed in these studies. Since this pioneering work, there have been numerous additional studies [9,k7,30,M,M, ai d from which two distinct types of intramolecular motion, chaotic and quasiperiodic [14], have been identified. Both are depicted in figure A3,12,7. Chaotic vibrational motion is not regular as predicted by tire nonnal-mode model and, instead, there is energy transfer between the modes. If all the modes of the molecule participate in the chaotic motion and energy flow is sufficiently rapid, an initial microcanonical ensemble is maintained as the molecule dissociates and RRKM behaviour is observed [9], For non-random excitation initial apparent non-RRKM behaviour is observed, but at longer times a microcanonical ensemble of states is fonned and the probability of decomposition becomes that of RRKM theory. 

By far the most common methods of studying aqueous interfaces by simulations are the Metropolis Monte Carlo (MC) technique and the classical molecular dynamics (MD) techniques. They will not be described here in detail, because several excellent textbooks and proceedings volumes (e.g., [2-8]) on the subject are available. In brief, the stochastic MC technique generates microscopic configurations of the system in the canonical (NYT) ensemble the deterministic MD method solves Newton s equations of motion and generates a time-correlated sequence of configurations in the microcanonical (NVE) ensemble. Structural and thermodynamic properties are accessible by both methods the MD method provides additional information about the microscopic dynamics of the system. 

In the canonical ensemble (P2) = 3kBTM and p M. In the microcanonical ensemble (P2) = 3kgT i = 3kBTMNm/(M + Nm) [49]. If the limit M —> oo is first taken in the calculation of the force autocorrelation function, then p = Nm and the projected and unprojected force correlations are the same in the thermodynamic limit. Since MD simulations are carried out at finite N, the study of the N (and M) dependence of (u(t) and the estimate of the friction coefficient from either the decay of the momentum or force correlation functions is of interest. Molecular dynamics simulations of the momentum and force autocorrelation functions as a function of N have been carried out [49, 50]. 

Calvo, F. Neirotti, J.P. Freeman, D.L. Doll, J.D., Phase changes in 38-atom Lennard-Jones clusters. II. A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical ensembles, J. Chem. Phys. 2000, 112, 10350-10357... 

Temperature effects are included explicitly in molecular dynamics simulations by including kinetic energy terms - the balls representing the atoms are now on the move The principles are simple. In the microcanonical ensemble (NVE) ... 

Molecular Dynamics Methods. In contrast to the MC method, both kinetic and structural properties of a molecular system can be evaluated from MD studies. These properties are evaluated as averages over configurations generated during time. In microcanonical ensemble studies with the MD method, the properties which are controlled... 

Once the boundary conditions have been implemented, the calculation of solution molecular dynamics proceeds in essentially the same manner as do vacuum calculations. While the total energy and volume in a microcanonical ensemble calculation remain constant, the temperature and pressure need not remain fixed. A variant of the periodic boundary condition calculation method keeps the system pressure constant by adjusting the box length of the primary box at each step by the amount necessary to keep the pressure calculated from the system second virial at a fixed value (46). Such a procedure may be necessary in simulations of processes which involve large volume changes or fluctuations. Techniques are also available, by coupling the system to a Brownian heat bath, for performing simulations directly in the canonical, or constant T,N, and V, ensemble (2,46). 

Figure 2). The calculations were done in the microcanonical ensemble at a temperature of 300K 5K. Energy was well conserved throughout the trajectories, and no overall drifts in molecular temperature were observed. Small ensembles of trajectories (12 for SI and 6 each for the other minima) were calculated for the averaging of system properties. Each trajectory was equilibrated by velocity reassignments during an initial period of 20ps, followed by another 20ps of dynamics used for data collection.

At variance from Xe, the presented properties for Kr require more computional efforts. In order to reach the small-4 range of S(q), large-scale molecular dynamics have been carried out in the microcanonical ensemble (NVE) with the usual periodic boundary conditions. The equations of motion are integrated in the same discrete form as for Xe. The time step At is the same as for Xe and g r) is extracted over a sample of 8000 time-independent configurations every lOAf. 

The molecular dynamics simulation was performed using the MOTECC suite of programs [54] in the context of a microcanonical statistical ensemble. The system considered is a cube, with periodic boundary conditions, which contains 343 water molecules. The molecular dynamic simulation of water performed at ambient conditions revealed good agreement with experimental measurements. The main contribution to the total potential energy comes from the two-body term, while the many-body polarisation term contribution amounts to 23% of the total potential energy. Some of the properties calculated during the simulation are reported in Table 3. 

In standard molecular dynamic simulations the temperature is not constant. The MD simulation samples the microcanonical ensemble, or NVE ensemble, as the volume (unit-cell size) is assumed to be constant. The control of temperature is on the other hand especially important in the simulation of chemical reactions, when the excess of heat dissipated or adsorbed during the reaction strongly influences the kinetic energy (temperature) of the system. 

There are many excellent reviews on the standard molecular dynamics method dealing with calculations in the microcanonical ensemble as well as on the Monte Carlo method involving calculations in the canonical, isothermal isobaric, and grand canonical ensemble (< ). In the present article, we shall limit ourselves exclusively to those developments that have taken place since the work of Andersen (4). In the molecular dynamics method, the developments are the constant-pressure, constant-temperature, constant-temperature-constant-pressure, variable shape simulation cell MD, and isostress calculations in the Monte Carlo method, it is the variable shape simulation cell calculation. 

A different approach to the structural problem is provided by the technique of computerized molecular dynamics. It is now possible to set up in a computer program a microcanonical ensemble of perhaps two... 

Temperature is defined in molecular dynamics simulations in the microcanonical ensemble (constant N, V, T) by making use of the fact that the average kinetic energy K per degree of freedom is kT. Consequently, one initially assigns velocities to the molecules that add up to this value and then solves the equations of motion for a period of time sufficient to reach equilibrium while rescaling the velocities to hold the system at the desired temperature. After equilibrium has been achieved, the velocity scaling is turned off and the system is allowed to evolve in... 

Constant-energy molecular dynamics simulations (NVE microcanonical ensemble) were performed assuming that the subsystems are rigid (quaternion formalism) the respective code [23] uses a fifth-order predictor-corrector formalism. 

The use of non-Hamiltonian dynamical systems has a long history in mechanics [8] and they have recently been used to study a wide variety of problems in molecular dynamics (MD). In equilibrium molecular dynamics we can exploit non-Hamiltonian systems in order to generate statistical ensembles other than the standard microcanonical ensemble NVE) that is generated by traditional Hamiltonian dynamics. These ensembles, such as the canonical (NVT) and isothermal-isobaric (NPT) ensembles, are much better than the microcanonical ensemble for representing the actual conditions under which experiments are carried out. 

The mean transition temperatures obtained in molecular dynamics simulations of a canonical (Can) and a microcanonical (Mic) ensembles differ slightly... 

---