Theory of solvation in polar nematics

We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phases. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates the experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra and makes several experimentally testable predictions. The equilibrium free energy of solvation, anisotropic in the nematic phase, is given by a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for the equilibrium Stokes shift and the Stokes shift time correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl and 4,4-n-heptyl-cyanopiphenyl solvents as a function of temperature in both the nematic and isotropic phases.     

Molecular density functional theory of solvation: from polar solvents to water

A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics.     

Towards a modellisation of the solvation energy in multi-component solvents. The interesting case of a charged solute imbedded in a polymer-containing electrolyte solution

In this paper we propose a mean-field theory to calculate the solvation free energy of a charged solute imbedded in a complex multi-component solvent. We considered a solvent made up of a mixture of small (electrolyte solution) and large (polymer) components. The presence of macromolecules ensures reduced mixing entropy among the different solvent components, an effect due to polymer connectivity. The reduced entropy favours strong preferential distribution of a particular solvent even in the presence of weak preferential solute–solvent interactions. In addition, two energy terms must be considered: (a) the interaction between the solute electrostatic potential and the electrolyte solution and (b) the formation of a polymer–solute interface. Because of the different dielectric permittivity of the solvent components, the electrolyte and polymer distribution functions are strongly coupled: ions, indeed, are more solvated in regions of higher local dielectric permittivity arising from the inhomogeneous mixing of solvent and polymer. We combined together the different energy terms in the framework of the de Gennes free energy functional for polymer solutions along with a generalised Poisson–Boltzmann equation developed for inhomogeneous dielectric media. Moreover, the preferential electrolyte solvation in regions of greater polarity was considered by an extension of the Born equation. Setting the polymer dielectric permittivity smaller than the solvent one and making null the specific polymer–solute interactions, we calculated enhanced electrolyte concentration and reduced polymer concentration near the solute surface on raising the solute surface charge density. The theory shows also the breakdown of the widely used separation between electrostatic and surface tension-dependent contributions to solvation energy when non-ideal mixed solvents are considered. In fact, according to the model, the surface tension of such mixed solvents strongly depends on the solute surface charge density: at high potentials the interfacial tension may increase rather than decrease on raising the polymer volume fraction. The theoretical results have been compared with experimental data on polymer+electrolyte solution surface tension and with solubility data of colloidal particles. The comparison evidences the complex behaviour of multi-component solvents going well beyond the trivial weighted average of the dielectric permittivity and surface tension of the isolated chemical components. Deviations from the simple behaviour predicted by an average picture of multi-component solvents could be understood by developing more sophisticated, but still simple, approaches like that proposed in this paper.Contribution to the Jacopo Tomasi Honorary Issue. This paper is dedicated to Jacopo Tomasi. I learned much of the difficult art of transforming complex problems into simple models after reading his early works on solvation energy.     

Activation entropy of electron transfer reactions

We report microscopic calculations of free energies and entropies for intramolecular electron transfer reactions. The calculation algorithm combines the atomistic geometry and charge distribution of a molecular solute obtained from quantum calculations with the microscopic polarization response of a polar solvent expressed in terms of its polarization structure factors. The procedure is tested on a donor–acceptor complex in which ruthenium donor and cobalt acceptor sites are linked by a four-proline polypeptide. The reorganization energies and reaction energy gaps are calculated as a function of temperature by using structure factors obtained from our analytical procedure and from computer simulations. Good agreement between two procedures and with direct computer simulations of the reorganization energy is achieved. The microscopic algorithm is compared to the dielectric continuum calculations. We found that the strong dependence of the reorganization energy on the solvent refractive index predicted by continuum models is not supported by the microscopic theory. Also, the reorganization and overall solvation entropies are substantially larger in the microscopic theory compared to continuum models.     

Solute rotation in polar liquids: microscopic basis for the Stokes-Einstein-Debye model

