RETRACTED ARTICLE: Comparing ion transport in ionic liquids and polymerized ionic liquids

Polymerized ionic liquids (polyILs) combine the unique properties of ionic liquids (ILs) with macromolecular polymers. But anion diffusivities in polyILs can be three orders of magnitude lower than that in ILs. Endeavors to improve ion transport in polyILs urgently need in-depth insights of ion transport in polyILs. As such in the work we compared ion transport in poly (1-butyl-3-vinylimidazolium-tetrafluoroborate) (poly ([BVIM]-[BF4])) polyIL and 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM]-[BF4]) IL. The diffusivities of ions in the polyIL and IL were measured and computed. According to the results of the molecular dynamics simulations performed, in the IL the coupling motion between an anion and the ions around determines the ion diffusivities, and the ion association lifetime gives the time scale of ion transport. But in the polyIL, the hopping of an anion among cages composed of cationic branch chains determines the diffusivity, and the associated anion transport time scale is the trap time, which is the time when an anion is caught inside a cage, not the ion association lifetime, as Mogurampelly et al. regarded. The calculation results of average displacements (ADs) of the polyIL chains show that, besides free volume fraction, average amplitudes of the oscillation of chains and chain translation speed lead to various diffusivities at various temperatures.

Polymerized ionic liquids (polyILs) integrate the unique properties of ionic liquids (ILs) with macromolecular polymers.PolyILs are attractive candidates for potential applications in energy storage and electrochemical devices, e.g.super capacitors and batteries [1][2][3][4][5][6][7] .Using these materials can avoid some disadvantages of liquid electrolytes, such as flammability, and furthermore utilize the polymeric nature of polyILs to produce materials with different shape, e.g.thin films.
But the anion diffusivities in polyILs can be three orders of magnitude lower than that in ILs 4 .Researches are being performed to promote anion transport in polyILs 6 , but such endeavors necessarily need in-depth insight of mechanisms of anion transport in polyILs.Zhang and Maginn 8 proposed that in ILs, the ion transport mechanism is the structural relaxation of ion association, but researchers have uncertainty and discrepancy in the ion transport mechanism in polyILs.Some researchers 9 observed that ion transport in polyILs is related to the segmental dynamics of polyILs, but some other researchers [10][11][12][13] reported a decoupling between ionic conductivity and structural dynamics in polyILs.Mogurampelly 14,15 et al. proposed that electrostatic interactions between anions and branch chains of polyILs and how anions move along polyIL chains are the key factors that influence the magnitude of ion diffusivities.They used the lifetime of ion pairs, the continuous ion pairing time, to determine the anion transport time scale in polyILs, but we did not think it is correct.The work shows that the associated anion transport time scale is the trap time, which is the time when an anion is caught inside a cage, not the lifetime of ion pairs.So we obtained totally novel results.Wrong theories would give totally misleading guidance to researches of other researchers, causing large waste of money and labor, both manual and mental.As a result our work is essential and important.
These three groups, MacFarlane and Forsyth [16][17][18][19] , Shaplov and Gerbaldi [20][21][22] , Mecerreyes and Brandell [23][24][25] , had done experimental researches to explore solid-state electrolytes for batteries, focusing on application research.Gerbaldi et al. 26 , and Ahmad and Deepa 27 developed composite materials containing IL used for betteries or electrochromic devices.But none above from a molecular level analyzed the mechanisms that lead to the differences between the conductivity between solid and liquid electrolytes.Such researches are essential, as they can be very important guidance for experimental researches.So if we only develop another solid electrolyte and do no theoretical exploration from a molecular level, then such research will surely lack novelty. 1 School of Resources and Chemical Engineering, Sanming University, Fujian, 365004, China. 2 Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China.✉ e-mail: In this work we did a novel research.From a molecular level, we compared the mechanisms of ion transport in poly(1-butyl-3-vinylimidazolium-tetrafluoroborate) (poly([BVIM]-[BF 4 ])) polymerized ionic liquid (polyIL) and 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM]-[BF 4 ]) ionic liquid (IL).The relationship between ion transport and structure dynamics in the polyIL and the IL was clarified.
The MD methods are extensively applied in scientific researches to obtain in-depth knowledge behind phenomena from a molecular level [37][38][39][40][41] .In the research, MD simulations were carried out to simulate ion transport in the polyIL and IL.

