Local reaction field method

Lee, F. S. and Warshel, A. (1992) A local-reaction-field method for fast evaluation of long-range electrostatic interactions in molecular simulations, J. Chem. Phys. 97, 3100-3107. 

Any of the methods used in classical Monte Carlo and molecular dynamics simulations may be borrowed in the combined QM/MM approach. However, the use of a finite system in condensed phase simulations is always a severe approximation, even when appropriate periodic or stochastic boundary conditions are employed. A further complication is the use of potential function truncation schemes, particular in ionic aqueous solutions where the long-range Coulombic interactions are significant beyond the cutoff distance.Thus, it is alluring to embed a continuum reaction field model in the quantum mechanical calculations in addition to the explicit solute—solvent interaaions to include the dielectric effect beyond the cutoff distance. - uch an onion shell arrangement has been used in spherical systems, whereas Lee and Warshel introduced an innovative local reaction field method for evaluation of long-... 

F. S. Lee and A. Warshel, J. Chem. Phys., 97, 3100 (1993). A Local Reaction Field Method for Fast Evaluation of Long-Range Electrostatic Interactions in Molecular Simulations. 

Moreau Y, Loss P-F, Assfeld X (2004) Solvent effects on the asymmetric Diels Alder reaction between cyclopentadiene and ( )-menthyl acrylate revisited with the three-layer hybrid local self-consistent field/molecular mechanics/self-consistent reaction field method. Theor Chem Ace 112 228-239... 

Basis Sets Correlation Consistent Sets Configuration Interaction Coupled-cluster Theory Density Functional Applications Density Functional Theory Applications to Transition Metal Problems G2 Theory Integrals of Electron Repulsion Integrals Overlap Linear Scaling Methods for Electronic Structure Calculations Localized MO SCF Methods Mpller-Plesset Perturbation Theory Monte Carlo Quantum Methods for Electronic Structure Numerical Hartree-Fock Methods for Molecules Pseudospectral Methods in Ab Initio Quantum Chemistry Self-consistent Reaction Field Methods Symmetry in Hartree-Fock Theory. 

However, another study concluded that the changes of the hydrogen-bond stability may be important in biological processes. For these, the influence of local electric fields created by Li+, Na+, and Mg2+ ions on the properties and reactivity of hydrogen bonds in HF and HC1 dimer has been carried out by means of ab initio self-consistent field (SCF) method [33]. A few years later, the effect of intensity and vector direction of the external electric field on activation barriers of unimole-cular reactions were studied using the semiempirical MINDO/3 method [34]. However, both semiempirical and ab initio calculations were performed to study the multiplicity change for carbene-like systems in external electric fields of different configurations (carbene and silylene) and the factor that determines the multiplicity and hence the reactivity of carbene-like structures is the nonuniformity of the field [35]. 

For molecular systems in the vacuum, exact analytical derivatives of the total energy with respect to the nuclear coordinates are available [22] and lead to very efficient local optimization methods [23], The situation is more involved for solvated systems modelled within the implicit solvent framework. The total energy indeed contains reaction field contributions of the form ER(p,p ), which are not calculated analytically, but are replaced by numerical approximations Efp(p,p ), as described in Section 1.2.5. We assume from now on that both the interface Y and the charge distributions p and p depend on n real parameters (A, , A ). In the geometry optimization problem, the A, are the cartesian coordinates of the nuclei. There are several nonequivalent ways to construct approximations of the derivatives of the reaction field energy with respect to the parameters (A1 , A ) ... 

It is also possible to combine the supermolecule and continuum approaches by using specific solvent molecules to capture the short-range effects (i.e., those involving specific noncovalent interactions between solute and solvent) and a reaction field to treat longer range effects.33-35 Alternatively, structures along the gas phase reaction coordinate can be immersed in a box of hundreds (or more) of explicit solvent molecules that are treated using force field approaches.36,37 Each type of method - the SCRF, solvent box, and supermolecule approaches - tests the importance of particular features of the solvent on the reactivity of the solute dielectric constant, multiple specific classical electrostatic interactions, and specific local directional noncovalent interactions, respectively. 

Qualitative trends for jS in simple model systems two-centre model system one-dimensional NLO-phores 143 Three-centre model system (C2 and D h symmetry) 145 Local field factors and effective polarizabilities the reaction field model Non-linear optical susceptibilities and experimental methods to evaluate and /3 153... 

Among the few determinations of of molecular crystals, the CPHF/ INDO smdy of Yamada et al. [25] is unique because, on the one hand, it concerns an open-shell molecule, the p-nitrophenyl-nitronyl-nitroxide radical (p-NPNN) and, on the other hand, it combines in a hybrid way the oriented gas model and the supermolecule approach. Another smdy is due to Luo et al. [26], who calculated the third-order nonlinear susceptibility of amorphous thinmultilayered films of fullerenes by combining the self-consistent reaction field (SCRF) theory with cavity field factors. The amorphous namre of the system justifies the choice of the SCRF method, the removal of the sums in Eq. (3), and the use of the average second hyperpolarizability. They emphasized the differences between the Lorentz Lorenz local field factors and the more general Onsager Bbttcher ones. For Ceo the results differ by 25% but are in similar... 

The effect of the solvent is usually modelled either by the use of the Onsager s self consistent reaction field (SCRF) [20] or by the polarizable continuum method (PCM) [21]. With regard to the relative stability of cytosine tautomers in aqueous solution, these methods provided results [14,15] which, in spite of some discrepancies, are in reasonable agreement with experimental data [3]. However, continuum-based methods do not explicitly take into consideration the local solvent-solute interaction which is instead important in the description of the proton transfer mechanism in hydrogen-bonded systems. A reasonable approach to the problem was recently proposed [22,23] in which the molecule of interest and few solvent molecules are treated as a supermolecule acting as solute, while the bulk of the solvent is represented as a polarizable dielectric. 

FIG U RE 3.17 The self-consistency problem in Pt electrocatalysis. The metal phase potential determines oxidation state and charging properties at the catalyst surface. These properties in turn determine the local reaction conditions at the Hehnholtz or reaction plane. At this point, structural design and transport properties of the catalyst layer come into play (as illustrated for conventional and ultrathin catalyst layers). Newly developed methods in the emerging field of first-principles electrochemistry attempt to find self-consistent solutions for this conpled problem. 

---