Origin of the Phase Transition in Lithium Garnets - The Journal of

Jan 12, 2018 - Molecular dynamic and density functional theory based simulations were performed to obtain a better understanding of the origin of phas...
0 downloads 5 Views 4MB Size
Article Cite This: J. Phys. Chem. C 2018, 122, 1963−1972

pubs.acs.org/JPCC

Origin of the Phase Transition in Lithium Garnets Fei Chen,*,† Junyang Li,† Zhifeng Huang,† Ying Yang,‡ Qiang Shen,† and Lianmeng Zhang† †

State Key Laboratory of Advanced Technology for Materials Synthesis and Processing and ‡Department of Engineering Structure and Mechanics, Wuhan University of Technology, Wuhan 430070, China S Supporting Information *

ABSTRACT: Molecular dynamic and density functional theory based simulations were performed to obtain a better understanding of the origin of phase transition in garnet-type Li7La3Zr2O12 solid electrolyte. The phase transition coincided with the lithium redistribution among all sites. With the investigation of lithium distribution and dynamics, we found one temperature-dependent lithium migration pathway in lithium garnets. Lithium ions exhibited uniformly 3-dimensional diffusion in cubic LLZO, while the lithium diffusion in tetragonal LLZO was mainly in the a and b direction. The constrained diffusion in the c direction in tetragonal LLZO could be ascribed to the blocking effect of 16f sites which were found to be thermodynamically more stable than tetragonal 8a and 32g sites through density functional theory based calculations. Besides, the stabilizing effect of supervalent doping on cubic phase was also studied through Ta-doped LLZO. Further site occupancy investigations indicated that supervalent doping introduced lithium vacancies and reduced octahedral sites (tetragonal 16f and 32g sites) occupancy due to site energy preference. The reduced octahedral sites occupancy weakened the blocking effect of tetragonal 16f sites, promoted the lithium redistribution, and eventually lowered the phase transition temperature.

1. INTRODUCTION With the introduction of solid electrolyte, all-solid-state batteries obtain great improvements in safety and energy density compared to flammable liquid organic electrolyte-based lithium ion batteries. Acting as a combination of electrolyte and separator, the solid electrolyte offers mechanical strength to prevent lithium dendrite growth, high lithium ionic conductivity, and chemical and electrochemical stability in contact with lithium anode and most cathode materials.1−3 Among the reported candidates, lithium garnets exhibit capable performance with ionic conductivity varying from 10−6 to 10−3 S cm−1 and excellent stability with metallic lithium.4,5 One of the lithium garnet families being widely studied is Li7La3Zr2O12 (LLZO) with cubic structure (space group Ia3̅d).5 Dodecahedral LaO8 and octahedral ZrO6 form the 3D framework by edge-sharing. Lithium ions are randomly and partially distributed in tetrahedral 24d sites and octahedral 96h sites as shown in Figure 1c. One LiO4 tetrahedron connects with four LiO6 octahedrons. One LiO6 octahedron is surrounded by two LiO4 tetrahedrons. The LiO4 tetrahedron and LiO6 octahedron are connected with each other by facesharing and finally form a 3D continuous lithium transport pathway as shown in Figure 1e. However, pure LLZO generally crystallize to thermodynamically more stable tetragonal phase (space group I41/acd) at room temperature.6 The transition from cubic to tetragonal phase is called tetragonal distortion, as shown in Figure 1f, which transforms the cubic 24d sites into fully occupied 8a sites and unoccupied 16e sites. Cubic 96h © 2018 American Chemical Society

sites are also transformed into two fully occupied 16f and 32g sites.7,8 This tetragonal distortion has little influence on the framework but greatly affects lithium distribution and lowers the ionic conductivity by 2 orders of magnitude.5,6 Several studies from both experimental and theoretical investigations had been focusing on the stabilization of the cubic phase. The most efficient way to obtain room temperature stable cubic phase is introducing lithium vacancy into a lithium sublattice by supervalent doping, i.e., Al3+ (Li site),9−11 Ga3+ (Li site),12−14 and Ta5+ (Zr site),15,16 as summarized in a recent review.17 This method has been widely used in the preparation of cubic lithium garnets. Jalem et al.18 performed molecular dynamic simulations to study the cubic phase stabilizing effects of gallium. Their results showed that gallium doping can reduce the phase transition temperature of garnets. Cubic phase could be stable in room temperature when phase transition temperature was lower than room temperature. A similar result was also obtained by Bernstein et al.7 in the determination of effect of composition on phase transition of LLZO. The decrease of phase transition temperature can be regarded as a sign of the stabilization of cubic phase. Thus, knowing the mechanism of the phase transition from tetragonal to cubic phase in lithium garnets is necessary for a Received: November 5, 2017 Revised: January 10, 2018 Published: January 12, 2018 1963

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C

lithium ordering by performing molecular dynamic simulations. Such closed cyclic diffusion did not make any contribution to ionic conductivity. By applying first-principles-based calculations, Meier et al.27 revealed two different Li-ion migration mechanisms. In tetragonal LLZO, high ordering of Li ion results in a fully collective motion, which requires higher activation energy. In cubic LLZO, single-ion jump mechanism brings much lower energetic cost and higher ionic conductivity compared with the fully collective mechanism in tetragonal LLZO. Similar concerted migration mechanisms were also obtained by Jalem et al.28 in their FPMD simulations. Chen et al.29 performed molecular dynamic simulation to elucidate the lithium diffusion mechanism in LLZO. Their results showed that lithium exhibits uncorrelated Poisson-like diffusion in the cubic LLZO and correlated diffusion in the tetragonal LLZO. They further confirmed that the uncorrelated diffusion is resulted from the weak site dependence of lithium diffusion in the cubic LLZO. Morgan30 studied the correlation effects for lithium garnets by lattice-gas Monte Carlo simulations and found that the correlation in lithium garnets is reproduced from both site energy differences and nearest-neighbor repulsion. Furthermore, the ionic conductivity can be maximized by tuning the mobile-ion stoichiometry. The above studies26−30 indicated that different lithium distributions could greatly influence the ionic conductivity by exhibiting different lithium migration mechanism. In general, the studies of the phase transition of LLZO were concentrated in the lithium redistribution during the phase transition from tetragonal to cubic phase and the effects of lithium distribution on ionic conductivity, but the investigations about the energetic origin of the lithium redistribution and how lithium distribution affect the migration mechanism were still empty. In this paper, we performed molecular dynamic simulations and density functional theory based calculations to reveal the origin of the phase transition in Li7La3Zr2−xTaxO12 (x = 0−1, LLZ-Tax). The relation between the lithium distribution and the lithium migration mechanism during the phase transition of LLZO was first studied. Then, the lithium sites energy differences were calculated to understand the energetic origin of the lithium redistribution during the phase transition. Besides, the effects of supervalent dopants on lithium distribution and phase transition temperature in doped-LLZO were also characterized.

