Solvent effects Monte Carlo simulation

In principle, simulation teclmiques can be used, and Monte Carlo simulations of the primitive model of electrolyte solutions have appeared since the 1960s. Results for the osmotic coefficients are given for comparison in table A2.4.4 together with results from the MSA, PY and HNC approaches. The primitive model is clearly deficient for values of r. close to the closest distance of approach of the ions. Many years ago, Gurney [H] noted that when two ions are close enough together for their solvation sheaths to overlap, some solvent molecules become freed from ionic attraction and are effectively returned to the bulk [12]. 

As mentioned above, for the simulation in dimethylformamide (DMF) of the same reaction [53], the parameters for the substrate were not changed from the parametrization in water. For DMF the parameters were adopted from the OPLS parametrization of the pure liquid. The transferability was tested in part by performing a Monte Carlo simulation for CT plus 128 DMF molecules and evaluating the heat of solution for the chloride ion. The obtained value compares favorably with the experimental estimate. It is important to remark here that when potentials are used to simulate different solutions to the ones used in the parametrization process, they no longer are "effective" potentials. This fact becomes more evident in the simulation of solutions of small ions with localized charge that polarizes the neighboring solvent molecules. In this case it is convenient to consider the n-body corrections. 

Tapia, O. and Lluch, J. M. Solvent effects on chemical reaction profiles. Monte Carlo simulation of hydration effects on quantum chemically calculated stationary structures, J. Chem.Phys., 83 (1983, 3970-3982... 

Xue, Y Kim, C. K. Effects of Substituents and Solvents on the Reactions of Iminophosphorane with Formaldehyde Ab Initio MO Calculation and Monte Carlo Simulation. . Phys. Chem. A 2003, 107, 7945-7951. 

To evaluate solvent effects, statistical mechanical Monte Carlo simulations have been carried out. An important quantity to be computed is the potential of mean force, or free energy profile, as a function of the reaction coordinate, X, for a chemical reaction in solution using free energy perturbation method. (44) A straightforward approach is to determine free energy differences for incremental changes of certain geometrical variables that characteristically reflect the... 

As a second example, we have determined the influence of solvation on the steric retardation of SN2 reactions of chloride with ethyl and neopentyl chlorides in water, which has recently been studied by Vayner and coworkers [91]. In their study solvent effects were examined by means of QM-MM Monte Carlo simulations as well as with the CPCM model. Solvation causes a large increase in the activation energies of these reactions, but has a very small differential effect on the ethyl and neopentyl substrates. Nevertheless, a quantitative difference was found between the stability of the transition states determined using discrete and continuum treatments of solvation, since the activation free energies for ethyl chloride and neopentyl chloride amount to 23.9 and 30.4kcalmoF1 according to MC-FEP simulations, but to 38.4 and 47.6 kcal moF1 from CPCM computations. 

J. Zeng et al., Solvent effects on molecular spectra. I. Normal pressure and temperature Monte Carlo simulations of the structure of dilute pyrimidine in water. J. Chem. Phys. 99,1482-1495 (1993)... 

Coutinho, K., Oliveira M.J.D. and Canuto S., Sampling configurations in Monte Carlo simulations for quantum mechanical studies of solvent effects. Int. J. Quantum Chem. (1998) 66 249-253. 

To express the collective solvent reaction coordinate as in equation (6), it is necessary to define the specific diabatic potential surface for the reactant and product state. This, however, is not a simple task, and there is no unique way of defining such diabatic states. What is needed is a method that allows the preservation of the formal charges of the fragments of reactant and product resonance states. At the same time, solvent effects can be incorporated into electronic structure calculations in molecular dynamics and Monte Carlo simulations. Recently, we developed a block-localized wave function (BLW) method for studying resonance stabilization, hyperconjugation effects, and interaction energy decomposition of organic molecules.20-23 The BLW method can be formulated to specify the effective VB states.14... 

Dynamic Monte Carlo simulations were first used by Verdier and Stockmayer (5) for lattice polymers. An alternative dynamical Monte Carlo method has been developed by Ceperley, Kalos and Lebowitz (6) and applied to the study of single, three dimensional polymers. In addition to the dynamic Monte Carlo studies, molecular dynamics methods have been used. Ryckaert and Bellemans (7) and Weber (8) have studied liquid n-butane. Solvent effects have been probed by Bishop, Kalos and Frisch (9), Rapaport (10), and Rebertus, Berne and Chandler (11). Multichain systems have been simulated by Curro (12), De Vos and Bellemans (13), Wall et al (14), Okamoto (15), Kranbu ehl and Schardt (16), and Mandel (17). Curro s study was the only one without a lattice but no dynamic properties were calculated because the standard Metropolis method was employed. De Vos and Belleman, Okamoto, and Kranbuehl and Schardt studies included dynamics by using the technique of Verdier and Stockmayer. Wall et al and Mandel introduced a novel mechanism for speeding relaxation to equilibrium but no dynamical properties were studied. These investigations indicated that the chain contracted and the chain dynamic processes slowed down in the presence of other polymers. 

The explicit modeling approach surrounds a solute molecule with solvent molecules and then examines each molecule in that solvated environment. Quantum chemical methods, both semiempiricaP and ab initio" have been used to do this however, molecular dynamics and Monte Carlo simulations using force fields are used most often.Calculations on ensembles of molecules are more complex than those on individual molecules. Dykstra et al. discuss calculations on ensembles of molecules in a chapter in this book series. Because of the many conformations accessible to both solute and solvent molecules, in addition to the great number of possible solute molecule-solvent molecule orientations, such direct QM calculations are very computer intensive. However, the information resulting from this type of calculation is comprehensive because it provides molecular structures of the solute and solvent, and takes into account the effect of the solvent on the solute. This is the method of choice for assessing specific bonding information. 

---