Two-dimensional B[Formula: see text]C as a potential anode material for Mg-ion batteries with extremely high theoretical capacity.

The development of new high-capacity anode materials using ions other than lithium as a charge carrier is one of the essential strategies in searching for next-generation high-performance rechargeable batteries. Herein, using first-principles computations, we explore a B[Formula: see text]C monolayer as a potential anode material for Mg-ion batteries. The high stability of the free-standing B[Formula: see text]C monolayer has been demonstrated via calculating the adsorption energy, phonon dispersion, and ab-initio molecular dynamics simulations. The metallic character of the B[Formula: see text]C monolayer, desirable from the point of view of energy storage, ensures good electronic conductivity during the battery charge/discharge process. The calculated migration energy barrier, open-circuit voltage, and theoretical specific capacity of the B[Formula: see text]C monolayer are much better than those of some other two-dimensional materials. These findings provide the B[Formula: see text]C monolayer as a potential candidate for Mg-ion battery anode material with a high theoretical specific capacity of 3187.55 mAh/g.

Introduction

Due to the rapid growth in global demand for energy, the development of high-performance energy storage devices with high energy density is much desirable[1-3]. Currently, rechargeable batteries have attracted great attention due to their advantages like portability and reusability. Li-ion batteries (LIBs) have been developing rapidly during the past decades, and they are use in many places from portable electronic devices, such as mobile phones, and laptop computers to electric cars[4, 5]. However, the limited lithium resources in the Earth’s crust and the associated continuously increasing cost of this material lead to the high production cost of LIBs[6]. Recently, ions of various alkali metals such as sodium (Na) and potassium (K) have attracted significant research attention because of their low-cost and similar chemistry to Li. On the other side, the multivalent ions of alkali earth metals such as magnesium-ion (Mg) and calcium-ion (Ca) are another appropriate choice that could replaced Li due to their relative abundance, good environmental compatibility, and greater charge densities that lead to higher theoretical capacities[7-10]. However, one of the most critical challenges for metal-ion batteries is the development of high-performance anode materials, which ensure adequate capacity, power density, charge and discharge rates, and cycle life of the battery[11].

Many two-dimensional (2D) materials like silicene, phosphorene, transition metal dichalcogenides (TMDs), and transition-metal carbides (MXenes) have shown a high potential to be used as anode materials for alkali and alkaline earth metal ion batteries due to the unique physical/chemical properties and high surface-to-volume ratios[12-16]. Thus they would provide rapid ions migration as well as sufficient channels for ions insertion, delivering high capacity and rate performances[3]. One of the most widely studied 2D materials is graphene[17]. However, its applications in batteries are severely limited due to its relatively low capacity (372 mAh/g for LIBs)[18]. Nevertheless, graphene-like nanostructures can be used as a potential negative electrode material for rechargeable ion batteries. For instance, Liu et al. showed that new graphene and nitrogen monolayer, CN, is a promising material for LIBs, which exhibits low open-circuit voltage (0.12 V), high reversible capacity (840.35 mAh g), fast charging/discharging rate, and good electronic conductivity[19, 20]. Jana et al. designed buckled graphene-like monolayer structure PC with remarkably high theoretical capacity (1200 mAh g) and an ultralow sodium diffusion barrier[21]. Good application performance in LiBs fields exhibits also PC with the storage capacity value of 478.61 mAh g and 717.09 mAh g reached for Li atom respectively adsorbed on single and double side of PC monolayer[22]. It was recently shown that boron and carbon 2D material named BC should find potential applications in future nanomechanics, electronics, and optoelectronics. For example, Dai et al., based on first-principles calculations, reported that BC is a 2D phonon-mediated superconductor with a relatively high transition temperature of 19.2 K[23]. Recent reports show also that BC has enormous potential to be applied as anode material for LIBs/NIBs. In this case, the results reveal that BC exhibits a very high theoretical capacity and a low diffusion barrier[24, 25].

Motivated by the above report on the excellent electrochemical performance of BC as an anode material for LIBs and NIBs, we performed first-principles simulations to investigate the properties of monolayer BC as anode material for rechargeable Mg-ion batteries.

Results and discussion