Figure 1. Crystal structure of (a) tetragonal LLZO, (b) transformed tetragonal LLZO, and (c) cubic LLZO. Lithium arrangement in (d) tetragonal and (e) cubic LLZO. (f) The loop structure of lithium arrangement in both tetragonal and cubic LLZO.

better understanding of the stabilization of cubic phase and the further design of lithium garnets with higher ionic conductivity. Previous reports19−24 showed that lithium distribution changed from order to disorder during the phase transition from tetragonal to cubic phase in lithium garnets. Adams’s group captured the whole transition process of LLZO from both neutron diffraction and molecular dynamic simulations with their own Morse-type force field.21,25 Their results showed an approximate agreement and revealed that tetrahedral site occupancies increased with increasing temperature, while octahedral site occupancies decreased in the meantime. Bernstein et al.7 performed first-principles-based molecular dynamic simulations to investigate the phase transition process of LLZO. Only the simple shift from 8a to 16e sites was observed during the phase transition. Wang and co-workers24 performed neutron diffraction and molecular dynamic calculations to study the local lithium distribution in the Li7−xLa3Zr2−xTaxO12 (x = 0−2) series. Their results showed an increasing tendency of lithium ordering as lithium content increased. The local lithium distribution suggested that tetrahedral sites lithium occupancy was generally lower than that obtained from neutron diffraction and the exclusion principle. The X-ray and neutron diffraction performed by Logéat et al.23 also showed that lithium ions were promoted from tetrahedral sites to octahedral sites during Ta doping which might reduce electrostatic repulsion. In the view of thermodynamics, the phase transition in LLZO is an entropydriven process,19,20 during which lithium ions were distributed among all sites. Supervalent doping reduced lithium content and disorder lithium sublattice,12,21 which promoted the shift of lithium distribution from order to disorder and eventually reduced the free energy.7 Efforts have been made to study the relationship between lithium distribution and ionic conductivity. Burbano et al.26 ascribed the low ionic conductivity of tetragonal LLZO to concerted ionic motion around a closed loop due to the strong

2. METHODS 2.1. Interatomic Potentials. Empirical interatomic potentials consisting of long-range Coulombic interaction and shortrange Buckingham interaction, as given by eq 1, were used to describe atomic interactions. ⎛ r ⎞ C qiqj ij ij Φij(rij) = Aij exp⎜⎜ − ⎟⎟ − 6 + 4π ϵ0rij ⎝ ρij ⎠ rij

(1)

where Aij, ρij, and Cij are empirical force-field parameters for interacting i and j ions, ϵ0 is the free space permittivity, rij is the distances between i and j ions, and qi represents the charges of i ion. All the force-field parameters except for the Buckinghamtype potential parameters of Ta−O pair were taken from the literature.31 Interactions parameters between Ta and O were fitted against the experimental structure of Ta-doped LLZO24 by the GULP package.32 The force-field parameters used in this work are listed in Table 1. The cutoff distances were set to 12 Å 1964

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C ⎡ 1 ⎤ D = lim ⎢ ⟨[ r (⃗ t )]2 ⟩⎥ ⎦ t →∞⎣ 2dt

for all short-range Buckingham interaction. The Ewald summation were used to summarize long-range Coulombic interaction.

where d is the dimension of diffusion. According to above equation, when it reaches diffusive regime, the slope of log⟨[r(⃗ t)]2⟩ versus log t should be 1. For these simulations in this study who reaches diffusive regime, D is calculated by the linear fitting to the dependence of ⟨[r(⃗ t)]2⟩ over 2dt. 2.4. Density Functional Theory Based Calculations. First-principles calculations were performed in 1 × 1 × 1 cell using the Perdew−Burke−Ernzerhof generalized gradient approximation (GGA-PBE)34 and the Projector Augmented Wave (PAW)35 method via the Vienna ab initio simulation package (VASP).36 The kinetic energy cutoff of 520 eV and 4 × 4 × 4 k-points mesh was used for all the calculations. All structures were relaxed in two steps. In the first step, the cell parameters were relaxed with fixed fractional coordinates. In the second step, the atoms were allowed to relax with fixed cell parameters. The energy of the relaxed structures were then calculated.

Table 1. Force-Field Parameters Buckingham potential parameters species

charge (e)

A (eV)

ρ (Å)

C (eV Å6)

Li−O La−O Zr−O Ta−O O−O

1.00 2.50 2.65 3.65 −1.65

1087.29 2075.26 1650.32 713.831 4870.00

0.260 0.326 0.311 0.359 0.267

0.00 3.25 5.10 0.36 77.00

