Geometry optimization of molecular structures in solution by the polarizable continuum model

A new implementation of analytical gradients for the polarizable continuum model is presented, which allows Hartree‐Fock and density functional calculations taking into account both electrostatic and nonelectrostatic contributions to energies and gradients for closed and open shell systems. Simplified procedures neglecting the derivatives of the cavity surface and/or using single spheres for XHn groups have also been implemented and tested. The solvent‐induced geometry relaxation has been studied for a number of representative systems in order to test the efficiency of the procedure and to investigate the role of different contributions. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 404–417, 1998

[1]  M. Räsänen.  A matrix infrared study of monomeric formamide , 1983 .

[2]  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 .

[3]  J. Tomasi,et al.  Solute—solvent electrostatic interactions with non-homogeneous radial dielectric functions , 1994 .

[4]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[5]  H. Schlegel,et al.  Optimization of equilibrium geometries and transition structures , 1982 .

[6]  J. Tomasi,et al.  Ab initio study of solvated molecules: A new implementation of the polarizable continuum model , 1996 .

[7]  A new general form of molecular force fields. Application to intra- and interresidue interactions in peptides , 1992 .

[8]  Iñaki Tuñón,et al.  GEPOL: An improved description of molecular surfaces II. Computing the molecular area and volume , 1991 .

[9]  Iñaki Tuñón,et al.  GEPOL: An improved description of molecular surfaces. III. A new algorithm for the computation of a solvent‐excluding surface , 1994, J. Comput. Chem..

[10]  V. Barone,et al.  Development and validation of reliable quantum mechanical approaches for the study of free radicals in solution , 1996 .

[11]  T. Yamazaki,et al.  Conformational studies of N‐acetyl‐N′‐methylamide derivatives of α‐aminobutyric acid, norvaline, and valine. I. Preferred conformations in solution as studied by 1H‐nmr spectroscopy , 1988 .

[12]  Vincenzo Barone,et al.  Inclusion of Hartree-Fock exchange in the density functional approach. Benchmark computations for diatomic molecules containing H, B, C, N, O, and F atoms , 1994 .

[13]  B. Honig,et al.  Accurate First Principles Calculation of Molecular Charge Distributions and Solvation Energies from Ab Initio Quantum Mechanics and Continuum Dielectric Theory , 1994 .

[14]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[15]  Alan R. Katritzky,et al.  Quantitative predictions of tautomeric equilibria for 2-, 3-, and 4-substituted pyridines in both the gas phase and aqueous solution: combination of AM1 with reaction field theory , 1989 .

[16]  H. Rzepa,et al.  An MCSCF study of the effect of substituents and solvent on the [2 + 2]cycloaddition of tert-butylcyanoketene to phenylethene , 1993 .

[17]  J. McDouall,et al.  Assessment of the Langevin dipoles solvation model for Hartree-Fock wavefunctions , 1996 .

[18]  Ian R. Gould,et al.  A quantum Mechanical Investigation of the Conformational Energetics of the Alanine and Glycine Dipeptides in the Gas Phase and in Aqueous Solution , 1994 .

[19]  Guntram Rauhut,et al.  Transferable Scaling Factors for Density Functional Derived Vibrational Force Fields , 1995 .

[20]  Manuel F. Ruiz-López,et al.  Computation of hydration free energies using a parameterized continuum model: Study of equilibrium geometries and reactive processes in water solution , 1996, J. Comput. Chem..

[21]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[22]  D. Case,et al.  Incorporating solvation effects into density functional theory: Calculation of absolute acidities , 1997 .

[23]  Kenneth B. Wiberg,et al.  Barriers to rotation adjacent to double bonds. 3. The carbon-oxygen barrier in formic acid, methyl formate, acetic acid, and methyl acetate. The origin of ester and amide resonance , 1987 .

[24]  Vincenzo Barone,et al.  Structure and ESR features of glycine radical in its zwitterionic form , 1995 .

[25]  J. Tomasi,et al.  Analytical expressions of the free energy derivatives for molecules in solution. Application to the geometry optimization , 1995 .

[26]  J. Tomasi,et al.  Analytical derivatives for molecular solutes. III. Hartree–Fock static polarizability and hyperpolarizabilities in the polarizable continuum model , 1996 .

[27]  Jean-Louis Rivail,et al.  Analytical energy derivatives for a realistic continuum model of solvation: Application to the analysis of solvent effects on reaction paths , 1996 .

[28]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[29]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

[30]  Benny G. Johnson,et al.  Kohn—Sham density-functional theory within a finite basis set , 1992 .

[31]  J. Tomasi,et al.  Cavitation and Electrostatic Free Energy for Molecular Solutes in Liquid Crystals , 1996 .

[32]  J. Tomasi,et al.  The implementation of density functional theory within the polarizable continuum model for solvation , 1994 .

[33]  Krishnan Raghavachari,et al.  Gaussian-2 theory using reduced Moller--Plesset orders , 1993 .

[34]  Andreas Klamt,et al.  Treatment of the outlying charge in continuum solvation models , 1996 .

[35]  J. Tomasi,et al.  Dispersion and repulsion contributions to the solvation energy: Refinements to a simple computational model in the continuum approximation , 1991 .

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

[37]  V. Barone,et al.  SOLVENT EFFECTS ON THE CONFORMATIONAL BEHAVIOR OF MODEL PEPTIDES. A COMPARISON BETWEEN DIFFERENT CONTINUUM MODELS , 1996 .

[38]  V. Barone.  CHARACTERIZATION OF THE POTENTIAL ENERGY SURFACE OF THE HO2 MOLECULAR SYSTEM BY A DENSITY FUNCTIONAL APPROACH , 1994 .

[39]  Jacopo Tomasi,et al.  A new definition of cavities for the computation of solvation free energies by the polarizable continuum model , 1997 .

