On the coupling between molecular diffusion and solvation shell exchange

The connection between diffusion and solvent exchanges between first and second solvation shells is studied by means of molecular dynamics simulations and analytic calculations, with detailed illustrations for water exchange for the Li(+) and Na(+) ions, and for liquid argon. First, two methods are proposed which allow, by means of simulation, to extract the quantitative speed-up in diffusion induced by the exchange events. Second, it is shown by simple kinematic considerations that the instantaneous velocity of the solute conditions to a considerable extent the character of the exchanges. Analytic formulas are derived which quantitatively estimate this effect, and which are of general applicability to molecular diffusion in any thermal fluid. Despite the simplicity of the kinematic considerations, they are shown to well describe many aspects of solvent exchange/diffusion coupling features for nontrivial systems.     

Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.     

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.     

Contributions of short-range and excluded-volume interactions to unperturbed polymer chain dimensions

A Monte Carlo study is made of the mean-square radius of gyration for the freely rotating chain with such fictitious excluded-volume interactions that the Lennard-Jones 6-12 potentials at the Theta temperature act only between the fourth-through (3+Delta)th-neighbor beads (Delta > or = 1) along the chain. The behavior of the asymptotic value (/n)infinity of the ratio /n as a function of the number n of bonds in the chain in the limit of n --> infinity is examined as a function of Delta. It is shown that the approach of (/n)infinity to its value for the real unperturbed chain with Delta = infinity is so slow that the interactions between even up to about 100th-neighbor beads should be taken into account in order to reproduce nearly its dimension. The result implies that the unperturbed polymer chain dimension as experimentally observed at the Theta temperature depends not only on short-range interactions but also to a considerable extent on the long-range excluded-volume interactions, and that the asymptotic value Cinfinity of the characteristic ratio Cn for the rotational isomeric state model in the limit of n --> infinity, which is determined only by the very local conformational energy, cannot be directly compared with the corresponding experimental value.     

Modeling the effect of solvation on solute retention in reversed-phase liquid chromatography

The retention of a homologous series of alkylbenzenes was determined on octyl and octadecyl reversed-phase columns in several polar organic liquids. Free energies of transfer were calculated by the SM5.0R classical solvation model for each organic liquid tested and for several alkanes. The relationships between the measured retention factors and the calculated free energies of transfer were then investigated. Although the natural logarithms of the retention factor and the calculated free energies of transfer were linearly correlated, the obtained free energies of transfer of the solutes did not completely explain the retention behavior of the solutes. Nonetheless, even in these pure organic liquids, the energetics of RPLC retention behaved very similarly to those of partitioning.     

Electrostatic interactions of a neutral dipolar solute with a fused salt: a new model for solvation in ionic liquids

Ionic liquids represent a novel and poorly understood class of solvents, and one challenge in understanding these systems is how one should view the electrostatic character of solute-solvent interactions. The highly structured nature of a fused salt makes a dielectric continuum approximation difficult to implement, and there is no obvious connection between the structure of an individual ion and the polarization character of the medium. We address this problem by making the ansatz that rather than polarizing the medium, the solute may be viewed as intercalating in the charge distribution of the neat liquid such that the solvent screens the electric field of the solute. This approach allows derivation of an analytical expression for the distribution of solvent charge about the solute, and this distribution is found to be a close match to simulation data. The theory also predicts that the electrostatic character of solute-solvent interactions should be determined primarily by the number density of solvent ions, a prediction proven correct by analysis of existing experimental data. The approach represents a new model for the interpretation of solvation phenomena in ionic liquids.     

Influence of the solvation process on solute adsorption in reversed phase liquid chromatography

The adsorption isotherms of phenol were acquired by frontal analysis on six different reversed phase adsorbents from five different organic solvent solutions. The end-capped octadecyl columns only differed in the bonding density of the C(18) ligands. The inverse method was used to confirm the estimated isotherm parameters derived from the frontal experiments. The effect of the bonding density of the end-capped octadecyl bonded phase on the adsorption properties of phenol from different mobile phase compositions was investigated. The adsorption behavior of phenol has changed from Langmuir type to BET type with the change of the organic modifier and the bonding density of the adsorbent.     

Effects of solute structure on local solvation and solvent interactions: results from UV/vis spectroscopy and molecular dynamics simulations