2.2. Lithium Arrangements. To balance the accuracy and computational efficiency, a 2 × 2 × 2 supercell containing 576 available Li sites were used in this study.31 In order to obtain the natural distribution of lithium ions, the disordered lithium distribution of the initial model was obtained from the ordered lithium distribution from thermal relaxation. Tetragonal and cubic LLZO shared the same framework structure, but their crystal structure was a little bit different. For the convenience of the simulation of the whole phase transition process, tetragonal LLZO (Figure 1a) was transformed into one cell (Figure 1b) whose framework had the same fractional coordinates with that of the cubic phase (Figure 1c). After this symmetric operation, Ta ions uniformly and symmetrically distributed in Zr site and replaced one neighbor octahedral site Li ion (32g or 16f) with vacancy simultaneously to maintain electrical neutrality. 2.3. Molecular Dynamic Simulations. Molecular dynamic simulations were performed from 300 to 1400 K with 2 fs time step and 1 atm pressure using the GULP package.32 In the first step, initial models were heated under constant number, pressure, and temperature (NPT) ensemble for 50 ps to determine cell sizes at target condition. In the second step, using cell sizes obtained from NPT simulations, we performed constant number, volume, and temperature (NVT) ensemble simulations for 1000 ps to obtain dynamic features. Thermostat and barostat were used with a relaxation time of 0.05 and 0.25 ps, respectively. Trajectories were sampled with a time step of 0.1 ps for NPT run and 1 ps for NVT run. Ionic conduction and local structure were analyzed from atomic trajectories with number density, pair distribution function (PDF33), and mean-squared displacement (MSD). Number density was calculated by dividing the cell into a 300 × 300 × 300 mesh grid. Then, the number of Li ions in each grid were counted and averaged over all time steps. To ensure that every Li was counted in the evaluation of the occupancies of lithium in different symmetrical position, the distances between Li and every stable Li sites in cubic LLZO (24d and 96h) were calculated in fractional coordinate. Li would be labeled as its nearest stable Li site. Occupancies were the number of Li averaged over all time steps and the number of stable sites. MSD was calculated over all simulation times from the following definition: ⟨[ r (⃗ t )]2 ⟩ =

1 N

3. RESULTS AND DISCUSSION 3.1. Phase Transition and Lithium Conduction Modeling. To better understand the phase transition from tetragonal to cubic phase and the stabilization of cubic phase, molecular dynamics simulations of LLZO with Ta dopants concentration varying from 0 to 2 were performed with empirical force-field parameters fitted against experimental structures as mentioned in section 2.1. The temperature dependence of lattice parameters of Ta-doped LLZO captured from NPT simulations are shown in Figure 2a. The phase transition process of LLZO can be captured from the changes of lattice parameter. When phase transition from tetragonal to

N

∑ ⟨[ ri(⃗ t + t0)]2 i=1

− [ ri(⃗ t0)]2 ⟩

(3)

Figure 2. (a) Temperature dependence of lattice parameter of Li7−xLa3Zr2−xTaxO12 (x = 0−2). (b) Vegard’s law. Experimental lattice parameter are taken from Wang,24,40 Thompson,38 Buschmann,39 Hamao,37 and Li15. (c) Mean-squared displacement of lithium ion at temperatures from 300 to 1400 K with an interval of 100 K. (d) A log−log plot of mean-squared displacement.

(2)

where N is the number of atoms. The self-diffusivity was defined as 1965

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C

Figure 3. (a) Li−Li pair distribution function as a function of temperature. (b) Physically reasonable lithium arrangements.41 Local lithium arrangements in (c) cubic and (d) tetragonal LLZO derived from 1400 and 300 K NPT run, respectively.

temperature range. Consequently, the MSD data of 300−500 K are not used to calculate the self-diffusivity. In general, the empirical force-field parameters used in this work are able to capture the phase transition process and obtain lithium ion diffusion properties. With Ta doping, the stabilization process of cubic phase is also observed. In the following sections, we will focus on the relationship between lithium distribution, diffusion pathway and the phase transition in LLZO first. Then, the effect of Ta doping on the stabilization of cubic phase are investigated on the basis of the understanding of the phase transition. 3.2. Lithium Distribution in LLZO. It can be confirmed that the phase transition of LLZO is accompanied by the lithium redistribution among all tetrahedral and octahedral sites. In order to understand lithium distribution in LLZO, Li− Li pair distribution function are calculated for all the simulations as shown in Figure 3a. Similar approaches have been made by Adams et al.21 with their own Morse-type force field and Wang and Klenk19,20 with a Buckingham-type interatomic potential. Li−Li pairs are separated with each other in a distance larger than 2 Å, which is the same as Adams et al.21 obtained. This distribution feature results from the electrostatic repulsion suggests that two neighboring sites in a distance smaller than 2 Å, i.e., two distorted octahedral sites in one octahedron, tetrahedral sites, and its nearest distorted octahedral sites, cannot be occupied simultaneously. Physically reasonable lithium arrangements, as described by O’Callaghan et al.,41 are shown in Figure 3b: I, two occupied tetrahedral sites and unoccupied octahedral site; II, one occupied tetrahedral site and an occupied distorted octahedral site away from the tetrahedral site; III, two unoccupied tetrahedral sites and an occupied octahedral site. Combining the above three basic arrangements, all possible local lithium arrangements that contribute to Li−Li pair distribution function are shown in Figure 3c,d. These local arrangements contain all the cases where lithium distance is less than 5 Å. It should be noted that in cubic structure at 1400 K (Figure 3c) arrangement I offers a Li−Li distance of 4.049 Å, and Li−Li PDF shows that there is a

