Composition Dependence of Dynamic Heterogeneity Time- and Length Scales in [Omim][BF4]/Water Binary Mixtures: Molecular Dynamics Simulation Study.
15 Dec 2015-Journal of Physical Chemistry B (American Chemical Society)-Vol. 119, Iss: 51, pp 15683-15695
TL;DR: Simulated structural aspects and DH time scales provide microscopic explanations to the existing experimental observations from time-resolved fluorescence and Kerr spectroscopic measurements.
Abstract: Composition dependence of four-point dynamic susceptibilities, overlap functions, and other dynamic heterogeneity (DH) parameters have been investigated by using all-atom molecular dynamics simulations for aqueous solutions of the ionic liquid (IL), 1-octyl-3-methyl imidazolium tetrafluoroborate ([Omim][BF4]) covering the pure-to-pure range. Upon addition of water in the IL, the DH time scales become faster and the four-point dynamic susceptibility time scale softens. Evidences for jump motions for both water and ions have been found from the simulated single particle displacements that show strong deviation from Gaussian distribution. Estimated dynamic correlation length for water reflects effects of IL, whereas those for ions remain largely insensitive to the mixture composition. Simulated structural aspects and DH time scales provide microscopic explanations to the existing experimental observations from time-resolved fluorescence and Kerr spectroscopic measurements.
Citations
More filters
TL;DR: This review comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.
Abstract: Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.
236 citations
TL;DR: This work discusses metal-free DESs, focusing on their preparation, characterization of physical properties, and considerations for their application, and highlights inconsistencies or omissions in literature reports, as well as factors to consider when investigating these systems.
Abstract: In the past decade, reports detailing the preparation, characterization, and application of deep eutectic solvents (DESs) have grown in number significantly, showing signs of increased interest and...
89 citations
TL;DR: The temporal heterogeneity aspects of these media have been investigated by examining the simulated particle motion characteristics and substantiated by estimating the dynamically correlated time scales and length-scales through simulations of four-point susceptibilities and density correlations.
Abstract: Dielectric relaxation (DR) measurements in the frequency window 0.2 ≤ ν(GHz) ≤ 50 for deep eutectic solvents (DESs) made of acetamide (CH3CONH2) and urea (NH2CONH2) with the general composition, [f CH3CONH2 + (1 − f) NH2CONH2] at f = 0.6 and 0.7, reveal three distinct relaxation time scales—τ1 ∼ 120 ps, τ2 ∼ 40 ps, and τ3 ∼ 5 ps. Qualitatively similar time scales have been observed for DR of neat molten urea, whereas the reported DR for neat molten acetamide in the same frequency window reflects two relaxation processes with no trace of ∼100 ps time scale. This slowest DR time scale (τ1) resembles closely to the long-time constant of the simulated structural H-bond relaxation (CHB(t)) involving urea pairs. Similarity in activation energies estimated from the temperature dependent DR measurements (335 ≤ T/K ≤ 363) and structural H-bond relaxations indicates that the structural H-bond relaxation overwhelmingly dominates the slowest DR relaxation in these DESs. Simulated collective reorientational correlation functions (Cl(t)), on the other hand, suggest that the second slower time scale (∼40 ps) derives contributions from both the single particle orientation dynamics and structural H-bond relaxation, leaving no role for hydrodynamic molecular rotations. The sub-10 ps DR time scale has been found to be connected to the fast reorientation dynamics of the component molecules (acetamide or urea). Fractional viscosity dependence for the longest DR times, τDR∝η/Tp, has been observed for these DESs with the fraction power p = 0.7. Subsequently, the temporal heterogeneity aspects of these media have been investigated by examining the simulated particle motion characteristics and substantiated by estimating the dynamically correlated time scales and length-scales through simulations of four-point susceptibilities and density correlations. These estimated dynamical time scales and length-scales assist in explaining the different inferences regarding solution heterogeneity drawn from different measurements on these DESs.
32 citations
TL;DR: Whether these cycloethers can induce stiffening and transition of water H-bond network structure and, if they do, whether such structural modification differentiates the chemical nature (dipolar or quadrupolar) of the cosolvent molecules is investigated.
Abstract: We have performed steady state UV-visible absorption and time-resolved fluorescence measurements and computer simulations to explore the cosolvent mole fraction induced changes in structural and dynamical properties of water/dioxane (Diox) and water/tetrahydrofuran (THF) binary mixtures. Diox is a quadrupolar solvent whereas THF is a dipolar one although both are cyclic molecules and represent cycloethers. The focus here is on whether these cycloethers can induce stiffening and transition of water H-bond network structure and, if they do, whether such structural modification differentiates the chemical nature (dipolar or quadrupolar) of the cosolvent molecules. Composition dependent measured fluorescence lifetimes and rotation times of a dissolved dipolar solute (Coumarin 153, C153) suggest cycloether mole-fraction (X(THF)/Diox) induced structural transition for both of these aqueous binary mixtures in the 0.1 ≤ X(THF)/Diox ≤ 0.2 regime with no specific dependence on the chemical nature. Interestingly, absorption measurements reveal stiffening of water H-bond structure in the presence of both the cycloethers at a nearly equal mole-fraction, X(THF)/Diox ∼ 0.05. Measurements near the critical solution temperature or concentration indicate no role for the solution criticality on the anomalous structural changes. Evidences for cycloether aggregation at very dilute concentrations have been found. Simulated radial distribution functions reflect abrupt changes in respective peak heights at those mixture compositions around which fluorescence measurements revealed structural transition. Simulated water coordination numbers (for a dissolved C153) and number of H-bonds also exhibit minima around these cosolvent concentrations. In addition, several dynamic heterogeneity parameters have been simulated for both the mixtures to explore the effects of structural transition and chemical nature of cosolvent on heterogeneous dynamics of these systems. Simulated four-point dynamic susceptibility suggests formation of clusters inducing local heterogeneity in the solution structure.
31 citations
TL;DR: The structural and dynamical properties of two imidazolium-based geminal dicationic ionic liquids (GDILs) have been studied to obtain a fundamental understanding of the molecular basis of the macroscopic and microscopic properties of the bulk liquid phase.
Abstract: In this work, the structural and dynamical properties of two imidazolium-based geminal dicationic ionic liquids (GDILs), i.e. [Cn(mim)2][NTf2]2 with n = 3 and 5, have been studied to obtain a fundamental understanding of the molecular basis of the macroscopic and microscopic properties of the bulk liquid phase. To achieve this purpose, molecular dynamics (MD) simulation, density functional theory (DFT) and atoms in molecule (AIM) methods were used. Interaction energies, charge transfers and hydrogen bonds between the cation and anions of each studied GDIL were investigated by DFT calculations and also AIM. The mean square displacement (MSD), self-diffusion coefficient, and transference number of the cation and anions, and also the density, viscosity and electrical conductivity of the studied GDILs, were computed at 333.15 K and at 1 atm. The simulated values were in good agreement with the experimental data. The effect of linkage alkyl chain length on the thermodynamic, transport and structural properties of these GDILs has been investigated. The structural features of these GDILs were characterized by calculating the partial site–site radial distribution functions (RDFs) and spatial distribution functions (SDFs). The heterogeneity order parameter (HOP) has been used to describe the spatial structures of these GDILs and the distribution of the angles formed between two cation heads and the middle carbon atom of the linkage alkyl chain was analyzed in these ILs. To investigate the temporal heterogeneity of the studied GDILs, the deviation of the self-part of the van Hove correlation function, Gs(,t), from the Gaussian distribution of particle displacement and also the second-order non-Gaussian parameter, α2(t), were used. Since, the transport and interfacial properties and ionic characteristics of these GDILs were studied experimentally in our previous studies as a function of linkage chain length and temperature, in this work, we try to give a better perspective of the structure and dynamics of these systems at a molecular level.
31 citations
References
More filters
TL;DR: The dynamical steady-state probability density is found in an extended phase space with variables x, p/sub x/, V, epsilon-dot, and zeta, where the x are reduced distances and the two variables epsilus-dot andZeta act as thermodynamic friction coefficients.
Abstract: Nos\'e has modified Newtonian dynamics so as to reproduce both the canonical and the isothermal-isobaric probability densities in the phase space of an N-body system. He did this by scaling time (with s) and distance (with ${V}^{1/D}$ in D dimensions) through Lagrangian equations of motion. The dynamical equations describe the evolution of these two scaling variables and their two conjugate momenta ${p}_{s}$ and ${p}_{v}$. Here we develop a slightly different set of equations, free of time scaling. We find the dynamical steady-state probability density in an extended phase space with variables x, ${p}_{x}$, V, \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{}, and \ensuremath{\zeta}, where the x are reduced distances and the two variables \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{} and \ensuremath{\zeta} act as thermodynamic friction coefficients. We find that these friction coefficients have Gaussian distributions. From the distributions the extent of small-system non-Newtonian behavior can be estimated. We illustrate the dynamical equations by considering their application to the simplest possible case, a one-dimensional classical harmonic oscillator.
17,939 citations
TL;DR: In this article, the authors compared the canonical distribution in both momentum and coordinate space with three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 77, 63 (1983); and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983).
Abstract: Three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 48, 1818 (1982)], and Evans and Morriss [Chem. Phys. 77, 63 (1983)]; and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983)] are examined analytically via calculating the equilibrium distribution functions and comparing them with that of the canonical ensemble. Except for effects due to momentum and angular momentum conservation, method (1) yields the rigorous canonical distribution in both momentum and coordinate space. Method (2) can be made rigorous in coordinate space, and can be derived from method (1) by imposing a specific constraint. Method (3) is not rigorous and gives a deviation of order N−1/2 from the canonical distribution (N the number of particles). The results for the constant temperature–constant pressure ensemble are similar to the canonical ensemble case.
13,921 citations
TL;DR: Weiner et al. as mentioned in this paper derived a new molecular mechanical force field for simulating the structures, conformational energies, and interaction energies of proteins, nucleic acids, and many related organic molecules in condensed phases.
Abstract: We present the derivation of a new molecular mechanical force field for simulating the structures, conformational energies, and interaction energies of proteins, nucleic acids, and many related organic molecules in condensed phases. This effective two-body force field is the successor to the Weiner et al. force field and was developed with some of the same philosophies, such as the use of a simple diagonal potential function and electrostatic potential fit atom centered charges. The need for a 10-12 function for representing hydrogen bonds is no longer necessary due to the improved performance of the new charge model and new van der Waals parameters. These new charges are determined using a 6-31G* basis set and restrained electrostatic potential (RESP) fitting and have been shown to reproduce interaction energies, free energies of solvation, and conformational energies of simple small molecules to a good degree of accuracy. Furthermore, the new RESP charges exhibit less variability as a function of the molecular conformation used in the charge determination. The new van der Waals parameters have been derived from liquid simulations and include hydrogen parameters which take into account the effects of any geminal electronegative atoms. The bonded parameters developed by Weiner et al. were modified as necessary to reproduce experimental vibrational frequencies and structures. Most of the simple dihedral parameters have been retained from Weiner et al., but a complex set of 4 and yj parameters which do a good job of reproducing the energies of the low-energy conformations of glycyl and alanyl dipeptides has been developed for the peptide backbone.
12,660 citations
TL;DR: On the other hand, in this paper, a superparamagnetically collapsed Mossbauer spectrum is obtained for carbon with fewer active sites, and these particles sinter and carburize in a manner more similar to that of Fe particles supported on graphite.
Abstract: on carbon. These particles can interact with the active sites of the carbon to form a species a t the metal-carbon interface which is stable against high-temperature reduction. This interaction may occur via oxygen atoms bonding with both the carbon surface and iron atoms contacting the carbon. This oxygen could be provided by dissociation of the C O ligands during decomposition. The remainder of the surface Fe can exist in either a reduced or carburized state, but the small size of these particles results in either case in a superparamagnetically collapsed Mossbauer spectrum. For carbons with fewer active sites, more of the particles are unbound, and these particles sinter and carburize in a manner more similar to that of Fe particles supported on graphite.
11,221 citations
TL;DR: In this paper, the effects of velocity rescaling on the self-diffusion coefficient D and radial distribution functions, gOO, gOH, and gHH for all five water models were determined and compared to experimental data.
Abstract: Molecular dynamics simulations of five water models, the TIP3P (original and modified), SPC (original and refined), and SPC/E (original), were performed using the CHARMM molecular mechanics program. All simulations were carried out in the microcanonical NVE ensemble, using 901 water molecules in a cubic simulation cell furnished with periodic boundary conditions at 298 K. The SHAKE algorithm was used to keep water molecules rigid. Nanosecond trajectories were calculated with all water models for high statistical accuracy. The characteristic self-diffusion coefficients D and radial distribution functions, gOO, gOH, and gHH for all five water models were determined and compared to experimental data. The effects of velocity rescaling on the self-diffusion coefficient D were examined. All these empirical water models used in this study are similar by having three interaction sites, but the small differences in their pair potentials composed of Lennard-Jones (LJ) and Coulombic terms give significant difference...
2,223 citations