Response properties of embedded molecules through the polarizable embedding model

The polarizable embedding (PE) model is a fragment-based quantum-classical approach aimed at accurate inclusion of environment effects in quantum-mechanical response property calculations. The aim of this tutorial review is to give insight into the practical use of the PE model. Starting from a set of molecular structures and until you arrive at the final property, there are many crucial details to consider in order to obtain trustworthy results in an efficient manner. To lower the threshold for new users wanting to explore the use of the PE model, we describe and discuss important aspects related to its practical use. This includes directions on how to generate input files and how to run a calculation.     

Liquid-vapor and liquid-liquid phase equilibria of the Brodholt-Sampoli-Vallauri polarizable water model

Liquid-vapor and liquid-liquid phase equilibria of the polarizable Brodholt-Sampoli-Vallauri water model have been investigated by Gibbs ensemble Monte Carlo computer simulations. The coexisting liquid and vapor densities and energy of vaporization of the model is found to be in a reasonable agreement with experimental data in the entire temperature range of liquid-vapor coexistence. The critical temperature and density of the model are found to be 615 K and 0.278 gcm(3), respectively, close to the experimental values of 647.1 K and 0.322 gcm(3). In the supercooled state two distinct liquid-liquid coexistence regions are observed. The existence of liquid-liquid phase separation of a polarizable water model is demonstrated for the first time.     

A classical polarizable model for simulations of water and ice

We develop a classical rigid polarizable model of water for molecular simulations of water and ice. The model uses the Rowlinson five-site geometry: oxygen bearing the Lennard-Jones interaction and linearly polarizable point dipole, two positively charged hydrogens, and two massless negative charges placed symmetrically off oxygen so that the experimental dipole moment is reproduced. The target properties are the densities of water and ice Ih, diffusivity, enthalpies of fusion and vaporization, and the ice Ih melting point. The surface tension at lower temperatures is by 7% underestimated whereas the dielectric constant by 6% overestimated. Diffusivity and viscosity worsen at higher temperatures, although the Stokes radius is overestimated only by 2-7%. The ice Ih melting temperature is 260 K and the temperature of maximum density is 269 K. Rescaling the charges by a factor of 1.01 and Lennard-Jones energy by 1.0201 improves the melting point and energy-related quantities but shifts the agreement of kinetic properties to higher temperatures. For the model we propose abbreviation POL4D.     

A polarizable,charge transfer model of water using the drude oscillator

The transfer of small amounts of charge between neighboring particles can be a significant part of interactions among particles. A model is developed for treating charge transfer (CT) combined with the Drude model for polarizability to create an efficient model for liquid water which includes both CT and polarizability. The model is shown to be accurate for a variety of liquid properties, including the density as a function of temperature and the dielectric constant. A new model for water with CT and polarization is developed and applied to the liquid. The inclusion of CT increases the accuracy of many properties, like the density as a function of temperature, indicating the importance of charge redistribution as induced by other particles. © 2016 Wiley Periodicals, Inc.     

Wavelet formulation of the polarizable continuum model

The first implementation of a wavelet discretization of the Integral Equation Formalism (IEF) for the Polarizable Continuum Model (PCM) is presented here. The method is based on the application of a general purpose wavelet solver on the cavity boundary to solve the integral equations of the IEF‐PCM problem. Wavelet methods provide attractive properties for the solution of the electrostatic problem at the cavity boundary: the system matrix is highly sparse and iterative solution schemes can be applied efficiently; the accuracy of the solver can be increased systematically and arbitrarily; for a given system, discretization error accuracy is achieved at a computational expense that scales linearly with the number of unknowns. The scaling of the computational time with the number of atoms N is formally quadratic but a N^1.5 scaling has been observed in practice. The current bottleneck is the evaluation of the potential integrals at the cavity boundary which scales linearly with the system size. To reduce this overhead, interpolation of the potential integrals on the cavity surface has been successfully used. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Solvation of N3- at the water surface: the polarizable continuum model approach

We present a new quantum mechanical model to introduce Pauli repulsion interaction between a molecular solute and the surrounding solvent in the framework of the Polarizable Continuum Model. The new expression is derived in a way to allow naturally for a position-dependent solvent density. This development makes it possible to employ the derived expression for the calculation of molecular properties at the interface between two different dielectrics. The new formulation has been tested on the azide anion (N3-) for which we have calculated the solvation energy, the dipole moment, and the static polarizability at the interface as a function of the ion position. The calculations have been carried out for different ion-surface orientations, and the results have also been compared with the parallel electrostatic-only solvation model.     

A simple, efficient polarizable coarse-grained water model for molecular dynamics simulations