experimental and computational methods
Experimental methods to measure diffusivities.To determine anion diffusion coefficients in the polyIL and IL, pulsed field gradient (PFG) NMR [42][43][44] and broadband dielectric spectroscopy (BDS) for lower temperature 45 were employed.
The PFG is a short, timed pulse with spatial-dependent field intensity.A PFG is identified with four parameters: axis, strength, shape and duration.Gradient enhanced NMR is a method for obtaining high resolution nuclear magnetic resonance spectra.PFG NMR techniques are employed widely for studies of diffusion via diffusion ordered nuclear magnetic resonance spectroscopy [42][43][44] .
BDS is a very powerful experimental tool for understanding the molecular level dynamics in ILs and PolyILs [46][47][48][49] .This technique allows the dielectric response to be measured over 12 decades of frequency at different temperatures, below as well as above its glass transition temperature, giving a valuable insight into molecular mobility of polyILs 50 .
BDS measures the complex dielectric function and the complex conductivity function.The real part of the complex conductivity function in ionic liquids is characterized on the intermediate frequency regime by a plateau, the value of which directly yields the DC conductivity, σ 0 .At lower frequencies, the real part of the complex conductivity function decreases from σ 0 due to electrode polarization.The diffusion coefficient can be determined from BDS spectra by a combination of the Einstein and Smoluchowski equations 45 .The DC conductivity σ 0 is given by where n denotes the effective number density (given by Nγ(T), where γ(T) is the fraction of mobile charge carriers at a given temperature T and N is the total number density of charge carriers available in the material), μ is the mobility of charge carriers, q is the elementary electric charge, D is the diffusion coefficient, k is the Boltzmann constant, l jump (0.5-1 nm) is the displacement anions travel in one jump and ω c is the characteristic (diffusion) rate.
preparation of the polyiL.) was confirmed using elemental analysis.Differential scanning calorimeter (DSC) technique was used to measure the glass transition temperature of the prepared polyIL at a heating and cooling rate of 10 K/min.The midpoint of the total heat flow curve in the thermal transition region was taken as the transition temperature.The obtained glass transition temperature of the prepared polyIL is 414 K. Differential scanning calorimeter (DSC) technique was also used to measure the melting temperature of the IL and the result is 202 K.  51 , the polyIL chain was modeled as H([BVIM]-[BF 4 ]) 31 H, 31 repeating units being adopted.The methyl and ethyl groups of [BMIM] + and poly([BVIM] + ) were treated as united units (a methyl group was regarded as an individual unit in a whole; a ethyl group is regarded as being composed of two individual units, a methyl unit and a methylene unit), while all the other atoms were treated as single units.According to our experience, treating all atoms of methyl or ethyl groups individually only increase the computation cost and relaxation difficulty of the system, while has no contribution to simulation results improvement.Atoms of [BF 4 ] − were treated individually.The simulation cell of the IL has 310 [BMIM]-[BF 4 ] ion pairs and that of the polyIL contains 10 polyIL chains.We applied the AMBER/OPLS force field [52][53][54][55] for parameters of bond, angle, dihedral, electrostatic and van der Waals (VDW) interactions among the units.We computed the electrostatic interaction and VDW interaction with Ewald sum algorithm and Lennard-Jones equation, respectively 55 .The Lorentz-Berthelot mixing rule was applied to compute the interactions between various types of units.

R E T R
VDW and electrostatic interactions are among those that matter most for MD simulation.Table 1 shows the VDW potential parameter and charge values [52][53][54] for the units of the polyIL and IL.The charges of B and F are correlated according to quantum chemical (QC) calculation results with the technique proposed by Singh and Kollman 56 and their VDW energy parameters were fitted through least square fit of the interaction potentials computed with QC method between two anions having different distance.

Potential parameters and control parameters of molecular dynamics (MD) simulation.
Figure 1 shows the geometries of [BF 4 ] − , [BMIM] + and poly([BVIM] + ).The bond and angle interactions were calculated with harmonic potential equation while dihedral interaction was estimated with OPLS dihedral angle potential equation (Eq.( 3)).The electrostatic interactions were estimated using Ewald sum algorithm, while the VDW interactions were computed using Lennard-Jones equation.
bonded unbonded As in MD simulation, VDW and electrostatic interactions are among those that matter most, Table 1 shows the VDW potential parameters and charge values [52][53][54] for units of the polyIL and IL.
The IL and the polyIL densities at 300 K, 1.21 g/cm 3 and 1.31 g/cm 3 , respectively, were employed to calculate the sizes of the simulation cells at 300 K.The simulation cells of the IL and the polyIL contains 310 [BMIM]-[BF 4 ] ion pairs and 10 polyIL chains, respectively.For the simulation of anion transport in the IL, the cell size is 4.00 nm   1. VDW Lennard Jones (LJ) potential parameters and charges for the units of the polyIL and IL (εenergy parameters for the LJ potential, σsize parameters for the LJ potential, q -charges of units).

