Solvent effects on infrared spectra of methyl 4-hydroxybenzoate in pure organic solvents and in ethanol/CCl4 binary solvents

Infrared spectroscopy studies of methyl 4-hydroxybenzoate (MHB) in 17 different organic solvents and in ethanol/CCl4 binary solvent were undertaken to investigate the solvent-solute interactions. The frequencies of carbonyl stretching vibration nu(C=O) of MHB in single solvents were correlated with the solvent acceptor number (AN) and the linear solvation energy relationships (LSER). The assignments of the two bands of nu(C=O) of MHB in alcohols and the single one of that in non-alcoholic solvents were discussed. The shifts of nu(C=O) of MHB in ethanol/CCl4 binary solvents showed that several kinds of solute-solvent hydrogen bonding interactions coexisted in the mixture solvents, with a change in the mole fraction of ethanol in the binary solvents.     

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.     

Solvent effects on infrared spectra of 2-acetylthiophene in organic solvents

Infrared spectroscopy studies of 2-acetylthiophene (ACTH) in 18 different organic solvents, both polar and non-polar, were undertaken to investigate the solvent-solute interactions. The frequencies of carbonyl stretching vibration upsilon(C=O) of ACTH were correlated with the properties such as the solvent acceptor number (AN) and the linear solvation energy relationships (LSER). The solvent-induced stretching vibration frequency shifts showed a better correlation with the LSER than the AN. A six-membered ring-like hydrogen bonding structure was presented and the solvent effects of ACTH in alcohol solvents were investigated in detail.     

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.     

Modeling of tracer diffusion in liquids when solute–solvent interactions are present

The linear solvation energy relationship (LSER) model is employed to correlate the tracer diffusion coefficients of 550 binary systems at 298.15 K. Among the selected solutes and solvents there exist apolar, polar and hydrogen-bonding substances that can interact with themselves (solvent polymerization) or with the other compound (solute–solvent complexes formation). The results of the proposed formulas are compared with those of other predictive equations.     

Solvatochromic study of 1,2-dihydroxyanthraquinone in neat and binary solvent mixtures

Preferential solvation of a solvatochromic probe has been studied in binary mixtures comprising of a non-protic and a protic solvent. The non-protic solvents employed are carbon tetrachloride (CCl(4)), acetonitrile (AcN) and N,N-dimethyl formamide (DMF) and the protic solvents are methanol (MeOH) and ethanol (EtOH). The probe molecule exhibits different spectroscopic characteristics depending upon the properties of the solubilizing media. The observed spectral features provide an indication of the microenvironment immediately surrounding the probe. Solvatochromic shifts of the ground and excited states of the probe were analysed by monitoring the charge transfer absorption band and the fluorescence emission spectra in terms of the solute-solvent and solvent-solvent interactions. Fluorescence emission spectra show the dual emission due to excited state proton transfer nature of the probe molecule. The effect of solvent and the excitation energy on dual emission are also studied. The observed magnitude of the Stokes shift in the above solvents has been used to deduce experimentally the dipole moment ratio of the probe molecule for the excited state to the ground state. The dipole moment of excited state is higher than the ground state.     

Solvent reorganization entropy of electron transfer in polar solvents

We report the results of molecular dynamics simulations of the solvent reorganization energy of intramolecular electron transfer in a charge-transfer molecule dissolved in water and acetonitrile at varying temperatures. The simulations confirm the prediction of microscopic solvation theories of a positive reorganization entropy in polar solvents. The results of simulations are analyzed in terms of the splitting of the reorganization entropy into the contributions from the solute-solvent interaction and from the alteration of the solvent structure induced by the solute. These two contributions mutually cancel each other, resulting in the reorganization entropy amounting to only a fraction of each component.     

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.     

Preferential solvation and solvation shell composition of free base and protonated 5, 10, 15, 20-tetrakis(4-sulfonatophenyl)porphyrin in aqueous organic mixed solvents

