Single particle and pair dynamics in water-formic acid mixtures containing ionic and neutral solutes: nonideality in dynamical properties

A series of molecular dynamics simulations of water-formic acid mixtures containing either an ionic solute or a neutral hydrophobic solute has been performed to study the extent of nonideality in the dynamics of these solutes for varying composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to show nonideal behavior with variation of composition, and similar nonideality is also observed for the diffusion and orientational relaxation of solvent molecules in these mixtures. The diffusion coefficient of a neutral hydrophobic solute, however, decreases monotonically with increase in water concentration. We have also investigated some of the pair dynamical properties such as water-water and water-formic acid hydrogen bond relaxation and residence dynamics of water molecules in water and formic acid hydration shells. The lifetimes of water-water hydrogen bonds are found to be longer than those between formic acid carbonyl oxygen-water hydrogen bonds, whereas the lifetimes of formic acid hydroxyl hydrogen-water hydrogen bonds are longer than those of water-water hydrogen bonds. In general, the hydrogen bond lifetimes for both water-water and water-formic acid hydrogen bonds are found to decrease with increase in water concentration. Residence times of water molecules also show the same trend with increase in formic acid concentration. Interestingly, these pair dynamical properties show a monotonic dependence on composition without any maximum or minimum and behave almost ideally with respect to changes in the composition of the mixtures. The present calculations are performed with fixed-charge nonpolarizable models of the solvent and solute molecules without taking into account many-body polarization effects in an explicit manner.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Dynamics of ionic and hydrophobic solutes in water-methanol mixtures of varying composition

We have carried out a series of molecular-dynamics simulations of water-methanol mixtures containing either an ionic or a neutral atomic solute to investigate the effects of composition of the mixture on the diffusion of these solutes. Altogether, we have considered 17 different systems of varying composition ranging from pure water to pure methanol. The diffusion coefficients of ionic solutes are found to show nonideal behavior with variation of composition of the solvent mixture. The extent of nonideality of the solute diffusion is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and methanol molecules in these mixtures and is attributed to the enhanced stability of the hydrogen bonds and formation of interspecies complexes in the mixtures. The neutral solute shows characteristics of hydrophobic solvation and its diffusion decreases monotonically with increase of methanol concentration. The present simulation results are compared with those of experiments wherever available.     

某些非极性和极性溶质在苯和甲醇及二者混合物中的偏摩尔体积V0, 298.15 K

密度法测定了298.15 K下乙醇、环己烷、三氯甲烷、甲苯、丙酮、四氯化碳、乙腈、二甲基甲酰胺、二甲基亚砜在甲醇或苯及两者混合物中的无限稀释偏摩尔体积. 密度测定所用溶液中溶质的浓度范围是0.2一1.5 m; 甲醇和苯混合物是全组成比范围. 溶质偏摩尔体积随甲醇-苯组成比的变化趋势反映了几种分子间相互作用结果即三种分子间物理型分子间相互作用; 溶质与甲醇分子氢键缔合相互作用; 溶质同甲醇或苯的弱络合作用。     

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.     

Solute-solvent complex switching dynamics of chloroform between acetone and dimethylsulfoxide-two-dimensional IR chemical exchange spectroscopy

Hydrogen bonds formed between C-H and various hydrogen bond acceptors play important roles in the structure of proteins and organic crystals, and the mechanisms of C-H bond cleavage reactions. Chloroform, a C-H hydrogen bond donor, can form weak hydrogen-bonded complexes with acetone and with dimethylsulfoxide (DMSO). When chloroform is dissolved in a mixed solvent consisting of acetone and DMSO, both types of hydrogen-bonded complexes exist. The two complexes, chloroform-acetone and chloroform-DMSO, are in equilibrium, and they rapidly interconvert by chloroform exchanging hydrogen bond acceptors. This fast hydrogen bond acceptor substitution reaction is probed using ultrafast two-dimensional infrared (2D-IR) vibrational echo chemical exchange spectroscopy. Deuterated chloroform is used in the experiments, and the 2D-IR spectrum of the C-D stretching mode is measured. The chemical exchange of the chloroform hydrogen bonding partners is tracked by observing the time-dependent growth of off-diagonal peaks in the 2D-IR spectra. The measured substitution rate is 1/30 ps for an acetone molecule to replace a DMSO molecule in a chloroform-DMSO complex and 1/45 ps for a DMSO molecule to replace an acetone molecule in a chloroform-acetone complex. Free chloroform exists in the mixed solvent, and it acts as a reactive intermediate in the substitution reaction, analogous to a SN1 type reaction. From the measured rates and the equilibrium concentrations of acetone and DMSO, the dissociation rates for the chloroform-DMSO and chloroform-acetone complexes are found to be 1/24 ps and 1/5.5 ps, respectively. The difference between the measured rate for the complete substitution reaction and the rate for complex dissociation corresponds to the diffusion limited rate. The estimated diffusion limited rate agrees well with the result from a Smoluchowski treatment of diffusive reactions.     

