Fast 3D-lithium-ion diffusion and high electronic conductivity of Li2MnSiO4 surfaces for rechargeable lithium-ion batteries.

High theoretical capacity, high thermal stability, the low cost of production, abundance, and environmental friendliness are among the potential attractiveness of Li2MnSiO4 as a positive electrode (cathode) material for rechargeable lithium-ion batteries. However, the experimental results indicated poor electrochemical performance in its bulk phase due to high intrinsic charge transfer resistance and capacity fading during cycling, which limit its large-scale commercial applications. Herein, we explore the surface stability and various lithium-ion diffusion pathways of Li2MnSiO4 surfaces using the density functional theory (DFT) framework. Results revealed that the stability of selected surfaces is in the following order: (210) > (001) > (010) > (100). Moreover, the Wulff-constructed equilibrium shape revealed that the Li2MnSiO4 (001) surface is the most predominant facet, and thus, preferentially exposed to electrochemical activities. The Hubbard-corrected DFT (DFT + U, with U = 3 eV) results indicated that the bulk insulator with a wide band gap (E g = 3.42 eV) changed into narrow electronic (E g = 0.6 eV) when it comes to the Li2MnSiO4 (001) surface. Moreover, the nudged elastic band analysis shows that surface diffusion along the (001) channel was found to be unlimited and fast in all three dimensions with more than 12-order-of-magnitude enhancements compared with the bulk system. These findings suggest that the capacity limitation and poor electrochemical performance that arise from limited electronic and ionic conductivity in the bulk system could be remarkably improved on the surfaces of the Li2MnSiO4 cathode material for rechargeable lithium-ion batteries. This journal is © The Royal Society of Chemistry.

Introduction

There is growing interest in the large-scale application of dilithium orthosilicate (Li2MSiO4, M = Fe, Mn, Co, and Ni) cathode materials[1-3] for rechargeable lithium-ion batteries (LIB) like that of recent developments in anode materials.[4] These are because of their cyclability, thermal stability and polymorphism possibilities over a wide range of temperatures just to mention a few. As a promising cathode material for the next-generation rechargeable LIB, Li2MSiO4 cathode materials are extensively investigated via both experimental and computational approaches. The designed LIBs from such cathode materials are widely used for portable electronic devices and also for sustainable transportation including electric vehicles (EVs) and plug-in hybrid electric vehicles (PHEVs). They also have significant potential to mitigate the intermittency of renewable energy sources such as solar, wind, and hydroelectric power.[5,6]

Arroyo Dambablo et al.[7] computationally studied the effect of the different polyoxyanions composed of (XO4), where X = Ge, Si, As, Sb, and P on Li2MXO4 materials based on the relation between the Li+ deintercalation voltage and their corresponding electronegativity. Their findings indicated the best fit for Li2MnSiO4 materials. However, Li+ deintercalation voltages for Co and Ni silicates were found to be too high for the current electrolyte window. Experimental studies[8,9] also confirmed that Li2CoSiO4 and Li2NiSiO4 have low capacities and poor reversibility. Li2MSiO4 (M = Mn and Fe) are also highly preferable LIB cathode materials compared to other polyanion compounds such as LiMPO4 (M = Fe, Mn, and Co), which are limited to only a single Li+ insertion and extraction per formula unit.[10,11] Moreover, other layered transition metal oxide cathode materials (for example, Li1/3Mn1/3Co1/3O2) are also toxic and costly.[12]

Previous studies have demonstrated that the method used to synthesize Li2MnSiO4 and other polyanion cathode materials also significantly influences their crystal structures, and so does their ionic conductivity.[1,3,13] One of the strategies to improve the electrochemical performance of lithium-ion batteries with Li2MnSiO4 is the optimization of the particle size and morphology. Particularly, the rate of electrochemical performance increases with the reduction in the particle size of Li2MnSiO4. This is because the porous structure significantly shortens the diffusion time of Li+ in Li2MnSiO4 cathode materials. Therefore, since the orthorhombic structure of Li2MnSiO4 has smaller particle size than that of the monoclinic structure, it will have better Li+ diffusion, resulting in better electrochemical performance.[1,3,13-15] Polymorphs of Li2MSiO4 tetrahedral structures are well known to be classified into low- and high-temperature forms, which differ in the distribution of cations within tetrahedral sites of a hexagonal close-packed (HCP)-based arrangement of oxygen. According to the literature, five different structures were observed with Li2FeSiO4,[16-19] three as-synthesized (two are orthorhombic, Pmnb and Pmn21; one is monoclinic, P21/n) and two cycled phases (Pmn21 cycled and P21/n cycled). Similarly, multiple phases have been reported for Li2MnSiO4 [2,13,20,21] with two orthorhombic (Pmn21 and Pmnb) and one monoclinic (P21/n) symmetries; Li2CoSiO4 [22-24] has two orthorhombic (Pnb21 and Pmn21) and one monoclinic (P21/n) symmetries, while Li2NiSiO4 [3,25,26] has one crystal symmetry.

