Acceleration of convergence in Monte Carlo simulations of aqueous solutions using the metropolis algorithm. Hydrophobic hydration of methane

Convergence problems encountered in the computer simulations of aqueous solutions are discussed. Solute–solvent radial distribution functions are shown to converge very poorly when the standard Metropolis Monte Carlo procedure is applied. To overcome this difficulty, several modifications are made in the Metropolis method. Optimum maximum step sizes are determined for simulations of liquid water. A scheme is employed for preferential sampling of both the solvent and the solute molecules. To test these modifications, a simulation is carried out for pure liquid water, treating a single water molecule as a “solute.” The convergence of the radial distribution functions is found to be accelerated significantly. A further test is made by simulating an aqueous solution of methane, consisting of one methane molecule (using the EPEN /2 potential for methane–water interactions) and 124 water molecules (using the MCY potential for water–water interactions). Again, the convergence of solute–solvent radial distribution functions is found to be accelerated. The computation of partial molar thermodynamic quantities, however, still suffers from convergence difficulties. This problem is discussed in detail. The EPEN /2 potential is found to yield structural and thermodynamic features of hydrophobic hydration that are consistent with available experimental and theoretical results for aqueous solutions of methane.     

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.     

Solvent Effects on the Structure-Activity Relationship of Phenytoin-like Anticonvulsant Drugs

The absorption spectra of nine compounds structurally related to phenytoin (5,5-diphenylhydantoin) were recorded in twelve solvents over the range of 200 to 400 nm. The effects of solvent dipolarity/polarizability and solvent/solute hydrogen bonding interactions were analyzed by means of the linear solvation energy relationship (LSER) concept proposed by Kamlet and Taft. The lipophilic activity of the investigated hydantoins was estimated by calculation of their log ₁₀ P values. The calculated values of log ₁₀ P were correlated with the ratio of the contributions of specific and non-specific solute/solvent interactions. The correlation equations were combined with the corresponding ED₅₀ values to generate new equations that demonstrate exact relationship between solute/solvent interactions and the structure-activity parameters.     

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.     

Solvatochromism of anthraquinone and symmetrical dihydroxy derivatives. Local interactions

The solvent effects on the electronic absorption spectra of 9,10-anthraquinone (AQ) and its symmetric dihydroxy derivatives namely 1,5-dihydroxyanthraquinone (1,5-DHAQ) and 2,6-dihydroxyanthraquinone (2,6-DHAQ) have been studied in pure solvents and some binary solvent mixtures. The frequencies of the absorption for AQ and 2,6-DHAQ are quite solvent sensitive while those for 1,5-DHAQ are not. Due to the intramolecular hydrogen bond between the CO and OH groups, no influence of solvent hydrogen bond acceptors is observed in 1,5-DHAQ. This hydrogen bond gives a stable six member cycle which is not broken even by the strongest hydrogen bond acceptor solvents used in this work, such as DMSO and DMF. The Taft and Kamlet's solvatochromic comparison method was applied for AQ and 2,6-DHAQ. Aromatic solvents and aliphatic amines were not included in the correlations since they strongly deviate suggesting another type of interactions. All the π→π^* bands of AQ and 2,6-DHAQ show strong influence of π^* despite the fact that their dipole moment is zero. Although it would be reasonable to expect that in the absence of a solute dipole moment there is not significant orientation of solvent molecules around the solute molecules, in this case dipolar interactions between solute and solvent due to local effects might be expected. AQ may be considered as formed by two carbonyl groups weakly interacting with the benzene rings; that means that the carbonyl group can behave as an isolated dipole and independently of the other. To detect possible specific interactions between the AQ and aliphatic amines and aromatic hydrocarbons, preferential solvation in mixed solvent was investigated. It is concluded that EDA interactions are important in the solvation of AQ with these compounds as solvents.     

Solvent and concentration effects on the visible spectra of tri-para-dialkylamino-substituted triarylmethane dyes in liquid solutions

