Yilong
Zhou
*ab,
S. O.
Kucheyev
a and
Liwen F.
Wan
*ab
aMaterials Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA. E-mail: zhou17@llnl.gov; wan6@llnl.gov
bLaboratory for Energy Applications for the Future (LEAF), Lawrence Livermore National Laboratory, Livermore, California 94550, USA
First published on 4th February 2025
Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B3N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.
Analogous to the structure of graphite, hBN is comprised of alternating layers of highly polarizable, sp2-hybridized B and N. The layers are stacked together via weak van der Waals (vdW) forces.17,18,22 As a result, pristine hBN is a wide bandgap insulator,27,28 which has historically hindered its practicality as an electrode material in energy storage and conversion applications. However, hBN has found applications as components in composite separators or electrolytes and as protective or artificial SEI for electrode materials in rechargeable batteries.23–26,29–31 For example, hBN-based composite separators in Li–S batteries can sufficiently block polysulfide migration, manage thermal distribution, and suppress dendrite growth.26 Composite electrolytes based on hBN are also found effective to extend the working electrochemical window, warranting their application in next-generation high-voltage Li-ion batteries.29,30 In addition, hBN has been used as a protective SEI in Li batteries to stabilize the Li anode, prevent dendrite growth, and enhance cycling stability.31 Furthermore, significant headway has been made through various engineering approaches to extend hBN functionality.32–35 For instance, hybridizing graphene oxide with hBN as a composite electrode material for Li-ion batteries can not only enhance thermal stability and increase the specific surface area of the electrode—both of which are beneficial for its electrochemical performance—but also help to maintain the structural integrity of the anode during cycling, especially at elevated temperatures.35
Defect engineering has been demonstrated as an effective approach to tune materials properties. For instance, atomic defects can be readily introduced with a variety of synthetic or post-processing techniques.36–41 In hBN, different types of vacancy defects, such as the isolated single B vacancy (VB), single N vacancy (VN), B–N pair vacancy (VBN) and B3N cluster vacancy (VB3N), have been commonly considered.42,43 Vacancy defects can effectively alter the electronic and physicochemical properties of hBN for advanced applications.44–46 When considering hBN as an active component in Li and Li–S batteries, a key question to address is the effect of lattice defects on Li-ion transport behavior in hBN. Most of previous work has been focused on the electronic properties of hBN,43,47–49 with fewer explorations of its interaction with Li.50,51
Here, we carry out first-principles evaluation of the impact of four distinct atomic defects, VB, VN, VBN, and VB3N (Fig. 1b), on Li-ion mobility in hBN. Both in-plane and out-of-plane Li-ion diffusion pathways (Fig. 1a) are considered. It is found that Li can be trapped at the defect sites due to local charge deficiency. However, once moving away from the defect sites, minimal impact is observed for in-plane Li-ion migration. In contrast, introducing these defects is found to have a major impact on out-of-plane Li-ion diffusion kinetics, where a significant decrease in activation energy is observed with increasing defect size. This finding is crucially important for the design and deployment of hBN as a separator or artificial SEI in next-generation high-performing Li batteries. Based on an exponential extrapolation of the activation energies of out-of-plane Li diffusion as a function of defect size, we speculate that defects with diameters larger than 4.4 Å, such as VBN, would allow sufficient Li-ion transport (with a barrier <0.5 eV) through the hBN layers.
Followed by detailed calculations and discussion of Li-ion transport behavior in defected hBN, we first present reference results for Li-ion diffusion in pristine hBN. As depicted in Fig. 2a, three non-equivalent, high-symmetry atomic sites for Li insertion are identified: the hollow (H) site, between the centers of two hexagonal rings at adjacent layers; the tip (T) site, between the B and N atoms sitting across the alternating layers; and the edge (E) site, between the midpoints of the B–N bonds at two adjacent layers. These three distinct Li insertion sites are associated with three different local coordination environments of Li; i.e., octahedral, tetrahedral, and linear with N atoms at the H, T and E site, respectively. The difference in local coordination environment effectively leads to variations in the Li intercalation energies as shown in Table S1 (ESI†). It is found that the thermodynamically most favorable potential intercalation site for Li is the H site, followed by E and T sites (also see Fig. S1, ESI†). The consistently positive intercalation energies obtained at these sites, regardless of the Li concentration (Fig. S2, ESI†), imply that minimal Li intercalation would be observed in pristine hBN.
Based on the layered geometry of hBN, two distinct Li-ion diffusion pathways are compared (Fig. 1a): in-plane diffusion, where Li migrates either within BN layers or between adjacent hBN layers, and out-of-plane diffusion, where Li passes across the hBN layers. Starting with Li at its thermodynamically most favorable intercalation site (H site), we calculate its diffusion barriers following two non-degenerated in-plane diffusion pathways, as shown in Fig. 2b. It is found that Li energetically prefers to migrate directly among the H sites, through the H–E–H pathway, as denoted in Fig. 2b, with an activation energy (Ea) of 0.34 eV. The slightly higher barrier observed along the H–T–E diffusion pathway (Ea = 0.4 eV) is due to the higher formation energy presented at the T site as shown in Table S1 (ESI†). Overall, these nudged elastic band (NEB) calculation results suggest sufficient Li-ion transport between the 2-dimensional BN layers. On the other hand, out-of-plane Li-ion diffusion, along the pathway depicted in Fig. 2c, would experience an extremely high barrier of 6.68 eV, implying impracticality of Li transport across the BN layers (also see Fig. S4, ESI†).
Intuitively, by introducing lattice defects, the hexagonal ring can be opened to permit Li-ion transport. Here, we compare effects of the following four types of hBN vacancies that are commonly discussed in the literature:42,43 single B vacancy (VB), single N vacancy (VN), BN pair vacancy (VBN) and B3N cluster vacancy (VB3N) in a single layer of BN. The final relaxed structures of hBN containing these vacancies are presented in Fig. 1b. As expected, larger vacancies tend to be less stable as indicated by their formation energies shown in Table S2 (ESI†). The formation energies of VB and VN are found similar, which agree with previous work.49 Considering these four types of vacancies, we evaluate their impact on Li intercalation energies. As shown in Fig. 3a, the introduction of these vacancies strongly modulates the local environment for Li intercalation. For example, Li is found energetically most favorable to reside at the center just above the single B vacancy site (T site), while donating its extra valence electron to stabilize the surrounding charge deficient N atoms (Fig. S3a, ESI†). This charge compensation effectively results in a low Li intercalation energy of ∼−3 eV as shown in Table S3 (ESI†). In contrast, in the case of a single N vacancy, Li remains equidistant between the two BN layers but is off-centered from the vacancy. This positioning arises from its tendency to minimize unfavorable electrostatic interactions with the four positively charged B atoms surrounding the vacancy—three from the vacancy layer and one from the adjacent layer (Fig. S3b, ESI†). For the scenarios of larger vacancies, such as BN or B3N, Li is attracted in closer proximity to the BN layer, similar to the case of a single B vacancy. Additionally, Li resides above the B vacant site in the case of a BN vacancy and between the two B vacant sites in the case of B3N vacancy. It is interesting to note that charge transfer is generally observed from Li to its neighboring N atoms in the presence of B, BN, and B3N vacancies, but the partial charges on the B atoms are largely unaffected (Fig. S3, ESI†). The enhanced electrostatic interactions resulting from this local charge redistribution effectively define favored Li intercalation energies at these defect sites, as shown in Table S3 (ESI†). In contrary, in the case of a N vacancy, charge transfer mainly occurs between Li and the two neighboring B atoms surrounding the vacancy while the charges on N atoms remain largely unaffected (Fig. S3b, ESI†), leading to unfavored Li incorporation (Table S3, ESI†).
We subsequently explore out-of-plane Li-ion diffusion from these intercalation sites (Fig. 3b and c). The distinct Li interaction with the single N vacancy, compared with other vacancies (i.e., B, BN, and B3N), also leads to its unique out-of-plane diffusion behavior. Unlike the cases of B, BN, and B3N vacancies, Li migration through a single N vacancy follows a winding path (see the inset in Fig. 3b), passing through the saddle point, where Li forms two local bonds with surrounding N atoms (see the inset in Fig. 3c), necessitating a high activation energy of 4.87 eV. Conversely, with B, BN, or B3N vacancies, Li migration follows almost straight paths through the vacancies (see insets in Fig. 3b), forming four local bonds with surrounding N atoms at the saddle points (see insets in Fig. 3c). The favorable local atomic environment of Li at the saddle point effectively leads to lowered migration barriers. For example, in the case of a single B vacancy, the energy barrier drops to 1.42 eV. When fitting the migration barriers as a function of the pore size, it follows a nearly exponential decay, as shown in Fig. 3c. Specifically, defects with radii larger than 2.2 Å are expected to allow sufficient Li-ion transport with a barrier less than 0.5 eV. These observations, along with the charge density difference analyses (Fig. S4, ESI†), collectively support the notion that vacancies can effectively promote out-of-plane Li-ion diffusion, thereby converting hBN from a 2D to 3D Li-ion conductor for enhanced versatility.
Next, we examine in-plane Li-ion diffusion in defective hBN, starting with the case of a single B vacancy. Fig. 4a and Fig. S5a (ESI†) present the calculated in-plane Li-ion diffusion barriers along three distinct pathways, originated from the most stable Li intercalation site (T) highlighted in Fig. 3a. Following the same notations as in Fig. 2, the paths 1, 2, and 3 considered here echo Li-ion migration pathways of T–E–T–H–E–H, T–H–T–H–T–H, and T–H–T–H–T–H, respectively, in pristine hBN. As shown in Fig. 4a, a strong trapping effect is observed at the B vacancy site (also see Fig. S6, ESI†), where a high activation energy of >1.5 eV is needed for Li to escape from the vacancy site. Even at the site above the defective BN layer, a noticeable trapping effect is evidenced (Fig. 4b) due to the strong electrostatic interactions between Li and the neighboring N atoms (Fig. S7, ESI†). These results are consistent with the formation energy calculation discussed above. As Li migrates away from the vacancy site (to a distance of ∼5 Å within the defected layer and 3 Å in adjacent layers), its diffusion resembles the bulk behavior (Fig. S5b, ESI†). Notably, with one Li trapped at the B vacancy site (Fig. S8b, ESI†), the trapping effect is significantly reduced for an additional Li, where the activation energy is lowered to 0.45 eV for the second Li to escape the trapping site. This implies that the vacancy does not necessarily hinder the overall in-plane Li-ion conduction at high Li loading.
![]() | ||
Fig. 4 In-plane diffusion of Li in hBN in the presence of B or N vacancies. (a) and (b) Energy profiles of Li migration in the presence of a B vacancy at the vacancy layer (a) along three distinct pathways: T–E–T–H–E–H (path 1), T–H–T–H–T–H (path 2), and T–H–T–H–T–H (path 3), and above the vacancy layer (b). The most stable intercalation site for the Li above the vacancy layer is identified as a T site (T1) adjacent to, rather than directly above (T0), the B vacancy (Table S4 and Fig. S7, ESI†). Li intercalation into the T1 site induces the formation of a BN bridge between two adjacent BN layers, where the B atom beneath the Li atom is pushed downward, and the N atom adjacent to the vacancy is pulled upward (side view in the inset). (c) and (d) Energy profile of Li migration in the presence of a N vacancy along a short H–T–H pathway adjacent to the vacancy (c) and along long diffusion pathways (d): H–E–H–E–H (path 1) and H–E–H–T–H (path 2). Associated activation energies (Ea) are labeled. Insets show the corresponding minimum energy diffusion pathways (yellow spheres indicate the initial, local-minimum, and final states; brown spheres indicate the saddle points), with the positions of vacancy highlighted by pink dashed circles. The activation energy of Li diffusion in pristine hBN is 0.34 eV. |
In the presence of a single N vacancy, Li-ion preferentially occupies the off-centered H site near the vacancy, as shown in Fig. 3a. Here, we explore in-plane Li-ion diffusion from this site. Fig. 4c reveals that Li-ion mobility around the vacancy site will not be strongly affected by the presence of the vacancy (also see path 2 in Fig. 4d), with a barrier of 0.28 eV compared to 0.34 eV of the pristine case (Fig. 2b). Note that for the single B vacancy case, no similar Li diffusion pathway (i.e., around the vacancy site) is identified due to the more pronounced trapping effect (see Fig. S6, ESI†). For Li-ion migrating away from the vacancy site (path 1 in Fig. 4d), significantly higher activation barriers are observed. This again indicates strong trapping effects at defect sites although less pronounced compared to the single B vacancy case. As Li moves further away from the vacancy site (to a distance of 5 Å), its diffusion behavior resembles the bulk one (Fig. S5c, ESI†), similar to the case of a single B vacancy described above.
We further explored in-plane Li-ion diffusion in the presence of BN and B3N vacancy clusters. Fig. S9 (ESI†) compares several in plane Li-ion diffusion pathways around the BN (Fig. S9a, ESI†) and B3N (Fig. S9b, ESI†) vacancies and associated energetics. These pathways generally share energy profiles and feature characteristics akin to those observed in the case of a single B vacancy (Fig. 4a). A strong trapping effect, observed near the vacancy site, is alleviated at a distance of ∼5–6 Å away from the vacancy. Additionally, the trapping effect is significantly reduced for a second Li when the vacancy is pre-filled with a Li (Fig. S8d, ESI†). Compared to smaller-size single vacancies, such as B, N, and B–N pair, larger vacancies, such as B3N, seem to show longer-range impact on Li-ion diffusion (Fig. S9b, ESI†), which can be attributed to the higher degree of lattice distortion experienced in the presence of these larger defects.
The above results underscore the current need for a better understanding of the impact of lattice defects on Li-ion diffusion going beyond the simplest scenario of a single Li atom interacting with an isolated point defect. For example, exploring Li transport behavior at higher Li concentrations (Fig. S10, ESI†) and defect densities is equally intriguing and important to address the pivotal role of lattice defects in devising the overall Li-ion storage mechanism and transport kinetics.52 Although this study primarily focuses on single vacancy defects, it would be interesting to extend this work to investigations of the impact of other defect types, such as interstitial atoms,47 dislocations,53 planar defects, and dopant atoms.50 For example, previous studies have shown that substituting three N atoms of the hBN hexagonal ring with larger P atoms can substantially enhance the performance of BN sheets as anode materials in Na-ion batteries.50 Furthermore, it is important for future research to evaluate the role of charged defects, introducing additional complex ion–electron interactions near or beyond the local defects sites and thereby altering Li-ion migration pathways and the associated energetics.
In summary, this study has investigated the influence of vacancy defects on Li intercalation and diffusion in hBN from first principles simulations. Four distinct vacancy defects (VB, VN, VBN, and VB3N) have been systematically examined, and their impact on in-plane and out-of-plane Li-ion diffusion kinetics have been compared to that in the pristine hBN lattice. Our results have revealed that pristine hBN is predominately a 2D Li-ion conductor, with vastly different activation energies for in-plane and out-of-plane diffusion of 0.34 and 6.68 eV, respectively. The introduction of atomic vacancies, such as VB, VBN, and VB3N, is found effective to promote out-of-plane Li-ion diffusion, with negligible energy barriers for Li diffusion through larger defects with effective pore sizes of >2.2 Å in hBN layers. Although these atomic defects exhibit a strong trapping effect for in-plane Li-ion conduction, they do not necessarily slow down the overall Li-ion diffusion kinetics as the thermodynamic driving force for Li trapping is highly localized. For Li ions at distances >5 Å from a defect site, a negligible change of its diffusion behavior is observed. Moreover, in realistic diffusion conditions, with vacancy sites decorated by trapped Li atoms, their impact on additional in-plane Li-ion diffusion is less significant. Overall, introducing vacancies can transform hBN from a 2D into a 3D Li conductor by realizing out-of-plane Li-ion diffusion. This transformation holds great promise for applications such as separators or artificial SEI for advanced Li batteries. By tailoring the type and densities of defects in hBN, Li-ion transport kinetics can be significantly improved while simultaneously mitigating Li dendrite growth. These findings underscore the potential of defect engineering as a powerful strategy to optimize hBN's performance for a wide range of energy storage applications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cp03655g |
This journal is © the Owner Societies 2025 |