cubic phase occurs, the lattice parameters transform from tetragonal feature (a = b ≠ c) to cubic feature (a = b = c). The rapidly merging of lattice parameters between 850 and 950 K can be regarded as a sign of the phase transition from tetragonal to cubic phase. Results show that the phase transition temperature decreases with the increasing Ta dopants concentration. When Ta dopants concentration is larger than 0.25, no tetragonal phase is observed. The phase transition temperature is so sensitive to lithium concentration that the stabilization process of the cubic phase is observed in a small concentration windows. It should be noted that due to the inherent limitation of the Buckingham-type interatomic potential, it is hard to get the right phase transition temperature and also the critical Ta dopants concentration for the stabilizing of cubic phase. Figure 2b shows the linear dependence of lattice parameters on dopants concentration, known as the Vegard’s law. Our simulations match well with the experimental lattice parameters.15,24,37−40 The occurrence of phase transition can also be seen from diffusion properties of lithium ions. Figure 2c shows the MSD of lithium ions from NVT simulations range from 300 to 1400 K with an interval of 100 K in LLZO. According to eq 3, the slope of MSD is 6D for 3-dimensional diffusion. The larger the slope, the greater the self-diffusivity. The negligible MSD of tetragonal phase represents its low diffusivity compared to that of cubic phase. The MSD of 900 K locates at the middle part of MSD of tetragonal and cubic phase. The data point that is clearly different from the point around and changes rapidly help us locate the phase transition process. Figure 2d shows the dependent on log⟨[r(⃗ t)]2⟩ on log t. As mentioned in the Methods section, when it reaches the diffusive regime, the slope of log⟨[r(⃗ t)]2⟩ versus log t should be 1. The simulation time used in this paper is long enough to reach diffusive regime for high temperature simulations. At lower temperature, diffusive regime can be reached in the last hundreds of picoseconds of the whole simulation time. For the simulations at 300−500 K, MSD keep constant values throughout the simulations, suggesting the oscillation behavior of lithium ions at this 1966

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C

Figure 4. Lithium density map derived from (a, e) 500 K, (b, f) 700 K, (c, g) 900 K, and (d, h) 1100 K NVT run. Lithium ions at 8a, 16f, and 32g sites are colored red, blue, and green, respectively. The observed migration pathway in tetragonal phase is marked in solid black line in (b) and (f).

Figure 5. (a) Temperature dependence of site occupancies in LLZO. (b) Local structure of lithium network. (c) Arrhenius plot of the total diffusivity and the diffusivity in ab plane and c direction.

valley. As a result, the first two combinations, I + I and I + II, should not exist or be relatively few; Li−Li PDF is mainly contributed from the last four combinations: (II + II)1, (II + II)2, II + III, and III + III. After that from a combination standpoint, there is no difference in local lithium arrangement between cubic and tetragonal phase: two structures both come from the combination of arrangements II and III. This result corresponds to statistical result for octahedral centered clusters performed by Klenk et al.20 To figure out the difference of Li−Li PDF between tetragonal and cubic phase, the distances between lithium ions in all local arrangement are measured and marked in Figure 3. In the distance less than 5 Å, cubic phase has two Li−Li PDF peaks while tetragonal has three peaks. The first peak of cubic PDF is contributed by a large amount of Li−Li pairs; the distances of these pairs are 2.389, 2.539, 2.704, 2.938, 3.396, 3.092, and 3.356 Å. During the tetragonal distortion of cubic phase, the first five distances change to 2.532, 2.683, and 2.645 Å, which contributes to the first peak of tetragonal PDF. 3.092 Å increases to 3.696 Å and contributes to the second peak of tetragonal PDF. Because of the splitting of angle between two arrangements II in combination (II + II)1, 3.566 Å splits into 3.794 and 4.473 Å, which contributes to the second and third peak of the tetragonal peak. Unlikely, Adams et al.21 observed two peaks (2.35 and 2.8 Å) in the Li−Li PDF in

tetragonal phase, between which there is a valley. However, in our measurements, no Li−Li pair with a distance of 2.35 and 2.8 Å is observed. In general, our simulations match well with experimental structure refined from neutron diffraction24 and suggest that (i) the splitting and widening of PDF peaks result from the distortion of lithium ring caused by tetragonal distortion and (ii) the first peak of Li−Li PDF moves to larger distance during temperature decrease from 1400 to 300 K. This result indicates that tetragonal structure can offer a distribution that allows lithium ions separate from each other in a larger distance, which may reduce Li−Li repulsion and make tetragonal structure more thermodynamically stable. To intuitively observe the lithium distribution in the crystal structure, lithium density maps in different temperatures are shown in Figure 4. While the differences of density map between adjacent temperatures are relatively small, only four temperature points with typical feature are shown, which are 500, 700, 900, and 1100 K. Lithium ion located at different symmetrical sites are marked with different colors. The unoccupied 16e site is surrounded by two 16f and two 32g sites and marked as pink hollow circle in Figure 4a,e. At 500 K (Figure 4a,e), lithium ions prefer to locate at their stable sites, and only oscillations are observed. It is worth noting that in this temperature 32g site lithium ions may jump to 16e sites, stay for a little while, and then jump back. When temperature 1967

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C increase to 700 K, the overlapping of some 8a (red part) and 32g sites (green part) is observed. The overlapping sites form a continuous pathway throughout the simulation box as marked by black solid line in Figure 4b,f. Such pathway is considered to be lithium migration pathway. Furthermore, these migration pathways exhibit spiral structure whose axis appears to be only in the a and b direction. As temperature increases to 900 K (Figure 4c,g), at which the phase transition begins, all stable sites including 16f sites are connected and preliminarily form a 3-dimensional lithium migration network. After the phase transition has occurred, lithium ions distribute more evenly over each site as mapped in Figure 4d,h. From above discussion, we can conclude that (i) the phase transition has little influence on lithium local arrangements but affects the Li−Li pair distance a lot and (ii) phase transition from tetragonal to cubic is accompanied by the lithium redistribution among all lithium sites rather than from tetragonal 8a sites to tetragonal 16e sites.7,20 Moreover, there may exist different lithium migration pathways before and after the phase transition. In the next section, the evolution of lithium migration pathway with phase transition from tetragonal to cubic will be analyzed and mapped by combining site occupancy and local structure. 3.3. Evolution of Lithium Migration Pathway in LLZO. Figure 5a shows the temperature dependence of occupancy of different sites. The local structure of lithium ions in tetragonal phase is shown in Figure 5b. Combining occupancy and local structure, the evolution of lithium migration pathway can be divided into four stages by temperature. In the first stage, from 300 to 500 K, 8a and 16f sites remain fully occupied, indicating their oscillation feature in this stage. This can also be visualized from the density map in Figure 4a,e. As mentioned in the last section, lithium ions shift from 32g sites to 16e sites. This shift forms a local migration between 32g-16e-32g as mapped in green in Figure 5b but does not contribute to ionic conductivity. In the second stage, there is a shift of lithium ions from 8a and 32g to 16e. The vacancies in 8a, 32g, and 16e sites are responsible for the forming of lithium migration pathway. The local structure of lithium network (Figure 5b) shows that these three sites can be connected into a continuous pathway, 8a-32g16e-32g-8a (the yellow part). Ab-initio molecular dynamics simulation by Meier et al.27 shows a similar migration pathway in tetragonal LLZO with which we described in Figure 6c. However, the migration pathway they proposed (16f-32g-8a32g-16f) cannot form a continuous pathway, which may be because of the neglect of 16e sites. The spiral structure, or helix, of the above pathway is shown in Figure 6d. Figure 6b,c shows the side view of two pathways in Figure 6a. In the direction of the axis of the helix, 8a, 32g, and 16e sites are arranged in circles as shown in Figure 6a and marked with blue solid line. One helix can be considered as 1-dimensional lithium migration pathway and contributes to the ionic conduction in the axis direction. The axis of the helix is parallel to a or b direction and is perpendicular to c direction. The helix-pathways in a and b directions are perpendicular to each other and intersect in 8a sites. The lithium ions in two helix pathways can be exchanged in 8a sites, which is the only contribution to ionic conduction in the c direction. To compare the differences of diffusivity of lithium ions between the c direction and ab plane, the diffusivities of them were calculated from the slope of MSD (as shown in Figure S1 of the Supporting Information) and plotted as Figure 5c. As shown, in this stage, ionic conduction is mainly

