PEACH 4 with ABINIT-MP: a general platform for classical and quantum simulations of biological molecules

A program package for molecular simulations of biological molecules was developed. The package, "PEACH version 4 with ABINIT-MP version 20021029," was constructed by incorporating ABINIT-MP, a program for the fragment molecular orbital (FMO) method [Chem.Phys. Lett. 313 (1999) 701], into PEACH, a program package for classical molecular dynamics simulations (MD). A few capabilities of the package were demonstrated. First, high parallel efficiency of FMO was demonstrated in a single point calculation of a protein. Second,FMO-MD simulations [Chem. Phys. Lett. 372 (2003) 342] of a peptide were performed with and without explicit solvent, and the simulations showed the influence of the solvent on the electronic state of the peptide.     

On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations

We have implemented the accelerated molecular dynamics approach (Hamelberg, D.; Mongan, J.; McCammon, J. A. J. Chem. Phys. 2004, 120 (24), 11919) in the framework of ab initio MD (AIMD). Using three simple examples, we demonstrate that accelerated AIMD (A-AIMD) can be used to accelerate solvent relaxation in AIMD simulations and facilitate the detection of reaction coordinates: (i) We show, for one cyclohexane molecule in the gas phase, that the method can be used to accelerate the rate of the chair-to-chair interconversion by a factor of ～1 × 10(5), while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H(2)O molecules, that A-AIMD can also be used in the condensed phase to accelerate the sampling of water conformations, without affecting the structural properties of the solvent; and (iii) The method is then used to compute the potential of mean force (PMF) for the dissociation of Na-Cl in water, accelerating the convergence by a factor of ～3-4 compared to conventional AIMD simulations.(2) These results suggest that A-AIMD is a useful addition to existing methods for enhanced conformational and phase-space sampling in solution. While the method does not make the use of collective variables superfluous, it also does not require the user to define a set of collective variables that can capture all the low-energy minima on the potential energy surface. This property may prove very useful when dealing with highly complex multidimensional systems that require a quantum mechanical treatment.     

An application of coupled reference interaction site model/molecular dynamics to the conformational analysis of the alanine dipeptide

We present an application of our recently proposed coupled reference interaction site model (RISM) molecular dynamics (MD) solvation free energy methodology [Freedman and Truong, Chem. Phys. Lett. 381, 362 (2003); J. Chem. Phys. 121, 2187 (2004)] to study the conformational stability of alanine dipeptide in aqueous solution. In this methodology, radial distribution functions obtained from a single MD simulation are substituted into a RISM expression for solvation free energy. Consequently, iterative solution of the RISM equation is not needed. The relative solvation free energies of seven different conformations of the alanine dipeptide in aqueous solution are calculated. Results from the coupled RISM/MD methodology are in good agreement with those from earlier simulations using the accurate free energy perturbation approach, showing that the alphaR conformation is most stabilized by solution. This study establishes a framework for applying this coupled RISM/MD method to larger biological systems.     

Overcoming free energy barriers using unconstrained molecular dynamics simulations

Association of unconstrained molecular dynamics (MD) and the formalisms of thermodynamic integration and average force [Darve and Pohorille, J. Chem. Phys. 115, 9169 (2001)] have been employed to determine potentials of mean force. When implemented in a general MD code, the additional computational effort, compared to other standard, unconstrained simulations, is marginal. The force acting along a chosen reaction coordinate xi is estimated from the individual forces exerted on the chemical system and accumulated as the simulation progresses. The estimated free energy derivative computed for small intervals of xi is canceled by an adaptive bias to overcome the barriers of the free energy landscape. Evolution of the system along the reaction coordinate is, thus, limited by its sole self-diffusion properties. The illustrative examples of the reversible unfolding of deca-L-alanine, the association of acetate and guanidinium ions in water, the dimerization of methane in water, and its transfer across the water liquid-vapor interface are examined to probe the efficiency of the method.     

Conformational landscape of (R,R)-pterocarpans with biological activity in vacuo and in aqueous solution (PCM and/or water clusters)

