Solvent effects on 17O nuclear magnetic shielding: N-methylformamide in polar and apolar solutions

The solvent effect on ¹⁷O isotropic shielding was computed by different methods: the polarizable continuum model and a mixed approach, including a few real solvent molecules treated as the solute. The experimental data show that the behavior of protic and aprotic solvents is markedly different: we found that the continuum approach describes well the observed shielding at various dielectric constants for aprotic solvents, while the mixed procedure is needed when hydrogen bonds to the magnetic centre are present.Contribution to the Jacopo Tomasi Honorary Issue     

Interactions of polarizable media in water: a molecular dynamics approach

We investigate the interactions of polarizable solutes in water as a function of the solute permittivity. A generic and computationally efficient simulation methodology for the investigation of systems involving dielectric discontinuities is introduced. We report results for interactions between two polarizable cylindrical solutes of nanometer dimensions, which demonstrate that the interactions between the solutes strongly depend on the solute permittivity epsilon. For low permittivity, epsilon approximately 1-2, the interactions are dominated by surface tension forces whose origin lies in the formation of a vapor cavity between the two hydrophobic solutes. This effect leads to a drying transition, where the intersolute force changes discontinuously at a specific solute-solute separation. We find that a moderate permittivity, epsilon approximately 20, enhances the solvation of the polarizable objects inhibiting this drying transition. In the limit of moderately high permittivity, the interactions are dominated by solvation forces. These forces are much larger than those calculated using macroscopic models of dielectrics, which consider water as a continuum dielectric medium. Our results emphasize the importance of including the solvent explicitly to investigate dielectric discontinuities and interactions between polarizable media in water.     

Effects of ion solvation on phase equilibrium and interfacial tension of liquid mixtures

We study the bulk thermodynamics and interfacial properties of electrolyte solution mixtures by accounting for electrostatic interaction, ion solvation, and inhomogeneity in the dielectric medium in the mean-field framework. Difference in the solvation energy between the cations and anions is shown to give rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The ion solvation affects the phase equilibrium of the solvent mixture, depending on the dielectric constants of the solvents, reflecting the competition between the solvation energy and translation entropy of the ions. Miscibility is decreased if both solvents have low dielectric constants and is enhanced if both solvents have high dielectric constant. At the mean-field level, the ion distribution near the interface is determined by two competing effects: accumulation in the electrostatic double layer and depletion in a diffuse interface. The interfacial tension shows a nonmonotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations and decreases approximately as the square root of the salt concentration for dilute solutions, reaching a minimum near 1 mM. We also find that, for a fixed cation type, the interfacial tension decreases as the size of anion increases. These results offer qualitative explanations within one unified framework for the long-known concentration and ion size effects on the interfacial tension of electrolyte solutions.     

Tautomeric conjugate acids of 2-aminopyrroles: effect of substituents, solvation and cosolute

The tautomeric preferences of the conjugated acids of 2-aminopyrrole derivatives have been examined both in the gas phase and in aqueous solution by using a combination of quantum mechanical, self-consistent reaction field and Monte Carlo–free-energy perturbation methods. The results show that the nature of substituents, the solvent and the presence of cosolute are relevant factors in modulating the relative stability between the tautomeric conjugate acids protonated at the heterocyclic ring and at the exocyclic amino nitrogen. Thus, attachment of electron-withdrawing groups to the ring, solvation in polar solvents, and the presence of negatively charged cosolutes tend to favor protonation at the exocyclic amino nitrogen. Nevertheless, none of these factors alone suffice to change the tautomeric preference for the ring-protonated forms. The results point out that the concerted occurrence of the three factors is necessary to shift the tautomeric preference towards the conjugated species protonated at the exocyclic nitrogen.Contribution to the Jacopo Tomasi Honorary Issue     

Solute–solvent interactions of acid–1,4-dioxane mixtures—By dielectric, FTIR, UV–vis and 13C NMR spectrometric methods

Results of the dielectric studies carried out on the binary mixture of n-butyric and caprylic acids with 1,4-dioxane over the entire composition range and at temperatures 303 K, 308 K, 313 K and 318 K, and FTIR, UV–vis and 13C NMR spectral studies are presented in this paper. The excess permittivity and excess free energy were fitted with the Redlich–Kister polynomial. The variation of Kirkwood correlation factors, excess permittivity and excess free energy of mixing with the concentration and temperature has been investigated in view of understanding the ordering of dipoles of solute and solvent molecules. The FTIR, UV–vis and 13C NMR spectral analysis reveals the formation of complex between solute and solvent molecules. The parallel alignment of electric dipoles of the complex predicted by dielectric studies is well supported by UV–vis spectral analysis. The structure of the complex molecule present in the clusters has been deduced.     

The static electric permittivity of solutions of organic molecules in dipolar solvents

The permittivity of solutions of pyrazine and of 1,4-diazabicyclo-[2,2,2]-octane (TED) in various dipolar solvents has been measured at 10 MHz by means of a sensitive bridge method. The solvents used were water, methanol, formamide, N-methylformamide, N,N-dimethylformamide, and dimethylsulfoxide. The results of the measurements are compared with those of theoretical mixture formulas. It is shown that — depending on structural properties — the permittivity attributed to the solvation shells around the organic solute molecules is changed with respect to the pure solvent value.     

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.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

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.     