The structural, electronic and electrochemical properties of Li2MSiO4 depend on the distributions of M. The density functional theory (DFT) predicted some parameter variations, for instance, cell, volume and voltage intercalations between intermediate phases (Li2MSiO4, LiMSiO4 and MSiO4) based on their crystal symmetry.[27] According to their report, Li2MnSiO4 in the orthorhombic symmetry with the Pmn21 space group was found to be the most stable structure. Due to the double reversible redox state of Mn ions (Mn2+/Mn3+ and Mn3+/Mn4+) within the electrochemical stability window of LIB common electrolytes,[28-31] its theoretical energy density reaches as high as 333 mA h g−1.[32] Li2MnSiO4 also has a higher operating voltage of 4.8 V against Li+/Li electrodes with stability of organic electrolytes.[28,30,31] Moreover, the Mn4+ is more accessible than Fe4+ and Co4+.[33,34] However, the electrochemical performance of Li2MnSiO4 has many drawbacks such as the very low Li+ diffusion coefficient estimated from 10−14 to 10−16 cm2 s−1 and electronic conductivity from 10−12 to 10−16 S cm−1 at room temperature[35,36] limiting its electrochemical performance. Unlike other polyanion cathode materials, for instance, LiMPO4,[10] bulk Li2MnSiO4 supports two-dimensional Li+ diffusion in its orthorhombic (Pmn21) symmetry.[20,37,38]

The development of battery nanoelectrodes plays a vital role in achieving higher ionic and electronic transport. These are possible in one way by maximizing fast lithium diffusion pathways and in another way by narrowing electronic bandgaps unlike their corresponding bulk insulator because of the quantum confinement effect.[39-41]

Previous studies have demonstrated the development of Li2MnSiO4 cathode materials at nanoscale[13,18,22] with improved charge/discharge kinetics. For instance, a nanostructured Li2MnSiO4 material was synthesized[42] by an easy and fast urea combustion method, which showed an extraordinarily intrinsic electronic conductivity of 4 ± 10−8 S cm−1. Ultrathin nanosheets of Li2MnSiO4 (M = Mn/Fe)[31] were also prepared by a rapid one-pot synthesis method using a mixed solvent of aqueous ethanol, and the results revealed a superior energy density of 340 mA h g−1. Very recently, Ding et al. have reported the (100), (010), (001), and (210) facets as stable surface structures of Li2MnSiO4 using DFT calculations combined with experimental characterizations.[43] Interestingly, they also managed to successfully synthesize elongated, hexagonal prism-shaped Li2MnSiO4@C nanoplates with the preferentially exposed (001) and (210) facets. The initial discharge capacity was estimated as high as 326 mA h g−1 that corresponds to 1.98 Li+ stored PFU of Li2MnSiO4, while the charge specific capacity was found to be 367.6 mA h g−1, which is higher than that of bulk Li2MnSiO4 (333 mA h g−1) determined by a galvanostatic charge–discharge method. Although the Li+ diffusion rate (1.0043 × 10−13 cm−2 s−1) for Li2MnSiO4@C nanoplates was calculated, the kinetic pathways, Li+ diffusion dimensionalities and electronic structures of preferentially exposed surfaces were not exhaustively emphasized. If the two-dimensional Li+ diffusion and electronic insulating properties of bulk Li2MnSiO4 enhanced to fast three-dimensional Li+ diffusion and electronic conductive properties, higher electrochemical performances could be achieved in Li2MnSiO4 cathode materials for LIBs. The DFT calculation is also a very powerful technique to predict metal–ion and metal–air batteries at bulk, surface, and interface structures emphasizing on the charge (ionic, electronic, and polaronic) transport mechanisms, thermodynamic stability, and their catalytic effect.[15,20,35,44-49]