All possible combinations of stable dihedral values have been considered in vacuo at the B3LYP/6-31G level for 3,9-dihydroxy-4,8-diprenylpterocarpan (erybraedin C), whose hydroxy out-out conformation had been examined earlier together with the conformational preferences of 3,9-dimethoxy-4-prenylpterocarpan (bitucarpin A) at the same level (Phys. Chem. Chem. Phys. 2004, 6, 2849). The structure with O5 trans with respect to H6a (O(t)) is about 2 kcal/mol less stable in vacuo than that with one of the H6 trans to it (H(t)); in aqueous solution its energy gap is nearly conserved. The in-in arrangement of the hydroxyl groups of erybraedin turns out to be preferred in vacuo (even considering zero point and thermal effects), where pseudo H-bonds are formed between hydroxy hydrogens and pi electron distributions of prenyl groups. The continuum solvent effect (water) at the IEF-PCM/B3LYP/6-31G level on the relative stability of the various rotamers is very limited both on bitucarpin and erybraedin. Considering the dihydrated derivatives, significant differences in the solvation energy are found between the distinct hydration sites, increasing in the order: methoxy O, ring O, hydroxy O, and hydroxy H. In hydroxy-water interactions, in fact, water prefers to behave as an H-bond acceptor unless nearby bulky groups prevent its approach. Interestingly enough, a bridging water molecule between the hydroxy H of erybraedin and the prenyl group can be found. The inclusion of BSSE corrections in hydroxy-water interactions decidedly favors out-out hydrated arrangements, followed by out-in and in-out ones. Bulk solvent effects with IEF-PCM about the dihydrated systems almost invert the stability order found in vacuo. When a four-water cluster is considered using QM methods, waters gather in H-bonded pairs around the solute OH groups. MD simulations, carried out on a pterocarpan solute (J. Phys. Chem. B 2005, 109, 16918), supply water adducts consistent with a liquid state that have also been embedded in the continuum solvent.     

Hybrid particle‐field molecular dynamics simulations: Parallelization and benchmarks

The parallel implementation of a recently developed hybrid scheme for molecular dynamics (MD) simulations (Milano and Kawakatsu, J Chem Phys 2009, 130, 214106) where self‐consistent field theory (SCF) and particle models are combined is described. Because of the peculiar formulation of the hybrid method, considering single particles interacting with density fields, the most computationally expensive part of the hybrid particle‐field MD simulation can be efficiently parallelized using a straightforward particle decomposition algorithm. Benchmarks of simulations, including comparisons of serial MD and MD‐SCF program profiles, serial MD‐SCF and parallel MD‐SCF program profiles, and parallel benchmarks compared with efficient MD program GROMACS 4.5.4 are tested and reported. The results of benchmarks indicate that the proposed parallelization scheme is very efficient and opens the way to molecular simulations of large scale systems with reasonable computational costs. © 2012 Wiley Periodicals, Inc.     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Molecular dynamics simulations of porphyrin-dendrimer systems: toward modeling electron transfer in solution

We have performed computational simulations of porphyrin-dendrimer systems--a cationic porphyrin electrostatically associated to a negatively charged dendrimer--using the method of classical molecular dynamics (MD) with an atomistic force field. Previous experimental studies have shown a strong quenching effect of the porphyrin fluorescence that was assigned to electron transfer (ET) from the dendrimer's tertiary amines (Paulo, P. M. R.; Costa, S. M. B. J. Phys. Chem. B 2005, 109, 13928). In the present contribution, we evaluate computationally the role of the porphyrin-dendrimer conformation in the development of a statistical distribution of ET rates through its dependence on the donor-acceptor distance. We started from simulations without explicit solvent to obtain trajectories of the donor-acceptor distance and the respective time-averaged distributions for two dendrimer sizes and different initial configurations of the porphyrin-dendrimer pair. By introducing explicit solvent (water) in our simulations, we were able to estimate the reorganization energy of the medium for the systems with the dendrimer of smaller size. The values obtained are in the range 0.6-1.5 eV and show a linear dependence with the inverse of the donor-acceptor distance, which can be explained by a two-phase dielectric continuum model taking into account the medium heterogeneity provided by the dendrimer organic core. Dielectric relaxation accompanying ET was evaluated from the simulations with explicit solvent showing fast decay times of some tens of femtoseconds and slow decay times in the range of hundreds of femtoseconds to a few picoseconds. The variations of the slow relaxation times reflect the heterogeneity of the dendrimer donor sites which add to the complexity of ET kinetics as inferred from the experimental fluorescence decays.     