Simulation of volume polarization for the influence of solvation on chemical shielding

The effect of bulk dielectric solvation on chemical shielding at nitrogen in CH₃CN is studied with reaction field theory. A previous work has demonstrated the strong influence on this property from volume polarization, which describes that part of the reaction field arising from solute charge density penetrating outside its cavity. The essentially exact treatment of volume polarization used in that work is computationally demanding, and a more facile method for simulation of the volume polarization has recently been proposed. It is found in the present work that this simulation of the volume polarization yields results in excellent agreement with the essentially exact treatment of the strong volume polarization effects on nitrogen shielding in CH₃CN.Contribution to the Jacopo Tomasi Honorary Issue     

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.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

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.     

Absorption spectral band width of charge transfer transition of E(T)(30) dye in homogeneous and heterogeneous media

Solvation characteristics in homogeneous and heterogeneous media have been probed by monitoring the band width of ICT band of 2,6-di-phenyl-4(2,4,6-triphenyl-1-pyridino) phenolate, the indicator solute for E(T)(30) scale, in pure, mixed binary solvents and aqueous micellar solution. Non-ideal solvation behaviour is observed in all the binary solvent mixtures. Index of preferential solvation has been calculated as a function of solvent composition. Study in micellar media indicates that the dye is located at the micelle-water interface. The effects of variation of micelle concentration, temperature and electrolyte concentration have also been studies.     

Study on Preferential Solvation of Water by Electrochemical Method

The preferential solvation of water plays an important role in ferrocene research which is a subject of current interest. Voltammetric investigations were carried out for Au electrode in acetonitrile/water, showing preferential solvation of water. In our work, the preferential solvation of water in acetonitrile/water was studied by electrochemical methods including cyclic volitammetry, electrochemical impedance spectra and double‐step chronoamperometry. Ferrocenemethanol (FcCH₂OH) molecules as a solute spontaneously adsorb on the electrode surface in anhydrous acetonitrile, resulting from acetonitrile molecules tend to form an acetonitrile solvent layer on the surface of the electrode and acetonitrile solvent layer has a lower energy barrier than the aqueous solvent layer, which has been obtained by modeling solvation. The solvent strongly influences electrochemical behavior of solute. Once there is an amount of water in acetonitrile solvent, FcCH₂OH that adsorbed on the electrode surface desorb. This is because water preferentially solvate with FcCH₂OH in term of intermolecular forces between solvent and solute. Moreover, hydrogen bond between water molecules and FcCH₂OH molecules is stronger than dipole‐dipole interaction between acetonitrile molecules and FcCH₂OH molecules in solvation effect. Through electrochemical behavior of FcCH₂OH changing, preferential solvation of water is analyzed by electrochemical methods.     

Compensation Between the Entropic and Enthalpic Changes Arising from Changes in Solvent–Solvent Interactions in Mixed Solvents

A theorem presented by Professor Ben-Naim (J Phys Chem 82:874–885, 1978) states that the standard state enthalpy and entropy changes arising from changes in the solvent structure that are induced by solvation of a solute cancel exactly in the standard state Gibbs energy. In this paper this is explored by consideration of the thermodynamics of transfer of electrolytes in mixed solvents, using previously developed models of the solvation process. Two cases are considered. One is random solvation, where curvatures in plots of the transfer enthalpies and entropies, which arise from changes in solvent–solvent interactions, exactly compensate in the transfer Gibbs (free) energies, which are sensibly linear with solvent composition. The second type of system are those with strong preferential solvation where it is found that the transfer Gibbs energies can be accounted for quantitatively in terms of changes in the solute–solvent interactions, with no contribution from changes in solvent–solvent interactions. The results are entirely consistent with the Ben-Naim theorem.     

Electronic spectroscopic study of solvation of a ketocyanine dye in ternary solvent mixtures

Solvation characteristics in ternary solvent mixtures have been studied by monitoring the solvent-sensitive electronic absorption band of a ketocyanine dye in two ternary solvent mixtures, water + ethanol + benzene and water + ethanol + cyclohexane, in which one of the pairs are partially miscible. Investigations have been done in a completely miscible region including the binodal curve. The maximum energy of absorption (E) of the solute in a ternary solvent mixture differs significantly from the mole fraction average of the E-values in the component solvents. Results in the corresponding binary solvent mixtures also show a deviation of the E-value from the mole fraction averaged E-values, indicating preferential solvation by a component solvent. The results in ternary solvent mixture have been explained in terms of a realistic model of solvation using the results on binary solvation.     

Use of a solvatochromic probe for study of solvation in ternary solvent mixture

Solvation characteristics of 2,6-diphenyl-4-(2,4,6-triphenyl-1-pyridino)phenolate in completely miscible ternary solvent mixtures (viz., methanol + acetone + water, methanol + acetone + benzene, and methanol + chloroform + benzene) have been studied by using an electronic spectroscopic procedure. The transition energy (E) corresponding to the charge-transfer band maximum of the solute in a ternary solvent mixture differs significantly from the average E-values in the component solvents weighted by the mole fraction of the solvents. A two-phase model of solvation has been invoked to explain the results. The excess or deficit of solvent components in the local region of the solute molecule over that in the bulk has been estimated using the knowledge of solvation in binary solvent mixtures.