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 [18]. 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 heterodynedetected Ramaninduced Kerr effect spectroscopy (OHDRIKES)[3, 9, 10] and THz timedomain spectroscopy (THzTDS)[11, 12] have been applied to investigate the intermolecular vibrational dynamics in ILs. In particular, with OHDRIKES 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[1518] and novel interionic dynamics under solvation dynamics [1921], dynamical properties [2226], 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 crosscorrelation terms between cations and anions, and the other is polarization effects due to manybody interactions caused by cations and anions in ILs. While a large number of experimental approaches have been applied to investigate these subjects, molecularlevel 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 wellparameterized potential functions and partial charges [29].
With the MD simulation procedures, it is possible to study the effects of the crosscorrelations 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 crosscorrelation 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 manybody 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 timedependent 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 offdiagonal 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 1butyl3methylimidazolium cation based ILs with hexafluorophosphate anion, [BMIm][PF_{6}], and bis(trifluoromethylsulfonyl)amide anion, [BMIm][NTf_{2}], as target systems. (See Figure 1 for all the molecular structures of [BMIm]^{ +}, [PF_{6}]^{–} and [NTf_{2}]^{–}.) Firstly, we focus on the collective properties of [BMIm][PF_{6}] with the crosscorrelation functions of ionic species in the IL and polarization effects on ionic motions. As a second, we investigate dynamical properties of [BMIm][NTf_{2}]. Below, we start from the explanation for the velocity crosscorrelations of a central atom with neighboring atoms. Then, we show how to evaluate both crosscorrelation 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 dipoledipole (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.
2. Tracking Ionic Motions Through Interionic Interactions: CrossCorrelation, 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 crosscorrelations 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, crosscorrelated ionic motions, and the influence of polarization coming from manybody 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 crosscorrelated 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 crosscorrelation functions. In particular, through the analysis of the longitudinal and nonlongitudinal contributions to the velocity crosscorrelation 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 nonpolarizable 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 cationanion is different from those of cationcation and anionanion 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 crosscorrelation analyses is important to investigate the interplay between cation and anion species in ILs, and it is shown that the crosscorrelation analyses provide contrastive features.
Different from a usual velocity autocorrelation functions (VACF), crosscorrelation 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 crosscorrelation 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 crosscorrelation function correlates with the alteration of the VACF. In particular, the crosscorrelation 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 crosscorrelations (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 crosscorrelation functions. Also, considering that these crosscorrelation functions could be modulated by the strength of interionic interactions and the coupling between cation and anion motions, crosscorrelation functions would be largely influenced by polarization effects.
On the other hand, as another type of the appearance of interplay between crosscorrelation 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 interactioninduced multipole effects. Obviously, it is expected that molecular polarizabilites are influenced by interionic crosscorrelations. 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 crosscorrelations.
3. Theoretical Background
Here, we introduce crosscorrelation functions, and then, give an explanation of a polarizable model and polarizability time correlation function. Computational details are also given.
3.1. CrossCorrelation Function
A velocity autocorrelation function (VACF) is defined for calculating the velocity correlation of a same particle as follows,
where v_{ i }(t) is the velocity of ith species. Different from the VACF, the time crosscorrelation 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],
where v_{
i
}(t) is the velocity of the species i and v_{
j
}(t) is the velocity of the species j.
where
where ρ is the number density of the system.
Here, we can define the crosscorrelation in ILs as that between a centered cation and the total contribution of other anions [28, 30],
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 crosscorrelation between a centered anion and other cations.
Also, with the computation of the crosscorrelation 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,
(6) 
where p _{+} and p _{‒} mean the momentums of cation and anion, respectively. The transference of momentum of an anion to distinct anions and cations, P _{ n } ^{A(total)} (t), is also given as follows,
3.2. Polarization Effects
3.2.1. A Polarizable Model
Here, we introduce a polarizable model considering induced dipole moments and explain the procedure for computing induced dipole moments and the manybody polarization energy [3335].
The total potential energy of the system under the resulting polarizable force field is defined as follows,
where the terms V _{ bond } and V _{ nonbond } are intra and intermolecular interaction energies. The polarization energy, V _{ pol }, is decomposed into the three terms as follows,
where the chargedipole contribution, V _{ chargedipole }, the dipoledipole contribution, V _{ dipoledipole }, and the selfpolarizability term, V _{ self }, are defined, respectively, as
In the above equations, E
_{
i
} is the electric field on atom i, produced by the partial charges of all other surrounding atoms, and
where r _{ij} = r _{i} – r _{j}. The induced dipole moment on atom i is given by
With Eqs. (10), (11), and (12), the polarization energy, E _{ pol }, is summarized as follows,
3.2.2. Polarizability TimeCorrelation Function (TCF)
The polarizability anisotropy of the system can be tractable by calculating the TCF of offdiagonal 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,
where t represents the time dependence, and subscripts M and II mean the molecular part and the interactioninduced part, respectively. The molecular part is given by the sum of the polarizability tensors of isolated gas phase molecular polarizability in the laboratory frame,
where N is the number of molecules, and α_{i} is the polarizability tensor of molecule i. For the formulation of the interactioninduced part, we employ the dipoleinduceddipole (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 interactioninduced polarizability in the DID approximation is given as follows,
where T
_{
ij
} means the dipole interaction tensor between molecules i and j.
Equation (19) can be solved by the calculation procedure that is called the allorders DID approximation [37]. It should be noted that the DID model employed here assumes a centercenter 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],
where the superscripts C and A represent cation and anion species, respectively. Referring to Equations. (16), (17), and (18), the
where the indices k and l go over all cations and anions, respectively. Then, the total polarizability of the system,
(23) 
Where
(24) 
where



The nuclear response function, R(t), is represented as the time derivative of the polarizability TCF,
3.3. Computational Details
Utilizing the sets of force field parameters for [BMIm]^{+}, [PF_{6}]^{–} and [NTf_{2}]^{–} [3840], 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][PF_{6}] and 5000 for [BMIm][NTf_{2}]) 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][PF_{6}] and [BMIm][NTf_{2}] at 298 K [13, 43]. 12 Å was set as cutoff length. The time step was 2 fs. The longrange 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][NTf_{2}]) and collected during 10 ns (20 ns for [BMIm][NTf_{2}]) 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,
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][NTf_{2}], we used molecular polarizabilities of 14.372 Å^{3} for [BMIm]^{+} and 11.259 Å^{3} for [NTf_{2}]^{–}.[48] The bodyfixed 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. CrossCorrelation, Momentum Correlation, and Polarization Effects
The computed velocity crosscorrelation functions [28] are shown in Figures 4 and 5. These figures show the comparison of velocity crosscorrelation functions observed around the cation, [BMIm]^{+}, and the anion, [PF_{6}]^{ –}, placed at the center, for both the nonpolarizable and polarizable models. From the result of the cationanion RDF in Figure 2, we selected 3.5, 8.4, 8.4, and 14.75 Å for the value of a_{1}, b_{1}, a_{2}, and b_{2} 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 (C_{1}(t)) and second (C_{2}(t)) coordination shells were specified. It should be noted that the C_{1}(t) between the center anion and distinct anions is not shown in Figure 5 because the C_{1}(t) is almost zero corresponding to the result that the anionanion 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, C_{n}(0), are negative since the system size is finite as have been pointed out [30].
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 crosscorrelation function, C_{1}(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 C_{1}(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 chargedipole and dipoledipole interactions by polarization effects in addition to chargecharge 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 crosscorrelation 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].
On the other hand, an interesting feature is seen in the C_{1}(t) between distinct cations (or the C_{2}(t) for distinct anions). As indicated in Figures 4 and 5, the crosscorrelation function for distinct ions reaches the maximum point earlier than the C_{1}(t) for the cationanion crosscorrelation, even though the peak height is smaller than that in the cationanion crosscorrelation 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 crosscorrelation 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 crosscorrelation 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 C_{2}(t) is shifted to the time region later than that of C_{1}(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 crosscorrelation functions. As shown in Figure 6, we consider the decomposition of a crosscorrelation function into the longitudinal (denoted as L in the figure), C_{n} ^{L}(t), and nonlogitudinal (denoted as NL in the figure), C_{n} ^{NL}(t), contributions as follows [28, 31, 49],
where CR_{n} ^{L}(t) is represented as the velocity crosscorrelation along the direction designated by the center of masses of distinct ions at t = 0. We can compute CR_{n} ^{L}(t) with the following equation [28, 31, 49],
where v_{
i
}
^{L}(t) = v_{
i
}(t)[r_{ij}(0)/r_{ij}(0)] and v_{
i
}(t) is the velocity of the ionic species i. Also, r_{
ij
}(0) means the direction vector between the center of masses of the distinct ions i and j.
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.
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, P_{1} ^{CC}(t), is smaller than that between the cation and anion molecules, P_{1} ^{CA}(t), the peak height of the P_{1} ^{CA}(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 chargedipole and dipoledipole interactions coming from polarization effects. As we could observe in the crosscorrelation 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.
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, P_{1} ^{AA}(t), indicates much smaller contributions to the total momentum correlations, P_{1}(t), in both the nonpolarizable and polarizable models, compared with the P_{1} ^{AC}(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 cationanion 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 P_{1} ^{AA}(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 P_{1}(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 chargecharge and dipoledipole interactions [28].
5. Relaxation Processes and Dynamical Properties
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 shows the computed polarizability TCF for each component and the time derivative of the total polarizability TCF for [BMIm][NTf_{2}]. 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 cationanion crosscorrelation indicates similar variation to the total polarizability TCF. Therefore, the relaxation behavior of the system strongly correlates with the cationanion crosscorrelation, and it is emphasized that tracking crosscorrelation terms is very important. On the other hand, corresponding to this result, the time derivative of the polarizability TCF shows longdecay 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 longtime 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 meansquared displacement (MSD) of the center of mass of ions,
Figure 12 displays the results calculated for the MSD and, in comparison, the nonGaussian parameter,
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 [NTf_{2}]^{–} indicate three typical dynamic ranges (regions), respectively. Also, the behavior of the short and longtime regimes of MSD for both the cation and the anion is similar to each other.
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 35 ns, and third regime is a sublinear time dependence (subdiffusive 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 shorttime maximum peak corresponds to the microscopic region of the MSD, while the longtime maximum is located at around the center of the crossover regime of the MSD. On the other hand, for the anion, [NTf_{2}]^{–}, α(t) has a double peak structure similar to that of the cation. In particular, the longtime maximum peak is slightly shifted to about 300 ps, while the shorttime maximum is locacted at the almost same time region as the cation. The time regime where the longtime 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 [NTf_{2}]^{–} as shown in Figure 12, it is expected that the relaxation behavior of the cation in perturbed [BMIm][NTf_{2}] 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 timedependence 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 glassforming 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 nonGaussian 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 subdiffusive, and the diffusive regimes have to be carefully examined, as one of criterions for reliable research. Thus, only with a longtime MD simulation, it is considered that we are able to investigate dynamical properties of ILs to which experimental procedures are not accessible.
6. Conclusion
In this chapter, we introduced crosscorrelation 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 crosscorrelation 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][PF_{6}]. With the MD simulation data for both the nonpolarizable and polarizable models, velocity crosscorrelation 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][NTf_{2}] and discussed their relaxation behavior. In addition, we investigated dynamical properties of ILs such MSDs and nonGaussian 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 crosscorrelation functions, it was shown that polarization effects could enhance the crosscorrelations between [BMIm]^{ +} and [PF_{6}]^{–} 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 chargedipole and dipoledipole interactions caused by polarization effects in addition to chargecharge Coulomb interactions. Based on the results of computed crosscorrelation 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 crosscorrelation function into the longitudinal and nonlongitudinal components, it was indicated that, at the short time region, the velocity crosscorrelation 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 [PF_{6}]^{–}, it was exhibited that the correlation between [BMIm]^{ +} and [PF_{6}]^{–} plays a important role. Also, it was shown that the contribution of the crosscorrelation 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 chargecharge and dipoledipole 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 manybody polarization effects.
Both the computed polarizability TCF for each component and the time derivative of the total polarizability TCF for [BMIm][NTf_{2}] showed longtime 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 longtime dynamics of ILs has to be studied with a longer MD simulation than about 15 ~ 20 ns for confirming reliable selfdiffusivities in ILs, carefully. Furthermore, we examined MSDs and nonGaussian parameters for [BMIm][NTf_{2}]. From our studies, it was exhibited that the magnitude of the deviation from the Gaussian behavior is remarkably larger for [BMIm]^{+} than for [NTf_{2}]^{–}, 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 nonGaussian parameter in addition to information on spatial relaxation such as intermediate scattering function. Lastly, it is concluded, only with a longtime MD simulation (> 15 ~ 20 ns), that we are able to investigate the dynamical properties of ILs to which experimental procedures are not accessible.