The development of coarse-grained (CG) models that correctly represent the important features of compounds is essential to overcome the limitations in time scale and system size currently encountered in atomistic molecular dynamics simulations. Most approaches reported in the literature model one or several molecules into a single uncharged CG bead. For water, this implicit treatment of the electrostatic interactions, however, fails to mimic important properties, e.g., the dielectric screening. Therefore, a coarse-grained model for water is proposed which treats the electrostatic interactions between clusters of water molecules explicitly. Five water molecules are embedded in a spherical CG bead consisting of two oppositely charged particles which represent a dipole. The bond connecting the two particles in a bead is unconstrained, which makes the model polarizable. Experimental and all-atom simulated data of liquid water at room temperature are used for parametrization of the model. The experimental density and the relative static dielectric permittivity were chosen as primary target properties. The model properties are compared with those obtained from experiment, from clusters of simple-point-charge water molecules of appropriate size in the liquid phase, and for other CG water models if available. The comparison shows that not all atomistic properties can be reproduced by a CG model, so properties of key importance have to be selected when coarse graining is applied. Yet, the CG model reproduces the key characteristics of liquid water while being computationally 1-2 orders of magnitude more efficient than standard fine-grained atomistic water models.     

From dimer to condensed phases at extreme conditions: accurate predictions of the properties of water by a Gaussian charge polarizable model

Water exhibits many unusual properties that are essential for the existence of life. Water completely changes its character from ambient to supercritical conditions in a way that makes it possible to sustain life at extreme conditions, leading to conjectures that life may have originated in deep-sea vents. Molecular simulation can be very useful in exploring biological and chemical systems, particularly at extreme conditions for which experiments are either difficult or impossible; however this scenario entails an accurate molecular model for water applicable over a wide range of state conditions. Here, we present a Gaussian charge polarizable model (GCPM) based on the model developed earlier by Chialvo and Cummings [Fluid Phase Equilib. 150, 73 (1998)] which is, to our knowledge, the first that satisfies the water monomer and dimer properties, and simultaneously yields very accurate predictions of dielectric, structural, vapor-liquid equilibria, and transport properties, over the entire fluid range. This model would be appropriate for simulating biological and chemical systems at both ambient and extreme conditions. The particularity of the GCPM model is the use of Gaussian distributions instead of points to represent the partial charges on the water molecules. These charge distributions combined with a dipole polarizability and a Buckingham exp-6 potential are found to play a crucial role for the successful and simultaneous predictions of a variety of water properties. This work not only aims at presenting an accurate model for water, but also at proposing strategies to develop classical accurate models for the predictions of structural, dynamic, and thermodynamic properties.     

On the calculations of interaction energies and induced electric properties within the polarizable continuum model

In this work we investigate the influence of a polarizable environment on the interaction energies and the interaction-induced (excess) static electric dipole properties for the selected model hydrogen-bonded complexes. The excess properties were estimated for water and hydrogen fluoride dimers using the supermolecular approach and assuming the polarizable continuum model (PCM) as a representation of the polarizable environment. We analyze in this context the performance of the counterpoise correction and the consequences of various possible monomer cavity choices. The polarizable environment reduces the absolute magnitudes of interaction energies and interaction-induced dipole moments, whereas an increase is observed for the absolute magnitudes of induced polarizabilities and first hyperpolarizabilities. Our results indicate that the use of either monomeric (MC) or dimeric (DC) cavities in calculations of monomer properties does not change qualitatively the resultant excess properties. We conclude that the DC scheme is more consistent with the definition of the interaction energy and consequently also the interaction-induced property, whereas the MC scheme corresponds to the definition of stabilization energy. Our results indicate also a good performance of the counterpoise correction scheme for the self-consistent methods in the case of all studied properties.     

Density maximum and polarizable models of water

To estimate accurately the density of water over a wide range of temperatures with a density maximum at 4?°C is one of the most stringent tests of molecular models. The shape of the curve influences the ability to describe critical properties and to predict the freezing temperature. While it was demonstrated that with a proper parameter fit nonpolarizable models can approximate this behavior accurately, it is much more difficult to do this for polarizable models. We provide a short overview of ρ-T diagrams for existing models, then we give an explanation of this difficulty. We present a version of the BK model [A. Baranyai and P. T. Kiss, J. Chem. Phys. 133, 144109 (2010); and ibid. 135, 234110 (2011)] which is capable to predict the density of water over a wide range of temperature. The BK model uses the charge-on-spring method with three Gaussian charges. Since the experimental dipole moment and the geometry is fixed, and the quadrupole moment is approximated by a least mean square procedure, parameters of the repulsion and dispersive attraction forces remained as free tools to match experimental properties. Relying on a simplified but plausible justification, the new version of the model uses repulsion and attraction as functions of the induced dipole moment of the molecule. The repulsive force increases, while the attractive force decreases with the size of the molecular dipole moment. At the same time dipole moment dependent dispersion forces are taking part in the polarization of the molecule. This scheme iterates well and, in addition to a reasonable density-temperature function, creates dipole distributions with accurate estimation of the dielectric constant of the liquid.     