Solvation of heterocyclic amines in CO(2)-expanded methanol (MeOH) has been explored with UV/vis spectroscopy and molecular dynamics (MD) simulations. A synergistic study of experiments and simulations allows exploration of solute and solvent effects on solvation and the molecular interactions that affect absorption. MeOH-nitrogen hydrogen bonds hinder the n-pi* transition; however, CO(2) addition causes a blue shift relative to MeOH because of Lewis acid/base interactions with nitrogen. Effects of solute structure are considered, and very different absorption spectra are obtained as nitrogen positions change. MD simulations provide detailed solvent clustering behavior around the solute molecules and show that the local solvent environment and ultimately the spectra are sensitive to the solute structure. This work demonstrates the importance of atomic-level information in determining the structure-property relationships between solute structure, local salvation, and solvatochromism.     

On the role of polarizability in chemical-biological interactions

This report considers the importance of electronic effects in their role in the QSAR of chemical-biological interactions. The problem of accounting for polarizability effects in ligand-substrate interactions is discussed in terms of molecular polarizability (MR) and NVE (number of valence electrons) using additive values for valence electrons. The two approaches give essentially the same result in examples of frog nerve toxicity and examples of nerve toxicity with rabbits and cockroaches. The point is made that no matter how one approaches QSAR, electronic interactions must be considered if we are to begin to develop a science of chemical-biological interactions.     

Effect of nonpolar solvents on the solute rotation and solvation dynamics in an imidazolium ionic liquid

Recognizing the potential of the mixed solvent systems comprising ionic liquid as one of the constituents in real applications, the steady-state and time-resolved fluorescence behavior of C153 has been studied in neat 1-butyl-3-methylimidazolium hexafluorophosphate and its mixtures with nonpolar solvents, namely, toluene and 1,4-dioxane. No significant effect of the cosolvent on the steady-state absorption or fluorescence spectra of C153 in ionic liquid has been observed. Time-resolved fluorescence anisotropy measurements show a decrease of the rotational correlation time of C153 with gradual addition of the cosolvent. Solvation dynamics in ionic liquid-cosolvent mixtures is found to be biphasic, and a decrease of the average solvation time is observed with increasing amount of the cosolvent in solution. The time-zero spectrum of C153 is found to shift toward higher energy with gradual addition of the nonpolar solvent, suggesting that the probe molecule experiences a more nonpolar environment at the early stage of the dynamics in mixed solvents. The blue shift of the time-zero spectrum caused by the addition of the nonpolar solvent results in a larger Stokes shift of the time-dependent spectra due to solvent relaxation in mixed solvents. A comparison of the time-dependent spectral data of the ionic liquid-toluene and ionic liquid-dioxane systems shows that, while a small amount of toluene can significantly affect the dynamics, comparatively, a larger amount of dioxane is required to bring about the same effect. This is explained in terms of favorable interactions between toluene and the imidazolium ring system leading to a more effective solubilization of toluene in the cybotactic region of the probe.     

The role of electrostatic interactions in protease surface diffusion and the consequence for interfacial biocatalysis

This study examines the influence of electrostatic interactions on enzyme surface diffusion and the contribution of diffusion to interfacial biocatalysis. Surface diffusion, adsorption, and reaction were investigated on an immobilized bovine serum albumin (BSA) multilayer substrate over a range of solution ionic strength values. Interfacial charge of the enzyme and substrate surface was maintained by performing the measurements at a fixed pH; therefore, electrostatic interactions were manipulated by changing the ionic strength. The interfacial processes were investigated using a combination of techniques: fluorescence recovery after photobleaching, surface plasmon resonance, and surface plasmon fluorescence spectroscopy. We used an enzyme charge ladder with a net charge ranging from -2 to +4 with respect to the parent to systematically probe the contribution of electrostatics in interfacial enzyme biocatalysis on a charged substrate. The correlation between reaction rate and adsorption was determined for each charge variant within the ladder, each of which displayed a maximum rate at an intermediate surface concentration. Both the maximum reaction rate and adsorption value at which this maximum rate occurs increased in magnitude for the more positive variants. In addition, the specific enzyme activity increased as the level of adsorption decreased, and for the lowest adsorption values, the specific enzyme activity was enhanced compared to the trend at higher surface concentrations. At a fixed level of adsorption, the specific enzyme activity increased with positive enzyme charge; however, this effect offers diminishing returns as the enzyme becomes more highly charged. We examined the effect of electrostatic interactions on surface diffusion. As the binding affinity was reduced by increasing the solution ionic strength, thus weakening electrostatic interaction, the rate of surface diffusion increased considerably. The enhancement in specific activity achieved at the lowest adsorption values is explained by the substantial rise in surface diffusion at high ionic strength due to decreased interactions with the surface. Overall, knowledge of the electrostatic interactions can be used to control surface parameters such as surface concentration and surface diffusion, which intimately correlate with surface biocatalysis. We propose that the maximum reaction rate results from a balance between adsorption and surface diffusion. The above finding suggests enzyme engineering and process design strategies for improving interfacial biocatalysis in industrial, pharmaceutical, and food applications.     