Molecular dynamics simulations of helix-forming, amine-functionalized m-poly(phenyleneethynylene)s

We present here the results of all-atom and united-atom molecular dynamics (MD) simulations that were used to examine the folding behavior of an amine-functionalized m-poly(phenyleneethynylene) (m-PPE) oligomer in aqueous environment. The parallelized GROMACS MD simulation code and OPLS force field were used for multiple MD simulations of m-PPE oligomers containing 24 phenyl rings in extended, coiled and helix conformations separately in water to determine the minimum energy conformation of the oligomer in aqueous solvent and what interactions are most important in determining this structure. Simulation results showed that the helix is the preferred minimum energy conformation of a single oligomer in water and that Lennard-Jones interactions are the dominant forces for the stabilization of the helix. In addition, these solvophobic interactions are strong enough to maintain the helix conformation at temperatures up to 523 K.     

FMO-MD simulations on the hydration of formaldehyde in water solution with constraint dynamics

Full-quantum mechanical fragment molecular orbital-based molecular dynamics (FMO-MD) simulations were applied to the hydration reaction of formaldehyde in water solution under neutral conditions. Two mechanisms, a concerted and a stepwise one, were considered with respect to the nucleophilic addition and the proton transfer. Preliminary molecular orbital calculations by means of polarized continuum model reaction field predicted that the hydration prefers a concerted mechanism. Because the calculated activation barriers were too high for free FMO-MD simulations to give reactive trajectories spontaneously, a More O'Ferrall-Jencks-type diagram was constructed from the statistical analysis of the FMO-MD simulations with constraint dynamics. The diagram showed that the hydration proceeds through a zwitterionic-like (ZW-like) structure. The free energy changes along the reaction coordinate calculated by means of the blue moon ensemble for the hydration and the amination of formaldehyde indicated that the hydration proceeds through a concerted process through the ZW-like structure, whereas the amination goes through a stepwise mechanism with a ZW intermediate. In inspection of the FMO-MD trajectories, water-mediated cyclic proton transfers were observed in both reactions on the way from the ZW-like structure to the product. These proton transfers also have an asynchronous character, in which deprotonation from the nucleophilic oxygen atom (or nitrogen atom for amination) precedes the protonation of the carbonyl oxygen atom. The results showed the strong advantage of the FMO-MD simulations to obtain detailed information at a molecular level for solution reactions.     

Application of the level-set method to the implicit solvation of nonpolar molecules