A polarizable model for ethylene oxide

A series of interaction models for ethylene oxide are developed for use in molecular simulation of the thermal properties of both the gas and liquid phases. While it is possible to develop nonpolarizable models to accurately generate either the gas or liquid properties separately, it was not possible to do so using a single model for both phases. A polarizable, rigid all-atom model was developed that reproduces the temperature dependence of the second virial coefficient B(T) and the pressure of the liquid at ambient conditions. The model consists of Lennard-Jones and Coulomb interactions between intermolecular atomic sites plus a scalar polarizability located at the midpoint of the line joining the carbon sites. The electrostatic charges and the polarizability are set to match the experimentally determined dipole and quadrupole moments and the molecular polarizability.     

Coupled cluster theory with the polarizable continuum model of solvation

The environment may significantly affect molecular properties. Thus, it is desirable to account explicitly for these effects on the wave function and its derivatives, especially when the latter are evaluated with accurate methods, such as those belonging to coupled cluster (CC) theory. In this tutorial review, we discuss how to combine CC methods with the polarizable continuum model of solvation (PCM). We describe useful approximations that include the solvent response to the correlation and excited state equations while maintaining the computational cost comparable to in vacuo calculations. Although applied to PCM, the theoretical framework presented in this review is general and can be used with any polarizable embedding model. Representative applications of the CC-PCM method to ground and excited state properties of solvated molecules are presented, and comparisons with experiment, and between the full and approximate schemes are discussed.     

Comparison of polarizable continuum model and quantum mechanics/molecular mechanics solute electronic polarization: study of the optical and magnetic properties of diazines in water

A combination of the polarizable continuum model (PCM) and the hybrid quantum mechanics/molecular mechanics (QM/MM) methodology, PCM-MM/QM, is used to include the solute electronic polarization and then study the solvent effects on the low-lying n→π(?) excitation energy and the (15)N nuclear magnetic shielding of pyrazine and pyridazine in aqueous environment. The results obtained with PCM-MM/QM are compared with two other procedures, i.e., the conventional PCM and the iterative and sequential QM/MM (I-QM/MM). The QM calculations are made using density functional theory in the three procedures. For the excitation energies, the time-dependent B3LYP/6-311+G(d) model is used. For the magnetic shielding, the B3LYP/aug-pcS2(N)/pcS2(C,O,H) is used with the gauge-including atomic orbitals. In both cases, i.e., PCM-MM/QM and I-QM/MM, that use a discrete model of the solvent, the solute is surrounded by a first shell of explicit water molecules embedded by an electrostatic field of point charges for the outer shells. The best results are obtained including 28 explicit water molecules for the spectral calculations and 9 explicit water molecules for the magnetic shielding. Using the PCM-MM/QM methodology the results for the n→π(?) excitation energies of pyridazine and pyrazine are 32,070 ± 80 cm(-1) and 32,675 ± 60 cm(-1), respectively, in good agreement with the corresponding I-MM/QM results of 32,540 ± 80 cm(-1) and 32,710 ± 60 cm(-1) and the experimental results of 33,450-33,580 cm(-1) and 32,700-33,300 cm(-1). For the (15)N magnetic shielding, the corresponding numbers for the gas-water shifts obtained with PCM-MM/QM are 47.4 ± 1.3 ppm for pyridazine and 19.7 ± 1.1 ppm for pyrazine, compared with the I-QM/MM values of 53.4?±?1.3 ppm and 19.5 ± 1.2 ppm and the experimental results of 42-54 ppm and 17-22 ppm, respectively. The agreement between the two procedures is found to be very good and both are in agreement with the experimental values. PCM-MM/QM approach gives a good solute polarization and could be considered in obtaining reliable results within the expected QM/MM accuracy. With this electronic polarization, the solvent effects on the electronic absorption spectra and the (15)N magnetic shielding of the diazines in water are well described by using only an electrostatic approximation. Finally, it is remarked that the experimental and theoretical results suggest that the (15)N nuclear magnetic shielding of any diazine has a clear dependence with the solvent polarity but not directly with the solute-solvent hydrogen bonds.     

A polarizable ion model for the structure of molten CuI

The results are reported of the molecular dynamics simulations of the coherent static structure factor of molten CuI at 938 K using a polarizable ion model. This model is based on a rigid ion potential to which the many body interactions due to the anions induced polarization are added. The calculated structure factor reproduces the clear sharp prepeak observed in neutron diffraction data. The corresponding partial structure factors and the related radial distribution functions calculated by molecular dynamics are compared with those found in the literature derived from a combination of neutron and x-ray diffraction data with the aid of the reverse Monte Carlo simulation technique, as well as those calculated by ab initio MD simulations.     