A first-principles prediction of a new two-dimensional inorganic material, namely, the BC was reported by Wu et al.[26]. As shown in Fig. 1a, different from graphene, BC is composed of hexagons and rhombi with one carbon atom and two boron atoms per unit cell (the black dashed framework). We find the calculated relaxed lattice parameters of freestanding BC to be Å and Å, and the bond lengths of B-B and B-C are 1.556 Å and 1.683 Å, respectively, which are in good agreement with the previous studies[23, 24, 27]. The corresponding electronic band structure along high symmetry lines in the Brillouin zone (see Fig. 1b) and density of states are given in Fig. 1c, from which one can see that supercell of BC exhibits metallic nature with several bands crossing the Fermi level. The metallic character ensures good electrical conductivity which is desired for efficient anode materials. Moreover, the calculated phonon dispersion curves for BC are shown in Fig. 1. A foremost important observation is the positive frequency of all phonon modes throughout the Brillouin zone, which indicates the dynamical stability of the BC structure.

Figure 1

(a) Optimized structure of supercell of BC with the rectangle indicating the unitcell. Green for B, brown for C. (b) The Brillouin zone for P2mm structure with special k-point paths. (c) The electronic band structure of BC supercell along the -Y-S-X--S high-symmetry line, together with the total and partial density of states. The Fermi level is set to zero. (d) Phonon-dispersion relations of pristine BC monolayer.

The Mg adatom can either adsorb on a bridge site (B) above a C–C bond, bridge site (B) above a C–B bond, on top (T) of a carbon site, on top (T) of a boron site, in the hole of a hexagonal ring (H) or in the hole of rhombus (H) - see Fig. 2. For each adsorption site, a geometry optimization was performed. As we can see the energetically more favorable is H position.

Figure 2

Possible adsorption sites on BC surface for Mg atom (left panel) and comparison of total energy for a single Mg atom adsorbed at different sites on BC monolayer (right panel). The red dotted line markers the values of the lowest energy.

Yu et al.[24] also reported that after full structural optimization the most stable adsorption site on BC surface for Li/Na is the sites above the center of the hexagonal BC ring. Starting from this most favorable position we considered a series of configurations with the chemical stoichiometry of () where Mg atoms were initially adsorbed onto the hollow site (H) on both-sides of BC. As a result of a structure optimization, the BC substrate is strongly modulated as shown in Fig. 3. In particular, for higher concentrations, the Mg atoms move significantly outwards from their initial hollow positions, due to their mutual Mg-Mg atoms interaction. Energetically more favorable is to maximize the distance between the Mg atoms. For this reason, the Mg atoms were moved, in the optimization process, on top of the atom or bond. However, no structural destruction was observed, and the Mg atoms were still tightly adsorbed on the BC surface. All of these observations indicate that the structures under our study are dynamically stable. From an application point o view only strong deformation of structure MgBC (Fig. 3d) is undesirable because can effects the operation of the battery under experimental conditions, especially the swelling effect (volumetric expansion) during multiple charge/discharge cycles can appear.

Figure 3

The top view and side view of relaxed structure of supercell of BC monolayer with (a) 2, (b) 6, (c) 12, (d) 16, (e) 18, and (f) 24 adsorbed Mg atoms. The corresponding electronic band structures are presented below.

In electrode materials, relatively large adsorption energy is important in the process of adsorbing metal ions. The average adsorption energy of Mg atoms we determined from:where is the total energy of pristine BC monolayer; is the energy of a single Mg atom in the bulk structure (hcp); n is the number of adsorbed Mg atoms, and means the total energy of BC with n adsorbed Mg atoms on both sides. From this definition, more negative adsorption energy means that the adsorption is more favorable and the adsorbing system is more stable. The results of the calculations of several Mg concentrations on double side configurations are collected in Table 1.

Table 1

The calculated average adsorption energy () and theoretical capacities (C) for BC with Mg adsorbed atoms.

Configuration	E_ads (eV)	C (mAh/g)
Mg_0.167B_2C	−0.9258	265.629
Mg_0.5B_2C	−1.0150	796.888
Mg_1B_2C	−1.4456	1593.78
Mg_1.333B_2C	−1.6179	2125.03
Mg_1.5B_2C	−1.6428	2390.66
Mg_2B_2C	−1.8376	3187.55