NMR Relaxation Studies in Aqueous Solutions of Amino Acids

Abstract

The proton magnetic resonance (PMR) spin-lattice and spin-spin relaxation times (T1 and T2) were measured in aqueous solutions of glycine and L-proline as a function of solute concentrations and at a temperature of 32°C. The relaxation times were measured using Bruker PC 120 NMR process analyser. The relaxation times were found to decrease with increase of solute concentrations. The results are interpreted on the basis of flickering cluster model and hydrogen bond formation between solute and solvent molecules.     

Properties and transport behavior of perfluorotripentylamine (FC-70)-doped amorphous teflon AF 2400 films

Teflon AF 2400 films are known to imbibe solvents, making films in the presence of solvents less fluorous than they might otherwise be. Herein, we demonstrate that doping films with perfluorotripentylamine (Fluorinert FC-70) maintains the fluorous nature of Teflon AF 2400 and improves transport selectivity for fluorine-containing organic compounds. Density measurements on the FC-70-doped films reveal that free volume decreases dramatically as the dopant concentration increases (0-12 wt %) and then increases to approach that of pure FC-70. Remarkably, films from 0 to 12 wt % FC-70 have the same w/v concentration of Teflon AF 2400, indicating that FC-70 fills the free volume of Teflon AF 2400. This is consistent with the observed increased storage modulus and significant decrease (compared to undoped films) of solute diffusion coefficients in the same range of FC-70 concentrations. In contrast, FC-70 at concentrations greater than 12 wt % dilutes Teflon AF 2400, leading to a decrease of storage modulus and dramatic increase in solute diffusion coefficients. Sorption of chloroform decreases from 11.8 g of chloroform/100 g of film (pure Teflon film) to 3.8 g of chloroform/100 g of film (27 wt % FC-70-doped Teflon film), less than the solubility of chloroform in pure FC-70 (4.06 g of chloroform/100 g of FC-70). Solute partition coefficients from chloroform to FC-70-doped films generally decrease with increased dopant concentration. However, within a series of toluenes and nitrobenzenes, selectivity for F-containing solutes over analogous H-containing solutes increases as dopant concentration increases if the substitution is on the aromatic ring but not if it is on the methyl group (toluene). Transport (partitioning × diffusion) rates, as they involve both thermodynamic and kinetic factors, are not simply related to composition.     

Enthalpies of solvation of alkylated uracils in water,nonaqueous and mixed solvents

The enthalpies of solution of uracil and its alkylated derivatives in water, methanol, N,N-dimethylformamide (DMF) and water+DMF mixtures were measured at 25°C. The enthalpies of solvation were determined. The enthalpies of cavity formation, corresponding to the enthalpies of solvent-solvent interactions were calculated and the enthalpies of solute-solvent interactions were obtained. The presence of the alkyl groups was found to have different effects on the enthalpy of interaction depending on the position and size of the substitution. The effect of alkylation at the nonpolar side of the uracil ring was found to arise mostly from the enhancement of the van der Waals interactions. The alkyl substitutions at the polar side resulted also in the removal of the solvent molecules interacting specifically with the polar groups of uracil. The enthalpy of those specific interactions was determined and found to be stronger in methanol and DMF than in water. Enthalpies of solvation in the binary water+DMF solvent were found to depend in a nonlinear way on the solvent composition. The nonlinearities in the water-rich region were found to arise from the decay of the hydrophobic hydration of the solutes with the increasing DMF content. The substitution of two methyl groups caused the uracil molecule to bahave as a predominantly hydrophobic solute. The nonlinearities in the DMF-rich region were found only for those solutes which can form hydrogen bonds with DMF.     