Here, we develop a framework for a molecular level understanding of the celebrated Stokes-Einstein-Debye (SED) formula. In particular, we explore reasons behind the surprising success of the SED model in describing dipolar solute rotation in complex polar media. Relative importance of solvent viscosity and solute-solvent dipolar interaction is quantified via a self-consistent treatment for the total friction on a rotating solute where the hydrodynamic contribution is modified by the friction arising from the longer ranged solute-solvent dipolar interaction. Although the solute-solvent dipolar coupling is obtained via the Mori-Zwanzig formalism, the inclusion of solvent structure via the wave vector dependent viscosity in the hydrodynamic contribution incorporates solvent molecularity in the present theory. This approach satisfactorily describes the experimental rotation times measured using a dipolar solute, coumarin 153 (C153), in protic and aprotic polar liquids, and more importantly, provides microscopic explanation for insignificant contribution of electrical interactions on solute rotation, in contrast to the substantial role played by the translational dielectric friction in the context of ionic mobility. It is also discussed on how the present theory can be suitably extended to study the rotation of a realistic solute in media other than dipolar solvents.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.     

Free energy of solvation from molecular dynamics simulation applying Voronoi-Delaunay triangulation to the cavity creation

The free energy of solvation for a large number of representative solutes in various solvents has been calculated from the polarizable continuum model coupled to molecular dynamics computer simulation. A new algorithm based on the Voronoi-Delaunay triangulation of atom-atom contact points between the solute and the solvent molecules is presented for the estimation of the solvent-accessible surface surrounding the solute. The volume of the inscribed cavity is used to rescale the cavitational contribution to the solvation free energy for each atom of the solute atom within scaled particle theory. The computation of the electrostatic free energy of solvation is performed using the Voronoi-Delaunay surface around the solute as the boundary for the polarizable continuum model. Additional short-range contributions to the solvation free energy are included directly from the solute-solvent force field for the van der Waals-type interactions. Calculated solvation free energies for neutral molecules dissolved in benzene, water, CCl4, and octanol are compared with experimental data. We found an excellent correlation between the experimental and computed free energies of solvation for all the solvents. In addition, the employed algorithm for the cavity creation by Voronoi-Delaunay triangulation is compared with the GEPOL algorithm and is shown to predict more accurate free energies of solvation, especially in solvents composed by molecules with nonspherical molecular shapes.     

Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies

Our previously developed approaches for integrating quantum mechanical molecular orbital methods with microscopic solvent models are refined and examined. These approaches consider the nonlinear solute–solvent coupling in a self-consistent way by incorporating the potential from the solvent dipoles in the solute Hamiltonian, while considering the polarization of the solvent by the potential from the solute charges. The solvent models used include the simplified Langevin Dipoles (LD) model and the much more expensive surface constrained All Atom Solvent (SCAAS) model, which is combined with a free energy pertubation (FEP) approach. Both methods are effectively integrated with the quantum mechanical AMPAC package and can be easily combined with other quantum mechanical programs. The advantages of the present approaches and their earlier versions over macroscopic reaction field models and supermolecular approaches are considered. A LD/MNDO study of solvated organic ions demonstrates that this model can yield reliable solvation energies, provided the quantum mechanical charges are scaled to have similar magnitudes to those obtained by high level ab initio methods. The incorporation of a field-dependent hydrophobic term in the LD free energy makes the present approach capable of evaluating the free energy of transfer of polar molecules from non polar solvents to aqueous solutions. The reliability of the LD approach is examined not only by evaluating a rather standard set of solvation energies of organic ions and polar molecules, but also by considering the stringent test case of sterically hindered hydrophobic ions. In this case, we compare the LD/MNDO solvation energies to the more rigorous FEP/SCAAS/MNDO solvation energies. Both methods are found to give similar results even in this challenging test case. The FEP/SCAAS/AMPAC method is incorporated into the current version of the program ENZYMIX. This option allows one to study chemical reactions in enzymes and in solutions using the MNDO and AM1 approximations. A special procedure that uses the EVB method as a reference potential for SCF MO calculations should help in improving the reliability of such studies.     

Dipole solvation in dielectrics