Figure 6. Lithium arrangement in half-cell and the schematic diagram of migration pathway in the a direction (solid black line) and b direction (solid blue line) from (a) (010) and (b, c) (001) crystal face in tetragonal LLZO. (b) and (c) differ in the choosing of half-cell. (d) Schematic diagram of the 8a-32g-16e-32g-8a helix migration pathway in tetragonal LLZO.

contributed by the conduction in ab plane. The ionic conduction in the c direction is relatively small but nonnegligible and is blocked by 16f sites. Thus, the lithium conduction in tetragonal can be considered as 2-dimensionallike conduction. In the third stage, when 16f sites finally participate in lithium conduction, a 3-dimensional lithium migration network is formed since there are no blocking ions. Density maps (Figure 4c,g) confirm the formation of a 3-dimensional lithium migration network. Occupancies show that phase transition and lithium redistribution from 16f sites to other sites begin simultaneously. Meanwhile, the gap of diffusivity between the ab plane and c direction is getting smaller as shown in Figure 5c, suggesting the blocking effect of 16f sites may be the key factor that prevents the phase transition from the tetragonal to cubic phase. In the final stage, the lithium redistribution among all lithium sites is completed. The occupancies of tetragonal 8a and 16e sites as well as 16f and 32g sites become the same in value. The first two sites merge into cubic 24d sites, and the last two sites merge into cubic 96h sites. This corresponds to the experimental refined structure. Figure 5c shows that lithium conduction in the ab plane and c direction has no difference in this stage. To obtain the thermodynamic origin of the differences of the temperature that lithium participates in conduction between different sites, DFT simulations of lithium defective LLZO model is performed. The lithium vacancy is obtained by substitute Zr4+ with Ta5+. Four lithium defective models, which are denoted as V8a, V16f, V32g‑near, and V32g‑far (Figure 7a), are built in consideration of the relative position of lithium ion and Ta5+. The energies of the defective models are shown in Figure 7b. The higher the energy of the defected structure, the harder it is to form. Since 16f vacancy has higher energy, lithium vacancy should prefer to locate at other sites, like 32g. Thus, the order of stability is 16f > 8a > 32g, which explains the temperature differences that lithium participates in conduction between different sites. Since the tetragonal phase and cubic phase differ only in the lithium distribution, the phase transition is actually one process of lithium redistribution from tetragonal distribution to cubic 1968

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C

parameter, Ta doping shifts the whole occupancy evolution to low temperature and keeps the order of participation of different sites in lithium conduction (32g > 8a > 16f), which confirms the four diffusion stages mechanism discussed above. Besides, the phase transition temperature is sensitive to lithium content and decreases rapidly with the decreasing lithium content. It also suggests that cubic phase can be obtained when Ta dopants concentration fall into a small range from 0.125 to 0.25, which corresponds to experimental fact that a small amount of dopants can stabilize the cubic phase.14,23,40 Figure 9a shows the change of occupancy during the stabilization of cubic phase. Ta doping can lower the occupancy of 16f sites and 32g sites simultaneously. The merging of occupancies is also observed. After the merging of occupancies, the redistribution among tetrahedral and octahedral sites are still going on. Figure 9b shows the occupancies of tetrahedral and octahedral sites as a function of Ta dopants concentration. Red and green area represent the permitted occupancy range21 of tetrahedral and octahedral sites derived from electrostatic repulsion,41 respectively. Both neutron diffraction5,6,24,42−46 and our molecular dynamic work show a strong trend that octahedral site occupancy exhibits the minimum value of permitted occupancy. It should be noted that the permitted minimum occupancies of octahedral sites are calculated with the assumption that lithium ions prefer to locate at tetrahedral sites. The increasing tetrahedral site occupancies can also push lithium ions from the undistorted octahedral site to distorted octahedral sites due to electrostatic repulsion. Thus, when the occupancy of tetrahedral sites is larger than 0.5, all lithium ions