The solvatochromic properties of the free base and the protonated 5, 10, 15, 20-tetrakis(4-sulfonatophenyl)porphyrin (TPPS) were studied in pure water, methanol, ethanol (protic solvents), dimethylsulfoxide, DMSO, (non-protic solvent), and their corresponding aqueous-organic binary mixed solvents. The correlation of the empirical solvent polarity scale (E(T)) values of TPPS with composition of the solvents was 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 and the synergistic 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 cybotactic region of the probe. The effective mole fraction variation may provide significant physico-chemical insights in the microscopic and molecular level of interactions between TPPS species and the solvent components and therefore, can be used to interpret the solvent effect on kinetics and thermodynamics of TPPS. The obtained results from the preferential solvation and solvent-solvent interactions have been successfully applied to explain the variation of equilibrium behavior of protonation of TPPS occurring in aqueous organic mixed solvents of methanol, ethanol and DMSO.     

Solvent Effects on the Structure-Property Relationship of Some Potentially Pharmacologically Active 3-(4-Substituted Benzyl)-5-Ethyl-5-Phenyl- and 3-(4-Substituted Benzyl)-5,5-Diphenylhydantoins

Two series of 3-(4-substituted benzyl)-5-ethyl-5-phenyl- and 3-(4-substituted benzyl)-5,5-diphenylhydantoins were synthesized and their UV absorption spectra were recorded in the region 200–400 nm in selected solvents of different polarity. The effects of solvent dipolarity/polarizability and solvent/solute hydrogen bonding interactions on the spectral shifts were analyzed by means of the linear solvation energy relationship (LSER) methodology of Kamlet and Taft. The quantitative relationships between hydrogen bonding interactions and the lipophilicity and blood-brain permeation of the studied compounds were discussed. Satisfactory linear dependences were obtained for moderate electron-donating and electron-withdrawing substituents at the benzyl moiety, while the strong electron-withdrawing substituent (NO₂) significantly modifies the solvation characteristics of the molecule. The paper clearly demonstrates how the solvatochromic comparison method may be applied to estimate the contributions of various modes of solvation to the pharmaceutically relevant properties of these newly synthetized hydantoin derivatives.     

The importance of excluded solvent volume effects in computing hydration free energies

Continuum dielectric methods such as the Born equation have been widely used to compute the electrostatic component of the solvation free energy, DeltaG(solv)(elec), because they do not need to include solvent molecules explicitly and are thus far less costly compared to molecular simulations. All of these methods can be derived from Gauss Law of Maxwell's equations, which yields an analytical solution for the solvation free energy, DeltaG(Born), when the solute is spherical. However, in Maxwell's equations, the solvent is assumed to be a structureless continuum, whereas in reality, the near-solute solvent molecules are highly structured unlike far-solute bulk solvent. Since we have recently reformulated Gauss Law of Maxwell's equations to incorporate the near-solute solvent structure by considering excluded solvent volume effects, we have used it in this work to derive an analytical solution for the hydration free energy of an ion. In contrast to continuum solvent models, which assume that the normalized induced solvent electric dipole density P(n) is constant, P(n) mimics that observed from simulations. The analytical formula for the ionic hydration free energy shows that the Born radius, which has been used as an adjustable parameter to fit experimental hydration free energies, is no longer ill defined but is related to the radius and polarizability of the water molecule, the hydration number, and the first peak position of the solute-solvent radial distribution function. The resulting DeltaG(solv)(elec) values are shown to be close to the respective experimental numbers.     

Solvent density inhomogeneities and solvation free energies in supercritical diatomic fluids: a density functional approach

Density functional theory is used to explore the solvation properties of a spherical solute immersed in a supercritical diatomic fluid. The solute is modeled as a hard core Yukawa particle surrounded by a diatomic Lennard-Jones fluid represented by two fused tangent spheres using an interaction site approximation. The authors' approach is particularly suitable for thoroughly exploring the effect of different interaction parameters, such as solute-solvent interaction strength and range, solvent-solvent long-range interactions, and particle size, on the local solvent structure and the solvation free energy under supercritical conditions. Their results indicate that the behavior of the local coordination number in homonuclear diatomic fluids follows trends similar to those reported in previous studies for monatomic fluids. The local density augmentation is particularly sensitive to changes in solute size and is affected to a lesser degree by variations in the solute-solvent interaction strength and range. The associated solvation free energies exhibit a nonmonotonous behavior as a function of density for systems with weak solute-solvent interactions. The authors' results suggest that solute-solvent interaction anisotropies have a major influence on the nature and extent of local solvent density inhomogeneities and on the value of the solvation free energies in supercritical solutions of heteronuclear molecules.