A first principles molecular dynamics study of lithium atom solvation in binary liquid mixture of water and ammonia: structural, electronic, and dynamical properties

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.     

Solvent-hindered intramolecular vibrational redistribution

Ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations of Mn(2)(CO)(10) in a series of linear alcohols reveal that the rate of intramolecular vibrational redistribution among the terminal carbonyl stretches is dictated by the average number of hydrogen bonds formed between the solute and solvent. The presence of hydrogen bonds was found to hinder vibrational redistribution between eigenstates, while leaving the overall T(1) relaxation rate unchanged.     

Enthalpies of dissolution of propylene carbonate,acetonitrile, and 1,4-dioxane in mixtures of water with acetone or dimethyl sulfoxide

Enthalpies of dissolution of acetonitrile, propylene carbonate, and 1,4-dioxane in mixtures of water with acetone or DMSO were measured in the whole concentration range of the mixed solvents. Standard enthalpies of dissolution and enthalpies of transfer of solutes from water to its mixtures with acetone or DMSO were determined. In the region of small proportions of the nonaqueous component, the enthalpy of cavity formation in the mixed solvent makes the main contribution to the variation of the enthalpy of dissolution. An increase in the proportion of the nonaqueous component leads to competition between the contributions of cavity formation and specific interaction between the solute and the solvent during solvation.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1747–1752, September, 1995.     

Orientational effects of hydrogen bonding in liquid-crystalline solutions containing Schiff bases

Intermolecular hydrogen bonding in binary mixtures containing nematogenic Schiff bases as solvents and proton-donating non-mesomorphic solutes has been considered. Reasons for the anomalous concentration dependences of solute order parameters are discussed. A solution structure model of acetic acid in nematic solvents is proposed; constants of complex-dimer equilibrium and coefficients of the orientational correlation of the non-mesogenic solute are calculated on the basis of this model. Hydrogen bonded complex structure using ¹³C NMR has been studied and stability constants in isotropic solutions in chloroform have been calculated. The influence of the solvent orientational ordering on the complex stability is discussed. Data on the solvation isotopic effects in the solutions investigated, which confirm the adequacy of the model are given.     

Conformational studies of new pseudotripeptide with pyrazine amidoxime motif and simplified analogs using IR,NMR spectroscopy,and molecular dynamic simulations

Solution structures of new pyrazine-based pseudotripeptide with amidoxime function and simplified pseudodipeptide analogs were determined by a combination of IR and NMR spectroscopic studies and molecular dynamic simulations using explicit chloroform as a solvent. It was found that proline-phenylalanine dipeptide residue and amidoxime moiety in o-position are essential for intramolecular hydrogen bonding including a seven-membered γ-turn formation. In addition, a cis/trans equilibrium study was present for prolyl amides in polar solvents (D₂O and DMSO). A phenylalanine substituent was found to exhibit profound effect on thermodynamic parameters in prolyl peptides. The presence of intramolecular hydrogen bonds dramatically increases the amount of trans isomer in non-hydrogen-bonding CHCl₃ and significantly favor cis isomer in hydrogen-bonding solvents such as DMSO and D₂O. All molecules are not cytotoxic therefore they can be further studied in relation to potent biological activities.

     

Observation of Unusual Molecular Diffusion Behaviour below the Lower Critical Solution Temperature of Water/2‐Butoxyethanol Mixtures by using Fluorescence Correlation Spectroscopy