We noticed that in all cases the adsorption energy keeps a negative value and systematically decreases by increasing the Mg concentration. The minimum value of adsorption energy ( eV) is achieved when both sides of BC are fully covered with 24 Mg atoms which correspond to the MgBC configuration. It should be noted, that usually the average adsorption energy gradually becomes less negative with increasing metal atoms coverage, indicating weaker and weaker binding. This is primarily caused due to the enhanced repulsions between the positively charged metal adatoms. However, in the case of Mg atoms, the situation is entirely different, similar to this previously observed for the borophene sheet explored as the Mg-ion anode material[28].

The maximum theoretical specific capacity (C) of BC can be computed via the following equation:where n is the number of adsorbed Mg atoms, z is the valence number ( for Mg), F is the Faraday constant (26801 mAh/mol), and is the molar mass of BC substrate. The results in Table 1 show that, for the maximum Mg concentration with both-side adsorption (MgBC), the BC monolayer could provide an extremely high capacity of 3187.55 mAh/g which is much larger compared to the storage capacities of Mg for other 2D anode materials like CN (588.4 mAh/g)[29], VO (815 mAh/g)[30], phosphorene (865 mAh/g)[31], arsenene (1430.9 mAh/g)[32], borophene (1960 mAh/g)[33], flat borophene films (2480 mAh/g)[28], BC monolayer (796 mAh/g)[34] and close to that of BSi (2749 mAh/g)[35].

Open circuit voltage (OCV) is another crucial parameter to characterize the performance of metal-ion batteries. In theory, the OCV curve can be obtained by calculating the energies in the concentration range of :where, , and are the energies of , and single Mg atom in the bulk structure (hcp), respectively. Symbol e denotes the fundamental charge and z is the electronic charge of Mg ion ( for Mg). Figure 4 shows the voltage profiles of the BC sheet for various Mg content. It is important to note that the voltage remains positive during the whole range of coverage values (the voltage potential ranges from 0.46 to 1.21 V which is desirable for anode materials). A negative value means that the metal ion prefers to form metallic states instead of adsorbing on the anode surface[28]. Thus, the lack of negative voltage indicates that BC monolayer is suitable to be the anode materials for Mg-ion batteries. The average potential value of 0.86 V is compared to value computed for VS monolayer (0.93 V)[36] and is higher than that of other 2D anode materials for Li-ions, such as MoS (0.26 V)[36], MoC (0.68 V)[37], TiC (0.43 V)[38]. From the above findings, we conclude that BC monolayer can be considered as a potential anode material for Mg-ion batteries.

Figure 4

Calculated voltage profile with respect to Mg content from 0 to 2.

Moreover, we have calculated the energy barrier for Mg atoms diffusion on monolayer BC. The mobility of metal-ions defines the rate performance, which is a significant character for electrode materials[39]. Considering the structural asymmetry of BC monolayer, metal ions can diffuse along three possible paths to the most stable adsorption sites (see insert in Fig. 5). The climbing-image nudge elastic (CI-NEB) calculations based on the diffusion energy profile of the Mg atom are presented in Fig. 5. The maximum barrier energy between the most stable adsorption site for path 1, path 2, and path 3 are calculated to be 0.30, 0.89, and 0.91 eV, respectively. Thus Mg diffusion in BC shows a strong directional anisotropy. The calculated diffusion barrier is comparable to that of commercially used anode materials based on graphite and TiO with a barrier of eV for Li[40, 41], indicating that Mg atoms can diffuse easily in BC monolayer along path 1. In contrast, for the diffusion in the other two cases (path 2 and path 3), rather large barriers mean an unfavorable Mg diffusion along these paths.

Figure 5

Migration barrier energies of Mg atom from the energetically stable absorption hole site to nearest hole site.