Prediction of solvation free energies from computed properties of solute molecular surfaces

We have shown that the solvation energies of a group of 12 solutes in 7 different solvents can be presented analytically in terms of quantities computed at the density functional B3P86/6‐31+G** level for the isolated solute molecules. These quantities include the molecular surface area and several properties of the electrostatic potential on the surface, e.g., the most positive and negative values, the average deviation of the potential, the positive and negative portions of the surface, and their average potentials. Overall, the average absolute deviation of the predicted from the experimental solvation free energies is 0.25 kcal/mol; the poorest results are obtained for the solute butanone, for which the average absolute deviation is 0.63 kcal/mol. The forms of the relationships reflect the natures of the solute–solvent interactions; for the solvents with low dielectric constants, these are primarily global, involving extended portions of the molecular surfaces, whereas for the more polar solvents, site‐specific interactions also play key roles. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 643–647, 2000     

The role of solvation in redox processes

In studying the role of solvation in redox processes we consider the influence of geometrical modifications in the discharge of monosolvated anion CN^?-H₂O. We define as “efficient” those modifications whose energy cost is less than the lowering which is produced in the ionization potential. we conclude that in the above-mentioned species the redox process is enhanced by desolvation. We try to generalise the results for the “stable” anions and cations.     

Theoretical representation of solvation in biochemical systems: From discrete solute-solvent interactions to bulk solvation

The interaction between solutes and simple solvents in dilute solutions is analyzed in a systematic way. Different theoretical methods to describe specific solute-solvent interactions and general solvation effects are presented. Finally, the importance of solvent-induced changes in solute properties is discussed through the use of mixed quantum mechanics/classical mechanics strategies. © 1996 John Wiley & Sons, Inc.     

Using molecular simulation to predict solute solvation and partition coefficients in solvents of different polarity

A methodology is proposed for the prediction of the Gibbs energy of solvation (Δ(Solv)G) based on MD simulations. The methodology is then used to predict Δ(Solv)G of four solutes (namely propane, benzene, ethanol and acetone) in several solvents of different polarities (including n-hexane, n-hexadecane, ethylbenzene, 1-octanol, acetone and water) while testing the validity of the TraPPE force field parameters. Excellent agreement with experimental data is obtained, with average deviations of 0.2, 1.1, 0.8 and 1.2 kJ mol(-1), for the four solutes respectively. Subsequently, partition coefficients (log P) for forty different solute/solvent systems are predicted. The a priori knowledge of partition coefficient values is of high importance in chemical and pharmaceutical separation process design or as a measure of the increasingly important environmental fate. Here again, the agreement between experimental data and simulation predictions is excellent, with an absolute average deviation of 0.28 log P units. However, this deviation can be decreased down to 0.14 log P units, just by optimizing partial atomic charges of acetone in the water phase. Consequently, molecular simulation is proven to be a tool with strong physical basis able to predict log P with competitive accuracy when compared to the popular statistical methods with weak physical basis.     

Molecular complexes and solvation interactions in the reaction of quinone imines with thiols

This study aims to investigate the role of complexation between reagents and the role of solvation of reagents by solvents in the kinetics of chain reactions of quinone imines with thiols. The thermodynamic characteristics of the complexation of quinone imines with thiophenol in CCl₄, chlorobenzene, and ethanol, as well as of the complexation of quinone imines and thiophenol with these solvents were calculated by quantum chemical methods (DFT calculations at the PBE/cc-pVDZ level of theory) and in terms of the additive-multiplicative model. Both approaches give consistent results. The formation of molecular complexes in quinone imine–thiphenol systems is accompanied by a 10–30 kJ mol^–1 decrease in enthalpy and has only a slight effect on the reaction mechanism.     

Characterizing ionic liquids on the basis of multiple solvation interactions

Room-temperature ionic liquids (RTILs) are useful in many chemical applications. Recent publications have attempted to determine the polarity of RTILs using empirical solvent polarity scales. The results have indicated that most RTILs have similar polarities. Nevertheless, RTILs are capable of behaving quite differently when used as solvents in organic synthesis, matrixes in matrix-assisted laser desorption/ionization (MALDI) mass spectrometry, liquid-liquid extraction, and as stationary phases in gas chromatography. The work presented in this study uses a linear free energy approach to characterize 17 RTILs on the basis of their distinct multiple solvation interactions with probe solute molecules. This model provides data that can be used to help identify the interactions and properties that are important for specific chemical applications.