Herein, we employed the DFT + U analysis to investigate the electronic properties, structural stability and lithium-ion diffusion pathways in the bulk and selected surface structures of the Li2MnSiO4 cathode material in the rechargeable lithium-ion batteries. Detail analysis will also be given to rationalize the reason behind the superior surface conduction observed on the Li2MnSiO4 (001) surface unlike the previously reported poor bulk Li2MnSiO4 conductivity.

Methodology

Bulk crystal structure

Li2MnSiO4 crystallizes in two main crystal systems: one is the orthorhombic symmetry of the Pmn21 and Pmnb[21,27,35,50] space groups at low temperatures and the other is the monoclinic symmetry with the Pn[51] and P21/n[52] space groups at high temperatures. These various polymorphs of Li2MnSiO4 showed different electrochemical performances due to the variety of Li+ diffusion pathways through their crystal lattices. Among the four polymorphs, the orthorhombic symmetry with the Pmn21 space group is widely studied due to its high thermal stability, although it supports two-dimensional lithium-ion diffusion whereas the monoclinic symmetry is less stable but supports three-dimensional lithium diffusion.[37,38] The unit cell of the Pmn21 crystal structure was retrieved from the Material Project open database source (see Fig. 1).[53] In the unit cell, the Pmn21 crystal structure has 16 atoms: Li = 4, Mn = 2, Si = 2 and O = 8 atoms. For the present calculations, we construct a 2 × 2 × 2 supercell, comprising 128 atoms: Li = 32, Mn = 16, Si = 16 and O = 64 atoms. Fig. 1(b) shows the supercell of the Li2MnSiO4 crystal structure, while Fig. 1(c) is the corresponding (001) crystal structure to be most emphasized in the subsequent sections. Observing Fig. 1(a) and (b), all the cations are arranged in the tetrahedral complexes of LiO4, MnO4, and SiO4, in which the structure may collapse as hexagonally closely packed oxygen arrays are partially distorted by the Jahn–Teller effect.[54]

Fig. 1

(a) Unit cell of DFT-optimized Li2MnSiO4 structure in the orthorhombic (Pmn21) symmetry as shown in polyhedral complexes in which all the cations Li, Mn, and Si residing in LiO4, MnO4 and SiO4 form a corner-sharing tetrahedral pyramid. (b) Constructed 2 × 2 × 2 supercell of the Li2MnSiO4 structure in the polyhedral symmetry showing all the cations Li, Mn, and Si residing in LiO4, MnO4 and SiO4 complexes forming a corner-sharing tetrahedral pyramid. (c) DFT-optimized Li2MnSiO4 (001) surface structure.

Surface crystal structures

Li2MnSiO4 surface structures were constructed from the bulk supercell with Li–O and Mn–Si–O terminations, especially those experimentally confirmed,[43] namely, (210), (001), (010), and (100) surfaces, were emphasized (Fig. 2). The number of atoms in all surfaces are the same (128) with Li = 32, Mn = 16, Si = 16 and O = 64 atoms. All the Li–O terminated Li2MnSiO4 surface structures, namely, (100), (010), and (001) facets are fixed with more than 14 Å vacuum space to get the ground-state energy. The Mn–Si–O terminated Li2MnSiO4 (210) surface structure is the thinnest with 5.506 Å among all the surface structures, which is indicative of the inconvenience as an electrochemically active surface, as it is predominantly covered by the other facets.

Fig. 2

Structures of low-energy Li2MnSiO4 surfaces.

Computational details