A level-set method is developed for numerically capturing the equilibrium solute-solvent interface that is defined by the recently proposed variational implicit solvent model [Dzubiella, Swanson, and McCammon, Phys. Rev. Lett. 104, 527 (2006); J. Chem. Phys. 124, 084905 (2006)]. In the level-set method, a possible solute-solvent interface is represented by the zero level set (i.e., the zero level surface) of a level-set function and is eventually evolved into the equilibrium solute-solvent interface. The evolution law is determined by minimization of a solvation free energy functional that couples both the interfacial energy and the van der Waals type solute-solvent interaction energy. The surface evolution is thus an energy minimizing process, and the equilibrium solute-solvent interface is an output of this process. The method is implemented and applied to the solvation of nonpolar molecules such as two xenon atoms, two parallel paraffin plates, helical alkane chains, and a single fullerence C(60). The level-set solutions show good agreement for the solvation energies when compared to available molecular dynamics simulations. In particular, the method captures solvent dewetting (nanobubble formation) and quantitatively describes the interaction in the strongly hydrophobic plate system.     

Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (～7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.     

Molecular level approaches for investigation of electron transfer in nonpolar solvents

The authors extend their previous work published in Leontyev and TachiyaJ. Chem. Phys. 123, 224502 (2005) and study not only forward but also reverse electron transfer between pyrene and dimethylaniline in a nonpolar solvent, n-hexane. The distribution function methodology and molecular dynamics technique adopted in their previous work are used. Two algorithms (I and II) are formulated for obtaining the reorganization energy and the solvation free energy difference in the linear response approximation. The two algorithms are combined with different cutoff schemes and tested for polarizable and nonpolarizable solvent models. Agreement between the results obtained by the two algorithms was achieved only for simulations employing the particle mesh Ewald treatment. It is concluded that algorithm I provides a reliable scheme for evaluation of the reorganization energy and the solvation free energy difference. Moreover, a new algorithm referred to as the G-function algorithm is formulated which does not assume the linear response approximation, and is tested on evaluation of the solvation free energy difference. Agreement between the results from the G-function algorithm and those from algorithms I and II is fairly good, although it depends on the degree of statistical consistency of the simulations. In the case of nonpolar solvents the G-function method has practical importance because, unlike the conventional thermodynamic integration approach, it requires equilibrium molecular configuration ensembles only for the initial and final states of the system.     

Incorporating the excluded solvent volume and surface charges for computing solvation free energy

Gauss's law or Poisson's equation is conventionally used to calculate solvation free energy. However, the near‐solute dielectric polarization from Gauss's law or Poisson's equation differs from that obtained from molecular dynamics (MD) simulations. To mimic the near‐solute dielectric polarization from MD simulations, the first‐shell water was treated as two layers of surface charges, the densities of which are proportional to the electric field at the solvent molecule that is modeled as a hard sphere. The intermediate water was treated as a bulk solvent. An equation describing the solvation free energy of ions using this solvent scheme was derived using the TIP3P water model. © 2013 Wiley Periodicals, Inc.     

Evidence in support of levitation effect as the reason for size dependence of ionic conductivity in water: a molecular dynamics simulation

We report extensive molecular dynamics simulations of (i) model ions in water at high concentrations as a function of the size and charge of the ion as well as (ii) realistic simulation of Cl- and Br- ions at low concentrations in water at room temperature. We also analyze existing experimental data in light of the results obtained here. The halide ion simulations have been carried out using the interaction potentials of Koneshan et al. (J. Phys. Chem. B 1998, 102, 4193). We compute structural and dynamical properties of ions in water and explore their variation with size and charge of the ion. We find that ions of certain intermediate sizes exhibit a maximum in self-diffusivity in agreement with previous experimental measurements and computer simulations. We analyze molecular dynamics trajectories in light of the previous understanding of the levitation effect (LE) and the recent suggestion that ionic conductivity has its origin in LE (J. Phys. Chem. B 2005, 109, 8120). We report the distribution of void and neck radii that exist amidst water. Our analysis suggests that the ion with maximum self-diffusivity is characterized by a lower activation energy and a single-exponential decay of F(s)(k,t). The behavior of these and other related quantities of the ion with maximum self-diffusivity are characteristic of the anomalous regime of the LE. The simulation results of Br- and Cl- ions in water also yield results in agreement with the predictions of LE. A plot of experimental conductivity data in the literature for alkali ions in water by Kay and Evans (J. Phys. Chem. 1966, 70, 2325) also yields a lower activation energy for the ion with maximum conductivity in excellent agreement with the LE. To the best of our knowledge, none of the existing theories predict a lower activation energy for the ion with maximum conductivity.     

Electrostatics of proteins in dielectric solvent continua. I. Newton's third law marries qE forces

The authors reformulate and revise an electrostatic theory treating proteins surrounded by dielectric solvent continua [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)] to make the resulting reaction field (RF) forces compatible with Newton's third law. Such a compatibility is required for their use in molecular dynamics (MD) simulations, in which the proteins are modeled by all-atom molecular mechanics force fields. According to the original theory the RF forces, which are due to the electric field generated by the solvent polarization and act on the partial charges of a protein, i.e., the so-called qE forces, can be quite accurately computed from Gaussian RF dipoles localized at the protein atoms. Using a slightly different approximation scheme also the RF energies of given protein configurations are obtained. However, because the qE forces do not account for the dielectric boundary pressure exerted by the solvent continuum on the protein, they do not obey the principle that actio equals reactio as required by Newton's third law. Therefore, their use in MD simulations is severely hampered. An analysis of the original theory has led the authors now to a reformulation removing the main difficulties. By considering the RF energy, which represents the dominant electrostatic contribution to the free energy of solvation for a given protein configuration, they show that its negative configurational gradient yields mean RF forces obeying the reactio principle. Because the evaluation of these mean forces is computationally much more demanding than that of the qE forces, they derive a suggestion how the qE forces can be modified to obey Newton's third law. Various properties of the thus established theory, particularly issues of accuracy and of computational efficiency, are discussed. A sample application to a MD simulation of a peptide in solution is described in the following paper [M. Stork and P. Tavan, J. Chem. Phys., 126, 165106 (2007).     