Gauge-origin-independent magnetizabilities of solvated molecules using the polarizable continuum model

We present an implementation of the polarizable continuum model in its integral equation formulation for the calculation of the magnetizabilities of solvated molecules. The gauge-origin independence of the calculated magnetizabilities and the fast basis set convergence are ensured through the use of London atomic orbitals. Our implementation can use Hartree-Fock and multiconfigurational self-consistent-field (MCSCF) wave functions as well as density-functional theory including hybrid functionals such as B3LYP. We present the results of dielectric continuum effects on water and pyridine using MCSCF wave functions, as well as dielectric medium effects on the magnetizability of the aromatic amino acids as a model for how a surrounding protein environment affects the magnetizability of these molecules. It is demonstrated that the dielectric medium effects on the magnetizability anisotropies of the aromatic amino acids may be substantial, being as large as 25% in the case of tyrosine.     

Electronic g-tensors of solvated molecules using the polarizable continuum model

We present the implementation of density functional response theory combined with the polarizable continuum model (PCM), enabling first principles calculations of molecular g-tensors of solvated molecules. The calculated g-tensor shifts are compared with experimental g-tensor shifts obtained from electron paramagnetic resonance spectra for a few solvated species. The results indicate qualitative agreement between the calculations and the experimental data for aprotic solvents, whereas PCM fails to reproduce the electronic g-tensor behavior for protic solvents. This failure of PCM for protic solvents can be resolved by including into the model those solvent molecules which are involved in hydrogen bonding with the solute. The results for the protic solvents show that the explicit inclusion of the solvent molecules of the first solvation sphere is not sufficient in order to reproduce the behavior of the electronic g-tensor in protic solvents, and that better agreement with experimental data can be obtained by including the long-range electrostatic effects accounted for by the PCM approach on top of the explicit hydrogen-bonded complexes.     

A polarizable ion model for the structure of molten AgI

The results are reported of the molecular dynamics simulations of the coherent static structure factor of molten AgI at 923 K using a polarizable ion model. This model is based on a rigid ion potential, to which the many body interactions due to the anions induced polarization are added. The calculated structure factor is in better agreement with recent neutron diffraction data than that obtained by using simple rigid ion pair potentials. The Voronoi-Delaunay method has been applied to study the relationship between voids in the spatial distribution of cations and the prepeak of the structure factor.     

Solvent effect on the stability and properties of platinum-substituted borirene and boryl isomers: The polarizable continuum model

The structure and properties of platinum borirene complex trans-[Cl(PMe₃)₂Pt{μ-BN(SiMe₃)₂C=C}Ph] and its isomer the platinum boryl complex trans–[Cl(PMe₃)₂PtBN(SiMe₃)₂C≡CPh] were investigated theoretically. The solvent effect on the stability, structural parameters, frontier orbital energies, HOMO–LUMO gaps, and hardness of isomers was investigated using the polarizable continuum model (PCM). It was found that borirene isomer is the most stable isomer in the gas phase and solvent. The calculated results show that the presence of solvent reduces the frontier orbital energy of the studied molecules. Geometries obtained from calculations were used to perform NBO analysis.     

The ice-vapor interface and the melting point of ice I(h) for the polarizable POL3 water model

We use molecular dynamics simulations to determine the melting point of ice I(h) for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys.1992, 97, 2659). Simulations are performed on a slab of ice I(h) with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields the melting point of the POL3 model to be T(m) = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.     

Charge-on-spring polarizable water models revisited: from water clusters to liquid water to ice

The properties of two improved versions of charge-on-spring (COS) polarizable water models (COS/G2 and COS/G3) that explicitly include nonadditive polarization effects are reported. In COS models, the polarization is represented via a self-consistently induced dipole moment consisting of a pair of separated charges. A previous polarizable water model (COS/B2), upon which the improved versions are based, was developed by Yu, Hansson, and van Gunsteren. To improve the COS/B2 model, which overestimated the dielectric permittivity, one additional virtual atomic site was used to reproduce the water monomer quadrupole moments besides the water monomer dipole moment in the gas phase. The molecular polarizability, residing on the virtual atomic site, and Lennard-Jones parameters for oxygen-oxygen interactions were varied to reproduce the experimental values for the heat of vaporization and the density of liquid water at room temperature and pressure. The improved models were used to study the properties of liquid water at various thermodynamic states as well as gaseous water clusters and ice. Overall, good agreement is obtained between simulated properties and those derived from experiments and ab initio calculations. The COS/G2 and COS/G3 models may serve as simple, classical, rigid, polarizable water models for the study of organic solutes and biopolymers. Due to its simplicity, COS type of polarization can straightforwardly be used to introduce explicit polarization into (bio)molecular force fields.