Figure 7. (a) Schematic diagram and (b) the relative energy of four lithium defective model.

distribution. The evolution of lithium distribution and migration pathway is derived from the site stability. With the fact that 32g and 8a site lithium does not affect the linear dependence of lattice parameter on temperature, while 16f site lithium promotes the merging of lattice parameters, we believe that 16f sites affect the phase transition from tetragonal to cubic phase by blocking the lithium redistribution in the c direction. 3.4. Effect of Ta Doping on Lithium Distribution. From above discussions, we known that the phase transition in garnets is actually the redistribution of lithium atoms. In this section, the effect of supervalent doping on lithium distribution is investigated. Figure 8 shows the evolution of occupancies of LLZ-Tax with different Ta content as a function of temperature. According to the separation of different diffusion stage as what has been done in section 3.3, the occupancies are colored differently. Similar with the temperature dependence of lattice

Figure 8. Temperature dependence of lithium site occupancy in Li7−xLa3Zr2−xTaxO12. (a) x = 0.03125, (b) x = 0.0625, (c) x = 0.125, (d) x = 0.25, (e) x = 0.375, and (f) x = 0.5. 1969

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C

4. CONCLUSION In summary, combining molecular dynamics and density function theory based calculations, we investigated the origin of phase transition from tetragonal to cubic phase and the stabilization of cubic phase in lithium garnets. First, temperature dependence of lattice constants and mean-squared displacement were employed to confirm that our simulations were able to capture the phase transition and model lithium conduction. Second, lithium distributions were explored via the Li−Li pair distribution function and lithium density map. The Li−Li pair distribution function showed that the phase transition did not change the lithium local arrangements. Moreover, tetragonal phase could provide a distribution scheme with enlarged Li−Li distance. The lithium density map described to us a 2-dimensional-like lithium migration pathway in tetragonal phase. Third, combining with lithium local structure, the structure and the formation of the different migration pathway were carefully analyzed from the temperature dependence of site occupancies. We divided the evolution of lithium migration pathway into four stages: oscillation, local migration, 2-dimensional-like migration, and 3-dimensional migration. DFT results showed that the evolution was the result of sites energy difference. The blocking effect of 16f site lithium ions was the main factor that hampered the transition from the tetragonal to cubic phase. Finally, the effects of supervalent doping on lithium distribution were investigated. The phase transition temperature was sensitive to lithium content and decreased rapidly with the decreasing lithium content. Ta doping also promoted the shift of lithium from octahedral sites to tetrahedral sites, which turned out to be the result of site energy preference. This work deepens the understanding about the phase transition and the stabilization of cubic phase in lithium garnets.

Figure 9. Lithium site occupancies at (a) 400 K and (b) 1400 K as a function of Ta dopant concentration. Red and green areas represent the permitted occupancy range of tetrahedral and octahedral sites derived from electrostatic repulsion, respectively. Open symbols are taken from neutron diffraction data in ref 5, 6, 24, and 42−46. Filled symbols represent the average occupancies of cubic tetrahedral and octahedral sites.

only locate at the tetrahedral sites and the distorted octahedral sites (see Figure 3b41). The assumption corresponds to our DFT results that 8a sites (tetrahedral sites) are more stable than 32g sites (distorted octahedral sites). Hence, 16f sites (undistorted octahedral sites) are responsible for the formation of tetragonal phase by blocking lithium diffusion in the c direction. Figure 10 shows the temperature dependence of lithium ion diffusivity in Ta-doped LLZO. Ta doping could improve



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.7b10911. Figure S1: the mean-square displacement of lithium ions in LLZO along a, b, and c directions (PDF)



AUTHOR INFORMATION

Corresponding Author

Figure 10. Arrhenius plot of the self-diffusivity of Ta-doped LLZO as a function of Ta dopant concentration.

*E-mail: [email protected] (F.C.). ORCID

Fei Chen: 0000-0001-9643-7191 Junyang Li: 0000-0002-1726-0656

lithium diffusivity by lowering phase transition temperature and keeping the linear trend of diffusivity in Arrhenius plot to room temperature. Further increasing the Ta dopants concentration will decrease the diffusivity. The ionic conductivity is the result of the combination of lithium concentration, lithium diffusivity, and the correlation effect. Our results can also be used to explain the stabilization of cubic phase by low-valent doping like Y3+.47,48 The increased lithium concentration derived from low-valent doping forces lithium ions to fill in the unoccupied 16e sites. This increases the occupancies of tetrahedral sites and pushes lithium ions from octahedral sites to distorted octahedral site, which weakens the blocking effect of octahedral sites.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work is financially supported by the National Natural Science Foundation of China (No. 51472188 and 51521001), Natural Research Funds of Hubei Province (No. 2016CFB583), Fundamental Research Funds for the Central Universities in China, State Key Laboratory of Advanced Electromagnetic Engineering and Technology (Huazhong University of Science and Technology), National Key Research and Development Program of China (No. 2017YFB0310400), and the “111” project (No. B13035). 1970

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C



