Tuning anisotropic ion transport in mesocrystalline lithium orthosilicate nanostructures with preferentially exposed facets

Li2TMSiO4 (TM = Mn, Fe, Co, etc.) is regarded as a new class of cathode materials for next generation Li-ion batteries because of the theoretical possibility of reversible deintercalation of two Li ions from the structure (ca. 330 mA hg−1). Nevertheless, the silicate cathode still suffers from low electronic conductivity, slow Li ion diffusion and structural instability upon deep cycling. To solve these problems, for the first time, we propose a rational design of mesocrystalline Li2FeSiO4 hollow discoids with an ordered single-crystal-like structure and highly exposed (001) facets. The Li2FeSiO4 mesocrystals display a near theoretical discharge capacity, superior rate capability and good cycling stability. The enhanced Li storage performance is ascribed to the unique structural features with a large surface area generated from the hollow mesocrystal structure and a shortened Li+ diffusion path along (001) exposed facets. This new facile, elegant synthesis method that enables the manipulation of crystal growth and subsequent improvements in the electronic and ionic kinetics and structural integrity should have a positive impact on the research and development of silicate materials as promising cathodes for next generation Li-ion batteries. A battery-electrode material that achieves the excellent performance predicted by theory has been synthesised by researchers in China and Canada. A lithium-ion battery is charged by removing lithium from one electrode and storing it in a second. The energy is released when a circuit is connected and the lithium ions flow back the other way, thus producing a current. Li2FeSiO4 is a new electrode material that is exciting battery researchers because it theoretically has a large capacity for the removal of two lithium ions. A team led by Weifeng Wei from the Central South University in Changsha has now developed a simple method for creating hollow Li2FeSiO4 discoids. They show that this material exhibits near theoretical discharge capacity plus good stability as the battery is repeatedly charged and discharged. Mesocrystalline Li2FeSiO4hollow discoids with an ordered single-crystal-like structure and highly exposed (001) facets have been successfully synthesized as a lithium-ion cathode. Li2FeSiO4@C mesocrystals display a near theoretical discharge capacity, superior rate capability and good cycling stability due to their unique structural features with a large surface area generated from the hollow mesocrystal structure and a shortened Li+ diffusion path along (001) exposed facets.

To address the kinetic issues, several efforts such as conductive coatings [15][16][17] , nanostructuring [18][19][20][21] , and heteroatom substitution [22][23][24][25][26] have been pursued to obtain high performance Li 2 FeSiO 4 -based cathodes. Among these strategies, nanoscale materials have been extensively studied because nanostructures generally offer a shorter diffusion length and more surface area for the insertion/ extraction process of charge carriers. In view of the high diffusion anisotropy of Li 2 FeSiO 4 , tailoring the crystallographic facets and hierarchical structures of electrodes is of great importance for the electrochemical performance 19,27 . Theoretical calculations indicate that for Li 2 FeSiO 4 with the Pmn2 1 space group, Li cations preferentially diffuse two-dimensionally along the a/c plane perpendicular to the b-axis in the [001] and [100] directions, with Li-vacancy migration energy barriers of 0.97-1.05eV and 0.85-0.90 eV in the lithiated and delithiated states, respectively 12,28 . Consequently, if the exposed facets of Li 2 FeSiO 4 materials are parallel to the (100) and (001) planes, the kinetics of the Li cation insertion/extraction process would be substantially enhanced compared to randomly oriented materials, as described in Fig. 1a.
In addition to the diffusion kinetics within a primary particle, the rate performance of the electrodes depends more on the interparticle transport of Li cations and electrons 29 . Hierarchical electrode materials, assembled with oriented primary nanoparticles, have therefore been developed to obtain improved rate performance [29][30][31] . For instance, Wu's group demonstrated hierarchically structured Li 1 . 2 Mn 0 . 6 Ni 0 . 2 O 2 cathodes with a high percentage of (110) exposed planes at the surface. This hierarchical cathode material exhibited excellent rate performance due to enhanced Li diffusion kinetics 32 . Qi's group prepared anatase TiO 2 mesocrystalline anode materials with a single crystal-like superstructure that consisted of periodic arrangements of numerous oriented TiO 2 nanocrystals 33 . This unique mesocrystal electrode demonstrated superior rate capability relative to a typical TiO 2 nanoparticle electrode under fast charge-discharge conditions. The ordered mesocrystalline structure offers efficient interparticle electronic conduction, and the uniformly porous structure of the mesocrystals guarantees perfect contact with the electrolyte and leads to fast Li ion transport between the mesocrystalline electrode and the electrolyte.
Motivated by this concept, we propose a facile and versatile solvothermal strategy (Fig. 1b) to fabricate Li 2 FeSiO 4 cathode materials with hollow mesocrystalline nanodiscoids as secondary particles, which are assembled by primary nanograins with enhanced growth of the electrochemically active (001) plane. A mesoscale assembly process involving the formation of metastable solid precursors to generate primary Li 2 FeSiO 4 nuclei, self-assembly of Li 2 FeSiO 4 nanosheets and Ostwald ripening growth is proposed for the formation of the unique, hollow discoid-like Li 2 FeSiO 4 mesocrystals. Combining experiments and density functional theory (DFT) calculations, it is found that ethylene glycol (EG) molecules act as a capping agent for the manipulation of crystal growth. To ensure efficient electron conduction to the current collector, a conformal, nanoscale carbon coating was also applied to the mesocrystalline discoids. Compared to randomly oriented nanoparticles, the unique structure of the hollow mesocrystals enables faster Li ion diffusion along the (001) exposed facets and higher interparticle electronic conduction. Hence, as a cathode material for LIBs, the obtained composite possesses high reversible capacity, outstanding rate capability, and excellent cycle stability.

