Ion solvation dynamics in water–methanol and water– dimethylsulfoxide mixtures

Selective Self-Assembly of 5-Fluorouracil through Nonlinear Solvent Response Modulates Membrane Dynamics

Controllable self-assembly and understanding of the interaction between single metabolite fibrils and live-cell membranes have paramount importance in providing minimal treatment in several neurodegenerative disorders. Here, utilizing the nonlinear nature and peculiar hydrogen bonding behavior of the dimethyl sulfoxide (DMSO)-water mixture, the selective self-assembly of a single metabolite 5-fluorouracil (5-FU) is achieved. A direct correlation between water availability and selective self-assembly of 5-FU is ratified from the excited-state dynamics. The specific fibrillar structures of 5-FU exhibit a great potential to modulate live cell membrane fluidity and model membrane lipid distribution. After 5-FU fibril addition, a disorder of H-bonded water molecules arises several layers beyond the first hydration shell of the polar headgroups, which essentially modifies interfacial water structure and dynamics. Overall, our results shed light on the role of solvent to govern specific self-assembly and also lay the foundation accounting for the earlier stage of several diseases and multidrug resistance.

195Pt NMR and Molecular Dynamics Simulation Study of the Solvation of [PtCl6]2- in Water-Methanol and Water-Dimethoxyethane Binary Mixtures

The experimental ¹⁹⁵Pt NMR chemical shift, δ(¹⁹⁵Pt), of the [PtCl₆]^2- anion dissolved in binary mixtures of water and a fully miscible organic solvent is extremely sensitive to the composition of the mixture at room temperature. Significantly nonlinear δ(¹⁹⁵Pt) trends as a function of solvent composition are observed in mixtures of water-methanol, or ethylene glycol, 2-methoxyethanol, and 1,2-dimethoxyethane (DME). The extent of the deviation from linearity of the δ(¹⁹⁵Pt) trend depends strongly on the nature of the organic component in these solutions, which broadly suggests preferential solvation of the [PtCl₆]^2- anion by the organic molecule. This simplistic interpretation is based on an accepted view pertaining to monovalent cations in similar binary solvent mixtures. To elucidate these phenomena in detail, classical molecular dynamics computer simulations were performed for [PtCl₆]^2- in water-methanol and water-DME mixtures using the anionic charge scaling approach to account for the effect of electronic dielectric screening. Our simulations suggest that the simplistic model of preferential solvation of [PtCl₆]^2- by the organic component as inferred from nonlinear δ(¹⁹⁵Pt) trends is not entirely accurate, particularly for water-DME mixtures. The δ(¹⁹⁵Pt) trend in these mixtures levels off for high DME mole fractions, which results from apparent preferential location of [PtCl₆]^2- anions at the borders of water-rich regions or clusters within these inherently micro-heterogeneous mixtures. By contrast in water-methanol mixtures, apparently less pronounced mixed solvent micro-heterogeneity is found, suggesting the experimental δ(¹⁹⁵Pt) trend is consistent with a more moderate preferential solvation of [PtCl₆]^2- anions. This finding underlines the important role of solvent-solvent interactions and micro-heterogeneity in determining the solvation environment of [PtCl₆]^2- anions in binary solvent mixtures, probed by highly sensitive ¹⁹⁵Pt NMR. The notion that preferential solvation of [PtCl₆]^2- results primarily from competing ion-solvent interactions as generally assumed for monatomic ions, may not be appropriate in general.

Preferential solvation and ion association properties in aqueous dimethyl sulfoxide solutions

We study the solvation and the association properties of ion pairs in aqueous dimethyl sulfoxide (DMSO) solution by atomistic molecular dynamics (MD) simulations. The ion pair is composed of two lithium and a single sulfonated diphenyl sulfone ion whose properties are studied under the influence of different DMSO concentrations. For increasing mole fractions of DMSO, we observe a non-ideal behavior of the solution as indicated by the derivatives of the chemical activity. Our findings are complemented by dielectric spectra, which also verify a complex DMSO-water mixing behavior. In agreement with these results, further simulation outcomes reveal an aqueous homoselective solvation of the ion species which fosters the occurrence of pronounced ion association constants at higher DMSO mole fractions. The consequences of this finding are demonstrated by lower ionic conductivities for increasing concentrations of DMSO.

A molecularly based theory for electron transfer reorganization energy

