Implementation of reaction field methods in quantum chemistry computer codes

The embedding of a quantum mechanically described subsystem by classical representations of its surroundings is reviewed. The choices for a distributed monopole representation and a distributed (group) polarizability representation, as well as the continuum approach to model bulk effects, are discussed. Focus is on the practical implementation of the classical description in quantum chemistry codes (in particular, HONDO8.1). Expressions are given for the self‐consistent coupling between the classical partitions (dipole polarizabilities and boundary surface dipoles and charges) and for the coupling between classical and quantum partitions. The latter is mediated through expanded, rather than exact, potentials and fields. In this way, the computation of only a limited number of formal interactions between unit charge distributions located at the expansion centers suffices to evaluate the reaction field contributions. The electronic part of the coupling can be included in the Hamiltonian via the Fock matrix. The field operators, as well as the one‐ and two‐electron matrix elements over the basis functions, are simple. The expressions for these are given explicitly.

[1]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[2]  J. Rullmann,et al.  Analysis of discrete and continuum dielectric models : Application to the calculation of protonation energies in solution , 1987 .

[3]  Luaña,et al.  Simulation of ionic transition-metal crystals: The cluster model and the cluster-lattice interaction in the light of the theory of electronic separability. , 1989, Physical review. B, Condensed matter.

[4]  Accurate non-local electron—argon pseudopotential for condensed phase simulation , 1991 .

[5]  H. Scheraga,et al.  Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids , 1975 .

[6]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[7]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[8]  Jean-Louis Rivail,et al.  LIQUID-STATE QUANTUM CHEMISTRY : AN IMPROVED CAVITY MODEL , 1994 .

[9]  J. Rullmann,et al.  The active site of papain. All-atom study of interactions with protein matrix and solvent. , 1989, Journal of molecular biology.

[10]  P. V. Duijnen,et al.  THE DIRECT REACTION FIELD HAMILTONIAN - ANALYSIS OF THE DISPERSION TERM AND APPLICATION TO THE WATER DIMER , 1982 .

[11]  Patrick W. Fowler,et al.  Theoretical studies of van der Waals molecules and intermolecular forces , 1988 .

[12]  J. Rullmann,et al.  Intermolecular interactions with the direct reaction field method , 1990 .

[13]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[14]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[15]  D. Rinaldi,et al.  Influence of dispersion forces on the electronic structure of a solvated molecule , 1986 .

[16]  P. V. van Duijnen,et al.  Do parallel β‐strands have dipole moments? An ab initio molecular‐orbital‐direct reaction field study , 1985 .

[17]  H. Berendsen,et al.  The electric potential of a macromolecule in a solvent: A fundamental approach , 1991 .

[18]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[19]  K. Morokuma,et al.  The effective fragment potential method. An approximate ab initio mo method for large molecules , 1983 .

[20]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[21]  Sarah L. Price,et al.  A DISTRIBUTED MULTIPOLE ANALYSIS OF THE CHARGE-DENSITIES OF THE AZABENZENE MOLECULES , 1983 .

[22]  K. Mikkelsen,et al.  Electron-transfer reactions in solution. An ab initio approach , 1987 .

[23]  J. Ángyán,et al.  Are direct reaction field methods appropriate for describing dispersion interactions , 1990 .

[24]  M. Zerner,et al.  Reaction field effects on the electronic structure of carbon radical and ionic centers , 1990 .

[25]  L. Seijo,et al.  Ab initio model potential study of local distortions around Cr+ and Cr3+ defects in fluorite , 1991 .

[26]  S. Lifson,et al.  Energy functions for peptides and proteins. I. Derivation of a consistent force field including the hydrogen bond from amide crystals. , 1974, Journal of the American Chemical Society.

[27]  A theoretical model for SN1 ionic dissociation in solution. 1. Activation free energetics and transition-state structure , 1992 .

[28]  K. Szalewicz,et al.  Symmetry-adapted double-perturbation analysis of intramolecular correlation effects in weak intermolecular interactions , 1979 .

[29]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[30]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[31]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[32]  J. A. C. Rullmann,et al.  Atomic charge models for polypeptides derived from ab initio calculations , 1989 .

[33]  C. Cramer,et al.  General parameterized SCF model for free energies of solvation in aqueous solution , 1991 .

[34]  D. Remler,et al.  Molecular dynamics without effective potentials via the Car-Parrinello approach , 1990 .

[35]  F.J.Olivares del Valle,et al.  A computation procedure for the dispersion component of the interaction energy in continuum solute-solvent models , 1989 .

[36]  A. H. Juffer,et al.  DYNAMIC SURFACE BOUNDARY-CONDITIONS - A SIMPLE BOUNDARY MODEL FOR MOLECULAR-DYNAMICS SIMULATIONS , 1993 .

[37]  B. Thole.  Molecular polarizabilities calculated with a modified dipole interaction , 1981 .

[38]  A. H. Vries,et al.  Theoretical calculation of tautomer equilibria of 4-substituted imidazoles in the gas phase and in solution , 1992 .

[39]  P. V. van Duijnen,et al.  Reaction field effects on proton transfer in the active site of actinidin. , 1983, Biophysical Chemistry.

[40]  M. L. Connolly.  Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

[41]  Kenneth B. Wiberg,et al.  Solvent effects. 1. The mediation of electrostatic effects by solvents , 1991 .

[42]  K. Szalewicz,et al.  On the convergence properties of the Rayleigh–Schrödinger and the Hirschfelder–Silbey perturbation expansions for molecular interaction energies , 1977 .

[43]  J. Bentley.  Atomic Multipole Expansions of Molecular Charge Densities. Electrostatic Potentials , 1981 .

[44]  János G. Ángyán,et al.  Common theoretical framework for quantum chemical solvent effect theories , 1992 .

[45]  J. A. C. Rullmann,et al.  A polarizable water model for calculation of hydration energies , 1988 .

[46]  Trygve Helgaker,et al.  A multiconfigurational self‐consistent reaction‐field method , 1988 .

[47]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[48]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[49]  Susumu Narita,et al.  Extended Mulliken electron population analysis , 1990 .

[50]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[51]  Antonio Rizzo,et al.  Selected topics in ab initio computational chemistry in both very small and very large chemical systems , 1991 .

[52]  Minoru Sakurai,et al.  Medium effects on the molecular electronic structure. I. The formulation of a theory for the estimation of a molecular electronic structure surrounded by an anisotropic medium , 1987 .

[53]  Mati Karelson,et al.  Theoretical treatment of solvent effects on electronic spectroscopy , 1992 .

[54]  R. Zauhar,et al.  The rigorous computation of the molecular electric potential , 1988 .