Materials synthesis
Mesocrystalline Li 2 FeSiO 4 hollow discoids (M-LFS) were synthesized via a solvothermal method where EG was used as the solvent, and the synthetic routes are shown in Fig. 1b 4 were dissolved in 40 mL EG, while 12 mmol FeCl 2 •H 2 O was dissolved in 20 mL EG. The two solutions were then mixed and stirred for 30 min to form a uniform suspension solution. The precursor solution was then transported to a 100 mL Teflon-lined autoclave and heated to 230°C for 48 h. The M-LFS products were collected, washed with deionized water and ethanol several times and then dried at 80°C. To improve the conductivity, a nanoscale carbon coating was applied on the M-LFS by dispersing M-LFS (0.9 g) in a sucrose solution (2 mmol sucrose in 20 mL distilled water) and then carbonizing it for 3 h at 600°C in an Ar/H 2 (95: 5) atmosphere. The carbon content in M-LFS@C was approximately 12.4 wt% as measured by thermogravimetric analysis (TGA) ( Figure S14a). For comparison, Li 2 FeSiO 4 @C, which is a nanocomposite with a similar grain size (Figures S15 and S16) and a similar carbon content (13.7% carbon content, Figure S14b), was prepared via a sol-gel method using EG and citric acid as chelating agents (see the Supporting Information for details); this sample is referred to as SG-LFS@C.
Characterization X-ray diffraction (XRD) patterns were obtained using an X-ray diffractometer with Cu Kɑ radiation (λ = 1.54056 Å) at steps of 0.02°with a counting time of 2 s (Bruker AXS D8 Advance, Germany). Rietveld refinement was performed using the Maud 2.33 software. Morphological and structural properties were observed using a scanning electron microscope (SEM, Nova NanoSEM 230, USA) and a transmission electron microscope (JEM-2100F, JEOL, Japan). A TEM equipped with an aberration corrector for the imaging lenses was applied to achieve high angle annular dark-field (HAADF) images and energy dispersive X-ray (EDX) mapping (FEI Titan3 G2 60-300, USA). X-ray photoelectron spectroscopy (XPS) was performed using an X-ray photoelectron spectrometer (ThermoFisher Scientific ESCALAB 250Xi, USA). XPS spectra were adjusted by the C 1s peak with the binding energy (BE) of 284.6 eV. Nitrogen adsorption-desorption isotherms were measured by a Quadrasorb SI-3MP surface area and pore size analytical instrument at 77 K. The specific surface area was calculated by the standard Brunauer-Emmett-Teller equation, and the pore size distribution was calculated on the basis of the Barret-Joyner-Halenda (BJH) method 34 . TGA was conducted in air at a heating rate of 10°C/min to 800°C.