Using field-theoretic techniques, we develop a molecularly based dipolar self-consistent-field theory (DSCFT) for charge solvation in pure solvents under equilibrium and nonequilibrium conditions and apply it to the reorganization energy of electron transfer reactions. The DSCFT uses a set of molecular parameters, such as the solvent molecule's permanent dipole moment and polarizability, thus avoiding approximations that are inherent in treating the solvent as a linear dielectric medium. A simple, analytical expression for the free energy is obtained in terms of the equilibrium and nonequilibrium electrostatic potential profiles and electric susceptibilities, which are obtained by solving a set of self-consistent equations. With no adjustable parameters, the DSCFT predicts activation energies and reorganization energies in good agreement with previous experiments and calculations for the electron transfer between metallic ions. Because the DSCFT is able to describe the properties of the solvent in the immediate vicinity of the charges, it is unnecessary to distinguish between the inner-sphere and outer-sphere solvent molecules in the calculation of the reorganization energy as in previous work. Furthermore, examining the nonequilibrium free energy surfaces of electron transfer, we find that the nonequilibrium free energy is well approximated by a double parabola for self-exchange reactions, but the curvature of the nonequilibrium free energy surface depends on the charges of the electron-transferring species, contrary to the prediction by the linear dielectric theory.

Ion specificity at a low salt concentration in water-methanol mixtures exemplified by a growth of polyelectrolyte multilayer

By use of a quartz crystal microbalance with dissipation (QCM-D), we have investigated the specific ion effect on the growth of poly(sodium 2-acrylamido-2-methylpropanesulfonate)/poly(diallyldimethylammonium chloride) multilayer at a salt concentration as low as 2.0 mM in water-methanol mixtures. QCM-D results demonstrate that specific ion effect can be observed in methanol and water-methanol mixtures though it is negligible in water. Moreover, the specific ion effect is amplified as the molar fraction of methanol (xM) increases from 0% to 75% but is weakened again with the further increase of xM from 75% to 100%. Nuclear magnetic resonance measurements reveal that the counterion-polyelectrolyte segment interactions may not account for the observed ion specificity. By extending the Collins' concept of matching water affinities to methanol and water-methanol mixtures, we suggest that the ion-solvent interactions and the resulted counterion-charged group interactions are responsible for the occurrence of the specific ion effect. The conductivity measurements indicate that water and methanol molecules may form complexes, and the change of relative proportion of complexes with the xM causes the amplification or weakening of the specific ion effect.

SINGLE-PARTICLE AND PAIR DYNAMICAL PROPERTIES OF ACETONE–METHANOL MIXTURES CONTAINING CHARGED AND NEUTRAL SOLUTES: A MOLECULAR DYNAMICS STUDY

The dynamical properties of acetone–methanol mixtures containing either an ionic or a neutral hydrophobic solute are investigated by means of a series of molecular dynamics simulations. The primary goal has been to study how the solute and solvent dynamical properties change with variation of composition of the mixture ranging from pure acetone to pure methanol. The variations of structure and energetics of the mixture with composition are also calculated. The diffusion coefficients of both ionic and neutral solutes are found to show nonlinear variation with composition of the mixture, although the extent of nonlinearity in the diffusion of the neutral solute is much weaker. Calculations of appropriate solute-solvent distribution functions reveal the extent and nature of selective solvation of these solute species which play a role in determining the nonideal dynamical characteristics of these solutes. The free energies of solvation of the ionic solutes are also calculated and the results are discussed in the context of their dynamical behavior. The hydrogen bond statistics and dynamics of these mixtures are also calculated over their entire composition range. The energies and lifetimes of hydrogen bonds between an acetone and a methanol molecule or between two methanol molecules are found to increase with increase of acetone mole fraction of the mixture. Residence times of methanol molecules in solvation shells of acetone and methanol are also found to follow the same trend as relaxation times. However, these pair dynamical properties show essentially linear dependence on composition, thus behave almost ideally with respect to changes in composition of the mixture.

Insights into ion specificity in water-methanol mixtures via the reentrant behavior of polymer

In the present work, we have for the first time systematically investigated the ion specific reentrant behavior of poly(N-isopropylacryamide) (PNIPAM) in water-methanol mixtures. Turbidity measurements demonstrate that SCN(-) and ClO(4)(-) depress the reentrant transition, whereas other anions enhance the transition. As the anion changes from chaotropic to kosmotropic, the minimum critical phase transition temperature (T(min)) decreases and the corresponding volume fraction of methanol (X(M)) shifts to a larger value. Our results demonstrate that anion specificity is due to the anionic structure making/breaking effect on water/methanol complexes. Cations are found to have a lesser but still significant effect on the reentrant transition, and as T(min) decreases the corresponding X(M) also shifts to larger values as with the anions. Our studies show that cation specificity is induced by specific interactions between cations and PNIPAM chains. Furthermore, both anion and cation specificities are amplified as X(M) is increased due to the formation of additional water/methanol complexes. Calorimetry measurements demonstrate that the ion specificity is dominated by changes in entropy.

TIME DEPENDENT DENSITY FUNCTIONAL METHODS AND THEIR APPLICATION TO CHEMICAL PHYSICS