[40]  J. Tomasi,et al.  A theoretical model of solvation in continuum anisotropic dielectrics , 1995 .

[41]  T. Truong,et al.  Analytical first and second energy derivatives of the generalized conductorlike screening model for free energy of solvation , 1995 .

[42]  Benny G. Johnson,et al.  An implementation of analytic second derivatives of the gradient‐corrected density functional energy , 1994 .

[43]  Estanislao Silla,et al.  GEPOL: An improved description of molecular surfaces. I. Building the spherical surface set , 1990 .

[44]  Vincenzo Barone,et al.  Density Functional Study of Intrinsic and Environmental Effects in the Tautomeric Equilibrium of 2-Pyridone , 1995 .

[45]  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.

[46]  V. Barone.  Structure, epr parameters, and reactivity of organic free radicals from a density functional approach , 1995 .

[47]  J. Tomasi,et al.  Energy and energy derivatives for molecular solutes: Perspectives of application to hybrid quantum and molecular methods , 1996 .

[48]  J. Tomasi,et al.  Analytical Hartree–Fock calculation of the dynamical polarizabilities α, β, and γ of molecules in solution , 1996 .

[49]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[50]  F. Javier Luque,et al.  Optimization of solute cavities and van der Waals parameters in ab initio MST‐SCRF calculations of neutral molecules , 1994, J. Comput. Chem..

[51]  Arieh Warshel,et al.  Langevin Dipoles Model for ab Initio Calculations of Chemical Processes in Solution: Parametrization and Application to Hydration Free Energies of Neutral and Ionic Solutes and Conformational Analysis in Aqueous Solution , 1997 .

[52]  Jacopo Tomasi,et al.  Analytical derivatives for molecular solutes. I. Hartree–Fock energy first derivatives with respect to external parameters in the polarizable continuum model , 1994 .

[53]  Martin J. Field,et al.  Free energy perturbation method for chemical reactions in the condensed phase: a dynamic approach based on a combined quantum and molecular mechanics potential , 1987 .

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

[55]  D. Rinaldi,et al.  Studies of solvent effects using density functional theory. Co-operative interactions in H3N…HBr proton transfer , 1994 .

[56]  János G. Ángyán,et al.  CHOOSING BETWEEN ALTERNATIVE MP2 ALGORITHMS IN THE SELF-CONSISTENT REACTION FIELD THEORY OF SOLVENT EFFECTS , 1995 .

[57]  P. Kollman,et al.  A Simulation of the Catalytic Mechanism of Aspartylglucosaminidase Using ab Initio Quantum Mechanics and Molecular Dynamics , 1997 .

[58]  V. Barone.  Electronic, vibrational and environmental effects on the hyperfine coupling constants of nitroside radicals. H2NO as a case study , 1996 .

[59]  Roberto Cammi,et al.  Analytical first derivatives of molecular surfaces with respect to nuclear coordinates , 1996, J. Comput. Chem..

[60]  Richard J. Hall,et al.  Combined Density Functional, Self-Consistent Reaction Field Model of Solvation , 1995 .

[61]  Jacopo Tomasi,et al.  Electrostatic interaction of a solute with a continuum. Improved description of the cavity and of the surface cavity bound charge distribution. , 1987 .

[62]  David A. Case,et al.  Incorporating solvation effects into density functional electronic structure calculations , 1994 .

[63]  Richard A. Friesner,et al.  Quantum mechanical geometry optimization in solution using a finite element continuum electrostatics method , 1996 .

[64]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[65]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[66]  Jacopo Tomasi,et al.  Electron correlation and solvation effects. I, Basic formulation and preliminary attempt to include the electron correlation in the quantum mechanical polarizable continuum model so as to study solvation phenomena , 1991 .

[67]  H. Ågren,et al.  Multiconfigurational self-consistent reaction field theory for nonequilibrium solvation , 1995 .

[68]  Donald G. Truhlar,et al.  MODEL FOR AQUEOUS SOLVATION BASED ON CLASS IV ATOMIC CHARGES AND FIRST SOLVATION SHELL EFFECTS , 1996 .

[69]  A. Fortunelli.  A study of the MBPT correlation corrections within the polarizable continuum model for solvation , 1995 .

[70]  Kenneth B. Wiberg,et al.  Solvent effects. 3. Tautomeric equilibria of formamide and 2-pyridone in the gas phase and solution: an ab initio SCRF study , 1992 .

[71]  Jacopo Tomasi,et al.  Evaluation of the dispersion contribution to the solvation energy. A simple computational model in the continuum approximation , 1989 .

[72]  R. Pierotti,et al.  A scaled particle theory of aqueous and nonaqueous solutions , 1976 .

[73]  K. Kopple,et al.  Solvent-dependent conformational distributions of some dipeptides , 1980 .

[74]  T. Russo,et al.  Solvation Free Energy Calculations Using a Continuum Dielectric Model for the Solvent and Gradient-Corrected Density Functional Theory for the Solute , 1996 .

[75]  M. Robb,et al.  Calculation of solvatochromic shifts using MC-SCF theory. The n—π* transition of acetone , 1993 .

[76]  K. Wiberg,et al.  Solvent Effects. 5. Influence of Cavity Shape, Truncation of Electrostatics, and Electron Correlation on ab Initio Reaction Field Calculations , 1996 .

[77]  V. Barone,et al.  STRUCTURE AND ESR FEATURES OF GLYCINE RADICAL , 1995 .

[78]  Axel D. Becke,et al.  Density‐functional thermochemistry. IV. A new dynamical correlation functional and implications for exact‐exchange mixing , 1996 .

[79]  Andreas Klamt,et al.  Incorporation of solvent effects into density functional calculations of molecular energies and geometries , 1995 .