Quantification of particle morphology features and electrochemical response during lithiation or sodiation
As illustrated in Fig. 2a, when a negative potential is applied, Li/Na ions in the electrolyte first accumulate on the (010) channel openings before leaping across the carbon coating into the interstitial vacancies present in the first layer of the crystal, while the electrons in the carbon coating tunnel to the adjacent iron site to reduce the Fe3+ ions. After the charge transfer reaction, adjacent Fe2+ and Li/Na ions form a neutral quasiparticle, or polaron, capable of migrating along the preferred [010] channels35. Notably, depending on the particle features and electrochemical response upon lithiation or sodiation, phase evolutions can manifest as two-phase separation or solid solution transition. It becomes evident that controlling particle morphology, encompassing the (010) facet, [010] channels, as well as particle volumes, is of utmost importance, given the anisotropic nature of ion transport in olivine FePO4 crystal. Specifically, the relative area of each facet of a particle depends on its surface energy36. According to our constructed Wulff shape of LiFePO4 from calculated surface energies (Fig. 2b, Supplementary Table 1, and Supplementary Note 1 for computation details), surfaces (201), (100), and (010) have the lowest energies, which are also consistent with reported results36. Here, six distinct well-crystallized LiFePO4 particles were prepared using solvothermal approaches followed by the surface carbon-coating treatment under calcination (Fig. 2c-h, Supplementary Figs. 1–6, and See Methods and Supplementary Note 2 for more synthesis details). Based on Rietveld refinement, the anti-site defect level is estimated to be low for all six particles (Supplementary Fig. 7 and Supplementary Table 2). The delivered capacities, which will be discussed later, further verify the low defect concentrations, such as less than 0.1% for the biggest particle (Cuboid-6000 nm). The facets of synthesized particles predominantly exhibit two orientations, either (010)-oriented (platelet particles) or (100)-oriented (cuboid particles), both of which have low surface energies (Fig. 2b, Supplementary Table 1). Specifically, one of the most crucial morphology features, the [010] channel length (1D migration direction) covers a wide range with average dimensions of 20 nm, 45 nm, 87 nm, 600 nm, 1200 nm, and 6000 nm, respectively (Fig. 2c-h and Supplementary Figs. 1–6). More morphology features are considered to provide a comprehensive quantification of the size and morphology. As shown in Fig. 2i, Supplementary Table 3, and Supplementary Note 3, the average particle length in the [100] and [001] directions are also quantified. Additionally, we determine the average exposure ratio of the (010) facet to the total surface area, a metric ranging from 12–70%. Furthermore, the (010) facet area to [010] channel length ratio is evaluated (ranging from 2.37 × 102 nm to 1.67 × 104 nm), which reflects the accessibility of storage sites and can influence the exchange current density of the particles28,33,34,37. We also estimate the average particle volume by calculating the product of the (010) area and the [010] channel length, which ranges between 2.5 × 10− 4 µm³ and 24 µm³. This collection of particle characteristics was used to identify critical features associated with high Li preference.
The electrochemical lithiation and sodiation behaviors were characterized next. The empty FePO4 hosts were prepared by chemical Li extraction (See Methods for chemical extraction and electrode preparation details), with the structure verification from Rietveld refinement (Supplementary Fig. 8 and Supplementary Table 4). We first compared the constant current (de)intercalation curves of each particle in 1M LiCl or 1M NaCl aqueous solutions (Figs. 3a-b, Supplementary Figs. 9–10, and Supplementary Table 5). The decent capacity delivered during delithiation at 0.1C (17 mA/g), ranged from 129 mAh/g for Cuboid-6000 nm particles to 159 mAh/g for Platelet-20 nm and Platelet-600 nm particles, indicating the low Li-Fe anti-site defects level, especially when considering the channel length (e.g., < 0.1% for Cuboid-6000 nm particles)24.
The intercalation voltage difference between Li+ and Na+ is a good indicator for Li selectivity. The six particles exhibited distinctive group behaviors based on their lithiation C-rate response and sodiation behavior, leading us to categorize them into two groups. As shown in Figs. 3a-b, Supplementary Figs. 9–10, one group displayed minimal and small voltage hysteresis during (de)lithiation and (de)sodiation, which consists of small particles with channel length < 100 nm (Group 1). Particles within this group demonstrated excellent rate capability during (de)lithiation at 0.5C as well, suggesting faster kinetics. Conversely, notable differences in voltage hysteresis emerged during (de)sodiation for Platelet-600 nm, Platelet-1200 nm, and Cuboid-6000 nm particles, leading to their categorization into a separate group (Group 2). Furthermore, all three particles in the larger particle group experienced a reduced capacity retention during their first desodiation (Supplementary Table 5). However, except for Platelet-600 nm, both Platelet-1200 nm and Cuboid-6000 nm particles displayed significant capacity decay during (de)lithiation at 0.5C (Supplementary Fig. 10). Particularly noteworthy is the potential difference at the halfway capacity point of the initial sodiation at 0.1C, which can reach up to 0.36 V between the two groups (Supplementary Table 6). This difference can be attributed to the more pronounced strain/nucleation penalty and slower kinetics experienced by the larger particle group. Additionally, the big particles demonstrate a high degree of non-topochemical Na+ intercalation, as indicated by the decreased capacity retention at the first charge (Supplementary Table 5). It is worth highlighting that the Platelet-600 nm particle in Group 2 has the largest difference in cycling features between (de)lithiation and (de)sodiation. Its size is relatively small to release strain penalty and facilitate fast Li+ (de)intercalation kinetics but is large enough to induce kinetic and chemo-mechanical barriers during Na+ (de)intercalation.
To better isolate the overpotential gain attributed to kinetics and mechanics from thermodynamic energy difference, we monitored the potential change throughout the constant current intercalation until reaching a certain depth of discharge, followed by 20 hours of relaxation in the original solution (Fig. 3c and Supplementary Fig. 11). The voltage difference was labeled in bars with bottom and top boundaries indicating voltages after intercalation and relaxation, respectively. Preconditioned FePO4 particles (Cycled once at 17 mAh/g in 1M LiCl aqueous solution to extract the accessible capacity delivered at the first charge) are used here to follow the steps for the later Li+-Na+ co-intercalation process. Specifically, the calculations of applied current and depth of lithiation or sodiation are based on the delivered capacity in the first de-lithiation rather than the theoretical capacity. For instance, 0.1C' for the Platelet-20 nm particle corresponds to 15.9 mA/g (Supplementary Table 5), and DOD_Li0.35'/Na0.35' represents 55.65 mAh/g capacity usage. As depicted in the left panel of Fig. 3c and summarized in Supplementary Table 7, all particles, except Platelet-1200 nm and Cuboid-6000 nm, have relatively small lithiation barriers. The operando lithiation potential and the equilibrium potential after relaxation both have small differences and are close to the thermodynamic voltage of LixFePO4 (0.259 V vs. Ag/AgCl; See Supplementary Note 1 for computation details), which indicates a low kinetic/chemo-mechanical lithiation barrier at 0.1C'. For larger particles (channel length > 1000 nm), kinetic barriers scale significantly with size. For example, at DOD_Li70', the end intercalation potential can be ~ 0.1 V lower for Platelet-1200 nm and Cuboid-6000 nm particles. During 0.1C' sodiation, as shown in the right panel of Fig. 3c, all particles exhibited higher overpotential relative to the thermodynamic sodiation voltage. This deviation can be attributed to the sluggish intercalation of Na+ ions, accompanied by increased nucleation or strain energy penalties. Additionally, the Group 2 particles exhibited larger overpotentials (summarized in Supplementary Table 7), further highlighting the more pronounced effects of kinetics and mechanics on big particles.
Particle morphology dependent phase evolutions during lithiation or sodiation
The intercalation pathways and the associated phase evolutions or ion storage mechanisms are critical in determining the energy barriers for both Li+ and Na+ intercalation, consequently influencing Li preference7. To investigate the host response upon lithiation or sodiation, we used in situ synchrotron X-ray diffraction (XRD) to track the structure changes. Figures 4a-b, and Supplementary Fig. 12 reveal that, during lithiation, particles in Group 1 (channel length < 100 nm) undergo SS transitions from the initial FePO4 phase, evidenced by a continuous change of peak positions and lattice parameters. These solid solution phases are thermodynamically stable and persist even at lower lithiation rates such as 0.1C or 0.01C (Supplementary Fig. 13). The SS crystallographic insertion pathway indicates the vanishing of the miscibility gap and improved kinetics, which is also consistent with the observed small voltage hysteresis (Fig. 3a, and Supplementary Figs. 9–10). For Group 1 particles, during sodiation, a two-stage phase evolution pathway was observed (Figs. 4c-d, and Supplementary Fig. 14). This pathway commenced with a two-phase equilibrium between FePO4 and Na2/3FePO4 phases, which then progressed into a SS transition from Na2/3FePO4 to NaFePO4 phase. It is worth noting that, at the first stage of sodiation, the smallest Platelet-20 nm particle shows broader peaks and more intensity contributions from the intermediate compositions compared with the Cuboid-87 nm particle, which suggests some degree of solid solution formation, but the transition is still dominated by the two-phase equilibrium. Interestingly, even when doubling the sodiation rate for the Cuboid-87 nm particle, we couldn’t realize the out-of-equilibrium SS transition between FePO4 and Na2/3FePO4 phases (Supplementary Fig. 14b). This underscores the benefit of the intermediate Na2/3FePO4 phase formation in mitigating the volumetric strain during sodiation.
The relatively worse rate capability of Group 2 particles requires us to employ a slower lithiation or sodiation rate. During lithiation, phase-separation-dominated evolutions are witnessed for all three big particles (Figs. 4e-f, Supplementary Figs. 15 and 16a-b). Specifically, in the case of Platelet-600 nm particles, there is a noticeable broadening and a slight left shift of the FePO4 (020) peak in the beginning of lithiation, followed by the co-existence of LiFePO4 and FePO4 phases. This is due to a certain level of intrinsic Li solubility18,38. Interestingly, for the largest Platelet-1200 nm particle, additional features become apparent. (200) and (020) peaks originating from the intermediate phase Li0.625FePO4 manifest before the formation of the LiFePO4 phase (Fig. 4f and Supplementary Fig. 16b). Li0.625FePO4 is situated at the eutectoid point of the phase diagram and has been observed as a preferred intermediate phase at high currents39. Notably, it is intriguing that this intermediate phase can persist even at a relatively slow current (0.077C), and we attribute this to the large particle size. The presence of the Li0.625FePO4 eutectoid composition will assist in releasing volumetric strain, analogous to the Na2/3FePO4 buffer phase. Furthermore, the emergence of the intermediate composition is consistent with the more slanted chronopotentiometry curve observed compared to the Platelet-600 nm particle (Figs. 4e-f). Additional unexpected features were observed for the larger particles (e.g., Platelet-600 nm and Platelet-1200 nm particles) during sodiation. As illustrated in Figs. 4g-i, and the snapshots in Supplementary Fig. 17, the sodiation process primarily involves peak intensity decrease, notable peak broadening, and a slight left shift, without distinct phase transformations. Moreover, the peaks exhibit reduced symmetry during in situ sodiation which indicates strong lattice distortions27. Interestingly, when the electrodes were allowed to relax in the open air overnight (~ 10 hours), subsequent ex situ synchrotron XRD revealed the presence of Na2/3FePO4, which represents the thermodynamic equilibrium phase. These observations suggest significant lattice distortion occurring during Na+ intercalation (Fig. 4j). The pronounced volumetric strain and nucleation energy penalty experienced by the big particles disrupt their structural equilibrium, thereby suppressing in situ phase separation. The phase response of the larger particles is also consistent with the previously observed high overpotential and considerable kinetic barriers (Fig. 3). The phase evolution confirms the rationale behind the grouping of particles based on their morphological form factors that particles in Group 1 have SS lithiation evolution pathway (in equilibrium) paired with two-stage phase evolution pathway during sodiation (in equilibrium), while particles in Group 2 have phase-separation-dominated lithiation evolutions (in equilibrium) paired with out-of-equilibrium sodiation transition.
Li extraction performance and non-faradic ion-exchange in 1D LixFePO4 hosts
The Li extraction performance of the six particles was examined using 1 mM: 1 M Li to Na molar ratio solutions unless specified. As shown in Fig. 5a, employing 70% accessible capacity and a 0.1C' extraction rate, particles in Group 2 exhibited better Li selectivity than those in Group 1. Particularly, the Platelet-600 nm particle showed the highest Li preference (recovered Li/(Li + Na) ratio = 0.95 ± 0.012) with a Li selectivity of 2.1×104, approximately 34-fold higher than that of the Platelet-20 nm particle. Besides, at a higher mass loading, ~ 15 mg/cm2, decreased Li/(Li + Na) recovery ratios were observed for all particles (Supplementary Fig. 18), but the same Li selectivity trend was still observed. We speculate that the decrease in performance could be attributed to the mass transfer limitations on the dense electrode side, highlighting the importance of electrode engineering for future practical applications. The small mass loading of ~ 2.5 mg/cm2 was used here to better reflect the intrinsic Li preference of the particles. We further evaluated the effects of co-intercalation rates in Fig. 5b. Platelet-600 nm particles showed a monotonic decrease in Li selectivity with elevated extraction rates. As shown in Supplementary Fig. 19, when applying a small current (≤ 0.5C'), the overpotential is small so that there is not enough energy to overcome the kinetic, chemo-mechanical, and nucleation barriers associated with sodiation, resulting in excellent Li selectivity ranging from 8.5×103 to 2.1×104. However, as the overpotential increases at higher currents (> 0.5C'), mass transfer limitations of Li+ ion on the electrolyte side become more pronounced. The larger overpotentials overcome the energy barriers for Na+ intercalation, leading to a decline in selectivity (< 9.6×102). The substantial overpotential eventually leads to less than 70% accessible capacity at voltage cutoff.
Interestingly, for the small particles in Group 1, we observed a non-monotonic trend in Li selectivity as the extraction rates increased. The drop in selectivity at a high extraction rate (e.g., 1.0C'/2.0C') was expected due to mass transfer limitations on the electrolyte side caused by low Li+ concentrations (1 mM), similar to what was observed in the larger particles. At lower rates below 0.5C', the differences between Li+ and Na+ intercalation enlarge. This is substantiated by examining the sodiation and lithiation chronopotentiometry curves obtained in pure 1 M NaCl(aq) and 1 M LiCl(aq) under varying rates (Supplementary Fig. 20). At 50% depth-of-lithiation/sodiation, more considerable sodiation barriers (~ 0.05 V) were witnessed from 0.1C' to 1.0C', while there was almost no increase in end potential at elevated lithiation rates, which can be attributed to the excellent rate capability of the small particles. Therefore, the increase in Li selectivity is attributed to the kinetic barrier gained from sluggish Na+ intercalation.
To take advantage of the rate capability of small particles, we can intentionally increase the kinetic barrier of Na+ intercalation by using super-fast extraction rates (e.g., 6C'). At elevated extraction rates, only brines with a higher Li+ concentration are applicable (e.g., 10 mM: 1 M Li: Na) to circumvent the mass transfer limitations on the electrolyte side. In Figs. 5c-d, we conducted a comparison of the extraction performance between Platelet-20 nm and Platelet-600 nm particles at high extraction rates in a 10 mM: 1 M Li: Na (1:100) mixed solution. Platelet-20 nm and Platelet-600 nm particles were previously identified as the worst and best performers in Fig. 5a, at a low extraction rate (0.1C') in a 1 mM: 1 M Li: Na (1:1000) mixed solution. Obviously, Platelet-20 nm particles outperformed Platelet-600 nm particles in all aspects at faster extraction rates (> 4C'). Lower energy input from the smaller overpotentials during extraction, as well as the higher Li selectivity, make small particles a better choice for brines with relatively high Li+ concentrations, such as the biggest Li brine source, Atacama, in Chile (0.22 M: 4 M: 0.4 M Li: Na: Mg)8. In other words, for small particles, the selectivity primarily arises from thermodynamic Li preference. Due to the minimal nucleation barrier and rapid diffusion for both ions, the kinetic preference for Li during co-intercalation will quickly diminish when we use low extraction rates. Specifically, our calculated lithiation voltage of the olivine FePO4 host is 0.259 V vs. Ag/AgCl, which is only 0.146 V higher than the sodiation voltage (0.113 V vs. Ag/AgCl) (See Supplementary Note 1 for more calculation details). The 0.146 V difference cannot bear the three orders of concentration difference between Li+ (1mM) and Na+ (1M) based on the Nernst equation, if there is no kinetic barrier gain.
Moreover, non-faradic ion exchange experiments further verify the thermodynamic-dominated Li preference of small particles. Initially, we pre-intercalated pure Li+ into the hosts at 0.1C' until different depth-of-discharge (DOD_Li0.1'/0.35'/0.5'). Subsequently, we soaked the electrodes in a mixed solution containing 1 mM LiCl and 1 M NaCl while concurrently measuring the OCV (See Methods for more details). Supplementary Fig. 21 shows significantly different OCV curves between the two particle groups, although they all reached equilibrium after around 7 hours. For big particles, an exponential decay of the potential was observed and stabilized within one hour. In contrast, the small particles exhibited a peculiar upward tilt of the curve after the initial decay, bringing it closer to the equilibrium voltage in pure Na solutions. We then measured the Na/Li contents in the particles. As shown in Fig. 5e, substantial non-faradic ion exchange was witnessed for the Group 1 small particles at all depths of pre-lithiation. Specifically, in the case of DOD_Li0.35', 75.6% ± 0.2% of the structure Li+ in Platelet-20 nm particles was replaced by the solution Na+, whereas less than 5% exchange was observed for the three big particles. We also monitored the composition evolutions overtime during the soaking process (Figs. 5f-g and Supplementary Fig. 22). The non-faradaic ion exchange behavior indicates that the kinetics of Li+ intercalation is faster than Na+ intercalation at certain C rate ranges; however, since the thermodynamic preference for Li alone cannot tolerate three-order-of-magnitude difference in Li+ and Na+ concentrations, a significant amount of Na+ will slowly replace structural Li+ via ion exchange. The much smaller nucleation barrier and more rapid ion diffusion in Group 1 small particles facilitate such a significant non-faradic ion exchange process within the 1D olivine FePO4 hosts, which has not been reported before. Consequently, at higher C rates, the shorter amount of contact time can also reduce the degree of ion exchange and promote Li selectivity. The phase evolution during ion exchange was also tracked by the in situ synchrotron XRD (Fig. 5h and Supplementary Fig. 23). Soaking Platelet-20 nm DOD_Li0.5' particles in 1 mM: 1 M Li: Na mixed solution, we observed a continuous broadening of the peaks and a decrease in intensity, particularly in the in-plane direction, such as (101) facet. This indicates a more disordered structure and a shorter coherence length in the in-plane directions after ion exchange. Given the highly confined 1D structure of olivine FePO4 hosts, the observed non-faradic ion exchange process is both intriguing and unexpected. Further investigations are necessary to fully comprehend the underlying mechanism behind this phenomenon.
Identification of critical features with high Li preference
To identify the relationships of Li competitiveness to particle morphology and electrochemical characteristics, we generated correlation coefficient (R) maps for both groups of particles in Fig. 6, Supplementary Figs. 24–25. These maps quantitatively represent the extent of linear relationships between any two variables (See Supplementary Tables 8–9, Supplementary Note 3 for the definition of each variable and summary of the values for each particle). The R values fall within the range of -1 to 1. A correlation coefficient of 1 or -1 means a perfect positive or negative correlation, respectively.
Most importantly, the analysis identifies critical morphology and electrochemical features that indicates Li preference. For electrochemical characteristics, the operando voltage difference (e.g., 0.1C'_E_(Li-Na)0.5') and kinetic barrier difference (e.g., K.B._(Li-Na)0.5') between Li and Na show a strong positive correlation with Li preference for both particle groups. However, the dependence of Li preference on resting voltage difference (e.g., 0.1C'_R.E._(Li-Na)0.5'), which reflects the thermodynamic properties, is relatively weak for Group 2 particles. This indicates that the thermodynamic preference provides the baseline for Li selectivity but the morphologies of FePO4 particles can be further designed to increase the kinetic barrier differences between Li and Na to promote the Li selectivity. For both Group 1 and Group 2 particles, kinetic barrier difference, K.B._(Li-Na)0.5’, and operando intercalation voltage difference, 0.1C'_E_(Li-Na)0.5', are identified as general descriptors for Li selectivity.
For morphology features, Group 1 small particles show very sensitive dimension dependence of Li preference that [010] channel length (L[010]) and (010) area (Area (010)) showed strong positive and negative correlation, respectively. Therefore, (010)/[010], (010)*[010], and (010) ratio all show a strong correlation with Li selectivity. The negative correlation of Area (010) could be due to the non-faradaic ion exchange behavior observed that a larger (010) area increases the reaction interface, allowing for a greater exchange of Na+ ions into the FePO4 host. For Group 2, strong correlations were identified for log[010] and log(010)*[010]. The logarithmic positive growth pattern of these morphological features indicates that excessively large particles are unnecessary to enhance Li favorability. Besides, the correlation between log(010)*[010] and kinetic barrier differences is even stronger than that of log[010]. However, Area (010) displays a weak correlation. This underscores the significance of mechanical deformation in influencing Li preference, as demonstrated by the robust correlation of (010)*[010] (particle volume). Moreover, the reverse correlation of channel length to Li preference for the two groups implies the existence of optimal dimension around 155–420 nm for the best Li selectivity (Supplementary Fig. 26). Below this dimension, the barrier of sodiation is still low; above this dimension, the barrier of lithiation starts to increase. To validate our prediction, a new platelet particle (Platelet-340 nm) with ~ 340 nm [010] channel length was synthesized (Supplementary Figs. 27–28 and Supplementary Note 2). This particle exhibited the highest Li preference (0.1C'_Li% = (95.8 ± 0.3)%) with Li to Na selectivity of 2.3 × 104 and the most significant kinetic barrier difference (K.B._(Li-Na)0.5' = 0.24 V) was witnessed (Supplementary Fig. 29 and Supplementary Table 10).
Additionally, reversible capacity is also important during lithiation or sodiation (e.g., QLi_1st _deint and QNa_1st _deint) for Li extraction performance. Figure 6 shows strong negative trends between the reversible capacity and the [010] channel length or particle size. This suggests that the existence of optimal dimensions of FePO4 particles will also benefit reversibility.