This article reviews microscopic development of time dependent functional method and its application to chemical physics. It begins with the formulation of density functional theory. The time dependent extension is discussed after the equilibrium formulation. Its application is explained by solvation dynamics. In addition, it reviews studies of nonlinear effects on polar liquids and simple mixtures.

A gas-phase study of the preferential solvation of Mn(2+) in mixed water/methanol clusters

The kinetic shift that exists between two competing unimolecular fragmentation processes has been used to establish whether or not gas-phase Mn(2+) exhibits preferential solvation when forming mixed clusters with water and methanol. Supported by molecular orbital calculations, these first results for a metal dication demonstrate that Mn(2+) prefers to be solvated by methanol in the primary solvation shell.

Absorption and emission lineshapes and solvation dynamics of NO in supercritical Ar

We perform a theoretical study of electronic spectroscopy of dilute NO in supercritical Ar fluid. Absorption and emission lineshapes for the A(2)Sigma(+)<--X(2)Pi Rydberg transition of NO in argon have been previously measured and simulated, which yielded results for the NO/Ar ground- and excited-state pair potentials [Larregaray et al., Chem. Phys. 308, 13 (2005)]. Using these potentials, we have performed molecular dynamics simulations and theoretical statistical mechanical calculations of absorption and emission lineshapes and nonequilibrium solvation correlation functions for a wide range of solvent densities and temperatures. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down at near-critical temperature due to dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.

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.

Vibrational Spectroscopy and Dynamics of Azide Ion in Ionic Liquid and Dimethyl Sulfoxide Water Mixtures

Steady-state and time-resolved infrared spectroscopy of the azide (N(3)-) anion has been used to characterize aqueous mixtures both with the ionic liquid (IL) 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF(4)]) and with dimethyl sulfoxide (DMSO). In the DMSO-water mixtures, two anion vibrational bands are observed for low water mole fractions (0 > X(w) > 0.25), which indicates a heterogeneous ion solvation environment. The band at 2000 cm(-1) observed for neat DMSO does not shift but decreases in amplitude as the amount of water is increased. Another band appears at slightly higher frequency at low X(w) (=0.05). As the amount of water is increased, this band shifts to higher frequency and becomes stronger and is attributed to azide with an increasing degree of hydration. At intermediate and high X(w), a single band is observed that shifts almost linearly with water mole fraction toward the bulk water value. The heterogeneity is evident from the infrared pump-probe studies in which the decay times depend on probe frequency at low mole fraction. For the azide spectra in IL-water mixtures, a single azide band is observed for each mole fraction mixture. The azide band shifts almost linearly with mole fraction, indicating nearly ideal mixing behavior. As with the DMSO-water mixtures, the time-resolved IR decay times are probe-frequency-dependent at low mole fraction, again indicating heterogeneous solvation. In both the DMSO and IL mixtures with water, the relaxation times are slower than would be expected from ideal mixing, suggesting that vibrational relaxation of azide is more sensitive than its vibrational frequency to the solvent structure. The results are discussed in terms of preferential solvation and the degree to which the azide shift and vibrational relaxation depend on the degree of water association in the mixtures.

Nonlinear effects on solvation dynamics in simple mixtures

The authors applied the time dependent density functional method (TDDFM) and a linear model to solvation dynamics in simple binary solvents. Changing the solute-solvent interactions at t=0, the authors calculated the time evolution of density fields for solvent particles after the change (t>0) by the TDDFM and linear model. First, the authors changed the interaction of only one component of solvents. In this case, the TDDFM showed that the solvation time decreased monotonically with a mole fraction of the solvent strongly interacting with the solute. The monotonical decreases agreed with experimental results, while the linear model did not reproduce these results. The authors also calculated the solvation time by changing the interaction of both components. The calculation showed that the mole fraction dependence had the peak. The TDDFM presented a much higher peak than the linear model. The difference between the TDDFM and the linear model was caused by a nonlinear effect on an exchange process of solvent particles.

Solute size effects on the solvation structure and diffusion of ions in liquid methanol under normal and cold conditions

We have performed a series of molecular dynamics simulations of alkali metal (Li+, Na+, K+, Rb+, and Cs+) and halide (F-, Cl-, Br-, and I-) ions in liquid methanol at two different temperatures to investigate the effects of ion size on the hydration structure and diffusion of ions in methanol under normal and cold conditions. Simulations are also carried out for some of the larger cations such as I+, (CH3)4N+, and (C2H5)4N+ and also neutral alkali metal atoms in methanol at both temperatures. With the increase of ion size, the diffusion coefficients of both positive and negative ions are found to show anomalous behavior. For cations, it is found that the maximum of the diffusion coefficient versus ion size curve occurs at the rather large cation of (CH3)4N+ unlike in water where the maximum occurs at the relatively smaller ion of Rb+. For halide ions, the anomalous behavior, i.e., the increase of diffusion with ion size, continues up to iodide ion and no maximum is observed. These results are in good agreement with experimental observations. The diffusion coefficients of neutral atoms are found to be greater in methanol than that in water and they decrease monotonically with solute size, whereas the diffusion coefficients of the corresponding ions are found to be smaller in methanol. Accordingly, an ion experiences a smaller Stokes friction and a higher dielectric friction in methanol than in water. These contrasting effects are believed to be responsible for the shift of the maximum of ion diffusion toward a larger ion size when compared with similar anomalous size dependence in liquid water.