We have characterized the spectroscopy properties of crystal violet (CV+) and ethyl violet (EV+) in liquid solutions as a function of the solvent type and dye concentration. The analysis of how solvent properties and dye concentration affects the electronic spectra of these tri-para-dialkylamino substituted tryarylmethane (TAM+) dyes was performed on the basis of two spectroscopic parameters, namely the difference in wavenumber (deltanu) between the maximum and the shoulder that appears in the short-wavelength side of the respective maximum visible band (deltanu = 1/lambda(shoulder)-1/lambda(max) cm(-1)), and the wavelength of the maximum absorption (lambda(max)). The solvent and the concentration effects on lambda(max) and deltanu have indicated that both solute/solute (ion-pairing and dye aggregation) and solute/solvent (H-bonding type) interactions modulate the shape of the visible electronic spectra of these dyes in solution. In solvent with small dieletric constant (epsilon < approximately 10), the formation of ion-pairs represents a major contribution to the shaping of these spectra. Upon increasing dye concentration the formation of ion-pairs was characterized by an increase in deltanu observed concomitantly with a red shift in lambda(max) In chloroform and chlorobenzene the ion-pair association constant of CV+ and EV+ with Cl- ions were found to be in the order of 10(6) and 10(5) M(-1), respectively. In trichloroethylene the association constant for the CV+Cl- pair was 10(8) M(-1). In water, dye aggregation instead of ion-pairing represents a major contribution to the shaping of the visible spectra of CV+ and EV+. Dye aggregation was indicated by an increase in deltanu observed concomitantly with a blue shift in lambda(max) upon increasing dye concentration. The distinct behavior of deltanu for dye aggregation and ion-pairing as a function of dye concentration can therefore assist in the characterization of these two distinct phenomena. The solute/solvent interactions were studied in a series of polar solvents in which solute/solute interactions do not occur in any detectable extent. The dependence found for deltanu as a function of the Kamlet-Tafts solvatochromic parameters (alpha, beta and pi*) is in keeping with previous inferences indicating that the splitting in the overlapped absorption band of CV+ and EV+ in hydroxilated solvents arises from a perturbation in the molecular symmetry induced by hydrogen bonding (donor-acceptor) type interactions with solvent molecules. A distinction between the effects of solute/solute and solute/solvent interactions on the visible spectra of these dyes is provided.     

Effet de Solvant en RMN du Carbone-13. III—Cas d'un Soluté Polaire: l'Acétone dans Divers Solvants et en Phase Gazeuse

The chemical shifts of acetone carbons are measured in the gas phase and in nineteen solvents, thus allowing the separation of the screening constant terms arising from the different kinds of solute–solvent interactions. It is shown that for methyl carbons the van der Waals term σ_w, interpreted using Rummens's theory, is more important than the specific solute–solvent interaction term σ_H. In contrast, for the carbonyl carbon, the term σ_H, which is eight times greater than for the methyl carbons, dominates when dipole–dipole interactions and hydrogen bonding occur. No evidence for an electric field term proportional to the Onsager reaction field is found. But when there are dipole–dipole interactions, σ_H is proportional to the electric field of the dipole of a solvent molecule interacting with the dipole of a solute molecule, the two dipoles being antiparallel. The variation of σ_H with the acetone concentration in a non-associating solvent is interpreted as a consequence of the displacement of the acetone in self-association equilibrium, leading to the determination of the equilibrium constant.     

Volumetric and compressibility studies for (L-arginine + D-maltose monohydrate + water) system in the temperature range of (298.15 to 308.15) K

