1/f spectrum and memory function analysis of solvation dynamics in a room-temperature ionic liquid

To understand the non-exponential relaxation associated with solvation dynamics in the ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate, we study power spectra of the fluctuating Franck-Condon energy gap of a diatomic probe solute via molecular dynamics simulations. Results show 1/f dependence in a wide frequency range over 2 to 3 decades, indicating distributed relaxation times. We analyze the memory function and solvation time in the framework of the generalized Langevin equation using a simple model description for the power spectrum. It is found that the crossover frequency toward the white noise plateau is directly related to the time scale for the memory function and thus the solvation time. Specifically, the low crossover frequency observed in the ionic liquid leads to a slowly-decaying tail in its memory function and long solvation time. By contrast, acetonitrile characterized by a high crossover frequency and (near) absence of 1/f behavior in its power spectra shows fast relaxation of the memory function and single-exponential decay of solvation dynamics in the long-time regime.Comment: 10 pages, 4 figure

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calculate solvation free energies of charged hard spheres using density functional theory interaction potentials with molecular dynamics simulation (DFT-MD) and isolate the effects of charge and cavitation, comparing to the Born (linear response) model. We show that using uncorrected Ewald summation leads to unphysical values for the single ion solvation free energy and that charging free energies for cations are approximately linear as a function of charge but that there is a small non-linearity for small anions. The charge hydration asymmetry (CHA) for hard spheres, determined with quantum mechanics, is much larger than for the analogous real ions. This suggests that real ions, particularly anions, are significantly more complex than simple charged hard spheres, a commonly employed representation.Comment: 28 pages, 5 figure

Ab initio molecular dynamics simulations of Aluminum solvation

The solvation of Al and its hydrolyzed species in water clusters has been studied by means of ab initio molecular dynamics simulations. The hexa-hydrate aluminum ion formed a stable complex in the finite temperature cluster simulation of one aluminum ion and 16 waters. The average dipole moment of strongly polarized hydrated water molecules in the first solvation shell of the hexa-hydrate aluminum ion was found to be 5.02 Debye. The deprotonated hexa-hydrate complex evolves into a tetra-coordinated aluminate ion with two water molecules in the second solvation shell forming hydrogen bonds to the hydroxyl groups in agreement with the observed coordination.Comment: 12 pages in Elsevier LaTeX, 5 figures in Postscript, 2 last figures are in color, submitted to Chemical Physics Letter

Competitive Solvation of the Imidazolium Cation by Water and Methanol

Imidazolium-based ionic liquids are widely used in conjunction with molecular liquids for various applications. Solvation, miscibility and similar properties are of fundamental importance for successful implementation of theoretical schemes. This work reports competitive solvation of the 1,3-dimethylimidazolium cation by water and methanol. Employing molecular dynamics simulations powered by semiempirical Hamiltonian (electronic structure level of description), the local structure nearly imidazolium cation is described in terms of radial distribution functions. Although water and methanol are chemically similar, water appears systematically more successful in solvating the 1,3-dimethylimidazolium cation. This result fosters construction of future applications of the ternary ion-molecular systems

Ab Initio Molecular Dynamics Study of Aqueous Solvation of Ethanol and Ethylene

The structure and dynamics of aqueous solvation of ethanol and ethylene are studied by DFT-based Car-Parrinello molecular dynamics. We did not find an enhancement of the structure of the hydrogen bonded network of hydrating water molecules. Both ethanol and ethylene can easily be accommodated in the hydrogen-bonded network of water molecules without altering its structure. This is supports the conclusion from recent neutron diffraction experiments that there is no hydrophobic hydration around small hydrophobic groups. Analysis of the electronic charge distribution using Wannier functions shows that the dipole moment of ethanol increases from 1.8 D to 3.1 D upon solvation, while the apolar ethylene molecule attains an average dipole moment of 0.5 D. For ethylene, we identified configurations with $\pi$-H bonded water molecules, that have rare four-fold hydrogen-bonded water coordination, yielding instantaneous dipole moments of ethylene of up to 1 D. The results provide valuable information for the improvement of empirical force fields, and point out that for an accurate description of the aqueous solvation of ethanol, and even of the apolar ethylene, polarizable force fields are required.Comment: 15 pages, 10 figures, 4 tables, revtex4, submitted to J. Chem. Phy

Dynamical coupling between protein conformational fluctuation and hydration water: Heterogeneous dynamics of biological water

We investigate dynamical coupling between water and amino acid side-chain residues in solvation dynamics by selecting residues often used as natural probes, namely tryptophan, tyrosine and histidine, located at different positions on protein surface and having various degrees of solvent exposure. Such differently placed residues are found to exhibit different timescales of relaxation. The total solvation response, as measured by the probe is decomposed in terms of its interactions with (i) protein core, (ii) side-chain atoms and (iii) water molecules. Significant anti cross-correlations among these contributions are observed as a result of side-chain assisted energy flow between protein core and hydration layer, which is important for the proper functionality of a protein. It is also observed that there are rotationally faster as well as slower water molecules than that of bulk solvent, which are considered to be responsible for the multitude of timescales that are observed in solvation dynamics. We also establish that slow solvation derives a significant contribution from protein side-chain fluctuations. When the motion of the protein side-chains is forcefully quenched, solvation either becomes faster or slower depending on the location of the probe.Comment: 12 pages and 6 figures(coloured

Li+ solvation in pure, binary and ternary mixtures of organic carbonate electrolytes

Classical molecular dynamics (MD) simulations and quantum chemical density functional theory (DFT) calculations have been employed in the present study to investigate the solvation of lithium cations in pure organic carbonate solvents (ethylene carbonate (EC), propylene carbonate (PC) and dimethyl carbonate (DMC)) and their binary (EC-DMC, 1:1 molar composition) and ternary (EC-DMC-PC, 1:1:3 molar composition) mixtures. The results obtained by both methods indicate that the formation of complexes with four solvent molecules around Li+, exhibiting a strong local tetrahedral order, is the most favorable. However, the molecular dynamics simulations have revealed the existence of significant structural heterogeneities, extending up to a length scale which is more than five times the size of the first coordination shell radius. Due to these significant structural fluctuations in the bulk liquid phases, the use of larger size clusters in DFT calculations has been suggested. Contrary to the findings of the DFT calculations on small isolated clusters, the MD simulations have predicted a preference of Li+ to interact with DMC molecules within its first solvation shell and not with the highly polar EC and PC ones, in the binary and ternary mixtures. This behavior has been attributed to the local tetrahedral packing of the solvent molecules in the first solvation shell of Li+, which causes a cancellation of the individual molecular dipole vectors, and this effect seems to be more important in the cases where molecules of the same type are present. Due to these cancellation effects, the total dipole in the first solvation shell of Li+ increases when the local mole fraction of DMC is high