This paper presents an exact solution for the free energy of linear solvation of a dipolar solute in an arbitrary dielectric material with a microscopic spectrum of polarization fluctuations. The solution is given in terms of wave vector-dependent longitudinal and transverse structure factors of the polarization fluctuations in the pure dielectric. Good agreement with computer simulations of dipole solvation in dipolar and dipolar--quadrupolar liquids is achieved.     

Coupling quantum Monte Carlo to a nonlinear polarizable continuum model for spherical solutes

Starting from the nonlinear dielectric response model of Sandberg and Edholm, we derive an analytical expression of the polarization contribution to the solvation free energy in terms of the electronic density of the solute and the dielectric properties of the solvent. The solvent inhomogeneity is taken into account with the use of a smooth switching function whose spacial variation is established on the basis of how the solvent is arranged around the solute. An explicit form of a local potential representing the solvent effect on the solute is thus obtained by functional analysis. This effective potential can be combined with density functional or quantum chemical methods for the quantum mechanical treatment of the solute. Here, we use quantum Monte Carlo techniques for the solute and apply the method to the hydration of atomic ions finding very good agreement with experimental data.     

Dynamic stokes shift and excitation wavelength dependent fluorescence of dipolar molecules in room temperature ionic liquids

Room temperature ionic liquids (RTILs) are viscous media consisting entirely of ions. Because of the complex nature of various interactions in these media, the solvent properties of the RTILs are very little understood. Since the fluorescence response of molecules comprising conjugated electron donor and acceptor groups, referred to as dipolar molecules, is one of the most frequently exploited sources of information on complex media, whose properties are largely unknown, it is possible to obtain insight into the structure and dynamics of the RTILs by studying the fluorescence behavior of dipolar solutes in these complex media. The most commonly exploited utility of a fluorescent dipolar system is in the estimation of the polarity of the media from its steady state fluorescence response. While several dipolar systems do provide estimates of the polarity of various RTILs, there can be circumstances when the steady state emission frequency of a dipolar system may not truly reflect the equilibrium solvation energy and, hence, the polarity of the medium. The fluorescence response of a dipolar system can be dependent on the excitation wavelength, an observation not commonly encountered in conventional solvents of similar polarities. On the other hand, the time-resolved fluorescence behavior of a dipolar solute in polar medium is one of the primary sources of information on the time-scale of reorganization of the solvent molecules around the photoexcited species. As the RTILs are sufficiently polar media, the time-dependent fluorescence data of the dipolar systems provide insight into the dynamics and mechanism of solvation in these media, which differ considerably from the conventional solvents. These aspects have been discussed taking into consideration the inherent absorption and fluorescence behavior of the imidazolium ionic liquids.     

Simulations of solvation free energies and solubilities in supercritical solvents

Computer simulations are used to study solvation free energies and solubilities in supercritical solvents. Solvation free energies are calculated using the particle insertion method. The equilibrium solvent configurations required for these calculations are based on molecular dynamics simulations employing model solvent potentials previously tuned to reproduce liquid-vapor coexistence properties of the fluids Xe, C(2)H(6), CO(2), and CHF(3). Solutes are represented by all-atom potentials based on ab initio calculations and the OPLS-AA parameter set. Without any tuning of the intermolecular potentials, such calculations are found to reproduce the solvation free energies of a variety of typical solid solutes with an average accuracy of +/-2 kJmol. Further calculations on simple model solutes are also used to explore general aspects of solvation free energies in supercritical solvents. Comparisons of solutes in Lennard-Jones and hard-sphere representations of Xe show that solvation free energies and thus solubilities are not significantly influenced by solvent density fluctuations near the critical point. The solvation enthalpy and entropy do couple to these fluctuations and diverge similarly to solute partial molar volumes. Solvation free energies are also found to be little affected by the local density augmentation characteristic of the compressible regime. In contrast to solute-solvent interaction energies, which often provide a direct measure of local solvent densities, solvation free energies are remarkably insensitive to the presence of local density augmentation.     