Electrochemical measurements
The electrochemical performance was tested using a coin-cell assembled in an Ar atmosphere glove box (MBRAUN LABstar, Germany). For the cathode, 80 wt% active materials, 10 wt% acetylene black and 10 wt% polyvinylidene fluoride were mixed in N-methyl-2pyrrolidone to form a slurry. Then, the slurry was coated onto an Al foil. After drying at 110°C in a vacuum oven, the foil was punched into plates with a diameter of 12 mm, and the active material mass loading was~2 mg cm −2 . Lithium metal was used as the anode, separated by a Celgard 2500 membrane and filled with 1 M LiPF 6 in ethylene carbonate and dimethylcarbonate (1:1v/v) as the electrolyte. The cells were galvanostatically charged-discharged between 1.5 V and 4.7 V vs. Li + /Li (Land CT2100A, China). Cyclic voltammetry (CV) curves and electrochemical impedance spectra (EIS) were obtained using an electrochemical workstation (PARSTAT 4000, Princeton Applied Research). The impedance spectra were recorded in the frequency range of 100 kHz to 0.1 Hz with an applied alternating voltage of 5 mV.

First-principles calculations
DFT calculations were performed using the Vienna ab initio simulation package (VASP) 35,36 . The electron-ion interactions were described through the projector augmented wave (PAW) 37 method, and the exchangecorrelation functional was constructed using the generalized gradient approximation (GGA) proposed by Perdew, Burke and Ernzerhof (PBE) 38 . An energy cutoff of 400 eV and appropriate k-point meshes were chosen to ensure that the total energies converged within 3 meV per formula unit of Li 2 FeSiO 4 [39][40][41] . The GGA + U method was applied to accurately calculate the electrochemical properties and partially correct for electron overdelocalization (and self-interaction errors) 42 . Within this approach, the U and J values were set to 5.3 eV and 1 eV, respectively, and a ferromagnetic high-spin Fe state was assumed.
The cell parameters of the relaxed bulk Li 2 FeSiO 4 structure are a = 6.3485 Å, b = 5.3619 Å and c = 4.9720 Å, which are in agreement with other experimental data 28 . The slab technique was used to model the surfaces of Li 2 FeSiO 4 , in which a set of infinite layers separated by vacuum layers were repeated regularly along the surface normal. The vacuum layers were set to 10 Å. When creating a surface, the SiO 4 tetrahedra were considered to be preserved due to their highly covalent Si-O bonds. The supercell lattice parameters (including slab and vacuum) were fixed, and the inner part of the slab was frozen at the bulk positions. Only atoms near the surface were able to relax until the forces were less than 0.03 eV/Å. Surface free energies (E surf ) were calculated from Eq. (1): where E bulk is the total energy per formula unit of bulk Li 2 FeSiO 4 . E slab is the total energy of the given supercell containing n formula units of Li 2 FeSiO 4 . S is the base area of the supercell.
EG was used as the solvent in the solvothermal conditions, and hence, the adsorption energy of the EG molecule ((CH 2 OH) 2 ) adsorbed on the (1 × 1) unit cell of the (001), (010), and (100) surfaces was calculated. The molecule was placed on top of the slab, and the dipole correction (VDW) was applied in the calculations. The (001), (010), and (100) surfaces had the same cutoff energy, and the corresponding (2 × 3 × 1), (2 × 1 × 3) and (1 × 2 × 3) k-point meshes were employed.The adsorption energy was calculated from Eq. (2), where E ads is the adsorption energy, E surf is the total energy of the relaxed slab model, E adsorbate is the energy of the adsorbed molecule and E total is the energy of the system after adsorption of EG.

