Piet Th. van Duijnen

A charge analysis derived from an atomic multipole expansion

A new charge analysis is presented that gives an accurate description of the electrostatic potential from the charge distribution in molecules. This is achieved in three steps: first, the total density is written as a sum of atomic densities; next, from these atomic densities a set of atomic multipoles is defined; finally, these atomic multipoles are reconstructed exactly by distributing charges over all atoms. The method is generally applicable to any method able to provide atomic multipole moments, but in this article we take advantage of the way the electrostatic potential is calculated within the Density Functional Theory framework. We investigated a set of 31 molecules as well as all amino acid residues to test the quality of the method, and found accurate results for the molecular multipole moments directly from the DFT calculations. The deviations from experimental values for the dipole/quadrupole moments are also small. Finally, our Multipole Derived Charges reproduce both the atomic and molecular multipole moments exactly.

A discrete solvent reaction field model within density functional theory

In this work we present theory and implementation for a discrete reaction field model within Density Functional Theory (DFT) for studying solvent effects on molecules. The model combines a quantum mechanical (QM) description of the solute and a classical description of the solvent molecules (MM). The solvent molecules are modeled by point charges representing the permanent electronic charge distribution, and distributed polarizabilities for describing the solvent polarization arising from many-body interactions. The QM/MM interactions are introduced into the Kohn–Sham equations, thereby allowing for the solute to be polarized by the solvent and vice versa. Here we present some initial results for water in aqueous solution. It is found that the inclusion of solvent polarization is essential for an accurate description of dipole and quadrupole moments in the liquid phase. We find a very good agreement between the liquid phase dipole and quadrupole moments obtained using the Local Density Approximation and r...

The first hyperpolarizability of p-nitroaniline in 1,4-dioxane: a quantum mechanical/molecular mechanics study.

In this work we have investigated the first hyperpolarizability of pNA in 1,4-dioxane solution using a quantum mechanics/molecular mechanics (QM/MM) model. The particular model adopted is the recently developed discrete solvent reaction field (DRF) model. The DRF model is a polarizable QM/MM model in which the QM part is treated using time-dependent density-functional theory and local-field effects are incorporated. This allows for direct computation of molecular effective properties which can be compared with experimental results. The solvation shift for the first hyperpolarizability is calculated to be 30% which is in good agreement with the experimental results. However, the calculated values, both in the gas phase and in solution, are by a factor of 2 larger than the experimental ones. This is in contrast to the calculation of the first hyperpolarizability for several small molecules in the gas phase where fair agreement is found with experimental. The inclusion of local-field effects in the calculations was found to be crucial and neglecting them led to results which are significantly larger. To test the DRF model the refractive index of liquid 1,4-dioxane was also calculated and found to be in good agreement with experiment.

A discrete solvent reaction field model for calculating molecular linear response properties in solution

A discrete solvent reaction field model for calculating frequency-dependent molecular linear response properties of molecules in solution is presented. The model combines a time-dependent density functional theory (QM) description of the solute molecule with a classical (MM) description of the discrete solvent molecules. The classical solvent molecules are represented using distributed atomic charges and polarizabilities. All the atomic parameters have been chosen so as to describe molecular gas phase properties of the solvent molecule, i.e., the atomic charges reproduce the molecular dipole moment and the atomic polarizabilities reproduce the molecular polarizability tensor using a modified dipole interaction model. The QM/MM interactions are introduced into the Kohn–Sham equations and all interactions are solved self-consistently, thereby allowing for the solute to be polarized by the solvent. Furthermore, the inclusion of polarizabilities in the MM part allows for the solvent molecules to be polarized ...

Iodine-benzene charge-transfer complex: Potential energy surface and transition probabilities studied at several levels of theory

We present the results of detailed studies of the potential energy surfaces of the iodine)benzene charge-transfer complex obtained from fully counterpoise . . corrected ab initio calculations at the second-order Moller)Plesset MP2 level and from .semi- classical calculations. The most stable conformations found were the above-bond and the above-carbon conformations. The axial conformation was found to be somewhat less stable. The remarkable difference in intermolecular distance for different orientations of the iodine is explained in terms of the polarization anisotropy. This feature makes the construction of an accurate classical force field rather difficult because of the marked dependence of the repulsion parameter}usually the radius}for iodine on both orientation and polarization of the iodine. Investigation of the oscillator strengths of different complex geometries shows that there are many conformations in which the charge-transfer excitation can take place. Q 1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 709)723, 1999