The density and speed of sound of L-arginine (0.025–0.2 mol kg^?1) in aqueous + D-maltose (0–6 mass% of maltose in water) were obtained at temperatures of (298.15, 303.15 and 308.15) K. The apparent molar volume, limiting apparent molar volume, transfer volume, as well as apparent molar compressibility, limiting apparent molar compressibility, transfer compressibility, pair and triple interaction coefficients, partial molar expansibilities, coefficient of thermal expansion and also the hydration number, were calculated using the experimental density and speed of sound values. The results have been discussed in terms of solute–solute and solute–solvent interactions in these systems. Solute–solvent (hydrophilic–ionic group and hydrophilic–hydrophilic group) interactions were found to be dominating over solute–solute (hydrophobic–hydrophilic group) interactions in the solution, which increases with increase in maltose concentration.     

An ab initio study of solvent polarity and hydrogen bonding effects on the nitrogen NMR shieldings of N,N‐dimethylacetamidine

Density functional theory combined with the polarizable continuum model (PCM) and continuous set of gauge transformations method is applied to investigate the effects of solvent polarity on the nitrogen NMR shieldings of N, N‐dimethylacetamidine. Hydrogen bonding effects on shielding are likewise calculated using a supermolecule approach, where the imino group of the solute is hydrogen bonded with solvent. Theoretical results are compared with published experimental data. The PCM shielding calculations utilizing PCM‐optimized solute geometries yield results comparable to those obtained with the supermolecule approach. Geometry optimization of the solute appears to be more important in PCM shielding calculations than in the supermolecule approach. The large solvent shifts observed in water can only be reproduced when the N·H distance used in the calculation indicates full proton transfer from water to the imino nitrogen of the solute. Copyright © 2002 John Wiley & Sons, Ltd.     

A New Equation for Solute Retention in Liquid Chromatography

In consideration of the adsorption of solvent, diluent and solute molecules on the surface of a stationaryphase, a new equation for solute retention in liquid chromatography is presented. This equation includesthree parameters: the displacement equilibrium constant (Ksd) between the solvent and diluent molecules onthe surface of the stationary phase, the total number(N) of the solvent and diluent molecules released fromthe stationary phase after one solute molecule being adsorbed, and the parameter (I) related to the thermody-namic equilibrium constant for the solute adsorption on the stationary phase. Over the whole concentrationrange of the solvent in the mobile phase, the experimental retention data can be well described by this equa-tion, parameters K~, N and I can be obtained by the regression analysis of the experimental retention data,and consequently the number of the solvent and the diluent molecules displaced by one solute molecule fromthe stationary phase can also be derived at different solvent concentrations in the mobile phase,     

Solvent dependent analysis of isotropic Raman band shape of C=O stretching vibration

The isotropic Raman band shape corresponding to C=O stretching vibration of some molecules has been studied in neat liquids and as a function of solvent concentration using both polar and non-polar solvents. The Raman band shape was analyzed on the basis of correlation with the Lorentzian line shape by employinga simple method of linear curve fitting. In neat liquids and in low solvent concentration region, the band shape was found to be non-Lorentzian. With the gradual increase in solvent concentration the band shape approaches a Lorentzian function. The plot of the correlation coefficient for a Lorentzian shape shows a discontinuity in the intermediate range of solvent concentration. The influence of the structural characteristics of the solute and the solvent systems on the reference mode and various multipolar interactions together with the time varying spatial distribution of solvent molecules with respect to the reference molecule are expected to govern the microenvironmental fluctuations. This may be responsible for the discontinuity in the intermediate solvent concentration region.     

D-甘露醇与D-山梨醇在氯化钠水溶液中的稀释焓

应用等温流动微量热法测定了298.15 K时互为旋光异构体的D-甘露醇与D-山梨醇在不同浓度的氯化钠水溶液中的稀释焓, 利用 McMillan-Mayer 理论计算了D-甘露醇与D-山梨醇在不同浓度的氯化钠水溶液中的焓对相互作用系数. 结果表明, D-甘露醇和D-山梨醇在氯化钠水溶液中的焓对相互作用系数h2均为正值, h2的值随着氯化钠浓度的增加皆逐渐增大, 但D-山梨醇的焓对相互作用系数h2增大的速率[dh2 /dm(NaCl)]比 D-甘露醇的要大. 根据两多元醇分子构象结构的差异， 溶质-溶质相互作用和溶质-溶剂相互作用对结果进行了解释.     

