Monte Carlo simulations solvation forces

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. 

Within the framework of Monte Carlo simulations, the relation between measurable quantities and the microscopic structure of confined phases can now be examined. An example of such a measurable quantity is the solvation force F h)/2 KR (see Sec. IIA 1). From a theoretical perspective and according to the discussion in Sec. IIA 3 its investigation requires the stress T zisz) exerted normally by a confined fluid on planar substrates [see Eqs. (19) and (22)]. Using Eqs. (11) and (53) one can derive a molecular expression for Tzz from... 

The atomic radii may be further refined to improve the agreement between experimental and theoretical solvation free energies. Work on this direction has been done by Luque and Orozco (see [66] and references cited therein) while Barone et al. [67] defined a set of rules to estimate atomic radii. Further discussion on this point can be found in the review by Tomasi and co-workers [15], It must be noted that the parameterization of atomic radii on the basis of a good experiment-theory agreement of solvation energies is problematic because of the difficulty to separate electrostatic and non-electrostatic terms. The comparison of continuum calculations with statistical simulations provides another way to check the validity of cavity definition. A comparison between continuum and classical Monte Carlo simulations was reported by Costa-Cabral et al. [68] in the early 1980s and more recently, molecular dynamics simulations using combined quantum mechanics and molecular mechanics (QM/MM) force-fields have been carried out to analyze the case of water molecule in liquid water [69],... 

Jorgensen has parameterized by fitting properties of bulk liquids to Monte Carlo simulations to give the AMBER/OPLS force field (26,157, 158). Conceptually, one is attracted ly the use of liquids and their observable properties as constraints during the derivation of a force field that is destined to study the properties of solvated molecules. 

In Figure 5.27 a curve calculated from Equation 5.216 is compared with the predictions of other studies. The dotted line is calculated by means of the Henderson theory. The theoretical curve calculated by Kjellander and Sarmatf for ( ) = 0.357 and h>2 by using the anisotropic Percus-Yevick approximation is shown by the dashed line the crosses represent grand canonical Monte Carlo simulation results due to Karlstrom. We proceed now with separate descriptions of solvation, depletion, and colloid structural forces. 

Soon the excitement to see for the first time graphical representations of computed solvation shells for solvated ions from Monte Carlo simulations [81]. The next step was to go to even more complex systems, like enzymes, proteins, and particularly nucleic acids, A-DNA [82], B-DNA [83], without and with counterions [84] and in solution. The quantum biology community was taken by surprise, but soon accepted the new path as a new but necessary computational standard. I was proud to have forcefully recalled that the correct dictionary of quantum biology must contain terms like temperature, volume and free energy eventually, I was elected president of the International Society of Quantum Biology. 

Once the force field is chosen, a proper simulation method needs to be selected. Molecular dynamics simulations are applied to determine the solvation behaviour of ionic liquids by means of solving the Newtonian equations of motion for all molecules in the presence of a gradient in potential energy. Ionic liquid phase equilibria are determined by using Monte Carlo simulations in the isothermal isobaric Gibbs ensemble, grand canonical ensemble or osmotic ensemble with clever sampling schemes. 

---