Structure and morphology characterization
The powder XRD pattern for the hollow discoid-like Li 2 FeSiO 4 mesocrystals (M-LFS) is shown in Fig. 2. All diffraction peaks can be indexed to the orthorhombic structure (space group Pmn2 1 ) with cell parameters a = 6.2564 Å, b = 5.3408 Å, and c = 4.9685 Å (see the Rietveld refinement results in Fig. 2a and Table S1), indicating that pure Li 2 FeSiO 4 phase was synthesized by the solvothermal method 43,44 . The diffraction peak broadening observed in Fig. 2a suggests that the materials are nanostructured, and the average grain size was estimated to be~19.4 nm on the basis of the Scherrer formula. Furthermore, it is interesting to note that, as shown in Fig. 2b, the peak intensities of the (210), (020), and (002) reflections deviate from those of the standard 44 . Generally, differences in a given peak intensity between the experimental and the reference data hint at a preferred orientation along a particular facet 45,46 . As shown in Fig. 2c and Table S2, the I (210) /I (010) , I (020) /I (010) , and I (002) /I (010) intensity ratios for the M-LFS sample are 2.67, 1.39 and 2.02, respectively, which are much higher than the standard values of 1.083, 0.621, and 1.037 for the randomly oriented standard. This is indicative of preferential growth of the (210), (020), and (002) planes in the M-LFS sample.
Representative SEM micrographs of the M-LFS material, shown in Fig. 3a, demonstrate the formation through the solvothermal process of sub-micron, discoid-shaped particles with diameters ranging from 300 to 500 nm and a thickness of~100 nm. The low magnification TEM images shown in Fig. 3b, c reveal that the discoid-shaped particles consist of numerous primary particles with diameters of approximately 20-30 nm, consistent with the abovementioned XRD analysis results. Moreover, as shown in Fig. 3b, the discoid-like particles possess a hollow structure, and the shell thickness coincides with the diameters of the primary particles. This unique feature indicates that the hollow discoids are assembled with a single layer of primary nanoparticles.
The corresponding nitrogen adsorption/desorption analysis results for the hollow discoids ( Figure S1a) display a typical type-II isotherm with an H3-type hysteresis loop in the relative pressure (P/P 0 ) range from 0.5 to 1.0. This is indicative of a porous structure and slit-shaped pores formed by nanoparticle aggregation 34 . The BET surface area and total pore volumes were estimated to be 65.91 m 2 g −1 and 0.198 cm 3 g −1 , respectively (Figure S1). A pore size distribution with two peaks centered at 3-6 nm and 6-80 nm was observed, suggesting that both microand macropores are present in the M-LFS material.
The corresponding selected area electron diffraction (SAED) pattern of a single hollow discoid (Fig. 3d) shows diffraction spots with a six-fold rotational symmetry, suggesting that the hollow discoid has a "single crystallike" structure. The diffraction spots can be indexed to planes along the [001] zone axis, but the pattern is different from any simulated single crystal diffraction patterns along the [001] zone axis (Fig. 3f). As shown in a typical high resolution (HR) TEM image of the M-LFS material (Fig. 3e), many primary nanoparticles with the same [001] zone axis but different orientations exist. This was confirmed by the corresponding fast Fourier transform (FFT) patterns taken from squares A, B, and C (Fig. 3d) that are indexed as the Pmn2 1 phase along the [001] zone axis. Interestingly, if the FFT pattern of square C is rotated 60°in a counter-clockwise direction, it will overlap with the FFT patterns of squares A and B. If the three simulated diffraction patterns from the same [001] zone axis with ± 60°rotations are merged, one can obtain a new diffraction pattern (Fig. 3f) that is in agreement with the experimental SAED pattern (Fig. 3d). Therefore, it is reasonable to conclude that each hollow discoid is mesocrystalline in nature and is composed of nanosized grains along the [001] zone axis with a six-fold rotational symmetry. The exposed surface of the hollow discoids is thus predominantly constructed with the (001) plane perpendicular to the [001] zone axis, as shown schematically in Fig. 3f, which is consistent with the results of the XRD analysis. Note that the preferentially grown (210) and (020) planes parallel to the [001] zone axis contribute more to the grain boundaries within the hollow discoids. Therefore, due to its porous mesocrystalline features and preferentially exposed active planes at the surface, the M-LFS material is expected to show superior Li + transport kinetics.
To improve the electronic conductivity, the M-LFS material was modified with a carbon coating; these materials are referred to as M-LFS@C. The XRD pattern and simulated pattern for the M-LFS@C composite are presented in Figures S2 and S3, and Rietveld refinement was carried out using the Maud software. The M-LFS@C material exhibits a mixed structure of the Pmn2 1 and P2 1 /n phases, and the weight percentages of Pmn2 1 and P2 1 /n were calculated as 96.35% and 3.65%, respectively, as shown in Figure S3 43,44 . Hence, after carbon coating, the M-LFS@C composite still mainly exhibits an orthorhombic Pmn2 1 structure with cell parameters a = 6.2753(1) Å, b = 5.3763(8) Å and c = 5.0201(4) Å, and the preferred growth of the (210), (020), and (002) planes remains intact, as shown in Figure S2b 43 . The single crystal-like, hollow discoid-shaped morphology is also retained in the M-LFS@C composite (Fig. 4a, b). Moreover, as shown in Fig. 4c, d,  Battery performance of mesocrystal Li 2 FeSiO 4 hollow discoids M-LFS@C with (001) exposed facets, a hollow mesocrystalline structure and a conformal carbon nanocoating is expected to provide excellent electrochemical performance as a LIB cathode. For comparison, the Li 2 Fe-SiO 4 @C composite with similar grain size and carbon content synthesized by sol-gel methods (SG-LFS@C) was also investigated as an LIB cathode. Figure S6 shows the electrochemical performance of the bare Li 2 FeSiO 4 hollow discoids without the carbon coating for the first five cycles at 0.2 C. As shown in Figure S6a, the initial discharge capacity is 171 mA hg −1 , corresponding to 1.07 Li ions per Li 2 FeSiO 4 formula unit. In contrast, the initial discharge capacity of SG-LFS@C is 155.6 mA hg −1 , corresponding to 0.972 Li ions per Li 2 FeSiO 4 formula unit ( Figure S7). After carbon coating, the initial discharge capacity of M-LFS@C is further increased to 314.9 mAhg −1 , with 1.97 Li ions stored per Li 2 FeSiO 4 formula unit (Fig. 5a). After 5 cycles at 0.  Two voltage plateaus at~3.2 and 4.6V are observed in the initial charge curve, corresponding to the stepwise redox reactions from Fe 2+ to Fe 4+ , which is consistent with the CV curves (Fig. 5b) 50,51 . Theoretically, half of the capacity (~165 mAhg −1 ) should be observed at the first charge plateau of~3.1V, corresponding to the oxidation of Fe 2+ to Fe 3+ , and another half capacity should be observed for the second charge plateau of~4.6V (Fe 3+ to Fe 4+ ) 13,47,[49][50][51][52][53] . However, the initial charge capacity at the 3.1 V plateau is less than 165 mA hg -1 (~70 mA hg −1 in this work), which has also been observed in previous work 17,18,20,21,24,27,47,50,51,[53][54][55][56][57][58][59][60][61] . This commonly observed non-ideal voltage behavior has been thoroughly investigated with theoretical calculations and experimental characterization 13,50,52,62 . Mustarelli et al. confirmed that voltage deviations from the ideal value are caused by the formation of undesired delithiated phases due to Li concentration gradients 62 . Moreover, Lv et al. 51 employed in situ XRD/XANES analyses to understand the oxidation states of Fe during initial charging of Li 2 FeSiO 4 , suggesting the formation of Fe 4+ for the second plateau. During the initial discharge process, a cathodic peak at 2.62 V and a shoulder at 2.90 V are observed, corresponding to the reduction of Fe 4+ to Fe 2+ 47,51 . The subsequent charge curve is fairly different from the primary charge curve (Fig. 5a), with the peak at 4.65 V disappearing and the peak at 3.26 V shifting to 3.05 V. This phenomenon (~0.3 V voltage drop) is ascribed to structural rearrangements involving Li and Fe atom position exchange during the first charge 18,44,49,63 . After Li/Fe cation exchange and structural rearrangement during the initial charge process, the Pmn2 1 orthorhombic structure changes to an inverse-β II orthorhombic structure as Li cations reinsert into the structure during the subsequent discharge process, in which Li ions diffuse mainly along the [010] and [001] directions 47,50,51,64 . This robust structural integrity is beneficial for maintaining Li  Fig. 5b, no deviation occurs in the anodic and cathodic peak potentials during the second and third scans, indicating that M-LFS@C exhibits excellent reversibility for Li + intercalation and deintercalation in the charge and discharge process.
To evaluate the advantages of the hollow mesocrystalline structure, the charge-discharge profiles of SG-LFS@C and M-LFS@C were collected at increasing current densities ranging from 0.2 C to 8 C (Fig. 5c, d). With increasing C rate, the charge plateau shifts to a higher voltage and the discharge plateau shifts to a lower voltage because of the polarization induced by kinetic limitations. However, it can be noted that in all situations, the discharge specific capacity of M-LFS@C is much larger than that of SG-LFS@C. For instance, the discharge capacity of M-LFS@C is 299.9, 281.6, 262.1, 231.5, 187.3, 157.5, and 135.5 mA hg −1 at the rates of 0.2, 0.5, 1, 2, 4, 6, and 8 C, respectively, while that of SG-LFS@C is 155.6, 131.8, 110.5, 94.8, 81.7, 63.8, and 52.9 mA hg −1 , respectively, for the same rates. Furthermore, as shown in Figure S9, the discharge capacities for the M-LFS@C materials achieved in this work are higher than those of most other previously reported Li 2 FeSiO 4 cathodes. Considering that both SG-LFS@C and M-LFS@C materials possess similar nanograins and carbon coatings, the significant improvement in rate capability can be attributed to the unique structure of mesocrystalline M-LFS@C, which provides short Li + ion diffusion pathways along the (001) exposed planes and a large reaction area in the hollow mesocrystal structure. Figure 5e compares the rate cycling performances of the SG-LFS@C and M-LFS@C materials. For M-LFS@C, the discharge specific capacity can recover to 303.7 mA hg −1 when the rate is set back to 0.2 C after 35 cycles, corresponding to a capacity retention of 100%. In contrast, the discharge capacity of SG-LFS@C is relatively low (155.6 mA hg −1 at 0.2 C), and the capacity retention is 83.6% when it is set back to 0.2 C (130.1 mA hg −1 at 0.2 C). Figure 5f shows the cycling performance of SG-LFS@C and M-LFS@C electrodes at 1 C for 100 cycles. The initial discharge capacity of M-LFS@C is 258.8 mA hg −1 , and the capacity of 201.1 mA hg −1 is retained after 100 cycles, giving a capacity retention of 77.7%. For SG-LFS@C, the discharge capacity decreases from 102.9 mA hg −1 to 72.1 mA hg −1 , corresponding to a capacity retention of 70.0%. In addition, the discharge curves for the subsequent cycles ( Figure S8) for M-LFS@C exhibit slower potential degradation profiles at different cycles compared to SG-LFS@C. Moreover, as shown in Fig. 5g, the M-LFS@C cathode demonstrates excellent cycling stability at a high rate of 8 C. The initial discharge capacity of M-LFS@C is 138.4 mA hg −1 , and the capacity then gradually increases to 146.2 mA hg −1 at the 40th cycle. M-LFS@C retains a discharge capacity of 130.3 mA hg −1 after 200 cycles, with a capacity retention of 94.1%. Additionally, the hollow mesocrystal morphology is well preserved after 200 cycles, as shown in Figure S10, suggesting robust structural stability. As discussed previously, Li 2 FeSiO 4 may be attacked by HF generated from the LiPF 6 -based electrolyte by the hydrolysis of LiPF 6 43,65 . When the lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) system is used as the electrolyte, the cycle stability may be improved 65 . Further studies on the effects of different electrolytes, such as LiTFSI-based and LiDFOB-based systems, on the longterm cyclability of silicate cathodes are essential for the advancement of silicate chemistry for Li-ion cathodes.
To understand the excellent rate capability and cycling stability, the Li ion diffusion coefficient (D Li+ ) for SG-LFS@C and M-LFS@C was determined using CV data collected at different scan rates ( Figure S11) and was calculated on the basis of Eqs. (1) and (2) given in the Supporting Information. As shown in Table S3, the Li ion diffusion coefficient (D Li+ ) of M-LFS@C is higher than that of SG-LFS@C, suggesting that Li + diffusion is facilitated by the hollow mesocrystals with the (001) exposed planes. Electrochemical impedance spectroscopy (EIS) was also performed for the SG-LFS@C, M-LFS, and M-LFS@C cathodes to evaluate their charge transfer resistance ( Figure S12 and Table S4). The charge transfer resistance values for the SG-LFS@C, M-LFS, and M-LFS@C cathodes are 505.7, 341.4 and 73.24 Ω, respectively (Table S4). The much lower charge transfer resistance suggests a more favorable charge transfer process in the M-LFS@C cathode, in accordance with the superior rate performance of the M-LFS@C (Fig. 5e). It has been established that the reversible volumetric changes for the extraction of the first Li ion and the second Li ion are 1.57% and 6.03%, respectively 13 . The unique hollow and mesocrystal structure could accommodate the large volume changes during the Li + intercalation and deintercalation processes, resulting in the improved cyclability observed in the M-LFS@C cathode with a long cycling life. This is also confirmed by the intact hollow mesocrystal structure after 200 cycles at the high rate of 8 C ( Figure S10).