Mass transport in Triton X series nonionic surfactant solutions: a new approach to solute-solvent interactions

Densities, viscosities and tracer diffusion coefficients for solutions of the nonionic surfactants Triton X-45, X-114, X-100 and X-102 in water (except for Triton X-45) and methanol, and for Triton X-100 in three water-methanol mixtures have been measured at 298.15 K and 308.15 K. The activation energy for viscous flow and the contributions to it from solute and solvent have been calculated. Comparison of the Gibbs energies of these systems shows the roles of polyoxyethylene chains of Triton X molecules in the interactions between solute and solvent, and also the effects of solvents on the inter-actions for different solutes. The viscosityB coefficients are positive for all of the surfactant solutions, and the temperature coefficients ofB are negative. In the non-micellar solutions in methanol, values ofB are small and the temperature coefficient ofB is nearly zero. In water-methanol mixtures the critical micelle concentration becomes larger as the methanol content increases up to 40% w/w; micelles are not formed at higher methanol concentrations. TheB coefficient decreases with increasing proportion of methanol in the solvent, and the temperature coefficient ofB changes from a fairly large positive value at low methanol contents to a small negative value at 80% w/w methanol.     

Interactions of Phenylalanine,Tyrosine and Histidine in Aqueous Caffeine Solutions at Different Temperatures

Densities, ρ, viscosities, η, and refractive indices, n_D of aqueous caffeine (0.5 M) and of solutions of amino acids, l‐phenylalanine (Phe), l‐tyrosine (Tyr) and l‐histidine (His), (0.01–0.05 M) in aqueous‐caffeine have been measured at 298.15, 303.15, 308.15 and 313.15 K. From these experimental data, apparent molar volume, ?_v, limiting partial molar volume, ?^º_ν and the slope, S_v, transfer volume, ?^º_ν,tr, Falkenhagen coefficient, A, Jones‐Dole coefficients, B, free energies of activation per mole of solvent, Δμ^o#₁ and per mole of solute, Δμ^o#₂, enthalpy, ΔH^* and entropy, ΔS^* of activation of viscous flow, and molar refraction, R_m were calculated. The results are interpreted from the point of view of solute‐solvent and solute‐solute interactions in these systems. It has been observed that there exist strong solute‐solvent and weak solute‐solute interactions in these systems. Further, the solute‐solvent interactions decrease, whereas solute‐solute interactions increase with rise in temperature. It is observed that these amino acids act as structure‐makers in aqueous‐caffeine solvent. The thermodynamics of viscous flow have also been discussed.     

Volumes and heat capacities of some watersurfactant-alcohol ternary systems

The densities and volumetric heat capacities of urea and alcohols were measured in aqueous solutions of octylammonium bromide (OABr) and of OABr in aqueous urea and alcohol solutions. The alcohols studies were methanol, ethanol, 1-propanol, 2-propanol, n-butanol, t-butanol, n-pentanol, n-hexanol and 2-butoxyethanol (BE). In most experiments, the concentration of the reference solute was kept low, and volumes and heat capacities of transfer from water to the mixed solvent were calculated. A more complete study was made with the system BE-OABr-H₂O where both solutes were systematically changed. The observed trends in the thermodynamic functions can be explained through three effects: interactions between the reference solute and the cosolvent in the premiceller region of the surfactant or pre-aggregation region of the alcohol, a distribution of the reference solute between water and the micelle or microphase and an equilibrium displacement of the system, monomer-aggregate, in the vicinity of the reference solute.     

Molecular rotation as a tool for exploring specific solute-solvent interactions.