Shape effect on nanoparticle solvation: a comparison of morphometric thermodynamics and microscopic theories

Conventional wisdom for controlling the nanoparticle size and shape during synthesis is that particle growth favors the direction of a facet with the highest surface energy. However, the particle solvation free energy, which dictates the particle stability and growth, depends not only on the surface area and surface free energy but also on other geometric measures such as the solvent excluded volume and the surface curvature and their affiliated thermodynamic properties. In this work, we study the geometrical effects on the solvation free energies of nonspherical nanoparticles using morphometric thermodynamics and density functional theories. For idealized systems that account for only molecular excluded-volume interactions, morphometric thermodynamics yields a reliable solvation free energy when the particle size is significantly larger than the solvent correlation length. However, noticeable deviations can be identified in comparison to the microscopic theories for predicting the solvation free energies of small nanoparticles. This conclusion also holds for predicting the potential of mean force underlying the colloidal "key-and-lock" interactions. Complementary to the microscopic theories, morphometric thermodynamics requires negligible computational cost, therefore making it very appealing for a broad range of practical applications.     

Spectral Investigations of Preferential Solvation and Solute–Solvent Interactions of Free Base and Protonated 5,10,15,20-Tetrakis(4-trimethyl-ammonio-phenyl)-porphine Tetratosylate in Aqueous Organic Mixed Solvents

The solvatochromic properties of the free base and the protonated 5,10,15,20-tetrakis(4-trimethyl-ammonio-phenyl)-porphine tetratosylate (TTMAPP) were studied in pure water, methanol, ethanol, 2-propanol, and their corresponding aqueous mixtures. The correlation of the empirical solvent polarity scale (E _T) values of TTMAPP with composition of the solvents were analyzed by the solvent exchange model of Bosch and Roses to clarify the preferential solvation of the probe dyes in the binary mixed solvents. The solvation shell composition effects in preferential solvation of the solute dyes were investigated in terms of both solvent–solvent and solute–solvent interactions and also the local mole fraction of each solvent composition was calculated in the cybotactic region of the probe. The effective mole fraction variation may provide significant physicochemical insights in the microscopic and molecular level of interactions between TTMAPP species and the solvent components and, therefore, can be used to interpret the solvent effect on kinetics and thermodynamics of TTMAPP.     

Treating entropy and conformational changes in implicit solvent simulations of small molecules

Implicit solvent models are increasingly popular for estimating aqueous solvation (hydration) free energies in molecular simulations and other applications. In many cases, parameters for these models are derived to reproduce experimental values for small molecule hydration free energies. Often, these hydration free energies are computed for a single solute conformation, neglecting solute conformational changes upon solvation. Here, we incorporate these effects using alchemical free energy methods. We find significant errors when hydration free energies are estimated using only a single solute conformation, even for relatively small, simple, rigid solutes. For example, we find conformational entropy (TDeltaS) changes of up to 2.3 kcal/mol upon hydration. Interestingly, these changes in conformational entropy correlate poorly (R2 = 0.03) with the number of rotatable bonds. The present study illustrates that implicit solvent modeling can be improved by eliminating the approximation that solutes are rigid.     

The influence of solvation and finite temperatures on theWittig reaction: A theoretical study

The effects of the solvent and finite temperature (entropy) on the Wittig reaction are studied by using density functional theory in combination with molecular dynamics and a continuum solvation model. Standard gas-phase zero-temperature calculations are found to give similar results to previous studies. Gas-phase dynamics simulations allow the free energy profile of the reaction to be calculated through thermodynamic integration. The free energy profile is found to have a significant entropic barrier to the addition step of the reaction where only a small barrier was present in the potential energy curve. The introduction of the solvent dimethyl sulfoxide causes a change in the structure of the intermediate from the oxaphosphetane structure to the dipolar betaine structure. The overall reaction energy is changed only slightly. When the effects of both entropy and the solvent are included a significant entropic barrier to the addition reaction is obtained and the predicted intermediate again has the betaine structure.