Finally, to validate the BC monolayer as anode material for Mg-ion batteries, ab-initio molecular dynamics (AIMD) simulations were conducted under an NVT (constant number of atoms, volume, and temperature) ensemble at the finite temperature (300 K) to evaluate and check the thermal stability of freestanding BC. Firstly, the MgBC was chosen as the highest Mg concentration to investigate the structural stability of the fully charged BC monolayer. Figure 6a illustrates an overview of the AIMD simulation results for MgBC. The slight fluctuations in the kinetic energy imply the stability of the investigated system at 300 K even after 5 ps. The Mg adatoms were found to be intact and unclustered. Afterward, taking into account the resulting structure from the simulations, the Mg atoms were removed from the BC monolayer (situation of the fully discharged system), and further AIMD calculations were carried out (see Fig. 6b). The obtained results confirm the stability of the initial pristine structure and prove the reversibility in the charge/discharge process of BC material.

Figure 6

Kinetic energy fluctuation as a function of time and corresponding snapshots before and after simulation for (a) fully charged BC monolayer (MgBC) and (b) fully discharged BC monolayer at 300 K.

The structure degeneration (the appearance of various kinds of defects) of anodes is one of the major causes of damage or deterioration of metal-ion batteries. Self-healing materials recently have been shown to improve the cycle life of metal-ion batteries[42]. To examine the susceptibility of BC material to degradation or self-regeneration, we have created a hypothetical situation in which we induce a hole and check the time-dependent behavior of the material at room temperature (300 K). Single-atom (B or C) was pulled out of the BC layer creating the defected region. The self-healing process of BC as a function of time is shown in Fig. 7a and b. The atoms move to the positions of the vacant and heal the system. The final bonds are slightly modified and BC layer is slightly bulged, nevertheless, the hexagons and rhombi arrangement of atoms is maintained. The self-healing is faster in the case of carbon-defect, however, in both cases, the layers can heal without the Stone-Wales type defects in the final structures of BC. A successful self-healing is another feature that demonstrates the possibility of practical industrial applications of BC monolayer to achieve a high-energy electrode for Mg-ion batteries.

Figure 7

Spontaneous self-healing of BC monolayer without any external stimulus at 300 K. The point defect was introduced by extracting the one boron atom (a) or one carbon atom (b) from the surface of the material to a distance of 2.7 Å.

Conclusions

In summary, by using systematic first-principles calculations, the BC monolayer has been explored as an anode material for the Mg-ion battery for the first time. Combining the binding energy and the degree of deformation of BC, we conclude that the most stable configuration of the adsorbed Mg ions on the BC surface is MgBC. Monolayer BC keeps the metallic properties after Mg adsorption, which is fundamental for its use as an electrode. Finally, the BC monolayer exhibits a high theoretical storage capacity of 3187.55 mAh/g. To ensure whether the monolayer BC will not be damaged during the charge/discharge process in rechargeable batteries, AIMD simulations including two layers of Mg that cover both sides of BC at a temperature of 300 K were conducted. The result confirmed that the integral BC structure was well maintained. Considering all these advantages, the BC monolayer can be regarded as a good candidate for Mg-ion batteries.

Computational methods

To study electronic properties of investigated materials first-principles calculations are performed within the framework of the density-functional theory (DFT)[43] as implemented in the Quantum Espresso package[44, 45]. The generalized gradient approximation of Perdew-Burke-Ernzerhof (GGA-PBE) is used for the exchange-correlation functional together with projector-augmented wave (PAW) method. After proper convergence tests we obtained well converged values for the kinetic energy cutoff of the wavefunction equal to 60 Ry and the kinetic energy cutoff for charge density equal to 600 Ry.

The pristine model consists of a supercell of BC (containing 36 atoms). To avoid the interaction between neighboring layers, a vacuum layer of 20 Å in the z-direction is introduced. The van der Waals (vdW) interaction with a DFT-D correction of Grimme was considered[46]. The optimized atomic structures, were obtained by fully relaxing of both atomic positions as well as cell parameters by using the Broyden-Fretcher-Goldfarb-Shanno (BFGS) quasi-Newton algorithm until all forces are smaller than 0.01 eV/Å. The Brillouin zone is sampled utilizing a k-mesh in the Monkhorst-Pack scheme. The thermal stability of monolayer BC was examined by performing ab-initio molecular dynamics (AIMD) simulations at 300 K with a time step of 1 fs using CP2K software[47].

To visualize the results, we have used the Gnuplot (version 5.4.2)[48] and Vesta (version 3.5.7) software[49].