R E T R
.00 nm × 6.00 nm in x, y, z direction respectively, while for the simulation of anion transport in the polyIL, the cell size is 4.00 nm × 4.00 nm × 5.87 nm in x, y, z direction respectively.The DLPOLY 55 software system was used.Under condition of constant pressure and temperature, the Verlet algorithm was used, in purpose of solving equations of motion with the time step of 2.5 fs.To obtain a reasonable starting configuration of the IL and polyIL systems, energy minimization was done first, using the steepest descent algorithm.Subsequently the structure was relaxed by a 5 ns MD equilibration run in the NPT hoover ensemble.The weak coupling technique was used to adjust the temperature and pressure with relaxation times of 0.1 ps and 0.5 ps, respectively 55 .After the equilibration run, production runs were performed in the NVT ensemble.In this ensemble the averaged kinetic energy E kinetic (T) of molecules was held constant.The temperature was set to 300 K to do simulation at the temperature first, the pressure was set to 1 bar, the VDW potentials cut off distance was set to 1.0 nm and the electrostatic potential cut off distance was set to 1.2 nm.To do simulation at other temperatures, equilibrium runs in the NPT hoover ensemble at those temperatures were performed again, after which production runs in the NVT ensemble at those temperatures were performed.
MD production run time was estimated based on Eq. ( 4): where l jump (~0.5-1 nm) is the jump length, the displacement an anion can make in one single jump.MD production run time is at least several times of the time Eq. ( 4) gives so anions have enough time to complete at least several jumps (5-10) and long enough trajectories can be obtained.Besides, the time should be also long enough for anions to reach the Fickian diffusion regime.Only in the Fickian diffusive regime can diffusivities be calculated through the slope of MSDs versus time plot.
Quantum chemical (Qc) calculations.QC calculations were performed for ) 3 H with Gaussian 09 program package.The density functional theory (DFT) level of calculation was done, using the B3LYP functional with the standard 6-31++G(2d,p) basis set.Vibrational frequencies calculation was performed to verify that the optimized structures are globally optimal.The optimized structures were used to build the models for MD simulations.
The results of QC optimization structures for an ion pair are presented in Fig. 2. The structures of [BMIM] + and [BF 4 ] − in Fig. 1 only differ slightly from that in Fig. 2. Such minor differences would have no influence on MD simulation results.The interaction potential between ion pair is −37.04 kcal/mol.The ion association interactions give the IL brilliant properties, e.g.high viscosity, low vapor pressure and good thermal and chemical stability.

Diffusivity calculation.
The mean squared displacements (MSD) which were averaged on various time origins were applied to calculate anion diffusivities according to Eq. ( 5) [57][58][59] where R(t) represents the anion location at t and <…> stands for an average of moving time and the average is calculated with respect to t 0 .The logarithmic plot of MSD versus time was plotted to determine if the diffusion was in the Fickian diffusive regime.The MSD were averaged over different time origin.Then diffusivities were calculated through the slope of the plot of MSD versus time.The slope was obtained through a least square fit of the plot.Long time of production run was done to ensure ion transport trajectories are long enough.The diffusivities obtained are the average of all the anions.The sampling size is very large.All these ensure the accuracy of the calculated diffusivities and repeated simulations show that the relative standard deviation (RSD) (standard deviation as a percentage of average diffusivity) in the calculated diffusivities is smaller than 5%.
For Δt, the average displacement (AD) of a chain is estimated according to Eq. ( 6).

R E T R
is to compute the AD of unit i of the chain for a time interval of Δt. <…> stands for a moving time average.