Grid-based projector augmented wave method (GPAW) code was implemented to solve electronic-structure problems with the DFT[55,56] framework. The exchange correlations were examined with generalized gradient approximations (GGA) by the Perdew–Burke–Ernzerhof (PBE) functional.[57] The core electrons were treated by the projected-augmented wave (PAW) method.[58] GPAW uses a real space basis set with grid points. The ground-state energies are calculated using grid space fixed at 0.18. With a Monkhorst–Pack grids 2 × 2 × 2 for bulk, 2 × 2 × 1 for the (001) and (210) surfaces, and 2 × 1 × 2 for the (010) surface and 1 × 2 × 2 for the (100) surface. Li2MnSiO4 is found to be antiferromagnetic in agreement with the previous reports,[13,27,32] and thus, the magnetic moment for Mn ions is set to be 4. All the structures were relaxed to their ground-state levels by employing spin-polarized calculations until the Hellmann–Feynman forces are less than 0.03 eV A−1. GGA (PBE) often leads to significant deviations from experimental results, in particular the band gap (Eg), due to the incomplete cancellation of the electronic self-interaction energy (SIE) for systems exhibiting a strong localization of the d and/or f-orbital electrons, for instance, Mn(3d) in Li2MnSiO4. Previous studies solved this problem[59-61] by employing the appropriate Hubbard correction. For example, U = 3.9 eV was used for Mn(3d), which is not always the case as it also depends on the chemical element incorporated within a particular compound.[59-62] With this consideration, Eg was calculated by varying U-values from 1 to 6 eV with an interval of 0.5. U = 3 eV was found to provide the appropriate Eg and used throughout the calculations. Surface free energies of the selected Li2MnSiO4 surface structures were determined using the surface energy equation (eqn (1)):where Ebulk is the total energy per formula unit of Li2MnSiO4, Eslab is the total energy of the slab constructed from the supercell containing n formula units of Li2MnSiO4 and A is the base area of the slab.

The minimum energy pathway (MEP) and the activation barriers for Li+ diffusion in Li2MnSiO4 were calculated by the climbing-image nudged elastic band (CI-NEB) method.[63] A 2 × 2 × 2 supercell was considered throughout the calculations. For each CI-NEB calculation, five intermediate images were constructed to interpolate the initial and final states along the Li+ diffusion path. To calculate the energy barrier (Eb), all images and other ions in the supercell were also allowed to relax to its equilibrium position. From the calculated Eb, the Li+ diffusion rate (r) and Li+ diffusion coefficient (D) were calculated using the Arrhenius equation, D = a2r and r = νe−,[45] where a is the minimum jump length of the Li atom from one site to another along the facile channel, e is the electronic charge, ν is the hopping rate, kB is the Boltzmann constant and T is the temperature. Throughout the present work, ν = 1013 s−1, kB = 8.6173326 × 10−5 eV K−1, and T = 298 K were used.

Results and discussion

Calculation of surface energies

The structural stability of the selected low index surfaces was identified from the surface free energy calculations by manipulating eqn (1), and the results are summarized in Table 1 and the respective crystal structure facets of Li2MnSiO4 are indicated in Fig. 2. The negative sign indicates that the selected surface structures are stable. Among the facets of Li2MnSiO4, the (210) surface orientation is the most stable surface with a free energy of about −55.6 J m−2.

Calculated surface free energies of selected Li2MnSiO4 surfaces

Orientation	(100)	(010)	(001)	(210)
γ (J m−2)	−42.58	−48.98	−48.82	−55.6

The results indicate that the identified surfaces are in the low index facet with their relative decreasing order of surface stability as (210) > (010) > (001) > (100), which is in agreement with the recent experiment.[43]

Wulff construction

The Wulff shape for Li2MnSiO4 based on the surface free energies in Table 1 is shown in Fig. 3. The green, red, blue, and yellow colors show that the contribution of 35.64, 24.9, 23.4, and 16.47% of a surface by the (001), (010), (100), and (210) facets is based on the Wulff shape-constructed areal surfaces, respectively. These results indicate that all the Li–O terminated surface structures of Li2MnSiO4 are better than that of Mn–Si–O terminated surface structures, which is the (210) facet.

Fig. 3

Wulff-constructed equilibrium shape of Li2MnSiO4 based on the calculated relative surface areas and surface free energies.

In contrast to the calculated surface free energy, the Wulff shape illustrates that the areal surface contributions of the (210) and (100) surfaces are relatively minimum with respect to the formation of facets. The majority of the Wulff shape is dominated by the (001) facet, designated by green color in Fig. 3, with intermediate surface stability, hence limiting wide electrochemical kinetics due to the coverage of an electrochemically active surface area of the (210) facet. Therefore, the (001) facet was selected for further analysis because of the higher surface coverage.

