Molecular dynamics simulations framework

Predicting the solvent or density dependence of rate constants by equation (A3.6.29) or equation (A3.6.31) requires the same ingredients as the calculation of TST rate constants plus an estimate of and a suitable model for the friction coefficient y and its density dependence. While in the framework of molecular dynamics simulations it may be worthwhile to numerically calculate friction coefficients from the average of the relevant time correlation fiinctions, for practical purposes in the analysis of kinetic data it is much more convenient and instructive to use experimentally detemiined macroscopic solvent parameters. 

Combined Quantum and Molecular Mechanical Simulations. A recentiy developed technique is one wherein a molecular dynamics simulation includes the treatment of some part of the system with a quantum mechanical technique. This approach, QM/MM, is similar to the coupled quantum and molecular mechanical methods introduced by Warshel and Karplus (45) and at the heart of the MMI, MMP2, and MM3 programs by AUinger (60). These latter programs use quantum mechanical methods to treat the TT-systems of the stmctures in question separately from the sigma framework. 

For the Tar—Tar kissing loops, the P—B calculations are unable to discern their propensity to accumulate counterions accumulation at the loop—loop interface (data not shown). This is because the fully hydrated ions as defined by the Stem layer cannot penetrate into the central cation binding pocket (data not shown). Similarly, the axial spine of counterion density observed in the A-RNA helix (Fig. 20.5) is not captured by the P—B calculation (Fig. 20.7). No noticeable sequence specificity is observed in the counterion accumulation patterns in the P—B calculations, even though the sequence effects are explicitly represented in the P—B calculation through the appropriate geometry and assignment of point-charges. This is because the sequence specificity observed in the molecular dynamics simulations usually involves first shell interactions of base moieties with partially dehydrated ions, which cannot be accurately represented in the P—B framework. 

In order to overcome the limitations of currently available empirical force field param-eterizations, we performed Car-Parrinello (CP) Molecular Dynamic simulations [36]. In the framework of DFT, the Car-Parrinello method is well recognized as a powerful tool to investigate the dynamical behaviour of chemical systems. This method is based on an extended Lagrangian MD scheme, where the potential energy surface is evaluated at the DFT level and both the electronic and nuclear degrees of freedom are propagated as dynamical variables. Moreover, the implementation of such MD scheme with localized basis sets for expanding the electronic wavefunctions has provided the chance to perform effective and reliable simulations of liquid systems with more accurate hybrid density functionals and nonperiodic boundary conditions [37]. Here we present the results of the CPMD/QM/PCM approach for the three nitroxide derivatives sketched above details on computational parameters can be found in specific papers [13]. 

A problem with this interpretation relates to electrostriction, a process in which the density of the solvent changes about a solute. Shim et al. [243] noted evidence of electrostriction in molecular dynamics simulations of a model chromophore in an IL, and the degree of electrostriction was sensitive to the charge distribution of the solute. This observation does not necessarily contradict the framework above, as some local disruption of solvent structure due to dispersive interactions is inevitable. However, it is desirable to obtain a clearer understanding of the competition between these local interactions and the need to maintain a uniform charge distribution in the liquid. 

An innovative approach due to Haider et al. [113] may help to sidestep the challenges involved in explicit molecular dynamics simulation and obtain information on these slow dynamics. The authors use the results of dielectric reflectance spectroscopy to model the IL as a dielectric continuum, and study the solvation response of the IL in this framework. The calculated response is not a good description of the subpicosecond dynamics, a problem the authors ascribe to limited data on the high frequency dielectric response, but may be qualitatively correct at longer times. We have already expressed concern regarding the use of the dielectric continuum model for ILs in Section IV. A, but believe that if the wavelength dependence of the dielectric constant can be adequately modeled, this approach may be the most productive theoretical analysis of these slow dynamics. 

Successful first-principles molecular dynamics simulations in the Car-Paxrinello framework requires low temperature for the annealed electronic parameters while maintaining approximate energy conservation of the nuclear motion, all without resorting to unduly small time steps. The most desirable situation is a finite gap between the frequency spectrum of the nuclear coordinates, as measured, say, by the velocity-velocity autocorrelation function. 

Section 2 of this chapter notes will be devoted to the framework for separation of the ionic and electronic dynamics through the Born-Oppenheimer approximation. Atomic motion, with forces on the ions at each timestep evaluated through an electronic structure calculation, can then be propagated by Molecular Dynamics simulations, as proposed by first-principle Molecular Dynamics. This allows for a description of the electronic reorganisation following the atomic motion, e.g. bond rearrangements in chemical reactions. 

On the other hand, it is not the goal of this book to present an exhaustive review of the measurement of physico-chemical properties carried out in the last two or three decades. This would certainly exceed the framework of this book. The electrochemistry of molten salts is not included in this book either, since in my opinion, this topic would need a separate book due to its complexity and wealth of applications. The quantum chemistry and molecular dynamics simulation methods are omitted as well. 

Although it turned out that a number of essential features concerning dynamics of molecular liquids can be well captured by the theory of the previous subsection (see Secs. 3 and 5), an intense investigation through experimental, theoretical, and molecular-dynamics simulation studies for simple liquids has revealed that the microscopic processes underlying various time-dependent phenomena cannot be fully accounted for by a simplified memory-function approach [18, 19, 20]. In particular, the assumption that the decay of memory kernels is ruled by a simple exponential-type relaxation must be significantly revised in view of the results of the kinetic framework developed for dense liquids (see Sec. 5.1.4). This motivated us to further improve the theory for dynamics of polyatomic fluids. 

A distinctive feature of our theory is that it offers a self-contained framework for solving dynamical problems without making reference to any dynamical information from outside, i.e., it requires only the knowledge of parameters for potential functions and molecular geometry (such as bond lengths) as in simulation studies. This is in contrast to previous theories, the SSSV and RMFA theories, in which diffusion coefficients or memory functions of some reference dynamical variables have to be determined in advance from molecular dynamics simulations. 

Membrane modelling has been considered from both the nano/microscopic and the macroscopic viewpoints, but little has been done to bridge these two limits. The breadth of microscopic modelling work for PEMs encompasses molecular dynamics simulations [17] and statistical mechanics modelling [18-21]. Most applications have focused on Nafion, and interestingly, some models even apply macroscopic transport relations to the microscopic transport within a pore of a membrane [41]. While the focus of this Chapter is on macroscopic models required for computational simulations of complete fuel cells [12,13,15], the proposed modelling framework is based on fundamental relations describing molecular transport phenomena. 

---