Polarizabilities in the condensed phase and the local fields problem: A direct reaction field formulation

A consistent derivation is given for local field factors to be used for correcting measured or calculated static (hyper)polarizabilities in the condensed phases. We show how local fields should be used in the coupled perturbative Hartree–Fock or finite field methods for calculating these properties, specifically for the direct reaction field (DRF) approach, in which a quantum chemically treated “solute” is embedded in a classical “solvent” mainly containing discrete molecules. The derivation of the local fields is based on a strictly linear response of the classical parts and they are independent of any quantum mechanical method to be used. In applications to two water dimers in two basis sets it is shown that DRF matches fully quantum mechanical results quite well. For acetone in eleven different solvents we find that if the solvent is modeled by only a dielectric continuum (hyper)polarizabilities increase with respect to their vacuum values, while with the discrete model they decrease. We show that the use of the Lorentz field factor for extracting (hyper)polarizabilities from experimental susceptibilities may lead to serious errors.

Medium perturbations on the molecular polarizability calculated within a localized dipole interaction model

We have studied the medium effects on the frequency-dependent polarizability of water by separating the total polarizability of water clusters into polarizabilities of the individual water molecules. A classical frequency-dependent dipole–dipole interaction model based on classical electrostatics and an Unsold dispersion formula has been used. It is shown that the model reproduces the polarizabilities of small water clusters calculated with time-dependent density functional theory. A comparison between supermolecular calculations and the localized interaction model illustrate the problems arising from using supermolecular calculations to predict the medium perturbations on the solute polarizability. It is also noted that the solute polarizability is more dependent on the local geometry of the cluster than on the size of the cluster.

Solvatochromism of the ?* ?n transition of acetone by combined quantum mechanical?classical mechanical calculations

The solvent shift of the π* ← n transition of acetone in water, acetonitrile, and tetrachloromethane was calculated in a combined quantum mechanical—classical mechanical approach, using both dielectric continuum and explicit, polarizable molecular solvent models. The explicit modeling of solvent polarizability allows for a separate analysis of electrostatic, induction, and dispersion contributions to the shifts. The calculations confirm the qualitative theories about the mechanisms behind the blue shift in polar solvents and the red shift in nonpolar solvents, the solvation of the ground state due to electrostatic interactions being preferential in the former, and favorable dispersion interaction with the excited state, in the latter case. Good quantitative agreement for the solvent shift between experiment (+1,700, +400, and −350 cm−1 in water, acetonitrile, and tetrachloromethane, respectively) and the explicit solvent model (+1,821, +922, and −381 cm−1) was reached through a modest Monte Carlo sampling of the solvent degrees of freedom. A consistent treatment of the solvent could only be realized in the molecular solvent model. The dielectric‐only model needs reparameterization for each solvent

A discrete solvent reaction field model for calculating frequency-dependent hyperpolarizabilities of molecules in solution

We present a discrete solvent reaction field (DRF) model for the calculation of frequency-dependent hyperpolarizabilities of molecules in solution. In this model the solute is described using density functional theory (DFT) and the discrete solvent molecules are described with a classical polarizable model. The first hyperpolarizability is obtained in an efficient way using time-dependent DFT and the (2n+1) rule. The method was tested for liquid water using a model in which a water molecule is embedded in a cluster of 127 classical water molecules. The frequency-dependent first and second hyperpolarizabilities related to the electric field induced second harmonic generation (EFISH) experiment, were calculated both in the gas phase and in the liquid phase. For water in the gas phase, results are obtained in good agreement with correlated wave function methods and experiments by using the so-called shape-corrected exchange correlation (xc)-potentials. In the liquid phase the effect of using asymptotically c...

Time-dependent density functional study of the static second hyperpolarizability of BB-, NN- and BN-substituted C60

In this work we have investigated the effects of substituting carbon atoms with B and N on the second hyperpolarizability of C60 using time-dependent density functional theory. We have calculated the second hyperpolarizability of the double substitute-doped fullerenes C58NN, C58BB and C58BN. For C60 only small changes in the second hyperpolarizability were found when doping with either 2B or 2N. However, by doping C60 with both B and N, creating an donor-acceptor system, an increase in the second hyperpolarizability with about 50% was found.