Bulk and surface electronic structures

Bulk and surface electronic structures were considered to examine the electronic properties and their influence on the overall electrochemical performance of Li2MnSiO4. After relaxations, the total magnetic moments of the bulk and (001) surface structures of anti-ferromagnetic Li2MnSiO4 were calculated to be 80 and 73.54, respectively. The average magnetic moment for the Mn ion in the bulk and (001) surface of Li2MnSiO4 structures were calculated to be 4.43 and 4.17, respectively, slightly decreasing in the surface.

The strongly correlated electrons in the case of d-orbital or f-orbitals of transition metal cathode materials need the Hubbard correction (U) due to electron SIE, otherwise underestimate some calculated parameters when compared with the experimental results, for instance, Eg calculated at the DFT (PBE) theory level. Thus, to calculate reasonable Eg values equivalent to experimental results, one needs to add the Hubbard correction, U (eV).[59-62] Although it is challenging to fix the specific U-value for a particular d- or f-orbital element, researchers used different U-values without clear-cut reasons especially when the experimental result is not found.

To fix this, we used different U-values from 1 to 6 eV with 0.5 intervals in the GPAW code, as clearly indicated in Table 2. For this work, the calculated Eg values were found to be 3.42 and 0.6 eV for the bulk and (001) surface structures of Li2MnSiO4, respectively. To the best of our knowledge, there is no experimental or calculated Eg value for Li2MnSiO4 using the GPAW code. According to the literature, some reported Eg (eV) values based on the Vienna Ab initio Simulation Package (VASP) code were found for Li2MnSiO4 in the range between 2.21 and 3.44 eV.[27,53,64-66] It is obvious that the values of the Hubbard correction used in the VASP code are quite different and higher than the values used in the GPAW code employed in this work. To select the appropriate Eg value for this material, we calculated different Eg values as a function of the Hubbard correction (U-value) about the average Eg reported in previous similar studies, and the results are summarized in Table 2. Accordingly, the U-value of 3 eV is chosen as the appropriate Hubbard correction for both the bulk and surface systems with Eg values of 3.42 and 0.6 eV, respectively.

Calculated Eg values with the U value (eV) for the bulk and surface of Li2MnSiO4. The corresponding Eg values for the surface are presented inside the parenthesis

U	0	1	1.5	2	2.5	3	3.5	4	4.5	5	5.5	6
E g	2.5	2.84	3	3.16	3.29	3.42	3.54	3.64	3.74	3.82	3.89	3.95
	(0)					(0.6)

Fig. 4 shows the detailed analysis of the electronic projected density of state (PDOS) based on PBE and PBE + U functional for the bulk and (001) surface structures of Li2MnSiO4. Even though Mn(d) and O(p) orbitals are dominant in both cases, as shown in Fig. 4(a) and (c), their contributions are mixed at the PBE functional, i.e. near the Fermi energy level, where the contribution of Mn(d) is more dominant, while O(p) becomes dominating as far away in the higher occupied molecular orbital (HOMO) levels.

Fig. 4

Calculated PDOS indicated for bulk Li2MnSiO4: (a) PBE and (b) PBE + U (at U = 3 eV). Calculated PDOS for the Li2MnSiO4 (001) surface using (c) PBE and (d) PBE + U (at U = 3 eV) functionals.

When the SIE is corrected with the PBE + U functional, the contribution O(p) is more or less uniformly dominant in the lower occupied molecular orbital (LUMO) level. The influence of SIE correction is observed in Fig. 4(b) and (d).

From Fig. 4, it is clearly indicated that O(2p) dominated all the orbitals due to Mn(3d) delocalization, when PBE + U penalty was incorporated in the calculations.

The comparison of dominant spectra PDOS for the bulk and surface is presented in Fig. 5(a) and (b), and it shows that the PDOS of the (001) surface at the Fermi level is more dominant and has exactly metallic properties than the one in the bulk, indicating a higher conductivity as a result of strongly lowering the Eg value of the (001) surface. The calculated Eg value of 0.6 eV in the (001) surface is due to the contribution of other orbitals, as shown in the total density of states (DOS) of Fig. 5(c). Moreover, it should be very clear here that DFT (GGA at PBE) underestimates electronic Eg in the bulk and surface structures of Li2MnSiO4, while DFT + U (PBE + U) calculated Eg in this work is in good agreement with the previous studies.[64,65]

