Dynamic continuum models

F. Ingrosso, B. Mennucci and J. Tomasi, Quantum mechanical calculations coupled with a dynamical continuum model for the description of dielectric relaxation time dependent Stokes shift of coumarin Cl53 in polar solvents, J. Mol. Liq., 108 (2003) 21 -6. 

Chapter 6 contains the most advanced mathematics in the book, and introduces a particular nonlinear static and dynamic continuum model for smectic C liquid crystals. This continuum theory is a natural extension of nematic theory, originating from ideas that are familiar from the continuum description of nematics. 

Van der Zwan G and Hynes J T 1982 Dynamical polar solvent effects on solution reactions A simple continuum model J. Chem. Phys. 76 2993-3001... 

The present discussion of continuum modeling of dynamic fracture is not an exhaustive review. Rather, it points out the variety of approaches which have been, and are still being, pursued to provide methods for calculating dynamic fracture phenomena. Such work is still quite active and considerable effort... 

The flat interface model employed by Marcus does not seem to be in agreement with the rough picture obtained from molecular dynamics simulations [19,21,64-66]. Benjamin examined the main assumptions of work terms [Eq. (19)] and the reorganization energy [Eq. (18)] by MD simulations of the water-DCE junction [8,19]. It was found that the electric field induced by both liquids underestimates the effect of water molecules and overestimates the effect of DCE molecules in the case of the continuum approach. However, the total field as a function of the charge of the reactants is consistent in both analyses. In conclusion, the continuum model remains as a good approximation despite the crude description of the liquid-liquid boundary. 

The theoretical aspects highlighted in this section show that describing heterogeneous ET dynamics by the continuum model can be regarded as reasonable. In the next section, we shall discuss in detail the experimental advances and how some results cast doubts in the picture outlined so far. 

Archontis, G. Simonson, T. Dielectric relaxation in an enzyme active site molecular dynamics simulations intepreted with a macroscopic continuum model, J. Am. Chem. Soc. 2001,123, 11047-11056. 

Abstract This chapter reviews the theoretical background for continuum models of solvation, recent advances in their implementation, and illustrative examples of their use. Continuum models are the most efficient way to include condensed-phase effects into quantum mechanical calculations, and this is typically accomplished by the using self-consistent reaction field (SCRF) approach for the electrostatic component. This approach does not automatically include the non-electrostatic component of solvation, and we review various approaches for including that aspect. The performance of various models is compared for a number of applications, with emphasis on heterocyclic tautomeric equilibria because they have been the subject of the widest variety of studies. For nonequilibrium applications, e.g., dynamics and spectroscopy, one must consider the various time scales of the solvation process and the dynamical process under consideration, and the final section of the review discusses these issues. 

Continuum models remove the difficulties associated with the statistical sampling of phase space, but they do so at the cost of losing molecular-level detail. In most continuum models, dynamical properties associated with the solvent and with solute-solvent interactions are replaced by equilibrium averages. Furthermore, the choice of where the primary subsystem ends and the dielectric continuum begins , i.e., the boundary and the shape of the cavity containing the primary subsystem, is ambiguous (since such a boundary is intrinsically nonphysical). Typically this boundary is placed on some sort of van der Waals envelope of either the solute or the solute plus a few key solvent molecules. 

The two avenues above recalled, namely ab-initio computations on clusters and Molecular Dynamics on one hand and continuum model on the other, are somewhat bridged by those techniques where the solvent is included in the hamiltonian at the electrostatic level with a discrete representation [13,17], It is important to stress that quantum-mechanical computations imply a temperature of zero K, whereas Molecular Dynamics computations do include temperature. As it is well known, this inclusion is of paramount importance and allows also the consideration of entropic effects and thus free-energy, essential parameters in any reaction. 

The dynamics of carbon-halogen bond reductive cleavage in alkyl halides was studied by MP3 ab initio calculations, using pseudopotentials for the halogens and semidiffuse functions for the heavy atoms [104], The effect of solvent was treated by means of the ellipsoidal cavity dielectric continuum model. Both a concerted (i.e., a one-step) and a stepwise mechanism (in which an anion radical is formed at first) were... 

In chapter 1, Profs. Cramer and Truhlar provide an overview of the current status of continuum models of solvation. They examine available continuum models and computational techniques implementing such models for both electrostatic and non-electrostatic components of the free energy of solvation. They then consider a number of case studies with particular focus on the prediction of heterocyclic tautomeric equilibria. In the discussion of the latter they focus attention on the subtleties of actual chemical systems and some of the danger in applying continuum models uncritically. They hope the reader will emerge with a balanced appreciation of the power and limitations of these methods. In the last section they offer a brief overview of methods to extend continuum solvation modeling to account for dynamic effects in spectroscopy and kinetics. Their conclusion is that there has been tremendous progress in the development and practical implementation of useful continuum models in the last five years. These techniques are now poised to allow quantum chemistry to have the same revolutionary impact on condensed-phase chemistry as the last 25 years have witnessed for gas-phase chemistry. 

Coarse-grained polymer models neglect the chemical detail of a specific polymer chain and include only excluded volume and topology (chain connectivity) as the properties determining universal behavior of polymers. They can be formulated for the continuum (off-lattice) as well as for a lattice. For all coarse-grained models, the repeat unit or monomer unit represents a section of a chemically realistic chain. MD techniques are employed to study dynamics with off-lattice models, whereas MC techniques are used for the lattice models and for efficient equilibration of the continuum models.36 2 A tutorial on coarse-grained modeling can be found in this book series.43... 

Quantitative models of solute-solvent systems are often divided into two broad classes, depending upon whether the solvent is treated as being composed of discrete molecules or as a continuum. Molecular dynamics and Monte Carlo simulations are examples of the former 8"11 the interaction of a solute molecule with each of hundreds or sometimes even thousands of solvent molecules is explicitly taken into account, over a lengthy series of steps. This clearly puts a considerable demand upon computer resources. The different continuum models,11"16 which have evolved from the work of Bom,17 Bell,18 Kirkwood,19 and Onsager20 in the pre-computer era, view the solvent as a continuous, polarizable isotropic medium in which the solute molecule is contained within a cavity. The division into discrete and continuum models is of course not a rigorous one there are many variants that combine elements of both. For example, the solute molecule might be surrounded by a first solvation shell with the constituents of which it interacts explicitly, while beyond this is the continuum solvent.16... 

An interesting combined use of discrete molecular and continuum techniques was demonstrated by Floris et al.181,182 They used the PCM to develop effective pair potentials and then applied these to molecular dynamics simulations of metal ion hydration. Another approach to such systems is to do an ab initio cluster calculation for the first hydration shell, which would typically involve four to eight water molecules, and then to depict the remainder of the solvent as a continuum. This was done by Sanchez Marcos et al. for a group of five cations 183 the continuum model was that developed by Rivail, Rinaldi et al.14,108-112 (Section III.2.ii). Their results are compared in Table 14 with those of Floris et al.,139 who used a similar procedure but PCM-based. In... 

This Chapter has outlined several different approaches to the computational determination of solution properties. Two of these address solute-solvent interactions directly, either treating the effects of individual solvent molecules upon the solute explicitly or by means of a reaction field due to a continuum model of the solvent. The other procedures establish correlations between properties of interest and certain features of the solute and/or solvent molecules. There are empirical elements in all of these methods, even the seemingly more rigorous ones, such as the parameters in the molecular dynamics/Monte Carlo intermolecular potentials, Eqs. (16) and (17), or in the continuum model s Gcavitation and Gvdw, Eqs. (40) and (41), etc. 

---