(21) Adams, S.; Rao, R. P. Ion Transport and Phase Transition in Li7−xLa3(Zr2−xMx)O12 (M = Ta5+, Nb5+, X = 0, 0.25). J. Mater. Chem. 2012, 22, 1426−1434. (22) Matsui, M.; Sakamoto, K.; Takahashi, K.; Hirano, A.; Takeda, Y.; Yamamoto, O.; Imanishi, N. Phase Transformation of the Garnet Structured Lithium Ion Conductor: Li7La3Zr2O12. Solid State Ionics 2014, 262, 155−159. (23) Logéat, A.; Köhler, T.; Eisele, U.; Stiaszny, B.; Harzer, A.; Tovar, M.; Senyshyn, A.; Ehrenberg, H.; Kozinsky, B. From Order to Disorder: The Structure of Lithium-Conducting Garnets Li7‑xLa3TaxZr2‑xO12 (X = 0−2). Solid State Ionics 2012, 206, 33−38. (24) Wang, Y.; Huq, A.; Lai, W. Insight into Lithium Distribution in Lithium-Stuffed Garnet Oxides through Neutron Diffraction and Atomistic Simulation: Li7‑xLa3Zr2‑xTaxO12 (X = 0−2) Series. Solid State Ionics 2014, 255, 39−49. (25) Rao, R. P.; Gu, W.; Sharma, N.; Peterson, V. K.; Avdeev, M.; Adams, S. In Situ Neutron Diffraction Monitoring of Li7La3Zr2O12 Formation: Toward a Rational Synthesis of Garnet Solid Electrolytes. Chem. Mater. 2015, 27, 2903−2910. (26) Burbano, M.; Carlier, D.; Boucher, F.; Morgan, B. J.; Salanne, M. Sparse Cyclic Excitations Explain the Low Ionic Conductivity of Stoichiometric Li7La3Zr2O12. Phys. Rev. Lett. 2016, 116, 1−6. (27) Meier, K.; Laino, T.; Curioni, A. Solid-State Electrolytes: Revealing the Mechanisms of Li-Ion Conduction in Tetragonal and Cubic LLZO by First-Principles Calculations. J. Phys. Chem. C 2014, 118, 6668−6679. (28) Jalem, R.; Yamamoto, Y.; Shiiba, H.; Nakayama, M.; Munakata, H.; Kasuga, T.; Kanamura, K. Concerted Migration Mechanism in the Li Ion Dynamics of Garnet- Type Li7La3Zr2O12. Z. Phys. Chem. 2012, 226, 525−537. (29) Chen, C.; Lu, Z.; Ciucci, F. Data mining of molecular dynamics data reveals Li diffusion characteristics in garnet Li7La3Zr2O12. Sci. Rep. 2017, 7, 40769. (30) Morgan, B. J. Lattice-Geometry Effects in Garnet Solid Electrolytes: A Lattice-Gas Monte Carlo Simulation Study. R. Soc. Open Sci. 2017, 4, 170824. (31) Klenk, M. J.; Lai, W. Finite-Size Effects on the Molecular Dynamics Simulation of Fast-Ion Conductors: A Case Study of Lithium Garnet Oxide Li7La3Zr2O12. Solid State Ionics 2016, 289, 143− 149. (32) Gale, J. D.; Rohl, A. L. The General Utility Lattice Program (GULP). Mol. Simul. 2003, 29, 291−341. (33) Kitamura, N. Atomic-Configuration Modeling of Ion-Conducting Crystalline Oxide by Diffraction Technique and Theoretical Calculation. J. Ceram. Soc. Jpn. 2015, 123, 637−642. (34) Perdew, J.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (35) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953−17979. (36) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (37) Hamao, N.; Kataoka, K.; Kijima, N.; Akimoto, J. Synthesis, Crystal Structure and Conductive Properties of Garnet-Type Lithium Ion Conductor Al-Free. J. Ceram. Soc. Jpn. 2016, 124, 678−683. (38) Thompson, T.; Sharafi, A.; Johannes, M. D.; Huq, A.; Allen, J. L.; Wolfenstine, J.; Sakamoto, J. A Tale of Two Sites: On Defining the Carrier Concentration in Garnet-Based Ionic Conductors for Advanced Li Batteries. Adv. Energy Mater. 2015, 5, 1500096. (39) Buschmann, H.; Berendts, S.; Mogwitz, B.; Janek, J. Lithium Metal Electrode Kinetics and Ionic Conductivity of the Solid Lithium Ion Conductors “Li7La3Zr2O12” and Li7−xLa3Zr2−xTaxO12 with GarnetType Structure. J. Power Sources 2012, 206, 236−244. (40) Wang, Y.; Lai, W. High Ionic Conductivity Lithium Garnet Oxides of Li7−xLa3Zr2−xTaxO12 Compositions. Electrochem. Solid-State Lett. 2012, 15, A68−A71. (41) O’Callaghan, M. P.; Cussen, E. J. Lithium Dimer Formation in the Li-Conducting Garnets Li(5+x)Ba(x)La(3‑x)Ta2O12 (0 < X < or = 1.6). Chem. Commun. (Cambridge, U. K.) 2007, 12, 2048−2050.

REFERENCES