Fig. 5

Calculated PDOS for the most dominant spectra for the bulk (a) and (001) surface (b) separately, while (c) presents a comparison of the calculated DOS for the bulk and (001) surface.

Bulk and surface diffusion

The ion diffusion paths and dimensionalities within the crystal structures of cathode materials are among the most important fundamental features and also critical for the charge–discharge rate capabilities.[15]Fig. 6(a), (b) and Table 3 demonstrate the proposed possible kinetic pathways and calculated parameters including the minimum possible lithium atom (Li) hopping distance (a), energy barriers (Eb), diffusion rate (r) and diffusion coefficient (D). The calculated results indicate that bulk lithium-ion diffusion in Li2MnSiO4 with the orthorhombic symmetry and Pmn21 space group is confined to two-dimensional diffusion along the [100] and [001] directions with an estimated diffusion coefficient from 3.12 × 10−15 to 6.09 × 10−17 cm2 s−1. The corresponding energy barriers are also found to be in the range between Eb = 0.73 and Eb = 0.88 eV, whilst lithium diffusion along the A–G–H pathway (y-axis) is completely limited mainly in the second half hopping. This is because of the tetrahedral units such as MnO4 and SiO4 blocking the Li+ diffusion channels, resulting in a high energy barrier as high as 1.88 eV that corresponds to very slow diffusion coefficient (3.12 × 10−32 cm2 s−1) of Li+ diffusion from the hopping site G to H pathway, as indicated in Fig. 6(a), (b) and Table 3. The calculated results also indicate that the microscopic diffusions following the A–B–C pathway (x-axis) are twice hopping from A to B hopping site (4.4 × 10−15 cm2 s−1) and from B to C hopping site (6.1 × 10−17 cm2 s−1) with their corresponding energy barriers (0.73 and 0.84 eV), respectively. The microscopic diffusions follow the A–G–H pathway from A to G (4.4 × 10−15 cm2 s−1) and G to H (3.1 × 10−34 cm2 s−1) hopping sites with the calculated energy barriers (0.73 and 1.88 eV), respectively.

Fig. 6

Proposed pathways for Li+ and microscopic Li+ diffusion in bulk (a and b) and (001) surface (c and d) structures of Li2MnSiO4, respectively.

Calculated E, r and D parameters for the bulk and (001) surface structures of Li2MnSiO4

Pathway		a (Å)	E b (eV)	r (s−1)	D (cm2 s−1)
Bulk	AB	3.13	0.73	4.50 × 10−7	4.40 × 10−15
	BC	3.13	0.84	6.0 × 10−9	6.10 × 10−17
	AE	3.19	0.74	3.50 × 10−7	3.10 × 10−15
	EF	3.18	0.88	1.00 × 10−9	1.30 × 10−17
	AG	3.13	0.73	4.51 × 10−7	4.40 × 10−15
	GH	4.41	1.88	1.60 × 10−19	3.10 × 10−34
(001) surface	AB (XY)	4.44	0.06	1.09 × 1012	2.14 × 10−3
	BC (X)	3.14	0.25	5.90 × 108	5.83 × 10−7
	BD (Y)	2.94	0.39	2.09 × 106	1.81 × 10−9
	AE (Z)	2.99	0.22	1.90 × 109	1.71 × 10−6

Similarly, microscopic diffusion follows the A–E–F pathway (z-axis) from hopping site A to E (3.1 × 10−15 cm2 s−1) and E to F (1.3 × 10−17 cm2 s−1) with calculated energy barriers (0.74 and 0.88 eV), respectively. These results indicate that the microscopic diffusion along the (010) direction from hopping site G to H (1.88 eV) is a rate-limiting step during charge/discharge kinetics, which proves that bulk Li2MnSiO4 with the orthorhombic symmetry and Pmn21 space group is truly confined to two-dimensional Li+ transport pathways, as summarized in Fig. 6 and Table 3.