Solute-solvent interactions play an important role in determining the physicochemical properties of liquids and solutions. As a consequence, understanding these interactions has been one of the long-standing problems in physical chemistry. This Minireview describes our approach towards attaining this goal, which is to investigate rotational relaxation of a pair of closely related, medium-sized nondipolar solutes in a set of appropriately chosen solvents. Our studies indicate that solute-solvent hydrogen bonding significantly hinders solute rotation. We have also examined the role of solvent size both in the absence and presence of specific interactions and it has been observed that the size of the solvent has a bearing on solute rotation especially in the absence of specific interactions. Our results point to the fact that only strong solute-solvent hydrogen bonds have the ability to impede the rotation of the solute molecule because, in such a scenario, hydrogen-bonding dynamics and rotational dynamics transpire on comparable time scales. This aspect has been substantiated by measuring the reorientation times of the chosen solutes in solvents such as ethanol and trifluoroethanol, which have distinct hydrogen-bond donating and accepting abilities, and correlating them with solute-solvent interaction strengths. As an alternative treatment, it has been shown that specific interactions between the solute and the solvent can be modeled as dielectric friction with the extended charge distribution model. This approach is not unrealistic considering the fact that specific as well as non-specific interactions are electrostatic by nature and the differences between them are subtle.     

Physicochemical Properties of Amino Acids in AqueousCaffeine Solution at 25, 30, 35 and 40 ℃

Density, viscosity, and refractive index, for glycine, DL-alanine, L-serine and DL-valine have been determined in aqueous solution of 0.05 mol/kg caffeine as a function of amino acid （AA） concentration at 25, 30, 35, and 40 ℃ The density data have been used to compute apparent molar volume. The partial molar volume （limiting apparent molar volume） was obtained by applying the Masson＇s equation. The viscosity data have been analyzed by means of Jones-Dole equation. The values of Falkenhagen coefficient and Jones-Dole coefficient thus obtained are used to interpret the solute-solute and solute-solvent interactions, respectively. Hydration number was also computed. The transition-state theory was applied to obtain the activation parameters of viscous flow, i.e., free energy of activation per mole of solvent, and solute. The enthalpy and entropy of activation of viscous flow were computed for the system. Refractive index was used to calculate molar refractivity of the mixtures. The results have been interpreted in the light of various interactions occurring between the components of the mixtures under applied experimental conditions.     

Solvent effects on electronic absorption spectra of nitrochlorobenzenes,nitrophenols and nitroanilines—II. Studies in polar solvents

A general method for estimating the specific solute—solvent interaction energies (E_s) from analysis of solvent effects on electronic absorption spectra has been outlined. The E_s values for a number of mono- and disubstituted benzenes in a variety of solvents such as water, alcohols, chloroalkanes, ether and acetonitrile have been estimated and the results are discussed in relation to the interacting groups in the solute and solvent molecules. These interactions, which can be classified as H bonding or EDA type, are mainly electrostatic in nature. These studies indicate that, in the case of disubstituted benzenes, the intramolecular electronic interactions are stronger than the solute—solvent interactions.     

Solubility diagrams in solvent-antisolvent systems by titration calorimetry

Isothermal titration calorimetry (ITC) has been used to develop a method to construct the solid-liquid equilibrium line in ternary systems containing the solute to precipitate and an aqueous mixed solvent. The method consists in measuring the heat of dissolution of a solid component (the solute) during successive additions of the liquid solvent. The cumulated heat, resulting from the successive heat peaks obtained for the different injections of known volumes of solvent, plotted vs. the ratio of the numbers of moles n _solvent/n _solute is represented by two nearly straight lines. The intercept of the two lines gives the solubility limit and the corresponding enthalpy of dissolution of the solute in the solvent. Solubility diagrams have been established at 303.15 K in binary mixed solvents ethanol-water over the whole concentration range for seven compounds of pharmaceutical interest, namely: urea, phenylurea, l-valine, dl-valine, l-valine ethyl ester hydrochloride, tris(hydroxymethyl)amino methane.