1774
Structural Relaxation and Cation (Li+, H+, Na+) Diffusion in Crystalline Polymer Electrolytes: Ab Initio Molecular Dynamics Simulations
In this work, ab initio molecular dynamics (AIMD) simulations have been undertaken to study the structural stability and cation transport in PEO6:XPF6(X = H+, Li+, Na+) systems with doped cation ions, including interstitials and vacancies. The simulations were performed using the VASP code.4 During the molecular dynamics simulation, the supercell is heated up to 300 K via repeated velocity rescaling for 1 ps; and then being maintained at 300 K for 2.5 ps in a NVT (canonical) ensemble. In the density functional (DFT) calculation, a generalized gradient approximation is applied within the projector-augmented-wave method in the VASP code. Ion diffusion barriers for cations were calculated using the nudged elastic band (NEB) method. The atomic positions and cell parameters of PEO6:LiPF6 were relaxed from the experimental structures, the other two systems are fitted based on the assumption that these three systems are isostructural and adopt the same space group.
In the systems without doped cations (Fig. 1(b)), only the PEO chains in PEO6:LiPF6 maintain the chain shape, and lithium ions are constrained by the PEO chains. Interestingly, the calculated Li-O bond distance shows the Li-O bonding environment differs remarkably from previous studies. For PEO6:HPF6 system, although the protons still reside in the chains, the PEO chains are stretched in the x-y plane. This indicates that the protons have serious impact on the chain structure, resulting in distortions. In PEO6:NaPF6 system, the polymer chains are severely damaged with sodium ions outside of the PEO chains. For systems doped with cation interstitials and vacancies, all three systems show more severe damages in PEO chains. This suggests the introduction of cation point defect cannot guarantee improvement of cation diffusion inside the PEO chains. Other measures, such as end group modification, have to be taken into consideration for enhancing mechanical strength of the PEO chains. Cation diffusion barriers are also calculated using NEB method at 0 K. A vacancy diffusion mechanism is considered in this study (Fig. 1(c)). The proton in PEO6:HPF6 has a small barrier (0.80 eV) when comparing with other two systems (Li+ in PEO6:LiPF6: 4.16 eV, Na+ in PEO6:NaPF6: 3.86 eV) (Fig. 1(d)).
In summary, our AIMD simulations indicate that the introduction of cation point defects alone cannot guarantee improvement of cation transport inside the PEO chains. The diffusion barrier calculations clearly indicates the substantially high mobility of protons in comparison to the other cations in these crystalline polymers.
Acknowledgement
The authors gratefully acknowledge the support by the U.S. Army Research Office under Contract No. W911NF-07-1-0085.
References
1. Z. Gadjourova, Y. G. Andreev, D. P. Tunstall, and P. G. Bruce, Nature 412, 520 (2001).
2. Z. Stoeva, I. Martin-Litas, E. Staunton, Y. G. Andreev, and P. G. Bruce, Journal of American Chemical Society 125(15), 4619 (2003).
3. E. Staunton, Y. G. Andreev, and P. G. Bruce, Faraday Discussions 134, 143 (2007).
4. G. Kresse, and J. Furthmüller, Physcial Review B 54, 11169 (1996).