(1) Kamaya, N.; Homma, K.; Yamakawa, Y.; Hirayama, M.; Kanno, R.; Yonemura, M.; Kamiyama, T.; Kato, Y.; Hama, S.; Kawamoto, K.; Mitsui, A. A Lithium Superionic Conductor. Nat. Mater. 2011, 10, 682−686. (2) Knauth, P. Inorganic Solid Li Ion Conductors: An Overview. Solid State Ionics 2009, 180, 911−916. (3) Thangadurai, V.; Narayanan, S.; Pinzaru, D. Garnet-Type SolidState Fast Li Ion Conductors for Li Batteries: Critical Review. Chem. Soc. Rev. 2014, 43, 4714−4727. (4) Kotobuki, M.; Kanamura, K.; Sato, Y.; Yoshida, T. Fabrication of All-Solid-State Lithium Battery with Lithium Metal Anode Using Al2O3-Added Li7La3Zr2O12 Solid Electrolyte. J. Power Sources 2011, 196, 7750−7754. (5) Murugan, R.; Thangadurai, V.; Weppner, W. Fast Lithium Ion Conduction in Garnet-Type Li7La3Zr2O12. Angew. Chem., Int. Ed. 2007, 46, 7778−7781. (6) Awaka, J.; Kijima, N.; Hayakawa, H.; Akimoto, J. Synthesis and Structure Analysis of Tetragonal Li7La3Zr2O12 with the Garnet-Related Type Structure. J. Solid State Chem. 2009, 182, 2046−2052. (7) Bernstein, N.; Johannes, M. D.; Hoang, K. Origin of the Structural Phase Transition in Li7La3Zr2O12. Phys. Rev. Lett. 2012, 109, 1−5. (8) Awaka, J.; Takashima, A.; Kataoka, K.; Kijima, N.; Idemoto, Y.; Akimoto, J. Crystal Structure of Fast Lithium-Ion-Conducting Cubic Li7La3Zr2O12. Chem. Lett. 2011, 40, 60−62. (9) Wang, D.; Zhong, G.; Pang, W. K.; Guo, Z.; Li, Y.; McDonald, M. J.; Fu, R.; Mi, J. X.; Yang, Y. Toward Understanding the Lithium Transport Mechanism in Garnet-Type Solid Electrolytes: Li+ Ion Exchanges and Their Mobility at Octahedral/Tetrahedral Sites. Chem. Mater. 2015, 27, 6650−6659. (10) Deviannapoorani, C.; Ramakumar, S.; Janani, N.; Murugan, R. Synthesis of Lithium Garnets from La2Zr2O7 Pyrochlore. Solid State Ionics 2015, 283, 123−130. (11) Zhang, Y.; Chen, F.; Tu, R.; Shen, Q.; Zhang, L. Field Assisted Sintering of Dense Al-Substituted Cubic Phase Li7La3Zr2O12 Solid Electrolytes. J. Power Sources 2014, 268, 960−964. (12) Howard, M. A.; Clemens, O.; Kendrick, E.; Knight, K. S.; Apperley, D. C.; Anderson, P. A.; Slater, P. R. Effect of Ga Incorporation on the Structure and Li Ion Conductivity of La3Zr2Li7O12. Dalt. Trans 2012, 41, 12048−12053. (13) Wu, J.-F.; Chen, E.-Y.; Yu, Y.; Liu, L.; Wu, Y.; Pang, W. K.; Peterson, V. K.; Guo, X. Gallium-Doped Li7La3Zr2O12 Garnet-Type Electrolytes with High Lithium-Ion Conductivity. ACS Appl. Mater. Interfaces 2017, 9, 1542−1552. (14) El Shinawi, H.; Janek, J. Stabilization of Cubic Lithium-Stuffed Garnets of the type “Li7La3Zr2O12” by Addition of Gallium. J. Power Sources 2013, 225, 13−19. (15) Li, Y.; Han, J.-T.; Wang, C.-A.; Xie, H.; Goodenough, J. B. Optimizing Li+ Conductivity in a Garnet Framework. J. Mater. Chem. 2012, 22, 15357−15361. (16) Ren, Y.; Deng, H.; Chen, R.; Shen, Y.; Lin, Y.; Nan, C. W. Effects of Li Source on Microstructure and Ionic Conductivity of AlContained Li6.75La3Zr1.75Ta0.25O12 Ceramics. J. Eur. Ceram. Soc. 2015, 35, 561−572. (17) Ramakumar, S.; Deviannapoorani, C.; Dhivya, L.; Shankar, L. S.; Murugan, R. Lithium Garnets: Synthesis, Structure, Li+ Conductivity, Li+ Dynamics and Applications. Prog. Mater. Sci. 2017, 88, 325−411. (18) Jalem, R.; Rushton, M. J. D.; Manalastas, W.; Nakayama, M.; Kasuga, T.; Kilner, J. A.; Grimes, R. W. Effects of Gallium Doping in Garnet-Type Li7La3Zr2O12 Solid Electrolytes. Chem. Mater. 2015, 27, 2821−2831. (19) Wang, Y.; Klenk, M.; Page, K.; Lai, W. Local Structure and Dynamics of Lithium Garnet Ionic Conductors: A Model Material Li5La3Ta2O12. Chem. Mater. 2014, 26, 5613−5624. (20) Klenk, M.; Lai, W. Local Structure and Dynamics of Lithium Garnet Ionic Conductors: Tetragonal and Cubic Li7La3Zr2O7. Phys. Chem. Chem. Phys. 2015, 17, 8758−8768. 1971

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972

Article

The Journal of Physical Chemistry C (42) Xie, H.; Alonso, J. A.; Li, Y.; Fernández-Díaz, M. T.; Goodenough, J. B. Lithium Distribution in Aluminum-Free Cubic Li7La3Zr2O12. Chem. Mater. 2011, 23, 3587−3589. (43) O’Callaghan, M. P.; Cussen, E. J. The Structure of the LithiumRich Garnets Li6SrLa2M2O12 and Li6.4Sr1.4La1.6M2O12 (M = Sb, Ta). Solid State Sci. 2008, 10, 390−395. (44) Cussen, E. J. Structure and Ionic Conductivity in Lithium Garnets. J. Mater. Chem. 2010, 20, 5167−5173. (45) Awaka, J.; Kijima, N.; Kataoka, K.; Hayakawa, H.; Ohshima, K.; Akimoto, J. Neutron Powder Diffraction Study of Tetragonal Li7La3Hf2O12 with the Garnet-Related Type Structure. J. Solid State Chem. 2010, 183, 180−185. (46) Rao, R. P.; Sharma, N.; Peterson, V. K.; Adams, S. Formation and Conductivity Studies of Lithium Argyrodite Solid Electrolytes Using in-Situ Neutron Diffraction. Solid State Ionics 2013, 230, 72−76. (47) Hitz, G. T.; Wachsman, E. D.; Thangadurai, V. Highly Li-Stuffed Garnet-Type Li7+xLa3Zr2‑xYxO12. J. Electrochem. Soc. 2013, 160, A1248−A1255. (48) Murugan, R.; Ramakumar, S.; Janani, N. High Conductive Yttrium Doped Li7La3Zr2O12 Cubic Lithium Garnet. Electrochem. Commun. 2011, 13, 1373−1375.

1972

DOI: 10.1021/acs.jpcc.7b10911 J. Phys. Chem. C 2018, 122, 1963−1972