Research Article
Research Article
Research Article
Research
Volume 2019, Article ID 2198906, 10 pages
https://doi.org/10.34133/2019/2198906
Research Article
Facilitating Lithium-Ion Diffusion in Layered Cathode
Materials by Introducing Li+/Ni2+ Antisite Defects for High-Rate
Li-Ion Batteries
Zhongfeng Tang,1 Sen Wang,2 Jiaying Liao,1 Shuo Wang,1 Xiaodong He,1 Bicai Pan,2
Haiyan He,2 and Chunhua Chen 1
1
CAS Key Laboratory of Materials for Energy Conversions, Department of Materials Science and Engineering & Collaborative
Innovation Center of Suzhou Nano Science and Technology, University of Science and Technology of China,
Correspondence should be addressed to Haiyan He; [email protected] and Chunhua Chen; [email protected]
Copyright © 2019 Zhongfeng Tang et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a
Creative Commons Attribution License (CC BY 4.0).
Li+/Ni2+ antisite defects mainly resulting from their similar ionic radii in the layered nickel-rich cathode materials belong to one of
cation disordering scenarios. They are commonly considered harmful to the electrochemical properties, so a minimum degree of
cation disordering is usually desired. However, this study indicates that LiNi0.8Co0.15Al0.05O2 as the key material for Tesla
batteries possesses the highest rate capability when there is a minor degree (2.3%) of Li+/Ni2+ antisite defects existing in its
layered structure. By combining a theoretical calculation, the improvement mechanism is attributed to two effects to decrease
the activation barrier for lithium migration: (1) the anchoring of a low fraction of high-valence Ni2+ ions in the Li slab pushes
uphill the nearest Li+ ions and (2) the same fraction of low-valence Li+ ions in the Ni slab weakens the repulsive interaction to
the Li+ ions at the saddle point.
(003)
(003)
(003)
(104)
(104)
(104)
(101)
(101)
(101)
(108/110)
(108/110)
(108/110)
(006/104)
(006/104)
(006/104)
(107)
(015)
(113)
(107)
(107)
(015)
(015)
(113)
(113)
10 20 30 40 50 60 70 10 20 30 40 50 60 70 10 20 30 40 50 60 70
2 theta (degree) 2 theta (degree) 2 theta (degree)
obs obs obs
cal cal cal
diff diff diff
(a) (b) (c)
4.0
(003)
(003)
3.5
3.0
(104)
(104)
2.5
2.0
(101)
(101)
(108/110)
(108/110)
1.45
(006/104)
(015)
(107)
(113)
(015)
(107)
(113)
1.40
1.35
1.30
1.25
10 20 30 40 50 60 70 10 20 30 40 50 60 70 720 735 750 765 780
2 theta (degree) 2 theta (degree) Sintering temperature (°C)
Figure 2: Rietveld refinement of XRD. (a–e) Rietveld analysis of NCA-720, NCA-735, NCA-750, NCA-765, and NCA-780, respectively. (f)
The refinement results of the ratio of Li/Ni exchange and the peak intensity ratio of (003)/(104).
unit cell [34]. The energy relaxation iterates until the forces For the five samples, the average valences of transition
acting on all the atoms were less than 10-2 eV/Å. The climb- metals (Ni and Co) are almost the same according to the
ing image nudged elastic band (CI-NEB) method was used semiquantitative analysis of X-ray photoelectron spectroscopy
to determine the minimum energy paths for diffusion of Li (XPS) (Figures S2 and S3). For the cation disordering, only
atoms in LiNi0.8Co0.2O2 [35–37]. With these settings, we Li+ and Ni2+ occupations are considered to be exchanged
obtained the optimal crystallographic parameters of LiNi0.8- because Co3+ (0.061 nm) and Al3+ (0.054 nm) in the lithium
Co0.2O2, a = 2 79 Å, b = 2 95 Å, and c = 14 09 Å, which are layer would cause a drastically increased system energy due
in good agreement with the experimental values to the big differences in ionic radii and valences with Li+.
(a = 2 87 Å, b = 2 87 Å, and c = 14 13 Å) [38]. The refinement results are given in Table 1 and Figure 2(f),
according to which two crystal diagrams of the TM layer
and the Li layer (seen from c-axis direction) for NCA-750
3. Results
are given in Figure S4 to show the atom arrangements. It is
3.1. Morphology and Structure Characterization. Figure S1 quite clear that the ratio of Li+/Ni2+ exchange is highly
shows the SEM images of the NCA samples sintered at related to the intensity ratio of I 003 /I 104 that has been
different temperatures, which consist of irregular particles used to qualitatively analyze the degree of cation mixing in
around 700 nm-1 μm and without significant differences layered materials. It is widely believed that the lower the
in the morphology and particle size. As indicated in degree of cation disordering, the better the electrochemical
Figures 2(a)–2(e), all their diffraction peaks can be well properties [12, 14–16, 42]. NCA-720 has a relatively high
indexed in a hexagonal structure of α-NaFeO2 type without degree of Li/Ni mixing because the temperature is not high
any impurity [39, 40]. In order to reveal the detailed crystal enough to remove all of the lattice defects. The rest of the
information, the XRD data have been refined by the samples show an approximate linear relationship between
Rietveld method using GSAS software. In the refinement, the ratio of Li/Ni exchange and sintering temperature,
Li1/Ni2 are set at 3b sites (0, 0, 0.5), Li2/Ni1/Co1/Al1 at 3a which may be attributed to the phase transition of LiNiO2
sites (0, 0, 0), and O at 6c sites (0, 0, z) where z = 0 25 [41]. with increasing sintering temperature. It seems that NCA-
4 Research
Table 1: Summary of the Rietveld result of the five NCA samples based on the space group of R-3m.
Samples a (Å) c (Å) c/a χ2 Rwp% Rp% I 003 /I 104 Li/Ni exchange
NCA-720 2.8615 (4) 14.1738 (6) 4.9533 6.11 9.6 7.6 1.37131 3.19%
NCA-735 2.8614 (8) 14.1756 (4) 4.9541 7.629 10.3 7.6 1.45642 1.69%
NCA-750 2.8619 (3) 14.1776 (5) 4.9539 7.394 10.4 7.9 1.44492 2.39%
NCA-765 2.8628 (3) 14.1776 (3) 4.9539 9.188 11.4 8.5 1.41527 3.20%
NCA-780 2.8632 (1) 14.1781 (2) 4.9518 8.351 11 8.1 1.31637 3.82%
1.0
4.4 NCA-720
0.8 NCA-735
4.2 NCA-765
4.0 0.6
NCA-780
Current (mA)
3.8
Voltage (V)
0.4 NCA-750
3.6
0.2
3.4
3.2 0.0
3.0 −0.2
3.8
−0.4
2.6
0.5C 0.1C
180 1C 150
2C
160 5C
120
−Z″ (𝛺)
140 10C
120 90
100 60
80 30
60
0
40
0 5 10 15 20 25 30 35 40 0 30 60 90 100 150 180
Cycle number Z′ (Ω)
NCA-720 °C NCA-765 °C NCA-720 °C NCA-765 °C
NCA-735 °C NCA-780 °C NCA-735 °C NCA-780 °C
NCA-750 °C NCA-750 °C
(c) (d)
Figure 3: Electrochemical measurements and characterization: (a) the initial charge-discharge profiles; (b) the initial CV curves; (c) rate
capabilities; and (d) Nyquist plots of Li/NCA half cells for the five samples.
735 which has the least Li/Ni exchange in the layered on the sintering temperature. NCA-750 achieves the highest
structure should have exhibited the best electrochemical initial discharge capacity of 187.9 mAh g-1 as well as the high-
properties. However, our electrochemical measurements est initial Coulombic efficiency of 89.5%. The CV curves with
give us some anomalous but interesting results as shown in a scan rate of 0.1 mV s-1 in Figure 3(b) show the redox behav-
the following part. iors of the transition metals as well as the phase transitions
during the charging/discharging process. All the samples
3.2. Electrochemistry. The electrochemical properties of NCA show the similar oxidation of Ni3+ to Ni4+ at around 3.8 V
electrode materials are shown in Figure 3. In Figure 3(a), the but with slight difference in the oxidation peaks which may
charge/discharge profiles are obtained at 0.1C in the range of be caused by the larger polarization for NCA-720 and
2.75-4.3 V. Obviously, the discharge capacity depends highly NCA-765. It should also be noted that NCA-780 has another
Research 5
4 5
Oxidation Oxidation 5 Oxidation
3 4
4
3
Current (mA)
Current (mA)
Current (mA)
2 3
2 2
1 1 1
0 0 0
−1 −1 −1
−2 −2 −2
Reduction −3 Reduction −3 Reduction
2.5 3.0 3.5 4.0 4.5 2.5 3.0 3.5 4.0 4.5 2.5 3.0 3.5 4.0 4.5
Voltage (V) Voltage (V) Voltage (V)
0.1 mVs−1 0.8 mVs−1 0.1 mVs−1 0.8 mVs−1 0.1 mVs−1 0.8 mVs−1
0.3 mVs−1 1.0 mVs−1 0.3 mVs−1 1.0 mVs−1 0.3 mVs−1 1.0 mVs−1
0.5 mVs−1 0.5 mVs−1 0.5 mVs−1
(a) (b) (c)
5
4.0 4 4
3.2 Oxidation Oxidation
3 3
2.4
Current (mA)
Current (mA)
2 2
Ip (mA)
1.6
0.8 1 1 Linear fit of oxidation peak
0.0 0 0
Linear fit of reduction peak
−0.8 −1
Figure 4: CV measurements. (a–e) CV curves of NCA-720, NCA-735, NCA-750, NCA-765, and NCA-780, respectively; (f) linear fit of the
peak currents and the square root of scanning rates for the five samples.
Table 2: The calculated diffusion coefficient of Li+ for antisite defects in the layered structure is beneficial to the fast
deintercalation and intercalation processes. diffusion of Li+ in the lithium layer.
The lithium-ion diffusion coefficient was calculated based
Samples O‐DLi (cm2 s-1) R‐DLi (cm2 s-1)
on the CV curves with different scan rates (Figure 4) [44]. In
NCA-720 5 79 ∗ 10‐11 2 49 ∗ 10‐11 a diffusion controlled process, the peak current is propor-
NCA-735 1 08 ∗ 10‐10 3 63 ∗ 10‐11 tional to the square root of the scanning rates, which is in line
‐10 with the Randles-Sevcik equation:
NCA-750 1 45 ∗ 10 3 66 ∗ 10‐11
NCA-765 6 54 ∗ 10‐11 3 26 ∗ 10‐11 ip = 2 69 ∗ 105 ∗ n3/2 CLi Av1/2 DLi 1/2 , 1
‐11 ‐11
NCA-780 5 51 ∗ 10 2 16 ∗ 10
where ip is the peak current, n is the number of the electrons
transferred, C Li is the concentration of Li+, A is the surface
area of the electrode, v is the scan rate, and DLi is the diffusion
oxidation peak at around 4.5 V which mainly is attributed to coefficient of lithium ions. As shown in Figures 4(a)–4(e), the
the phase transition from H1 to O1 [43]. As mentioned deintercalation O‐DLi and intercalation R‐DLi are calculated
before, the LiNiO2-based materials decompose at high tem- from the oxidation peak and reduction peak, respectively,
peratures and form the nonstoichiometric Li1‐x NiO2 phase. and the data are summarized in Table 2. The slope of the lin-
And that is the main reason for its lowest capacity in the ear fit of each peak, as shown in Figure 4(f), is positively pro-
voltage range of 2.75-4.3 V. As shown in Figure 3(c), the rate portional to the diffusion coefficient of lithium ions. Clearly,
performances differ significantly with different sintering the diffusion coefficient of the deintercalation process has a
temperatures. NCA-750 shows the best rate capability even higher dependence on temperature, which is consistent with
with a capacity of 129.8 mAh g-1 at 10C rate while it is only the rate performance, than that of the intercalation process.
119.5 mAh g-1 for NCA-735. And the EIS also shows the
results consistent with the rate performance (Figure 3(d)). 3.3. Calculation Results. In order to understand the anoma-
These results and the diffusion coefficient of lithium ions lous experimental result that NCA-750 with higher Li/Ni
(see below) indicate that an appropriate amount of Li+/Ni2+ antisite defects show better electrochemical properties,
6 Research
(b)
(c)
(d)
(e)
O
TM ions
Li+
Figure 5: Structure and diffusion path of Li+ ions of the Ni-rich cathode materials. (a) Layered structure of the Ni-rich cathode materials; the
diffusion of Li+ ions under different conditions: (b) without Li/Ni antisite defects, (c) Li ion occupies TM ion sites, (d) Ni2+ occupies Li+ sites;
(e) simultaneously with two defects near the diffusion pathway.
especially the rate capability, than that of NCA-735 which [23]. And that is the main reason for the activation barrier
contains the lowest antisite defects, we perform the calcu- in layered electrode materials. If the high-valence TM ions
lations from the perspective of activation barrier which are replaced by the Li+ ions, the activation barrier for Li+ dif-
strongly affects the diffusion of lithium ions in electrode fusion decreases due to the weaker electrostatic repulsion
materials. However, due to the limited calculation abilities, (Figure 5(c)). On the other hand, when a high-valence TM
we can only build a 120-atom supercell consisting of 2 × ion occupies Li+ sites, i.e., Li/Ni exchange (Figures 5(d)
5 × 1 LiNi0.8Co0.2O2 primitive unit cells as mentioned and 5(e)), it should also have some effect on the activation
before but not a LiNi0.8Co0.15Al0.05O2 supercell in which barrier for the six nearest lithium ions’ diffusion, which will
the five percent aluminum ions would need a huge amount be discussed in the following part in detail.
of computations. We then calculate the activation barriers in lithium
In a typical layered structure as shown in Figure 5(a), the vacancy environments, i.e., a lithium ion migrates to occupy
lithium ions and transition metal ions alternately occupy a its nearest lithium vacancy in Li slab, as shown in Figure 6.
layer of octahedron sites, which means that there is a two- This is similar to the Li diffusion in LiCoO2 [45]. Firstly,
dimensional transfer pathway for lithium ions. When a lith- we explore Li diffusion in LiNi0.8Co0.2O2 without Li-Ni
ium ion migrates to the nearest lithium vacancy, it has to exchange and find that the activation barriers are approxi-
pass through the adjacent tetrahedron sites, in which the mately 0.69 eV, as shown in Figure 6(a). Then, we consider
neighboring transition metal ion hinders the Li+ diffusion the Li diffusion in LiNi0.8Co0.2O2 with 4.17% Li-Ni exchange,
because of the strong electrostatic repulsion (Figure 5(b)) i.e., one Ni for one Li in the 120-atom supercell. The initial
Research 7
0.80
0.20 a
c
0.00 b
0.00 1.00 2.00 3.00 4.00
Reaction coordinate
(a)
0.80 Ni Li
Relative energy (eV)
0.60 Li vacancy
0.45 eV
0.40 0.64 eV
c
0.20
Li
0.80 Li vacancy Ni
Relative energy (eV)
0.60
0.25 eV 0.60 eV
0.40
c
0.20
b
0.00 a
0.00 1.00 2.00 3.00
Reaction coordinate
(c)
Figure 6: Calculation results of the activation barriers in lithium vacancy environments. (a) Energy profiles for Li diffusion and diffusion path
in LiNi0.8Co0.2O2, (b) energy profiles for Li diffusion in LiNi0.8Co0.2O2 with 4.17% Li-Ni exchange and diffusion path under Ni ion, (c) energy
profiles for Li diffusion in LiNi0.8Co0.2O2 with 4.17% Li-Ni exchange and diffusion path under Li ion. Squares are lithium vacancies, black
arrows represent diffusion path, and the blue double-sided arrows represent Li-Ni exchange.
site of a Li ion is adjacent to a Ni ion in the Li slab, and the Then, we focus on Li diffusion in LiNi0.8Co0.2O2 with dif-
diffusion path is under the Ni ion in Ni(Co) slab, as shown ferent ratios of Li-Ni exchange. Our calculations reveal that
in Figure 6(b). The calculated activation barriers decrease to the activation barrier further decreases with the increase of
be 0.45 eV, indicating that Li-Ni exchange favors the diffu- the ratio of Li-Ni exchange, as depicted in Figure 7(a). From
sion of Li. For this ratio of Li-Ni exchange, we also concern Figure 7(a), we can see that for each concerned ratio of Li-Ni
another case of Li diffusion path; that is, a diffusing Li ion exchange, the activation barrier decreases by 0.2 eV at least
is under an exchanged Li ion in the Ni(Co) slab, as shown when the diffusion path is under a Ni ion, as compared
in Figure 6(c). Surprisingly, the activation barrier decreases to that under a Li ion. Therefore, we can conclude that
to 0.25 eV. Moreover, we further investigate diffusion of Li Li-Ni exchange promotes the diffusion possibilities of Li
at the next neighboring Ni ion in the Li slab and find that in LiNi0.8Co0.2O2.
the activation barriers are almost the same as those in LiNi0.8- To understand this interesting phenomenon above, we
Co0.2O2 without Li-Ni exchange. This means that the Ni ions investigate the electronic properties of the system. The charge
in the Li slab just influence the diffusion of Li ions within the distribution for the states at the saddle points of the diffusion
scope of the neighboring region. paths for the typical cases is displayed in Figures 7(b)–7(d).
8 Research
0.70
0.50
0.40
0.30
Ni Li
0.20
0.10
0.00
0.00% 4.17% 8.33% 12.50% 16.67%
Li-Ni exchange
Path under Ni ion Saddle point
Path under Li ion
(a) (b)
Figure 7: Calculation results of the activation barriers in Li-Ni exchange environments. (a) Calculated activation barrier for Li migration with
different Li-Ni exchanges under Li or Ni ion, (b) charge distribution for Li diffusion in LiNi0.8Co0.2O2 without Li-Ni exchange, (c) charge
distribution for Li diffusion in LiNi0.8Co0.2O2 with 4.17% Li-Ni exchange and diffusion path under Ni ion, and (d) charge distribution for
Li diffusion in LiNi0.8Co0.2O2 with 4.17% Li-Ni exchange and diffusion path under Li ion.
3.0
cies, and it reduces the activation barriers of Li diffusion
2.5 around the Ni ion. This is why the Ni ions in Li slab can only
Relative energy (eV)
2.0 alter the diffusion of Li ions within the scope of the neighbor-
ing region. For the same reason, when the Li diffusion path is
1.5
under a Li ion, the above Li ion will reduce the interactions
1.0 with surrounding O ions and then the distance between O
0.5 ions and Li diffusion path, which results in the further
decrease of activation barriers.
0.0
However, as shown in Figure 8, if there are too many Li-
0.00% 4.17% 8.33% 12.50% 16.67% Ni antisite defects in Ni-rich layered materials, the total sys-
Ratio of Li/Ni exchange tem energy drastically increases, which means that the lay-
ered structure becomes quite unstable with the increasing
Figure 8: Change of total energy with different amounts of Li-Ni Li/Ni antisite defects. Therefore, it seems impossible to syn-
exchange in the supercell.
thesize the layered structure materials with so many defects
under the mentioned experimental conditions. On the other
We can see clearly from Figure 7(b) that the interactions hand, the excessive inactive TM ions in Li slab would inevita-
between Li and O ions are weak while the interactions bly cause a longer diffusion path for lithium ions. Therefore,
between Ni and O ions are strong in LiNi0.8Co0.2O2. When the highest rate of Li+ diffusion can be achieved when the
Li-Ni exchange occurs with a Ni ion appearing in the Li slab, benefit from the decrease of activation barriers and hin-
stronger interactions between Ni ion in the Li slab and O ions drance reach a balance, which corresponds to the experimen-
drive Li ions near this Ni ion diffuse to its neighboring vacan- tal result of 2.39% for NCA-750 in our study.
Research 9
[17] J. Zheng, T. Liu, Z. Hu et al., “Tuning of thermal stability in [34] H. J. Monkhorst and J. D. Pack, “Special points for brillouin-
layered Li(NixMnyCoz)O2,” Journal of the American Chemical zone integrations,” Physical Review B, vol. 13, no. 12,
Society, vol. 138, no. 40, pp. 13326–13334, 2016. pp. 5188–5192, 1976.
[18] Y. Xiao, T. Liu, J. Liu et al., “Insight into the origin of lithium/- [35] G. Mills, H. Jónsson, and G. K. Schenter, “Reversible work
nickel ions exchange in layered Li(NixMnyCoz)O2 cathode transition state theory: application to dissociative adsorption
materials,” Nano Energy, vol. 49, pp. 77–85, 2018. of hydrogen,” Surface Science, vol. 324, no. 2-3, pp. 305–337,
[19] W. Lee, S. Muhammad, T. Kim et al., “New insight into Ni- 1995.
rich layered structure for next-generation Li rechargeable [36] G. Henkelman, B. P. Uberuaga, and H. Jonsson, “A climbing
batteries,” Advanced Energy Materials, vol. 8, no. 4, article image nudged elastic band method for finding saddle points
1701788, 2018. and minimum energy paths,” The Journal of Chemical Physics,
[20] H. Xie, K. Du, G. Hu, Z. Peng, and Y. Cao, “The role of sodium vol. 113, no. 22, pp. 9901–9904, 2000.
in LiNi0.8Co0.15Al0.05O2 cathode material and its electrochem- [37] G. Henkelman and H. Jonsson, “Improved tangent estimate in
ical behaviors,” The Journal of Physical Chemistry C, vol. 120, the nudged elastic band method for finding minimum energy
no. 6, pp. 3235–3241, 2016. paths and saddle points,” The Journal of Chemical Physics,
[21] Z. Y. Li, R. Gao, J. Zhang, X. Zhang, Z. Hu, and X. Liu, “New vol. 113, no. 22, pp. 9978–9985, 2000.
insights into designing high-rate performance cathode mate- [38] L. Xiao, Y. Yang, Y. Zhao, X. Ai, H. Yang, and Y. Cao, “Synthe-
rials for sodium ion batteries by enlarging the slab-spacing of sis and electrochemical properties of submicron LiNi0.8-
the Na-ion diffusion layer,” Journal of Materials Chemistry A, Co0.2O2 by a polymer-pyrolysis method,” Electrochimica
vol. 4, no. 9, pp. 3453–3461, 2016. Acta, vol. 53, no. 6, pp. 3007–3012, 2008.
[22] Z. Huang, Z. Wang, Q. Jing, H. Guo, X. Li, and Z. Yang, [39] Y. Makimura, T. Sasaki, T. Nonaka et al., “Factors affecting
“Investigation on the effect of Na doping on structure and cycling life of LiNi0.8Co0.15Al0.05O2 for lithium-ion batteries,”
Journal of Materials Chemistry A, vol. 4, no. 21, pp. 8350–