Formation mechanism of mesocrystalline Li 2 FeSiO 4 hollow discoids
The excellent electrochemical performance of M-LFS@C is closely related to the hollow mesocrystalline structure with exposed (001) planes, so it is crucial to determine the formation mechanism of this unique structure. Precipitates collected at different reaction times without calcination were characterized using SEM, TEM, and XRD, as shown in Fig. 6a-g and S13. During the initial stage (1 h), the products are large irregular-shaped sediments and are composed of Li 2 Si 2 O 5 and Fe 3 O 4 phases, as confirmed by XRD analysis (Fig. 6g). When the reaction time was extended to 6 h, the large sediments gradually decomposed into sub-micron discoids, which are composed of nanosheets (Fig. 6a, b and S13). Moreover, a diffraction peak corresponding to the (210) plane of Li 2 FeSiO 4 appeared in the XRD pattern (Fig. 6g). As the reaction time was continually increased to 12 h, the nanosheets self-assembled into hollow discoids (Fig. 6c, d  and S13). After aging for 24 h or longer, most of the nanosheets were converted into equiaxed nanograins and formed hollow discoid-like Li 2 FeSiO 4 mesocrystals by self-assembly (Fig. 6e, f). The diffraction peaks corresponding to the orthorhombic Pmn2 1 structure became stronger and sharper with increasing reaction time (Fig. 6g). Taking all the structural and morphological characterization results into account, a possible formation mechanism of the mesocrystalline Li 2 FeSiO 4 hollow discoids is proposed and illustrated in Fig. 6h. Initially, large sediments composed of amorphous precursors are formed based on Reaction (1). With increasing reaction time, the precursors decompose into nanosheets under solvothermal conditions. EG, used as the solvent in the experiment, can easily adsorb to the surface of the Li 2 FeSiO 4 nanosheets, leading to selfassembly of the nanosheets into a discoid-like structure. Meanwhile, the metastable nanosheets tend to dissolve and recrystallize into nanoparticles via an Ostwald ripening process. Moreover, under the influence of EG, the recrystallized nanoparticles primarily grow along the same direction, resulting in mesocrystalline discoids. In addition, the inner nanosheets in the discoids are dissolved and redeposit on the shell, leading to a hollow discoid-like structure. The formation mechanism of the preferentially exposed (001) facets in hollow discoid-like Li 2 FeSiO 4 mesocrystals was further investigated by first-principles theoretical calculations. From a thermodynamic perspective, it is well-known that the growth of a certain crystallographic plane depends on its surface energy 65,66 . Previous studies have shown that solvents with hydroxyl groups can easily adsorb on a certain crystallographic plane, alter its surface free energy and, in turn, influence the crystal growth process 65 . In this work, EG with two hydroxyl groups in each molecule was used as the solvent in the solvothermal process. The EG molecules tend to absorb on the Li 2 Fe-SiO 4 surfaces via hydrogen bonds between the hydroxyl groups and the oxygen atoms in different crystal facets. As shown in Fig. 6i and Table S5, the (001) facet of orthorhombic Li 2 FeSiO 4 possesses four oxygen atoms, twice the number for the (001) and (010) facets. Therefore, EG molecules show a stronger tendency to be adsorbed on the (001) surface of orthorhombic Li 2 FeSiO 4 . The adsorbed EG molecules will block the formation of new nuclei on the (001) surface and act as a capping agent to promote the preferential growth of the (001) facets, favoring the formation of the ordered mesocrystalline structure.

Conclusions
In summary, hollow discoid-like Li 2 FeSiO 4 mesocrystals with a single crystal-like structure and highly exposed (001) facets have been successfully prepared via a solvothermal method. A mesoscale decomposition/selfassembly process, involving the decomposition of solid precursors to nanosheets, self-assembly into discoid-like products from the nanosheets and Ostwald ripening crystal growth, is proposed for the formation of the hollow discoid-like Li 2 FeSiO 4 mesocrystals. Combining experimental data and DFT calculations, the EG solvent is shown to play an important role as a shape-directing agent for the manipulation of crystal growth. The mesocrystal Li 2 Fe-SiO 4 @C (M-LFS@C) realizes near theoretical reversible capacity, superior rate capability, and improved cycling stability. The enhanced electrochemical properties can be ascribed to the unique structural features with a uniform carbon coating, a large surface area generated from the hollow mesocrystal structure and a shorter path for fast Li + diffusion along the (001) exposed facets. These results should shed light on the research and development of silicates as promising cathode materials for next generation LIBs that require higher energy and power densities.