Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Japan
*Address all correspondence to: ishida@ims.ac.jp
1. Introduction
Ionic liquids (ILs) have attracted many researchers in the areas including physics and chemistry because of their characteristics that are different from conventional molecular liquids and, today, ILs have been one of interesting subjects of scientific study. ILs are consisted of ions and liquids at or near room temperature, and show negligible vapor pressure, thermal and chemical properties, and so on [1-8]. Also, ILs have been widely used as solvents for organic reactions with the expectation of high yields [1]. Most interesting features of ILs can be attributed to remarkable interionic interactions, and these can be an important key factor to study the characteristics of ILs at molecular level. From the results of both experimental and theoretical investigations, it has been recognized that the interionic interaction of ILs could determine physical and chemical properties.
From an experimental side, both femtosecond optically heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES)[3, 9, 10] and THz time-domain spectroscopy (THz-TDS)[11, 12] have been applied to investigate the intermolecular vibrational dynamics in ILs. In particular, with OHD-RIKES studies, the possibility to control a property such as shear viscosity by substituting an atomic element in an ionic unit has been reported [13, 14]. On the other hand, from theoretical and computational viewpoints in recent years, ILs have been chosen to study static properties such as structural and thermophysical properties[15-18] and novel interionic dynamics under solvation dynamics [19-21], dynamical properties [22-26], and Kerr spectra analyses [14, 27]. It has been suggested from the simulation studies by Ishida et al. [14] that interionic properties in ILs could be effectively adjustable by substituting an atomic unit in an ion unit in addition to a combination of cations and anions. Also, it has been pointed out that the interplay of motions between cation and anion species could play an important role in specific interionic interactions of ILs.[28]
Now, we can consider two factors important to understand specific interionic interactions in ILs. One is the interionic interaction depending on specific correlations such as cross-correlation terms between cations and anions, and the other is polarization effects due to many-body interactions caused by cations and anions in ILs. While a large number of experimental approaches have been applied to investigate these subjects, molecular-level understanding of many specific properties of ILs has been left unresolved. Obviously, theoretical researches are suitable to tackle these problems to which experimental procedures are not accessible. Thus, it is expected that such computational method as molecular dynamics (MD) simulations enables us to obtain significant information of ILs, utilizing the force field with well-parameterized potential functions and partial charges [29].
With the MD simulation procedures, it is possible to study the effects of the cross-correlations on dynamic collective motions of ions in ILs which are considered to govern the strength and behavior of couplings between ionic motions through interionic interactions [28]. The computation of the time correlations of velocity and momentum between a tagged ion and other unlike ions at different distances provides cross-correlation and momentum correlation functions [28, 30, 31]. Utilizing these calculated functions, we can investigate not only interionic interactions at molecular level but also how collective motions in ILs can proceed accompanied by the momentum transfer between ions in the target IL system.
On the other hand, electrostatic interactions between ions in ILs could be modulated due to many-body interactions and, then, it could emerge as polarization effects caused by the distortion of electron densities under anisotropic environment in ILs. It has been pointed out that the inclusion of polarization effects is significant to investigate the characteristics of ILs in structural and dynamical properties [23, 24]. Therefore, it is required for us to carry out the MD simulation, introducing such procedures as a point dipole model and a polarization energy term into the total potential energy representation of the system [32]. In addition to electrostatic interactions, describing the variation and relaxation of the polarizability anisotropy of ILs is important to investigate dynamics in ILs. For achieving this, we need to compute time-dependent polarization effects on a target system due to environmental effects in ILs. Theoretically, when we would take polarizability anisotropy into account, to track the change of molecular polarizability tensor partly dependent on molecular orientations such as rotational motions would be required. With the calculation of the time correlation function (TCF) of off-diagonal elements of the total polarizability of the system, we can investigate collective properties with the result of the polarizability anisotropy relaxation of the system.
In this chapter, we choose a 1-butyl-3-methylimidazolium cation based ILs with hexafluorophosphate anion, [BMIm][PF6], and bis(trifluoromethylsulfonyl)amide anion, [BMIm][NTf2], as target systems. (See Figure 1 for all the molecular structures of [BMIm] +, [PF6]– and [NTf2]–.) Firstly, we focus on the collective properties of [BMIm][PF6] with the cross-correlation functions of ionic species in the IL and polarization effects on ionic motions. As a second, we investigate dynamical properties of [BMIm][NTf2]. Below, we start from the explanation for the velocity cross-correlations of a central atom with neighboring atoms. Then, we show how to evaluate both cross-correlation and momentum correlation functions of the target IL. In addition, we describe the introduction of a polarization energy term and the procedure to calculate induced dipole moments and the polarization energy. Following those, the theoretical background of polarization TCFs is given. Computational details are also summarized. In later sections, we discuss polarization effects on interionic interactions and specific properties related to those based on the MD simulation results. We show computation results obtained from performing MD simulation and computing the polarization TCF with the dipole-dipole (DID) approximation [4,5]. Also, we examine the relation between the anisotropic polarizability relaxation and collective motions of ionic species. Lastly, we discuss relaxation processes of ILs including an explanation of important points in studying dynamical properties on ILs.
Figure 1.
Molecular structures and definitions of body-fixed coordinate axes: (a) [BMIm]+, (b) [PF6]– and (c) [NTf2]–. (See text.) In [BMIm]+, the Y direction is along the line connecting two nitrogen (blue colored) atoms in the ring, and the Z direction is set perpendicular to the ring plane and the Y direction axis. The X directionis set in the ring plane orthogonal to both the Y and the Z axes. In [PF6] –, the X, Y, and Z direction axes are set equivalently. In [NTf2]–, the Y axis is along the line connecting two sulfurs (dark yellow colored), and the Z axis is set perpendicular the S-N-S plane and the X axis. The X axis is set in the S-N-S plane orthogonal to both the Y and Z axes.
2. Tracking Ionic Motions Through Interionic Interactions: Cross-Correlation, Polarization Effects, and Dynamical Properties
As mentioned in the previous section, for many problems of ILs to which experimental procedures are not accessible, theoretical investigations with computer simulation procedures are promising and suitable to obtain detailed information at molecular level. In particular, specific correlations such as cross-correlations between cations and anions seem not to be feasible to detect experimentally, but the MD simulations enable us to evaluate those. Also, for tracing dynamics in ILs such as librational and reorientational dynamics of ionic species, cross-correlated ionic motions, and the influence of polarization coming from many-body effects caused by cations and anions in ILs from a microscopic point of view, the MD simulations are considered to be one of useful and powerful tools. In this chapter, we will show and explain the dynamical properties on ILs based on MD simulation results taking into consideration mainly following points:
how can cross-correlated ionic motions in ILs be modified by interionic dynamics and electronic polarizability effects ?
how can the collective dynamics through interionic interactions in ILs be tracked by computer simulation procedures ?
what kinds of properties with simulation data do we have to check and care in the investigation of dynamical properties of ILs ?
what kinds of subjects of ILs can be or should be investigated theoretically, considering important properties of ILs experimentally observed ?
The items 1 and 2 are main parts in this chapter, including the key results mentioned above. The item 3 is considered to be important for suggesting (or notifying) attention to researchers who are working on analyzing the dynamical properties on ILs. The item 4 includes future perspectives. In particular, we will show and discuss new aspects of ILs in each section related to the items 1, 2, and 3. Also, some of theoretical backgrounds and computational procedures for studying the dynamical properties of ILs are given in each related section, below.
At first, let us consider following points:
how can the collective dynamics through interionic interactions cause the unique physical and chemical properties of ILs ?
how can interionic dynamics be modified by electronic polarizability effects ?
The former includes the investigation of the contribution of ionic motions due to Coulombic interactions to velocity cross-correlation functions. In particular, through the analysis of the longitudinal and nonlongitudinal contributions to the velocity cross-correlation function in ILs, we will be able to investigate interionic interactions in detail. Also, important properties for physical and chemical interests such as case effects seems to be within the scope of unique collective dynamics in ILs. The later covers the relation between polarizability correlation functions and interionic interactions for ILs. To investigate these points, we can utilize useful information for static properties obtained from computer simulations, but those are not often enough to extract the details of specific interactions.
Here, we give an example that it is difficult for us to find the importance of the interactions between cation and anion species only from static properties. Figure 2 displays the computed radial distribution functions, g(r), comparing non-polarizable model with polarizable model [28]. As shown in the figure, we can observe only small difference in RDFs between two models, except that the g(r) of cation-anion is different from those of cation-cation and anion-anion as easily deduced. But, obviously, it is not feasible for us to observe information other than strong spatial correlations and sequential ordering of cation and anion pairs. In addition, the analysis of averaged static structures such as RDFs is not enough to investigate remarkable polarization effects such as screening effects influenced by polarization. Below, we give a schematic explanation and discuss that to consider interionic interactions with cross-correlation analyses is important to investigate the interplay between cation and anion species in ILs, and it is shown that the cross-correlation analyses provide contrastive features.
Figure 2.
Radial distribution functions (rdfs) for the center of mass of [BMIm][PF6] for both nonpolarizable (np) and polarizable (p) models: (a) [BMIm]+ - [PF6]– and (b) [BMIm]+ - [BMIm]+ and [PF6]– - [PF6]–.
Different from a usual velocity autocorrelation functions (VACF), cross-correlation functions describe interactions between unlike (ionic) species (that is, between a cation and an anion in ILs) and show opposite features to VACFs. In Figure 3, these features are explained schematically. Here, it should be emphasized that cross-correlation functions provide more information on interionic interactions than that static properties such as RDFs include. As seen in Figure 3 ((a) and (b)), the variation of cross-correlation function correlates with the alteration of the VACF. In particular, the cross-correlation function shows the increasing toward the maximum peak where the VACF approaches a minimum point. These features indicate that it is possible for us to track the time evolution of cross-correlations (between cation and anion species in ILs) as collective dynamics. In addition, it is observed in Figure 3 ((c)) that a particle (ionic species) is bouncing back and forth between like and unlike particles (ions). This implies collective (ionic) motions between coordination shells, thus, it is indicated that not only interionic interactions but also momentum transferring among ionic species could be extracted by analyzing cross-correlation functions. Also, considering that these cross-correlation functions could be modulated by the strength of interionic interactions and the coupling between cation and anion motions, cross-correlation functions would be largely influenced by polarization effects.
Figure 3.
Schematic explanation for a cross-correlation and a VACF: (a) short-time region: an ion approaches a counter ion, and, at the same time, the VACF decreases, (b) the cross-correlation functions start decreasing after the maximum (see text) and (c) bounced ions again approaches counter ions.
On the other hand, as another type of the appearance of interplay between cross-correlation and polarization effects, we consider the polarizability anisotropy and its relaxation of an IL system. These correspond to the variation of the sum of molecular polarizability depending on time. Therefore, it is required to compute the change of molecular polarizablity on each molecule due to interionic interactions and interaction-induced multipole effects. Obviously, it is expected that molecular polarizabilites are influenced by interionic cross-correlations. In later sections, we introduce the theoretical background of polarizability TCF and its application to the study of ILs, and show how degree these are effective and discuss the importance of considering cross-correlations.
Here, we introduce cross-correlation functions, and then, give an explanation of a polarizable model and polarizability time correlation function. Computational details are also given.
3.1. Cross-Correlation Function
A velocity autocorrelation function (VACF) is defined for calculating the velocity correlation of a same particle as follows,
C(t)=〈vi(t)⋅vi(0)〉〈vi(0)⋅vi(0)〉−1E1
where vi(t) is the velocity of i-th species. Different from the VACF, the time cross-correlation functions between the velocity of a central species i and the velocity of a neighboring species j different from the particle i is defined as follows [30, 31],
CRnij(t)=〈vi(0)⋅vj(t)〉n(〈vi2〉〈vj2〉)−1/2E2
where vi(t) is the velocity of the species i and vj(t) is the velocity of the species j. 〈vi2〉and 〈vj2〉represent the mean square velocities. Here, it should be noted that 〈〉n represents a restricted statistical average [30, 31],
where θ represents the step function. an and bn are set as the positions of the nth minima of the radial distribution function of the system, g(r). The Nn in Eqs. (2) and (3) is the coordination number in the region between r = an and bn. as follows,
Nn=4πρ∫anbnr2g(r)drE4
where ρ is the number density of the system.
Here, we can define the cross-correlation in ILs as that between a centered cation and the total contribution of other anions [28, 30],
CRnCA(t)=Nn〈vC(0)⋅vA(t)〉n(〈vC2〉〈vA2〉)−1/2E5
where “C” and “A” in super- and subscript represent cation and anion, respectively. By interchanging “C” and “A” in super- and subscript in the above equation, we can easily derive the formulation for the cross-correlation between a centered anion and other cations.
Also, with the computation of the cross-correlation functions between the velocity of a central ion and velocities of neighboring distinct ions, we can analyze the momentum transfer between distinct ions in ILs. Introducing the momentum correlation function [28, 31], the transfer of the momentum of a cation to distinct cations and anions is defined as follows,
where p+ and p‒ mean the momentums of cation and anion, respectively. The transference of momentum of an anion to distinct anions and cations, PnA(total) (t), is also given as follows,
Here, we introduce a polarizable model considering induced dipole moments and explain the procedure for computing induced dipole moments and the many-body polarization energy [33-35].
The total potential energy of the system under the resulting polarizable force field is defined as follows,
Vtot=Vbond+Vnonbond+VpolE8
where the terms Vbond and Vnonbond are intra- and intermolecular interaction energies. The polarization energy, Vpol, is decomposed into the three terms as follows,
Vpol=Vcharge−dipole+Vdipole−dipole+VselfE9
where the charge-dipole contribution, Vcharge-dipole, the dipole-dipole contribution, Vdipole-dipole, and the self-polarizability term, Vself, are defined, respectively, as
Vcharge−dipole=−∑iμi⋅EiE10
Vdipole−dipole=∑ijμi⋅Tij⋅μjE11
Vself=∑i12αiμi⋅μiE12
In the above equations, Ei is the electric field on atom i, produced by the partial charges of all other surrounding atoms, and μi represents the induced dipole moment on atom i. Also, αimeans the isotropic atomic polarizability of atom i. In Eq. (11), the Tij, is the dipole field tensor element defined as follows,
Tij=1rij3[1−3rijrijrij2]E13
where rij = ri – rj. The induced dipole moment on atom i is given by
μi=αi(Ei−∑j≠iTijμj)E14
With Eqs. (10), (11), and (12), the polarization energy, Epol, is summarized as follows,
Vpol=−12∑iμi⋅EiE15
3.2.2. Polarizability Time-Correlation Function (TCF)
The polarizability anisotropy of the system can be tractable by calculating the TCF of off-diagonal elements of the total polarizability. Here, the theoretical background of the polarizability TCF is summarized briefly.
At first, we define the total polarizability of the system, Π(t), that is the sum of the molecular polarizability, ΠM(t), and the interaction-induced polarizability, ΠII(t), as follows,
Π(t)=ΠΜ(t)+ΠΙΙ(t)E16
where t represents the time dependence, and subscripts M and II mean the molecular part and the interaction-induced part, respectively. The molecular part is given by the sum of the polarizability tensors of isolated gas phase molecular polarizability in the laboratory frame,
ΠΜ(t)=∑i=1Nαi(t)E17
where N is the number of molecules, and αi is the polarizability tensor of molecule i. For the formulation of the interaction-induced part, we employ the dipole-induced-dipole (DID) model approximation [36, 37], which assumes that the molecular polarizabilities are modified due to a dipolar coupling with the influence of higher order unconsidered. The interaction-induced polarizability in the DID approximation is given as follows,
ΠII(t)=∑i=1N∑j≠iNαi(t)⋅Tij(t)⋅α˜j(t)E18
where Tij means the dipole interaction tensor between molecules i and j. α˜i(t)is the effective polarizability for molecule i defined by the following equation including the interaction-induced effects,
α˜i(t)=αi(t)+∑j≠iNαi(t)⋅Tij(t)⋅α˜j(t)E19
Equation (19) can be solved by the calculation procedure that is called the all-orders DID approximation [37]. It should be noted that the DID model employed here assumes a center-center DID model, where it recognizes the polarizability as concentrated in the center of mass of the molecule.
Here, we give the representation of the total system polarizability including the cationic and anionic components [14, 27],
Π(t)=ΠC(t)+ΠA(t)E20
where the superscripts C and A represent cation and anion species, respectively. Referring to Equations. (16), (17), and (18), the ΠC(t) and ΠA(t) defined as follows,
ΠC(t)=∑i=1(αiC(M)(t)+αiC(II)(t))E21
ΠA(t)=∑j=1(αjA(M)(t)+αjA(II)(t))E22
where the indices k and l go over all cations and anions, respectively. Then, the total polarizability of the system, Π(t), can be rewritten, as follows,
Where βC(t)={ΠC(t)−13Tr(ΠC(t))I}and βA(t)={ΠA(t)−13Tr(ΠA(t))I}and Iis the unit tensor. Finally, the polarizability TCF, ϕ(t), can be rewritten as follows,
The nuclear response function, R(t), is represented as the time derivative of the polarizability TCF,
R(t)=−1kBT∂∂tϕ(t)E25
3.3. Computational Details
Utilizing the sets of force field parameters for [BMIm]+, [PF6]– and [NTf2]– [38-40], MD simulations have been performed with the DL_POLY molecular dynamics suite [41]. In all the simulations, all the stretching bonds were constrained with the SHAKE algorithm [42]. 125 ion pairs (4000 atoms for [BMIm][PF6] and 5000 for [BMIm][NTf2]) were set in a cubic box under the periodic boundary condition. The lengths of cubic box size were set to be 35.47 and 38.46 Å to reproduce the experimental data of the densities of [BMIm][PF6] and [BMIm][NTf2] at 298 K [13, 43]. 12 Å was set as cutoff length. The time step was 2 fs. The long-range Coulomb and polarization terms were computed with the Ewald’s summation technique [44]. Firstly, each system was equilibrated at 600 K for 15 ns in the NVE run and then successively cooled down to 298 K in several stages using velocity scaling. Then, 20 ns NVT run for equilibration at 298 K was carried out. After these equilibration runs, trajectories were recorded every 20 fs (50 fs for [BMIm][NTf2]) and collected during 10 ns (20 ns for [BMIm][NTf2]) production runs. For atomic polarizabilities, from the literature [45, 46], we adopted 1.152, 0.705 and 0.0885 Å3 for the C, N, and H atoms of the cation, respectively, and 0.121 and 3.630 Å3 for the F and P atoms of the anion, respectively. We considered the distance and vector between the atomic sites of distinct molecules in the dipole interaction tensor, T=(I−3r^r^)/r3, wherer^=r/rThen, the range of the attenuation of dipolar interactions at short distances, s, was evaluated with the Thole’s definition [47], s = 1.662(αiαj)1/6 with atomic polarizabilities, αiandαj. Also, the dipole interaction tensor we used is given by [47]
T=(I−3r^r^)/r3(rs)E28
T=[(4a3−3a4)I−3a4r^r^]/r3(r≤s)E29
where a = r/s. Equation (14) is solved with an iterative procedure at each time step, and then, the criterion value of iterative solution for induced dipole moment was set to 0.001 D. For the computation of the polarizability TCF of [BMIm][NTf2], we used molecular polarizabilities of 14.372 Å3 for [BMIm]+ and 11.259 Å3 for [NTf2]–.[48] The body-fixed coordinate axes set in the cation and the anion are shown in Figure 2. For more computational procedures in detail, interested readers should refer to references [14] and [28].
4. Cross-Correlation, Momentum Correlation, and Polarization Effects
The computed velocity cross-correlation functions [28] are shown in Figures 4 and 5. These figures show the comparison of velocity cross-correlation functions observed around the cation, [BMIm]+, and the anion, [PF6] –, placed at the center, for both the nonpolarizable and polarizable models. From the result of the cation-anion RDF in Figure 2, we selected 3.5, 8.4, 8.4, and 14.75 Å for the value of a1, b1, a2, and b2 in Equation (4) for the nonpolarizable model, respectively, and 3.5, 8.5, 8.5, and 15.0 Å for the polarizable model, respectively. With these values, the first (C1(t)) and second (C2(t)) coordination shells were specified. It should be noted that the C1(t) between the center anion and distinct anions is not shown in Figure 5 because the C1(t) is almost zero corresponding to the result that the anion-anion RDF result is almost zero at the region of the specified first coordination shell (see RDFs in Figure 2(b)). Also, it should be noted that the initial values, Cn(0), are negative since the system size is finite as have been pointed out [30].
Figure 4.
The Calculated velocity cross-correlation functions between the [BMIm]+ set at the center in the first coordination shell and [PF6] –, and between the center cation and other cations, including VACF: (a) in the nonpolarizable model and (b) in the polarizable model.
As shown in Figures 4 and 5, the VACFs are also displayed for comparison. As briefly explained in Section 2 with Figure 3, the initial rise of the cross-correlation function, C1(t), appears toward the maximum for both models, corresponding to the decay of the VACF to the minimum (see the (a) and (b) in Figure 3). These results indicate that the initial momentum of the central ion is gained by neighboring ions immediately after t = 0. In addition, following decay profiles (see the (c) in Figure 3) are seen. These are ascribed to the spread of transferred momentum to the outer coordination shells. The peak height of the C1(t) in the polarizable model is larger than that in the nonpolarizable model. These significant results indicate that the momentum transfer can be intensified from the more distant coordination shells to the first ones due to both the charge-dipole and dipole-dipole interactions by polarization effects in addition to charge-charge Coulomb interactions. As the characteristics of the VACF profile, it is noted that the minima of the VACF are located at the positions later than those of the maxima of the velocity cross-correlation functions for the first coordination shell. These features imply that the momentum transmitted to the neighbors at the first shell is regained partly by the central cation or anion, bouncing back and forth for some time, as have been pointed out in the literature on the computer simulations of simple liquid binary mixtures [28, 31].
Figure 5.
The Calculated velocity cross-correlation functions between the [PF6] – set at the center in the first and second coordination shells and [BMIm]+, and between the center anion and anions, including VACF: (a) in the nonpolarizable model and (b) in the polarizable model.
On the other hand, an interesting feature is seen in the C1(t) between distinct cations (or the C2(t) for distinct anions). As indicated in Figures 4 and 5, the cross-correlation function for distinct ions reaches the maximum point earlier than the C1(t) for the cation-anion cross-correlation, even though the peak height is smaller than that in the cation-anion cross-correlation function. On these results, it is considered that the alteration of cage effects increases the probability that a cation (or an anion) meets a distinct cation (or anion) at early time region.
Also, in Figures 4 and 5, it is clearly observed that the cross-correlation function between distinct like ions is influenced by the modulation of cage effects due to polarization effects and that its peak height is enhanced in the polarizable model [28]. As shown in Figures 4 and 5, when VACFs pass the value of zero to negative, the cross-correlation functions reach the maximum point. Thus, this implies that the central cation (or anion) is likely to lose its initial momentum. In addition, as seen in Figure 5, the peak position of C2(t) is shifted to the time region later than that of C1(t). This is consistent with the consideration that the momentum of the central ion is transferred from the first coordination shell to the second.
Here, we examine the characteristics of the cross-correlation functions. As shown in Figure 6, we consider the decomposition of a cross-correlation function into the longitudinal (denoted as L in the figure), CnL(t), and nonlogitudinal (denoted as NL in the figure), CnNL(t), contributions as follows [28, 31, 49],
CRn(t)=CRnL(t)+CRnNL(t)E30
where CRnL(t) is represented as the velocity cross-correlation along the direction designated by the center of masses of distinct ions at t = 0. We can compute CRnL(t) with the following equation [28, 31, 49],
CRnL(t)=Nn〈viL(0)⋅vjL(t)〉n(〈vi2〉〈vj2〉)−1/2E31
where viL(t) = vi(t)[rij(0)/rij(0)] and vi(t) is the velocity of the ionic species i. Also, rij(0) means the direction vector between the center of masses of the distinct ions i and j. 〈vi2〉and 〈vj2〉represent the mean square velocities. CRnNL(t) can be computed with CRn(t) and CRnL(t). Figure 7 shows the computed CnL(t) and CnNL(t) functions for a centered cation at the first coordination shell (n = 1) both for the nonpolarizable and polarizable models [28]. All the results indicate that the velocity cross-correlations at the short time region up to 0.4 ps are predominantly governed by the longitudinal function, CnL(t).
Figure 6.
Schematic explanation of the longitudinal (L) and nonlongitudinal (NL) of a cross-correlation function.
These results clearly indicate that the Coulomb interactions between neighboring ions are mainly effective on the ionic motions as “driving forces” at the short time. In addition, the magnitude of the longitudinal correlation of ions in the polarizable model is larger than that in the nonpolarizable model. Therefore, it is considered that the longitudinal motions are strongly affected by polarization effects. On the decay behavior, similar features are obtained in both models. These indicate that polarization effects on the longitudinal motion complete mostly at the short time region.
Figure 7.
The simulated longitudinal and nonlongitudinal contributions to the velocity cross-correlation function: (a) in the nonpolarizable model and (b) in the polarizable model.
Next, the computed results of momentum correlation functions [28] are displayed in Figures 8 and 9. As shown in Figures 8 and 9 in common, the initial momentum of the cation (or anion) is mainly transferred to the close anions (or cations) while the momentum correlation with distinct cations (or anions) is smaller. These results are consistent with the consideration that strong Coulomb interactions between the cations and the anions enhance the possibility of approaching or attracting each other. Also, while the contribution of the momentum correlation between distinct cations, P1CC(t), is smaller than that between the cation and anion molecules, P1CA(t), the peak height of the P1CA(t) in the polarizable model is larger than that in the nonpolarizable model. These results indicate that the transference of the momentum between distinct cations could be intensified by both charge-dipole and dipole-dipole interactions coming from polarization effects. As we could observe in the cross-correlation functions (see Figures 4 and 5), the maxima of the momentum correlation functions between distinct cations emerge at earlier time region in comparison with those of the momentum correlations between cation and anion species.
Figure 8.
The calculated momentum correlation functions for the [BMIm]+ at the center in the first coordination shell with the VACF: (a) in the nonpolarizable model and (b) in the polarizable model.
Based on these results, obviously, it is indicated that the variation of cage effects promotes the transfer of the initial momentum of the central cation to distinct cations rather than anions. Then, the cage effects is considered to be weakened by polarization effects, though the degree of this effect is likely to be relatively small as deduced from the figures.
On the other hand, the momentum correlations for the center anion in Figure 9 show distinct features from those for the center cation. The momentum correlation function between distinct anions, P1AA(t), indicates much smaller contributions to the total momentum correlations, P1(t), in both the nonpolarizable and polarizable models, compared with the P1AC(t). Therefore, these indicate that the initial momentum of the center anion is mainly transferred to neighboring cations, and the transference between distinct anions is not enhanced. These are consistent with the consideration that the cation-anion interactions promote the propagation of the momentum one after another, as mentioned in the cation case. In addition, Figure 9 shows more interesting features. Compared with the results for the central cation in Figure 8, while the P1AA(t) does not contribute negligibly to the total momentum correlation in the nonpolarizable model, it indicates characteristic oscillating behavior in the short time region up to 0.4 ps in the polarizable model, and the latter has a relatively larger contribution to the P1(t) up to 0.1 ps than in the nonpolarizable model. This clearly indicates the decrease of cage effects through polarization effects and implies that the interionic interactions between distinct anions could become effective by the charge-charge and dipole-dipole interactions [28].
Figure 9.
The calculated momentum correlation functions for the [PF6] – at the center in the first coordination shell with the VACF: (a) in the nonpolarizable model and (b) in the polarizable model.
Firstly, the relaxation feature of IL system is examined with the computed polarizability TCF and its time derivative, and then, we discuss dynamical behavior of cations and anions in ILs, and consider what kinds of properties with simulation data we have to check and care in the study of dynamical properties of ILs.
Figure 10.
The computed polarizability TCF for each component (a) and the time derivative of the total polarizability TCF (b) for [BMIm][NTf2] up to about 3 ns.
Figure 10 shows the computed polarizability TCF for each component and the time derivative of the total polarizability TCF for [BMIm][NTf2]. As seen in Figure 10(a), the relaxation of the polarizability TCFs seems to be slower than in usual liquids. In particular, the relaxation of the anion shows faster decay than that of the cation. Also, for each component, it takes about 2~3 ns to approach toward zero. Thus, this indicates that it is required for us to perform a production MD simulation run for, at least, a few ns to study the relaxation behavior of the polarizability TCF. Also, as seen in Figure 10(a), the cation-anion cross-correlation indicates similar variation to the total polarizability TCF. Therefore, the relaxation behavior of the system strongly correlates with the cation-anion cross-correlation, and it is emphasized that tracking cross-correlation terms is very important. On the other hand, corresponding to this result, the time derivative of the polarizability TCF shows long-decay feature extending up to the nanosecond times, as seen in Figure 10(b). Considering that the time derivative of the polarizability TCF is directly related to optical Kerr effect (OKE) response,[14] these results are indicative of how long a MD simulation has to be carried out to examine the relaxation of the total system polarizability and the OKE response.
The relaxation process of ILs implies long-time dynamics of ILs. In particular, it has been known that diffusive motion of each ion in ILs at room temperature is usually much slower than in usual liquids.[1, 16] Therefore, a MD simulation would need to be run for a few ns or more. Then, we are able to use a procedure to check whether a system is in the diffusive regime. The procedure is to compute β(t) = dlog(MSD)/dlog(t),[16] where MSD means the mean-squared displacement (MSD) of the center of mass of ions,〈∑i=1N[ri(t)−r(0)]2〉. In the case that β(t) = 1, the system is in the diffusive regime, while when b(t) < 1, the system is in the sub-diffusive regime. In Figure 11, computed β(t)s for [BMIm][NTf2] are shown. These results indicate that both the cation and the anion reach the β(t) region between 0.8 and 1 (the shadowed area in Figure 11) after about 16 ns. This β(t) region is considered to be almost in the diffusive regime. Therefore, Figure 11 obviously indicates that, for confirming reliable self-diffusivities carefully, we need to perform a longer MD simulation than about 15 ~ 20 ns.
Figure 11.
β(t) (see text) for [BMIm] + and [NTf2]–. Also, see text for the shadowed area.
Figure 12 displays the results calculated for the MSD and, in comparison, the non-Gaussian parameter, α(t), [50, 51]
α(t)=3〈r4(t)〉/5〈r2(t)〉2−1E32
where r(t) is the displacement of an ion at time t with respect to its position at t = 0. In Figure 12, the MSDs of [PMIm]+ and [NTf2]– indicate three typical dynamic ranges (regions), respectively. Also, the behavior of the short- and long-time regimes of MSD for both the cation and the anion is similar to each other.
Figure 12.
MSDs and non-Gaussian parameters for [BMIm] + and [NTf2]–. See text for the shadowed area.
For [BMIm]+, first region is a microscopic regime until about 10 ps. Second is a crossover regime (the shadowed area in the Figure 12) until about 3-5 ns, and third regime is a sub-linear time dependence (sub-diffusive region) toward 20 ns. Also, α(t) has a double peak structure. One is at about 0.5 ps and the other at about 800 ps. In addition, a small shoulder structure emerges at between two peaks. As clearly seen in Figure 12, the short-time maximum peak corresponds to the microscopic region of the MSD, while the long-time maximum is located at around the center of the crossover regime of the MSD. On the other hand, for the anion, [NTf2]–, α(t) has a double peak structure similar to that of the cation. In particular, the long-time maximum peak is slightly shifted to about 300 ps, while the short-time maximum is locacted at the almost same time region as the cation. The time regime where the long-time maximum appear, is in good accord with those of [BMIm]+, as shown in Figure 12. Therefore, considering that the magnitude of the deviation from the Gaussian behavior is remarkably larger for [BMIm]+ than for [NTf2]– as shown in Figure 12, it is expected that the relaxation behavior of the cation in perturbed [BMIm][NTf2] IL might be largely different from that of the anion. This consideration implies the possibility for us to find a distinct relaxation process in the OKE response (compare Figure 10 with Figure 12). Also, it should be noticed that the diffusive regime (linear time-dependence region) after about 10 ns in Figure 12 coincides with the diffusive region in Figure 11 (after about 16 ns) for both the cation and the anion.
In particular, the shadowed area in Figure 12 covers the time range of the full width at half maximum (FWHM) of the second peak of α(t). As pointed out previously,[52] this FWHM range is likely to include a slow decay region usually known as the α relaxation which is characteristic of glass-forming supercooled liquids in general. Also, recently, it has been reported that similar behavior could appear for ILs at room temperature.[51, 53, 54] Therefore, we had better study the α relaxation in ILs with the observations of the non-Gaussian parameter in addition to information on spatial relaxation such as intermediate scattering function. Also, it is suggested that, when we investigate dynamical properties of ILs, the crossover, the sub-diffusive, and the diffusive regimes have to be carefully examined, as one of criterions for reliable research. Thus, only with a long-time MD simulation, it is considered that we are able to investigate dynamical properties of ILs to which experimental procedures are not accessible.
In this chapter, we introduced cross-correlation function analyses with a polarizable model, polarizability TCFs for investigating the relaxation behavior of ILs, and some of procedures for studying dynamical properties. We showed the importance of considering cross-correlation functions and related properties, showing some of examples. Firstly, we employed the polarizable model based on point dipole treatment to the investigation of the polarization effect on the target IL, [BMIm][PF6]. With the MD simulation data for both the nonpolarizable and polarizable models, velocity cross-correlation analyses were shown, and we presented the momentum correlation functions between the cation and anion species in the IL. Next, we computed polarizability TCFs of [BMIm][NTf2] and discussed their relaxation behavior. In addition, we investigated dynamical properties of ILs such MSDs and non-Gaussian parameters and considered what kinds of properties with simulation data we have to check and care in the study of dynamical properties of ILs. This chapter is summarized as follows:
In the study and discussion of velocity cross-correlation functions, it was shown that polarization effects could enhance the cross-correlations between [BMIm] + and [PF6]– in the polarizable model in comparison with that in the nonpolarizable model. These features are ascribed to interionic interactions through attractive forces coming from the charge-dipole and dipole-dipole interactions caused by polarization effects in addition to charge-charge Coulomb interactions. Based on the results of computed cross-correlation between distinct cations (or anions), it was shown that, at early time region, the modulation of cage effects through polarization effects could improve the probability of approach between like ions. Also, by decomposing the cross-correlation function into the longitudinal and nonlongitudinal components, it was indicated that, at the short time region, the velocity cross-correlation is predominantly controlled by the longitudinal contribution. In addition, it was indicated that, compared with the longitudinal correlation in the nonpolarizable model, the longitudinal component is further modified in the polarizable model.
On the momentum correlation functions between [BMIm] + and [PF6]–, it was exhibited that the correlation between [BMIm] + and [PF6]– plays a important role. Also, it was shown that the contribution of the cross-correlation between distinct anions could be enhanced in the polarizable model. This result indicates that cage effects could be diminished with polarization effects, implying that the interionic interactions between distinct anions could be intensified by the charge-charge and dipole-dipole interactions related to polarization effects. Therefore, as has been pointed out in the literature [23, 24, 29], it is considered that the cage effect in ILs could be reduced by many-body polarization effects.
Both the computed polarizability TCF for each component and the time derivative of the total polarizability TCF for [BMIm][NTf2] showed long-time decay behavior. It was indicated that those took about 2~3 ns to approach toward zero. Therefore, this indicates that it is required for us to perform a production MD simulation run. Also, the relaxation behavior of ILs was investigated with the calculation of β(t) = dlog(MSD)/dlog(t) as an indicator of the diffusive region. Our results suggested that long-time dynamics of ILs has to be studied with a longer MD simulation than about 15 ~ 20 ns for confirming reliable self-diffusivities in ILs, carefully. Furthermore, we examined MSDs and non-Gaussian parameters for [BMIm][NTf2]. From our studies, it was exhibited that the magnitude of the deviation from the Gaussian behavior is remarkably larger for [BMIm]+ than for [NTf2]–, and that, comparing with corresponding MSDs, it is possible to study diffusive motion of the cation and the anion. In addition, it was suggested that these studies imply the possibility for us to find a distinct relaxation process in the OKE response. Therefore, it is suggested that we had better study the α relaxation in ILs with the observations of the non-Gaussian parameter in addition to information on spatial relaxation such as intermediate scattering function. Lastly, it is concluded, only with a long-time MD simulation (> 15 ~ 20 ns), that we are able to investigate the dynamical properties of ILs to which experimental procedures are not accessible.
This work was supported in part by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan (Grant-in Aid Scientific Research (C): 23550029).
References
1.WasserscheidP.WeltonT.2008Ionic Liquids in SynthesisWeinheimWiley-VCH
2.WasserscheidP.KeimW.2000Ionic liquids- New "solutions" for transition metal catalysisAngew. Chem., Int. Ed.3937733789
3.CastnerE. W. JrWishartJ. F.ShirotaH.2007Intermolecular Dynamics, Interactions, and Solvation in Ionic LiquidsAcc. Chem. Res.4012171227
4.WeingaertnerH.2007Understanding ionic liquids at the molecular level: facts, problems, and controversiesAngew. Chem., Int. Ed.47654670
5.OhnoH.2005Electrochemical Aspects of Ionic LiquidsHobokenWiley-Interscience
6.RogersR. D.VothG. A.2007Special Issue on Ionic LiquidsSpecial Issue on Ionic Liquids. Acc. Chem. Res.4011
7.WishartJ. F.CastnerE. W. Jr2007Special Issue on Physical Chemistry of Ionic Liquids.Special Issue on Physical Chemistry of Ionic LiquidsSpecial Issue on Physical Chemistry of Ionic Liquids. J. Phys. Chem. B11118
8.OhnoH.FukumotoK.2008Progress in ionic liquids for electrochemical reaction matricesElectrochemistry761623
9.ShirotaH.WishartJ. F.CastnerE. W.Jr 2007Intermolecular Interactions and Dynamics of Room Temperature Ionic Liquids That Have Silyl- and Siloxy-Substituted Imidazolium CationsJ. Phys. Chem. B11148194829
10.XiaoD.RajianJ. R.HinesJ.L. G.LiS.BartschR. A.QuitevisE. L.2008Nanostructural Organization and anion effects in the optical Kerr effect spectra of binary ionic liquids mixturesJ. Phys. Chem. B1121331613325
11.KoebergaM.WuC.-C.KimcD.BonnM.2007THz dielectric relaxation of ionic liquid:water mixturesChem. Phys. Lett.4396064
12.YamamotoK.TaniM.HangyoM.2007Terahertz Time-Domain Spectroscopy of Imidazolium Ionic LiquidsJ. Phys. Chem. B11148544859
13.ShirotaH.NishikawaK.IshidaT.2009Atom substitution effects of [XF6]- in ionic liquids. 1. Experimental studyJ. Phys. Chem. B11398319839
14.IshidaT.NishikawaK.ShirotaH.2009Atom substitution effects of [XF6]- in ionic liquids. 2. Theoretical studyJ. Phys. Chem. B11398409851
15.Refer to the sec.4 in ref.1 and references therein.
16.MaginnE. J.2007Atomistic Simulation of the Thermodynamic and Transport Properties of Ionic LiquidsAcc. Chem. Res.4012001207
17.LopesJ. N. A. C.PaduaA. A. H.2006Nanostructural Organization in Ionic LiquidsJ. Phys. Chem. B11033303335
18.BhargavaB. L.BalasubramanianS.2007Refined potential model for atomistic simulations of ionic liquid [bmim][PF6]J. Chem. Phys.12711451016
19.ZnamenskiyV.KobrakM. N.2004Molecular Dynamics Study of Polarity in Room-Temperature Ionic LiquidsJ. Phys. Chem. B10810721079
20.ShimY.DuanJ.ChoiM. Y.KimH. J.2003Solvation in molecular ionic liquidsJ. Chem. Phys.11964116414
21.KobrakM. N.2006Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulationJ. Chem. Phys.125064502111
22.HuZ. H.MargulisC. J.2006Heterogeneity in a room-temperature ionic liquid: Persistent local environments and the red-edge effectProc. Natl. Acad. Sci.U.S.A.103831836
23.YanT.BurnhamC. J.Del PopoloM. G.VothG. A.2004Molecular Dynamics Simulation of Ionic Liquids: The Effect of Electronic PolarizabilityJ. Phys. Chem. B1081187711881
24.JiangW.YanT.WangY.VothG. A.2008Molecular Dynamics Simulation of the Energetic Room-Temperature Ionic Liquid, 1-Hydroxyethyl-4-amino-1,2,4-triazolium Nitrate (HEATN)J. Phys. Chem. B11231213131
25.UrahataS. M.RibeiroM. C. C.2005Single particle dynamics in ionic liquids of 1-alkyl-3-methylimidazolium cationsJ. Chem. Phys.12202451119
26.UrahataS. M.RibeiroM. C. C.2006Collective excitations in an ionic liquidJ. Chem. Phys.12407451318
27.HuZ.HuangX.AnnapureddyH. V. R.MargulisC. J.2008Molecular Dynamics Study of the Temperature-Dependent Optical Kerr Effect Spectra and Intermolecular Dynamics of Room Temperature Ionic Liquid 1-Methoxyethylpyridinium DicyanoamideJ. Phys. Chem. B11278377849
28.IshidaT.2011Molecular dynamics study of the dynamical behavior in ionic liquids through interionic interactionsJ. Non-Cryst. Solids357454462
29.LeachA. R.2001Molecular Modeling, Principles and ApplicationsHarlowPearson Education
30.VerdaguerA.PadróJ. A.TrullàsJ.1998Molecular dynamics study of the velocity cross-correlations in liquidsJ. Chem. Phys.109228234
31.VerdaguerA.PadróJ. A.2001Computer simulation study of the velocity cross correlations between neighboring atoms in simple liquid binary mixturesJ. Chem. Phys.11427382744
32.FrenkelD.SmitB.2002Understanding Molecular Simulation, From Algorithms to ApplicationsLondonAcademic Press
33.BernardoD. N.DingY.Krogh-JespersenK.LevyR. M.1994An Anisotropic Polarizable Water Model: Incorporation of All-Atom Polarizabilities into Molecular Mechanics Force FieldsJ. Phys. Chem.9841804187
34.PředotaM.CummingsP. T.ChialvoA. A.2002Pair approximation for polarization interaction and adiabatic nuclear and electronic sampling method for fluids with dipole polarizabilityMol. Phys.10027032717
35.AhlströmP.WallqvistA.EngströmS.JönssonB.1989A molecular dynamics study of polarizable waterMol. Phys.68563581
36.FrenkelD.Mc TagueJ. P.1980Molecular dynamics studies of orientational and collision-induced light scattering in molecular fluidsJ. Chem. Phys.7228012818
37.GeigerL. C.LadanyiB. M.1987Higher order interaction-induced effects on Rayleigh light scattering by molecular liquidsJ. Chem. Phys.87191202
38.LopesJ. N. C.DeschampsJPaduaA. A. H.2004Modeling Ionic Liquids Using a Systematic All-Atom Force FieldJ. Phys. Chem. B10820382047
39.LopesJ. N. C.DeschampsJ.PaduaA. A. H.2004Additions and Corrections: Modeling Ionic Liquids Using a Systematic All-Atom Force FieldJ. Phys. Chem. B10811250
40.KöddermannTPaschekDLudwigR2007Molecular Dynamic Simulations of Ionic Liquids: A Reliable Description of Structure, Thermodynamics and DynamicsChemPhysChem824642470
41.SmithW.ForsterT. R.2001The DL_POLY_2 User ManualDaresbury, United KingdomDaresbury Laboratory
42.AllenM. P.TildesleyD. J.1987Computer Simulation of LiquidsClarendonOxford
43.ShirotaHMandaiTFukazawaHKatoT2011Comparison between Dicationic and Monocationic Ionic Liquids: Liquid Density, Thermal Properties, Surface Tension, and Shear ViscosityJ. Chem. Eng. Data5624532459
44.NymandT. M.LinseP.2000Ewald summation and reaction field methods for potentials with atomic charges, dipoles, and polarizabilitiesJ. Chem. Phys.11261526160
45.Van DuijnenP. T.SwartM.1998Molecular and Atomic Polarizabilities: Thole’s Model RevisitedJ. Phys. Chem. A10223992407
46.LideD. R.2006CRC Handbook of Chemistry and PhysicsBoca RatonCRC Press
47.TholeB. T.1981Molecular Polarizabilities Calculated with a Modified Dipole InteractionChem. Phys.59341350
48.Molecular polarizabilities of [BMIm]+ and [NTf2]- were calculated with the same method in [14]. For molecular polarizability tensor elements of [BMIm]+, the reported values were used [14]. With the same procedures as in [14], molecular polarizability tensor elements of [NTf2]- were calculated, and the values of 9.8855, 1.4440 x 10-3, 13.577-5x 10-4,-1.0539, and 10.315 Å3 were used for αXX, αXY, αYY, αXZ, αYZ, and αZZ, respectively.
49.VerdaguerA.PadróJ. A.2000Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential coresPhys. Rev. E62532537
50.RahmanA.1964Correlations in the Motion of Atoms in Liquid ArgonPhys. Rev.136A405A411
51.Del PópoloM. G.VothG. A.2004On the Structure and Dynamics of Ionic LiquidsJ. Phys. Chem. B10817441752
52.ColmeneroJ.AlvarezF.ArbeA.2002Self-motion and the α relaxation in a simulated glass-forming polymer: Crossover from Gaussian to non-Gaussian dynamic behaviorPhys. Rev. E65041804112
53.NgaiK. L.2011Relaxation and Diffusion in Complex SystemsNew YorkSpringer
54.BhargavaB. L.KleinM. L.BalasubramanianS.2008Structural Correlations and Charge Ordering in a Room-Temperature Ionic LiquidChemPhysChem96770