Results and discussion
Figure 3 shows the snapshots of the inside of the simulation cells of the IL and the polyIL.It is observed that in the polyIL, the cations are connected together by chains, so small or large cavities among cations are formed and anions are constrained inside these cavities before they get chances to jump to cavities nearby.On the contrary, in the IL, the cations are not linked, so they can move freely as anions do.
Figure 4 shows the diffusivities of anions of the polyIL and IL at various temperatures obtained via experiments and MD simulations.The values obtained via both methods are in agreement.At the same temperatures, anions in the polyIL have diffusivities 3 orders of magnitude lower than that of anions in the IL.In the calculated diffusivities the relative standard deviation (RSD) is smaller than 5%.
To clarify the underlying mechanisms, the exploration of the three diffusion regimes of anion transport, ballistic regime, the subdiffusive regime and Fickian diffusive regime, especially the subdiffusive regime is necessary.As presented in Fig. 5(a), in the ballistic regime, the subdiffusive regime and the Fickian diffusive regime, MSDs are proportional to t 2 , t 0.5 and t, respectively.
The logarithmic plots of MSD versus time for anion transport in the polyIL and IL at 300 K are presented in Fig. 5(b), which shows the three diffusion regimes above.The ballistic regime is shorter than the other two regimes (t < 2.0 ps).Compared to the IL, the polyIL has short ballistic regime.In the ballistic regime, anions fly freely, having no bounces with other ions 60,61 .Hence the results show that in the polyIL, the average cavity size is smaller than that in the IL and an anion is constrained more stringently in the polyIL than in the IL.
The subdiffusive regime is longer than the ballistic regime according to Fig. 5.The subdiffusive regime of the polyIL is longer than the IL.The hindrance influence of surrounding units on anions lead to the subdiffusion phenomena 60 , so studying the subdiffusive regime is essential for clarify the mechanisms behind.
Different models for anomalous diffusion models, e.g.random walk on a fractal (RWF), fractional Langevin equation (FLE) and continuous time random walk (CTRW) have been presented 62 .The schemes in Fig. 6 show these models.
A rather good example of where anomalous CTRW have been successfully used to model a physical process is DNA gel electrophoresis, as presented in Fig. 6(a).In the process, an electric field is applied to force DNA molecules through a porous gel.The rate at which molecules pass through the gel is determined by their mobility, which is a function of size and charge.For small sections of DNA, separation can be achieved successfully.But large molecules are often trapped in the gel, causing a long waiting time between displacements.The work of Weiss 63 shows that this process can be described successfully by subdiffusive CTRW model, with finite drift and zero variance in the step size, but a diverging mean waiting time In contrast, FLE is a Gaussian process in which anomalous diffusion comes about via statistical correlations between successive steps.Such FLE processes are widely associated with modeling of stochastic transport in viscoelastic environments.FLE equations can be postulated but may also arise as an effective description of more complex multiparticle Markovian dynamics 64 .