The effects of solute-solvent electrostatic interactions on solvation dynamics in supercritical CO2

We present here the results of molecular-dynamics simulation of solvation dynamics in supercritical CO(2) at a temperature of about 1.05T(c), where T(c) is the critical temperature, and at a series of densities ranging from 0.4 to 2.0 of the critical density rho(c). We focus on electrostatic solvation dynamics, representing the electronic excitation of the chromophore as a change in its charge distribution from a quadrupolar-symmetry ground state to a dipolar excited state. Two perturbations are considered, corresponding to different magnitudes of solute excited-state dipoles, denoted as d5 and d8. The d8 solute is more attractive, leading to a larger enhancement in CO(2) clustering upon solute electronic excitation. This has a large impact on solvation dynamics, especially at densities below rho(c). At these densities, solvation dynamics is much slower for the d8 than for the d5 solute. For both solutes, solvation dynamics becomes faster at densities above rho(c) at which solvent clustering diminishes. We show that the slowest solvation time scale is associated with solvent clustering and we relate it to solute-solvent mutual translational diffusion and the extent of change in effective local density resulting from solute electronic excitation.

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.

Molecular-dynamics simulations of dimethylsulfoxide-methanol mixtures

We present molecular-dynamics (MD) computer simulation results for the local structures, hydrogen (H)-bond distribution, and dynamical properties of methanol (MeOH) and dimethylsulfoxide (DMSO) binary mixtures at ambient conditions over the entire composition range. The simulated heat of mixing and site-site pair distribution functions suggest that the intermolecular structures of the pure liquids are not markedly altered upon mixing. Nevertheless, H-bonding statistics show that aggregates of the type 1DMSO:1MeOH are formed and represent the predominant form of molecular association in these mixtures. Only a small fraction (10%) of DMSO molecules in MeOH-rich mixtures (85% in mole) forms H-bonding trimers of type 1DMSO:2MeOH. No evidence of other types of interspecies association is found. The self-diffusion coefficient for DMSO (MeOH) increases (decreases) upon mixing. The characteristic reorientation time tau1 of both species increases in the mixture, but the composition dependence is weak. The frequency spectrum of MeOH reorientational time-correlation function shows significant redshifts of the principal librational band as DMSO is added to the system, whereas the librational band of DMSO shows small alterations upon mixing. Our results are discussed in the light of previous simulation analyses for a similar system, DMSO-water mixtures, and compared with available experimental results.

Bandwidth analysis of solvation dynamics in a simple liquid mixture

The time-dependent energy distribution of solvation dynamics is studied by molecular dynamics simulations of a Lennard-Jones mixture. We calculate the response functions of the average and the variance which correspond to the spectral peak shift and bandwidth. Our calculation shows that the variance relaxation is slower than that of the average. The result agrees qualitatively with the experimental results. Dividing the obtained response functions into subcomponents caused by each solvent, we find that the relaxation is dominated by that solvent which strongly interacts with the solute. Extracting the redistribution component from the response functions, we find that it causes the slower relaxation of the response function. Thus, we conclude that the difference of the slower relaxations between the average and variance is caused by the redistribution process.

Solvation dynamics in supercritical fluids: Equilibrium versus nonequilibrium solvent response functions

We present a theoretical study of solvation dynamics in supercritical fluids. Molecular dynamics simulations show a significant difference between equilibrium and nonequilibrium solvent response functions, especially pronounced at medium and low solvent densities. We propose an analytical theory for the nonequilibrium solvation function based on the generalized nonlinear Smoluchowski-Vlasov equation. The theory is shown to be in good agreement with simulation, providing an accurate description of the nonequilibrium time-dependent solvent density profile around the solute over a wide range of supercritical solvent densities. The nonequilibrium solvent response function is shown to reflect gradual solvent clustering around the excited solute.

Ion solvation dynamics in supercritical fluids

We present a theoretical study of ion solvation dynamics in a supercritical solvent. Molecular dynamics simulations show a significant difference between equilibrium and nonequilibrium solvent response functions, especially pronounced at medium and low solvent densities. We propose a simple analytical theory for the nonequilibrium solvation function based on the generalized nonlinear Smoluchowski-Vlasov equation. The theory is shown to be in excellent agreement with simulation over a wide range of supercritical solvent densities.