The effect of solute affinity on solute diffusion in binary liquids well below the lower critical solution temperature (LCST) was studied by using fluorescence correlation spectroscopy. We measured the hydrodynamic radii of a hydrophobic and an amphiphilic fluorescent dye under systematic variation of the relative molar fractions of water/2‐butoxyethanol and, for comparison, of water/methanol mixtures, which do not show phase separation. We found that the apparent hydrodynamic radius of the hydrophobic dye almost doubled in water/2‐butoxyethanol, whereas it remained largely unchanged for the amphiphilic dye and in water/methanol mixtures. Our results indicate that the translational diffusion of solutes is influenced by transient local solution structures, even at temperatures well below the LCST. We conclude that, even far below LCST, different solutes can experience different environments in binary liquid mixtures depending on both the solute and solvent properties, all of which impact their reactivity.     

Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H...Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.     

How fluorous is poly(2,2-bis(trifluoromethyl)-4,5-difluoro-1,3-dioxide-co-tetrafluoroethylene) (Teflon AF)?

We investigated poly(2,2-bis(trifluoromethyl)-4,5-difluoro-1,3-dioxide-co-tetrafluoroethylene (Teflon AF 2400, Tg = 240 degrees C), a stable and permeable fluorous polymer, as a transport/extraction medium for solutes for the first time. From the study of transport behavior of a series of solutes (in chloroform solution) through the film, and detailed measurement of the partitioning and diffusion of benzene in the film, we showed that the Teflon AF film is influenced by the solvent to which it is exposed. In particular, the solvent chloroform is sorbed in the film at a high concentration of 1.13 M. This plasticizes the film, so that the diffusion coefficient of benzene is about 3 orders of magnitude larger in a chloroform-equilibrated film compared to benzene coming from the gas phase into a "dry" film. The partition ratio of the polar solute 3-hydroxypyridine is dramatically higher in the case of partitioning between chloroform and a film (0.02) compared to partitioning between chloroform and the fluorous solvent FC-72 (6.7 x 10-5). Krytox FSH, a carboxylic acid terminated perfluoropolyether, plasticizes films. Tg in a 50% (w/w) film decreases to -40 degrees C. This carboxylic acid is capable of molecular recognition in the film. The noncovalent association between Krytox FSH (0.13 M in the film) and 3-hydroxypyridine increases the distribution ratio of the polar solute into the film by 41 times. In comparison, the partition ratio into a fluorocarbon solvent (FC-72) increases 15 000 times under the same conditions. As a result of imbibing organic solvent, the films of Teflon AF 2400 are not as fluorous as a fluorous liquid.     

Structure and dynamics of water confined in dimethyl sulfoxide.

We study the structure and dynamics of hydrogen-bonded complexes of H2O/D2O and dimethyl sulfoxide (DMSO) by infrared spectroscopy, NMR spectroscopy and ab initio calculations. We find that single water molecules occur in two configurations. For one half of the water monomers both OH/OD groups form strong hydrogen bonds to DMSO molecules, whereas for the other half only one of the two OH/OD groups is hydrogen-bonded to a solvent molecule. The H-bond strength between water and DMSO is in the order of that in bulk water. NMR deuteron relaxation rates and calculated deuteron quadrupole coupling constants yield rotational correlation times of water. The molecular reorientation of water monomers in DMSO is two-and-a-half times slower than in bulk water. This result can be explained by local structure behavior.     

Non‐Ideal Behaviour and Solution Interactions in Binary DMSO Solutions

The structure and dynamics of hydrogen‐bonded structures are of significant importance in understanding many binary mixtures. Since self‐diffusion is very sensitive to changes in the molecular weight and shape of the diffusing species, hydrogen‐bonded associated structures in dimethylsulfoxide–methanol (DMSO–MeOH) and DMSO–ethanol (DMSO–EtOH) mixtures are investigated using nuclear magnetic resonance (NMR) diffusion experiments and molecular dynamics (MD) simulations over the entire composition range at 298 K. The self‐diffusion coefficients of DMSO–MeOH and DMSO–EtOH mixtures decrease by up to 15% and 10%, respectively, with DMSO concentration, indicating weaker association as compared to DMSO–water mixtures. The calculated heat of mixing and radial distribution functions reveal that the intermolecular structures of DMSO–MeOH and DMSO–EtOH mixtures do not change on mixing. DMSO–alcohol hydrogen‐bonded dimers are the dominant species in mixtures. Direct comparison of the simulated and experimental data afford greater insights into the structural properties of binary mixtures.