6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

A combined quantum mechanical (QM) and molecular mechanical (MM) potential has been developed for the study of reactions in condensed phases. For the quantum mechanical calculations semiempirical methods of the MNDO and AM1 type are used, while the molecular mechanics part is treated with the CHARMM force field. Specific prescriptions are given for the interactions between the QM and MM portions of the system; cases in which the QM and MM methodology is applied to parts of the same molecule or to different molecules are considered. The details of the method and a range of test calculations, including comparisons with ab initio and experimental results, are given. It is found that in many cases satisfactory results are obtained. However, there are limitations to this type of approach, some of which arise from the AM1 or MNDO methods themselves and others from the present QM/MM implementation. This suggests that it is important to test the applicability of the method to each particular case prior to its use. Possible areas of improvement in the methodology are discussed.

A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations

Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics.

The integral equation formalism (IEF) is a recent method (the grounds have been elaborated at the beginning of 1997) addressed to solve the electrostatic solvation problem at the QM level with the aid of apparent surface charges (ASC). IEF uses a new formalism of this problem, based on integral operators never used before in the chemical community and it manages to treat on the same footing linear isotropic solvent models, as well as anisotropic liquid crystals and ionic solutions. In this overview we emphasize the good performances of IEF at the lowest level of its potentialities, i.e. for isotropic solvents, as a new approach to compute solvation free energies and properties (dipole hyperpolarizabilities) of molecular solutes, as well as energy gradients for geometry optimization procedures. Finally we present a new IEF implementation of the nonequilibrium problem for electronic spectra which appears to be decidedly competitive with the previous more standard ASC formulations.

The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level

Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models toward more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date, it has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with...

/pdf/current-status-of-the-amoeba-polarizable-force-field-5adk1teb7t.pdf

Current Status of the AMOEBA Polarizable Force Field

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. Nonequilibrium potentials and Monte Carlo sampling over classical degrees of freedom have been added to better mimic experimental conditions. (C) 1995 by John Wiley and Sons, Inc.

Implementation of reaction field methods in quantum-chemistry computer codes

The induction and dispersion terms obtained from quantum-mechanical calculations with a direct reaction field hamiltonian are compared to second order perturbation theory expressions. The dispersion term is shown to give an upper bound which is a generalization of Alexander's upper bound. The model is illustrated by a calculation on the interactions in the water dimer. The long range Coulomb, induction and dispersion interactions are reasonably reproduced.

The direct reaction field hamiltonian: Analysis of the dispersion term and application to the water dimer

In this paper a combined experimental and quantum chemical study of the geometry and opto-electronic properties of unsubstituted and dialkoxy-sustituted phenylene–vinylene oligomers (PV’s) is presented. The optical absorption spectra for PV cations with different chain lengths and substitution patterns were measured using pulse radiolysis with time-resolved spectrophotometric detection from 1380 to 500 nm (0.9 to 2.5 eV). The geometries of the PV’s studied were optimized using density functional theory (DFT) for both the neutral and singly charged molecule. The spectra for the PV radical cations were then calculated using singly excited configuration interaction with an intermediate neglect of differential overlap reference wave function method together with the DFT geometry. The agreement between experimental and theoretical absorption energies is excellent; most of the calculated radical cation absorption energies are within 0.15 eV of the experimental values. The pattern of dialkoxy-substitution is found to have a large effect on the optical absorption spectrum of the cation. Using the calculated charge distribution it is shown that the degree of delocalization of the charge correlates with the energy of the lowest absorption band. If alkoxy side chains are present on some of the rings the positive charge tends to localize at those sites.

/pdf/theoretical-and-experimental-studies-of-the-opto-electronic-36uks6vmzi.pdf

Theoretical and experimental studies of the opto-electronic properties of positively charged oligo (phenylene vinylene)s: Effects of chain length and alkoxy substitution

The direct reaction field (DRF) approach has proven to be a useful tool to investigate the influence of solvents on the quantum/classical behaviour of solute molecules. In this paper, we report the latest extension of this DRF approach, which consists of the gradient of the completely classical energy expressions of this otherwise QM/MM method. They can be used in (completely classical) molecular dynamics (MD) simulations and geometry optimizations, that can be followed by a number of single point QM/MM calculations on configurations obtained in these simulations/optimizations. We report all energy and gradient expressions, and results for a number of interesting (model) systems. They include geometry optimization of the benzene dimer as well as MD simulations of some solvents. The most stable configuration for the benzene dimer is shown to be the parallel-displaced form, which is slightly more stable (0.3 kcal/mol) than the T-shaped dimer.

DRF90: a polarizable force field

Ab initio MO calculations, using both minimal (STO-3G) and extended (Roos-Siegbahn) basis sets are reported for the systems methanethiol-imidazole, methanethiol-imidazole-formaldehyde, and methanethiol-imidazole-formamide, which, together with a point-change representation of a long α-helix, form models for the active site of papain. It is shown that the large electric field exerted by the helix in the active-site region is responsible for the presence of the essential residues Cys 25 and His 159 in the form of an ion pair RS- · · · ImH+, which is crucial for a recently proposed mechanism for the catalytic action of the enzyme. Also, an explanation is given for the anomalies in measured pK values for these residues. Detailed studies on the (sub)systems show that minimal basis sets lack the flexibility necessary for describing the type of proton transfer involved. We conclude that α-helices are essential parts of enzymes and that they play a significant role in the catalytic process.

P.Th. van Duijnen

Papers

Implementation of reaction field methods in quantum-chemistry computer codes

The direct reaction field hamiltonian: Analysis of the dispersion term and application to the water dimer

Theoretical and experimental studies of the opto-electronic properties of positively charged oligo (phenylene vinylene)s: Effects of chain length and alkoxy substitution

DRF90: a polarizable force field

Active-site α-helix in papain and the stability of the ion pair RS− ··· ImH+.Ab initio molecular orbital study