High resolution Dopplerimetry of correlated angular and quantum state-resolved CO(2) scattering dynamics at the gas-liquid interface

Full three dimensional (3D) translational distributions for quantum state-resolved scattering dynamics at the gas-liquid interface are presented for experimental and theoretical studies of CO(2) + perfluorinated surfaces. Experimentally, high resolution absorption profiles are measured as a function of incident (θ(inc)) and scattering (θ(scat)) angles for CO(2) that has been scattered from a 300 K perfluorinated polyether surface (PFPE) with an incident energy of E(inc) = 10.6(8) kcal mol(-1). Line shape analysis of the absorption profiles reveals non-equilibrium dynamics that are characterized by trapping-desorption (TD) and impulsive scattering (IS) components, with each channel simply characterized by an effective "temperature" that compares very well with previous results from rotational state analysis [Perkins and Nesbitt, J. Phys. Chem. A, 2008, 112, 9324]. From a theoretical perspective, molecular dynamics (MD) simulations of CO(2) + fluorinated self-assembled monolayer surface (F-SAMs) yield translational probability distributions that are also compared with experimental results. Trajectories are parsed by θ(scat) and J, with the results rigorously corrected by flux-to-density transformation and providing comparisons in near quantitative agreement with experiment. 3D flux and velocity distributions obtained from MD simulations are also presented to illustrate the role of in- and out-of-plane scattering.     

Can the anomalous aqueous solubility of beta-cyclodextrin be explained by its hydration free energy alone?

The well-documented anomalous solubility of beta-cyclodextrin (beta-CD), relative to alpha- and gamma-CD, has been examined by Naidoo et al. (J. Phys. Chem. B, 2004, 108, 4236-4238.) from the perspective of water organization and internal motion of the macrocyclic rings. Whether modulation in the hydration patterns and in the rigidity of the molecular scaffold can be reconciled with the hydration free energy of beta-CD to rationalize its notorious low solubility remains open to further investigation. In this contribution, multi-nanosecond molecular dynamics (MD) simulations have been carried out to investigate the hydration process of alpha-, beta- and gamma-CD. The distribution of water molecules involved in this process and the linearity of intramolecular hydrogen bonds have been analyzed. The results reported here demonstrate that the anomalous solubility for beta-CD can be essentially rationalized by its greater rigidity conferred by the participating intramolecular hydrogen bonds and the higher density of water molecules of lesser mobility. The hydration free energy of alpha-, beta- and gamma-CD was computed using the free energy perturbation method. This quantity is shown to increase with the number of glucose units, thereby suggesting that the anomalous solubility of beta-CD cannot be explained by its free energy of hydration alone.     

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. I. Theory and phase-space analysis

Current constant pressure molecular-dynamics (MD) algorithms are not consistent with the recent reformulation of the isothermal-isobaric (NpT) ensemble. The NpT ensemble partition function requires the use of a "shell" molecule to identify uniquely the volume of the system, thereby avoiding the redundant counting of configurations [e.g., G. J. M. Koper and H. Reiss, J. Phys. Chem. 100, 422 (1996); D. S. Corti, Phys. Rev. E, 64, 016128 (2001)]. So far, only the NpT Monte Carlo method has been updated to allow the system volume to be defined by a shell particle [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. A shell particle has yet to be incorporated into MD simulations. The proper modification of the NpT MD algorithm is therefore the subject of this paper. Unlike Andersen's method [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] where a piston of unknown mass serves to control the response time of volume fluctuations, the newly proposed equations of motion impose a constant external pressure via the introduction of a shell particle of known mass. Hence, the system itself sets the time scales for pressure and volume fluctuations. The new algorithm is subject to a number of fundamentally rigorous tests to ensure that the equations of motion sample phase space correctly. We also show that the Hoover NpT algorithm [W. G. Hoover, Phys. Rev. A. 31, 1695 (1985); 34, 2499 (1986)] does sample phase correctly, but only when periodic boundary conditions are employed.     

QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010