The calculated Eb in this work is better than the previously reported DFT and atomistic simulation studies,[30,36,38] proving that Li2MnSiO4 with the orthorhombic symmetry and Pmn21 space group supports two-dimensional Li+ transport. Despite the calculated high energy barriers (Eb ≈ 0.88 eV) along the (001) and (1000) channels as also reported in other studies, third-dimensional Li+ diffusion, specifically along the (010) channel following AG and GH paths, has not been reported yet. To optimize unlimited and fast Li+ transport paths, all the channels along the 3D direction are of paramount importance that we present by surface designing, particularly, on the Li2MnSiO4 (001) surface.

The Li2MnSiO4 (001) surface was explored to be fast and unlimited Li+ diffusion, and better electronic conductive surface, which is essential to the higher-power capability of rechargeable Li-ion batteries. For detailed analysis, identified Li+ diffusion pathways from bulk to the surface (AE) and on the surface (AB, BC, and BD) were emphasized. The DFT calculations show that surface diffusion along the [100], [010], and [001] directions with a dramatically lowered energy barrier (0.057–0.395 eV) corresponds to a higher estimated diffusion coefficient (2.14 × 10−3 to 1.81 × 10−9 cm2 s−1). Li+ diffusion paths follow similar mechanisms as discussed in the case of bulk diffusion. In Fig. 6(c), (d) and Table 3, the four microscopic Li+ diffusion pathways and estimated energy barriers Eb were identified within the surface and from the surface to bulk (see Fig. 6 and Table 3 for details). The calculated Eb value along the AB (XY), BC (X), BD (Y), and AE (Z) pathways were found to be 0.06, 0.25, 0.40, and 0.22 eV, respectively, attributing to the rapid diffusion coefficient estimated to be 1.90 × 10−3 cm2 s−1 (AB), 9.6 × 10−7 cm2 s−1 (BC), 1.6 × 10−9 cm2 s−1 (BD), and 1.8 × 10−6 cm2 s−1 (AE), respectively, enabling higher charge–discharge rates on the Li2MnSiO4 (001) surface. The surface conductivity of this work is because of the remarkable results of Fig. 4(d), 5(c), 6(c) and Table 3. On the Li2MnSiO4 (001) surface, Li+ diffusion channels were found to be open and Li+ diffuses with very minimum energy costs, inducing surface ionic conductors for this material. With regard to the electron conductivity, the Mn cations with d-orbital and O anions with p-orbital strongly overlapped, indicating the Mn(d) and O(p) hybridizations, which facilitates fast electron transfer, inducing a very narrow bandgap energy (0.6 eV) unlike the bulk insulating property (3.42 eV). The results suggest that Li+ diffusion on the Li2MnSiO4 (001) surface is three dimensional and very fast with a diffusion coefficient estimated to be 12 to 17 orders of magnitude higher than that of bulk Li2MnSiO4. This is also quite different from other polyanion compounds of LIB cathode materials. For instance, LiFePO4 is the well-known material confined to one-dimensional Li+ diffusion along the (010) direction in its bulk and (010) surface structures. However, the (010) surface Li+ diffusion was found to be lower than its bulk diffusion because of the higher energy barrier along this channel, revealing the surface architecture of the Li2MnSiO4 cathode material for high-performance lithium-ion batteries.[67]

Conclusion

Herein, we investigated the Li+ transport paths, electronic structures, and structural stabilities of the Li2MnSiO4 bulk and surfaces using both DFT and DFT + U analysis. The results of surface free energies predicted the following order of stability for low index surface: (210) > (010) > (001) > (100) facet. Based on the calculated free surface energies, an equilibrium shape of the Wulff construction was predicted and the relative percentage of exposure was found to be 35.64% for (001), 24.9% for (010), 23.4% for (100), and 16.47% for (210). Unlike the bulk insulator, electronic analysis results revealed that there are superior (001) surface conductivity pathways in Li2MnSiO4. Moreover, the NEB calculations confirmed unlimited and fast three-dimensional lithium-ion diffusion on the Li2MnSiO4 (001) surface, unlike the corresponding bulk phase due to the large energy barrier. The presence of conductive surface pathways for ion and electron transports along the Li2MnSiO4 surfaces could substantially improve the performance of the Li2MnSiO4-based rechargeable lithium-ion batteries, paving the way for designing and developing high-performance energy storage devices.

Data availability

The data that support the findings of this study are available within the article.

Conflicts of interest

There are no conflicts to declare.