R E T R
shows an example of the FLE process, the subdiffusion of fullerene in molten polymers 64 , which is incurred by the coupling between the motions of fullerene molecules and the molten polymers or fullerene molecules nearby.
RWF model can describe the effect of geometric constraints arising in complex geometries, is not space-filling, and displays non-Gaussian propagators 64 .Figure 6(c) shows one example of the RWF process, where penetrants can only "walk" along certain passages.When gas penetrants diffuse in polymers, like poly (4-methyl-2-pentyne) (PMP), such processes are observed as well, as presented in Fig. 6(d) 65 From Fig. 6(d), it is observed that there are small vacancies among the long chains of PMP, penetrants, like methane and n-butane can only diffuse among vacancies through passages connecting these vacancies.Figure 6(e) shows the snapshot of the PMP-gas system after the penetrants are placed in the vacuum portion at the left and right sides and the simulation production run has been performed for 2 ns.Most of the n-butane molecules penetrate deeply into the bulk, whereas most of the methane molecules placed at the same time as n-butane molecules only wander around the interface.Processes that can be described by the CTRW model have random penetrant jump time and penetrant jump length.In the research motion of anions is not in a random way, but is influenced by cavities nearby and chain oscillation.As a result, we cannot apply the CTRW model to describe the subdiffusive behavior of anions in the research.Furthermore, processes describable by the FLE model are Gaussian, while that describable by the RWF model are non-Gaussian 41 .
The minimal excursion ratio (MME) method, given by Meroz and Sokolov 57 , was used to determine if a transport process is non-Gaussian or Gaussian: where R(t) is the position of gas molecules at t, <…> denotes a moving time average and the sub-index t 0 indicates the average it taken with respect to t 0 .From Fig. 7, the MME ratio is approaching one for anion transport in the IL while that for anion transport in the polyIL is not tending to one.Hence the subdiffusive behavior of anions in the IL is Gaussian while that in the polyIL is non-Gaussian.The FLE model and the RWF model should be applied to explain the subdiffusion phenomena in the IL and the polyIL, respectively.
According to the FLE model, in the IL the coupling between the motions of an anion and the ions nearby causes the subdiffusive behavior of the anion.The same explanation is given by the scaling theory 66 , which explains that relaxation processes of the ions lead to subdiffusive phenomena and when the motions of an anion and the ions around complete coupling, the relaxation processes cease.For such type of subdiffusion phenomena, the ion association lifetime gives the time scale.
According to the RWF model, the trap of anions in cages consisting of chains lead to the subdiffusive behavior of anions in the polyIL.Cages in the polyIL are stable and anions can be caught for long time.The associated anion transport time scale is hence the trap time, the time when an anion is caught inside a cage, not the ion association lifetime.The following equation can be applied to give an estimate of the trap time: where l jump (~0.5-1 nm) is the jump length, the displacement an anion can make in one single jump.From Fig. 5, longer subdiffusive regime is observed for the polyIL than the IL.As presented above, it is because of the anion trap in the stable cages, which consist of polyIL chains.On the contrary, there are no such cages in the IL because cations are free to move and not constrained by the main chains.The motion of main and cationic branch chains that compose cages need to be investigated in order to have detailed knowledge of the cages.
Figure 8 shows the trajectory Scheme of an anion in [BMIM]-[BF 4 ] for a time interval of 2.5 ns and in poly([B-VIM]-[BF 4 ]) for a time interval of 550 ns.It is observed that the anion in the IL have remarkably higher diffusivity than that in poly([BVIM]-[BF 4 ]) according to Eq. ( 5).The anion in poly([BVIM]-[BF 4 ]) made a jump from a cage to another cage and the trajectories in the two cages are discrete.
Figure 9(a) is the plot showing the ADs of the polyIL chain units versus time intervals.The ADs increase sharply to a stair value in short time intervals (~3.5 ps).We know polymer chains swing around their balance locations rapidly.The stair value is an estimate of the average swing amplitudes.The ADs of chains are composed of the swing and the translation portions.The ADs nearly remain constant, after reaching a stair value, so the speed of the chain translation is tiny.The rates of the chain translation are around 0.00161, 0.00188, 0.00219 and 0.00251 Å/ps corresponding to 300 K, 400 K and 500 K, respectively.With the translation rate and the average swing amplitude shown in Fig. 9(b), ADs for any time intervals may be estimated.The promotion of both the average swing amplitude and the chain translation rate is observed with the increase of temperatures.
The free volume fraction (FVF) of a polymer also affects the magnitude of diffusivities 37,38 .In order to shun the FVF effects, for 400 K and 500 K we set the volume of the simulation cells the same.According to the simulation results, the calculated diffusivities of anions and the average swing amplitudes of the polyIL chains at 500 K are larger than that at 400 K.As a result, the average swing amplitude the polyIL chains is doing contributions to the various diffusivities at various temperatures.Anions will be caught in cages and the anion diffusivity will be null if the polyIL chains remain static.With the promotion of the average swing amplitude, anions achieve more chances to jump to cages nearby and the diffusivities get promotion accordingly.
With the increase of temperatures, the increase of both the average swing amplitude and the chain translation rate is observed.The cages are composed of branch and main chains.Higher chain swing amplitudes cause more opportunities in the same time periods for the emergence of large enough holes in the cages which anions can pass.Due to the promotion of the FVF of the polyIL with temperatures, larger spaces between the chains emerge.These factors mentioned above take effects, leading to increased opportunities for an anion to flee from a cage and hop to a cage around.
An anion should overcome the electrostatic and VDW interactions with the units of the chains composing the cage, in order to flee from a cage.When the anion approaches the cage wall, if there are no VDW interactions, the attractive electrostatic interactions between the anion and cationic branch chains that the anion is getting closer will become higher; the anion can hence flee from the cage when it makes the first try.Therefore electrostatic interactions alone cannot catch an anion inside a cage.When the anion approaches the cage wall, the repulsive VDW interactions between the anion and cage units that the anion is getting closer get higher, larger that the corresponding electrostatic interactions.The anion is pushed back to the cage inside.Hence it is VDW interactions that have the largest contribution in trapping anions.From Eq. ( 4), the length of the trap time determines the diffusivity.As a result VDW interactions exert a more important influence on anion diffusivities than electrostatic interactions, including ion association interactions.Analyzing only electrostatic interactions between ions cannot give the correct ion transport mechanisms.On the contrary, such analysis can only give wrong and misleading mechanisms, as Mogurampelly et al. did 14 .We cannot deny that ion associations affect anion diffusivities, but VDW interactions exert the most important influence, as presented above.
Mogurampelly et al. 14 proposed that if the cation-anion distance is below a benchmark value an ion pair is established and if the distance is above that value the ion pair is separated.They applied the continuous ion pairing time, the time interval when the ion pair remains paired continuously, to describe the time scales underlying ion transport.
Anions moves up and down, left and right rapidly inside cages composed of chains of the polyIL.Therefore at one moment an anion may be closer than a benchmark value to a cationic branch chain, but after a while the distance between the anion and cationic branch chain can be larger than that value.Then with time going on, the distance can vary to values larger or smaller than that value.However, in the whole process, the anion is inside the cage, i.e. whether the distance between the ion pair is larger or smaller than the benchmark value making no fundamental differences.Hence to use the continuous ion pairing time above to estimate the anion transport time scale is totally wrong.Definitely the trap time of an anion inside a cage should be applied instead.So we obtained totally novel results.Our work is essential and important, as wrong theories would give totally misleading guidance to researches of other researchers, causing large manual and mental waste of money and labor.
The diffusivities and trap times of anions in the polyIL at various temperatures are shown in Fig. 10.The trap times in the figure are the averaged trap times of all the anions.The promotion in diffusivities is observed with the increase of temperatures.Diffusivities go down with the increase of trap times and are inversely proportional to trap times, matching Eq. ( 4).The trap times decrease with temperatures.This can be explained by higher temperatures causing higher chain swing amplitudes, larger chain translation rate and larger spaces between chains and anions obtaining more hopping opportunities to flee from the cage.

conclusions
In conclusion, in the work we explored and compared the mechanisms and time scales of ion transport in the polyIL and the IL.In the IL the coupling between the motions of an anion the ions nearby cause the subdiffusive behavior of the anion, while in the polyIL, the trap of anions in cages consisting of chains lead to the subdiffusive behavior of anions.
Mogurampelly et al. stated that electrostatic interactions between anions and branch chains of polyILs and how anions move along polyIL chains are the key factors that influence the magnitude of ion diffusivities.They used the lifetime of ion pairs, the continuous ion pairing time, to determine the anion transport time scale in polyILs, but we did not think it is correct.The work shows that the associated anion transport time scale is the trap time, which is the time when an anion is caught inside a cage, not the lifetime of ion pairs.So we obtained totally novel results.Our work is essential and important, as wrong theories would give totally misleading guidance to researches of other researchers, causing large manual and mental waste of money and labor.
Higher temperature bring about higher chain swing amplitudes, leading to more opportunities in the same time periods for the emergence of large enough holes in the cages which anions can pass.Due to the promotion of the FVF of the polyIL with temperatures, larger spaces between the chains emerge.These factors mentioned above take effects, leading to increased opportunities for an anion to flee from a cage and hop to a cage around.
In the polyIL, VDW interactions exert a more important influence on anion diffusivities than electrostatic interactions, including ion association interactions.
The diffusivities go down with the increase of trap times and are inversely proportional to trap times.The trap times decrease with temperatures.This can be explained by higher temperatures causing higher chain swing amplitudes, larger chain translation rate and larger spaces between chains, so anions obtaining more hopping opportunities to flee from the cage.

Figure 2 .
Figure 2. The optimal geometry of an ion pair.

Figure 3 .
Figure 3. Snapshots of the inside specifics of the simulation cells of (a) the IL and (b) the polyIL (blue: [BF 4 ] − , red: [BMIM] + in (a) or chains in (b)).

Figure 4 .
Figure 4.The diffusivities of the anions of the polyIL and IL at various temperatures.

Figure 5 .
Figure 5. (a) The scheme of three diffusion regimes.(b) The logarithm plot of MSD versus time for the anion transport in the IL and polyIL.R E T R A C T E D A R T I C L E

Figure 7 .
Figure 7.The MME for the anion transport trajectories in the IL and polyIL.

Figure 8 .
Figure 8. Trajectories of an anion (a) in the IL for a time interval of 2.5 ns and (b) in the polyIL for a time interval of 550 ns.The red, blue and green points on XY, YZ and ZX planes are projections on the corresponding planes.

Figure 9 .
Figure 9. (a) The ADs of polyIL chain units for various time intervals.(b) The average oscillation amplitudes of the main and branch chains of the polyIL as a function of temperatures.

Figure 10 .
Figure 10.The anion diffusivities and trap times in the polyIL at different temperatures.R E T R A C T E D A R T I C L E