Li-Decorated β12-Borophene as Potential Candidates for Hydrogen Storage: A First-Principle Study

The hydrogen storage properties of pristine β12-borophene and Li-decorated β12-borophene are systemically investigated by means of first-principles calculations based on density functional theory. The adsorption sites, adsorption energies, electronic structures, and hydrogen storage performance of pristine β12-borophene/H2 and Li-β12-borophene/H2 systems are discussed in detail. The results show that H2 is dissociated into Two H atoms that are then chemisorbed on β12-borophene via strong covalent bonds. Then, we use Li atom to improve the hydrogen storage performance and modify the hydrogen storage capacity of β12-borophene. Our numerical calculation shows that Li-β12-borophene system can adsorb up to 7 H2 molecules; while 2Li-β12-borophene system can adsorb up to 14 H2 molecules and the hydrogen storage capacity up to 10.85 wt %.


Introduction
As the gap between energy supply and demand has become increasingly prominent, sources of renewable energy has been investigated urgently. Hydrogen is an inexhaustible source of clean energy, making it important for society to develop and utilize this energy [1,2]. Hydrogen storage is one of the most critical technical problems in the development of hydrogen energy sources. The average adsorption energy of the ideal physical hydrogen storage method should be between chemical and physical adsorption energy (0.1~0.8 eV) [3,4]. The US Department of Energy (DOE) and the International Energy Agency (IEA) reported that the ideal hydrogen storage capacity should be greater than 5.5 wt % [5]. At present, one of the best types of hydrogen storage methods involves physical adsorption, which results in low adsorption heat, small activation energy, fast hydrogen adsorption and desorption, and reversible cyclization performance. Carbon nanomaterials have become a hotspot of physical hydrogen storage materials due to their characteristics of a large specific surface area, good adsorption kinetic properties and reversible hydrogen storage [6,7]. However, clean carbon nanomaterials adsorb H 2 molecules with weak binding capacity, which means that they have low hydrogen storage capacity and are not ideal. Therefore, it is essential to find a suitable physical adsorbent.
Recently, 2D (two-dimensional) borophene created from Boron elements was artificially synthesized [8]. Although there are many theoretical studies about the possible 2D borophene structure [9], only three types of stable structures have been synthesized for borophene so far [8,10]. Borophene's unique metal properties, mechanical properties, and optical properties have been extensively studied [11][12][13][14][15], but only a few studies have considered its hydrogen storage properties. Borophene and graphene [16] have a similar 2D planar structure with a large specific surface area. Moreover, the relative atomic mass of B atom is smaller than the relative atomic mass of C atom. Therefore, we suspect that borophene has better hydrogen storage properties than graphene (it exhibits a triangular lattice with different periodic arrangements and is flat without obvious vertical undulation). Feng et al. [10] reported that β 12 -borophene is more stable than the other two types of borophene. Chen et al. [17] used the first-principles method to study the hydrogen storage properties of Ca-β 12 -borophene and found that it has a larger adsorption energy compared to other types of borophene. Therefore, we selected β 12 -borophene as the research focus. In this work, we performed theoretical calculations for the hydrogen storage properties of pure β 12 -borophene and Li-β 12 -borophene based on the first-principle study. We found that H 2 molecules were completely dissociated into two H atoms that were adsorbed on the B-B bridge sites to form H-B covalent bonds, thus making it difficult to dissociate. Comparison of the improvement in hydrogen storage properties of graphene found that the graphene surface was modified by alkali metal (Li, Na, K) [18], alkali-earth metal (Ca) [19], light metal (Al) [20] and transition metals (Cu, Pd, Y) [21][22][23][24], which can change the chemical activity of the graphene surface and could effectively change the hydrogen storagecapability. The quality of alkali metal (Li atoms) is very light, which helps to enhance the hydrogen storage density [25]. The transition metal atom-modified nanostructures are highly reactive and can easily cause the dissociation of H 2 molecules, which is detrimental to the reversible storage of hydrogen [26]. Therefore, we selected the lightest Li atom to modify the β 12 -borophenen. H 2 adsorbed on Li-β 12 -borophene by physical adsorption, which improved the reversible hydrogen storage performance and significantly increased the amount of hydrogen storage. It is expected that this work can provide theoretical support for β 12 -borophene being used as hydrogen materials.

Computational Methods
All density functional theory (DFT) calculations are carried out using the Cambridge Sequential Total Energy Package (CASTEP) [27], and the DFT evaluation is based on the plane-wave expansion. We use the Generalized Gradient Approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional [28] to describe exchange and correlation effects. The van der waals forces of H 2 adsorption on Li-β 12 -borophenen is modified by DFT-D methods. While the DFT-D perform poorly for energetics in layered materials [29], it is important to deal with the molecules adsorption system. We select the Ultrasoft Pseudopotential [30] to describe the interaction of electron-ion, and the electron wave functions are expanded by plane wave. The convergence tolerance energy, the force on each atoms and displacement convergence criterions are set to 5.0 × 10 −6 eV/atom, 0.01 eV/Å and 0.001 Å, respectively. All atoms are relaxed in our calculations. In order to eliminate the interaction of the interlayer we select the vacuum thickness 20 Å. Considering the calculation accuracy and computational efficiency, all calculations are using a cutoff energy of 600 eV and 9 × 16 × 5 k-point mesh in the Brillouin zone.
The adsorption energy (E ads ) and average adsorption energy (E ads ) of H 2 adsorption on Li-β 12 -borophene are calculated by the following formulas [31]: The average adsorption energy of Li atom on β 12 -borophenec [32] is defined as: where E iH2+nLi+β12−borophene , E (i−1)H2+nLi+β12−borophene and E nLi+β12−borophene are the total energy of the n Li-β 12 -borophene with i, i − 1 H 2 molecules and β 12 -borophene with n Li atoms, respectively. E β12−borophene , E Li and E H2 are the total energy of the β 12 -borophene, free Li atom and an isolated H 2 , respectively. n is the number of adsorbed Li atoms.

H2 Adsorption on β12-Borophene
The optimized lattice parameters of the primitive cell of β12-borophenen are a = 5.069 Å, b = 2.929 Å, agree well with the experimental result (a = 5 Å and b = 3 Å) [10] and other theoretical calculation results [33][34][35]. In our follow calculations, we choose a 2 × 2 unit cell (see Figure 1) of the β12-borophenen containing 20 B atoms (in Figure 1) to investigate the hydrogen storage adsorbed on β12-borophene. We first investigated the adsorption behavior of one H2 molecule on β12-borophene. The H2 molecule is initially placed in a parallel or vertical direction at different positions of the β12-borophene plane. We found that there are five stable adsorbed configurations in total, as illustrated in Figure 2. In all cases, the H2 molecule is dissociated into two separate H atoms after adsorption and the distance between the H atoms will change from 0.753 to 2.366 Å. Furthermore, the distance between H and its nearest B atom (rH-B) greatly increased from 1.217 to 1.358 Å. The most stable case among all the five adsorption configurations is shown in Figure 2a. In this case, the H2 molecule dissociated into two H atoms that are adsorbed on the B1-B3 and B2-B4 bridge sites with an value of −0.536 eV, which is related to chemical adsorption. Mulliken analysis demonstrates that there is 0.23 e − transferred from B to H, which occurs mainly in the H 1s orbital and B 2p orbital. The B-H bond population is 0.4, indicating it is a covalent bond, with difficult desorption of the β12-borophene/H2 system. In addition, we further studied the transition states of the stable adsorption configurations by combining linear synchronous and quadratic synchronous transits [21,36]. We found that the most stable adsorption configurations of the activation energy barrier from the reactant to transition state was 1.584 eV, which is smaller than the activation energy barrier of other adsorption methods, indicating it was difficult for the reaction of the H2 molecules adsorbed on the surface to take place. We first investigated the adsorption behavior of one H 2 molecule on β 12 -borophene. The H 2 molecule is initially placed in a parallel or vertical direction at different positions of the β 12 -borophene plane. We found that there are five stable adsorbed configurations in total, as illustrated in Figure 2. In all cases, the H 2 molecule is dissociated into two separate H atoms after adsorption and the distance between the H atoms will change from 0.753 to 2.366 Å. Furthermore, the distance between H and its nearest B atom (r H-B ) greatly increased from 1.217 to 1.358 Å. The most stable case among all the five adsorption configurations is shown in Figure 2a. In this case, the H 2 molecule dissociated into two H atoms that are adsorbed on the B1-B3 and B2-B4 bridge sites with an E ads value of −0.536 eV, which is related to chemical adsorption. Mulliken analysis demonstrates that there is 0.23 e − transferred from B to H, which occurs mainly in the H 1s orbital and B 2p orbital. The B-H bond population is 0.4, indicating it is a covalent bond, with difficult desorption of the β 12 -borophene/H 2 system. In addition, we further studied the transition states of the stable adsorption configurations by combining linear synchronous and quadratic synchronous transits [21,36]. We found that the most stable adsorption configurations of the activation energy barrier from the reactant to transition state was 1.584 eV, which is smaller than the activation energy barrier of other adsorption methods, indicating it was difficult for the reaction of the H 2 molecules adsorbed on the surface to take place.

H2 Adsorption on β12-Borophene
The optimized lattice parameters of the primitive cell of β12-borophenen are a = 5.069 Å, b = 2.929 Å, agree well with the experimental result (a = 5 Å and b = 3 Å) [10] and other theoretical calculation results [33][34][35]. In our follow calculations, we choose a 2 × 2 unit cell (see Figure 1) of the β12-borophenen containing 20 B atoms (in Figure 1) to investigate the hydrogen storage adsorbed on β12-borophene. We first investigated the adsorption behavior of one H2 molecule on β12-borophene. The H2 molecule is initially placed in a parallel or vertical direction at different positions of the β12-borophene plane. We found that there are five stable adsorbed configurations in total, as illustrated in Figure 2. In all cases, the H2 molecule is dissociated into two separate H atoms after adsorption and the distance between the H atoms will change from 0.753 to 2.366 Å. Furthermore, the distance between H and its nearest B atom (rH-B) greatly increased from 1.217 to 1.358 Å. The most stable case among all the five adsorption configurations is shown in Figure 2a. In this case, the H2 molecule dissociated into two H atoms that are adsorbed on the B1-B3 and B2-B4 bridge sites with an value of −0.536 eV, which is related to chemical adsorption. Mulliken analysis demonstrates that there is 0.23 e − transferred from B to H, which occurs mainly in the H 1s orbital and B 2p orbital. The B-H bond population is 0.4, indicating it is a covalent bond, with difficult desorption of the β12-borophene/H2 system. In addition, we further studied the transition states of the stable adsorption configurations by combining linear synchronous and quadratic synchronous transits [21,36]. We found that the most stable adsorption configurations of the activation energy barrier from the reactant to transition state was 1.584 eV, which is smaller than the activation energy barrier of other adsorption methods, indicating it was difficult for the reaction of the H2 molecules adsorbed on the surface to take place.

The Adsorption Structure of Li-β 12 -Borophene
It is well known that doping alkali metal atoms to modify hydrogen storage materials may can greatly improve the hydrogen storage properties and increase the hydrogen storage capacity. Specially, lithium (Li) has been widely employed to functionalize 2D materials and improve the hydrogen storage ability. Therefore, in the following section, we chose to add Li atoms to modify the hydrogen storage properties of β 12 -borophene.
We examined the adsorption of Li atoms on pure β 12 -borophene. After optimization, we obtained three different stable adsorption structures, as shown in Figure 3a-c. Similar to Li-decorating graphene [37], the most favorable Li adsorption site on β 12 -borophene is the hollow center of B ring (Figure 3a).

The Adsorption Structure of Li-β12-Borophene
It is well known that doping alkali metal atoms to modify hydrogen storage materials may can greatly improve the hydrogen storage properties and increase the hydrogen storage capacity. Specially, lithium (Li) has been widely employed to functionalize 2D materials and improve the hydrogen storage ability. Therefore, in the following section, we chose to add Li atoms to modify the hydrogen storage properties of β12-borophene.
We examined the adsorption of Li atoms on pure β12-borophene. After optimization, we obtained three different stable adsorption structures, as shown in Figure 3a-c. Similar to Li-decorating graphene [37], the most favorable Li adsorption site on β12-borophene is the hollow center of B ring (Figure 3a). Doping alkali metal atoms to modify hydrogen storage materials requires the average adsorption energy of the metal atoms on the substrate to be greater than the cohesive energy of the metal atoms in the solid form [38]. The average adsorption energy of Li atom on the β12-borophene is −3.088 eV, which is significantly greater than the cohesive energy of −1.795 eV of Li [39]. This indicates that Li atoms can be dispersed uniformly on β12-borophene, instead of forming metal clusters.
There are three stable adsorption structures of two Li atoms after adsorption on the β12borophene as shown in Figure 3d-f, respectively. One of the most stable adsorption sites involves the two Li atoms being located on both sides of the same B ring. The distance between Li and the nearest B is 2.298 Å. The average adsorption energy is −3.106 eV, which is larger than the cohesive energy of Doping alkali metal atoms to modify hydrogen storage materials requires the average adsorption energy of the metal atoms on the substrate to be greater than the cohesive energy of the metal atoms in the solid form [38]. The average adsorption energy of Li atom on the β 12 -borophene is −3.088 eV, which is significantly greater than the cohesive energy of −1.795 eV of Li [39]. This indicates that Li atoms can be dispersed uniformly on β 12 -borophene, instead of forming metal clusters.
There are three stable adsorption structures of two Li atoms after adsorption on the β 12 -borophene as shown in Figure 3d-f, respectively. One of the most stable adsorption sites involves the two Li atoms being located on both sides of the same B ring. The distance between Li and the nearest B is 2.298 Å. The average adsorption energy is −3.106 eV, which is larger than the cohesive energy of Li atoms. After optimization, the relaxation of β 12 -borophene is very small. Each Li atom in the Li-β 12 -borophene system is an active adsorption site, allowing a large number of H 2 molecules to be adsorbed around the Li atom in order to significantly increase the hydrogen storage capacity.
The charge transfer between atoms can be analyzed by Mulliken analysis [40], which shows that the charge was transferred from Li to B. From the Partial Densities of States (PDOS) of the Li-β 12 -borophene structure in Figure 4, we found the peak of B atom's 2p orbital overlaps with the peak of the Li atom's 1s orbital. This suggests a strong hybridization between B and Li atoms. A similar binding mechanism has also been confirmed in other metal-modified nanostructures [41]. In addition, it can be seen from the PDOS that the metal properties of β 12 -borophene did not change after modification of Li atom. Li atoms. After optimization, the relaxation of β12-borophene is very small. Each Li atom in the Li-β12borophene system is an active adsorption site, allowing a large number of H2 molecules to be adsorbed around the Li atom in order to significantly increase the hydrogen storage capacity. The charge transfer between atoms can be analyzed by Mulliken analysis [40], which shows that the charge was transferred from Li to B. From the Partial Densities of States (PDOS) of the Li-β12-borophene structure in Figure 4, we found the peak of B atom's 2p orbital overlaps with the peak of the Li atom's 1s orbital. This suggests a strong hybridization between B and Li atoms. A similar binding mechanism has also been confirmed in other metal-modified nanostructures [41]. In addition, it can be seen from the PDOS that the metal properties of β12-borophene did not change after modification of Li atom.

Adsorption of H2 Molecules on Li-β12-Borophene
We investigated the adsorption properties of H2 molecules on Li-β12-borophene. Figure 5 shows the optimized geometries of 1-7 H2 molecules adsorbed on the Li-modified β12-borophene. Table 1 lists the adsorption energy and average adsorption energy calculated by the GGA PBE functional and DFT-D methods. First, we investigate the adsorption sites of H2 molecules on Li-β12-borophene. For the first adsorbed H2 molecules, many adsorption sites were considered in order to find the most stable site. The most stable site involves H2 being parallel to the β12-borophene plane, which is opposite to the H2 vertical adsorption on Ca-β12-borophene [17]. After adsorption, the corresponding rH-H of the adsorbed H2 is 0.756 Å, which is larger than the distance of free H2 (0.753 Å). To investigate the maximum storage capacity of single Li atom-modified β12-borophene, more H2 was added around Li gradually. The minimum distance between the H and Li atom are range of 2.164 to 6.368 Å. The first four H2 molecules were parallel to the β12-borophene and were around the Li atom at the same level. When the fifth H2 molecule was added to the system, two H2 molecules moved to an upper layer after relaxation. This may be due to the limited space around the Li atom and the repulsive interactions between the adsorbed H2. The average adsorption energy slowly reduced from −0.385 to −0.210 eV/H2 due to the strong steric interactions between the adsorbed H2. Interestingly, the adsorption energy suddenly rose to -0.388 eV after the second H2 molecule was added to the system. With an increase in the number of H2 molecules, the H2 molecules becomes further away from the Li atom and the adsorption weakens. The average adsorption energy was at its minimum (−0.210 eV/H2) when the seventh H2 molecule was adsorbed. At this time, the hydrogen storage capacity reached 5.90 wt %, which exceeded the ideal hydrogen storage capacity (over 5.5 wt %). In order to further increase the hydrogen storage capacity, we added two Li atoms to decorate the β12-borophene to adsorb H2 molecules. 2Li-β12-borophene can adsorb up to 14 H2 molecules and the minimum average adsorption energy is −0.220 eV. The hydrogen storage capacity can reach up to 10.85 wt %, which is larger than the hydrogen storage capacity with a gravimetric hydrogen density of 9.5 wt % of the

Adsorption of H 2 Molecules on Li-β 12 -Borophene
We investigated the adsorption properties of H 2 molecules on Li-β 12 -borophene. Figure 5 shows the optimized geometries of 1-7 H 2 molecules adsorbed on the Li-modified β 12 -borophene. Table 1 lists the adsorption energy and average adsorption energy calculated by the GGA PBE functional and DFT-D methods. First, we investigate the adsorption sites of H 2 molecules on Li-β 12 -borophene. For the first adsorbed H 2 molecules, many adsorption sites were considered in order to find the most stable site. The most stable site involves H 2 being parallel to the β 12 -borophene plane, which is opposite to the H 2 vertical adsorption on Ca-β 12 -borophene [17]. After adsorption, the corresponding r H-H of the adsorbed H 2 is 0.756 Å, which is larger than the distance of free H 2 (0.753 Å). To investigate the maximum storage capacity of single Li atom-modified β 12 -borophene, more H 2 was added around Li gradually. The minimum distance between the H and Li atom are range of 2.164 to 6.368 Å. The first four H 2 molecules were parallel to the β 12 -borophene and were around the Li atom at the same level. When the fifth H 2 molecule was added to the system, two H 2 molecules moved to an upper layer after relaxation. This may be due to the limited space around the Li atom and the repulsive interactions between the adsorbed H 2 . The average adsorption energy slowly reduced from −0.385 to −0.210 eV/H 2 due to the strong steric interactions between the adsorbed H 2 . Interestingly, the adsorption energy suddenly rose to -0.388 eV after the second H 2 molecule was added to the system. With an increase in the number of H 2 molecules, the H 2 molecules becomes further away from the Li atom and the adsorption weakens. The average adsorption energy was at its minimum (−0.210 eV/H 2 ) when the seventh H 2 molecule was adsorbed. At this time, the hydrogen storage capacity reached 5.90 wt %, which exceeded the ideal hydrogen storage capacity (over 5.5 wt %). In order to further increase the hydrogen storage capacity, we added two Li atoms to decorate the β 12 -borophene to adsorb H 2 molecules. 2Li-β 12 -borophene can adsorb up to 14 H 2 molecules and the minimum average adsorption energy is −0.220 eV. The hydrogen storage capacity can reach up to 10.85 wt %, which is larger than the hydrogen storage capacity with a gravimetric hydrogen density of 9.5 wt % of the Ca-β 12 -borophene/H 2 system [17]. The optimized structure is shown in Figure 5h-n. The average adsorption energy (E ads ) is in the range of −0.381 to −0.220 eV/H 2 , which is necessary for practical application [3,4]. In addition, the calculated adsorption energy and average energy of H 2 adsorption by DFT-D method are larger than those calculated by the GGA PBE functional, but the overall adsorption method has not changed.
Materials 2017, 10, 1399 6 of 10 adsorption energy ( ) is in the range of −0.381 to −0.220 eV/H2, which is necessary for practical application [3,4]. In addition, the calculated adsorption energy and average energy of H2 adsorption by DFT-D method are larger than those calculated by the GGA PBE functional, but the overall adsorption method has not changed.   The density of states (DOS) reflects the number of states of the unit energy, which is important in further understanding the interaction between H 2 and Li-β 12 -borophene. The partial density of states (PDOS) of Li-β 12 -borophene/H 2 is shown in Figure 6. Obvious hybridizations between the Li-s orbit and H-s orbit can be found in 2.0 eV~5.0 eV, which demonstrates a strong interaction between H 2 and Li atoms. With an increase in the number of H 2 molecules, the peak values of H 2 molecules become smaller and further away from the fermi level. This indicates that the interaction between H 2 molecules and Li-β 12 -borophene weakens, which is consistent with the average adsorption energy becoming smaller. Another overlap between the B-p and H-s orbits was found at -10.0 eV~5.0 eV. Upon the adsorption of the second H 2 molecule, the H 1s orbit peaks move to the lift, implying an increased stability in the system. This is consistent with the increase in the Eads value after the second H 2 addition. With an increase in the number of H 2 molecules (an expected in the second H 2 ), the H-s orbits move to the right and the peak values become smaller, which indicates that the interaction between H 2 and β 12 -borophene becomes increasingly weaker. This conclusion is consistent with the decrease in the average adsorption energy (E ads ). The B-p and Li-s orbits also have hybridization, which implies an interaction between the Li and B atoms. The peaks near and below the fermi surface are mostly contributed by the B-s orbits, which means that the H 2 molecules and Li atom have less influence on the β 12 -borophene. The comparison of the PDOS of a single Li-β 12 -borophene show that the interaction between Li atom and β 12 -borophene is weakened due to the adsorption of H 2 molecules. The PDOS of two Li-β 12 -borophene/H 2 consistent with this analysis.

Electronic Properties of Li-β12-Borophene/H2
The density of states (DOS) reflects the number of states of the unit energy, which is important in further understanding the interaction between H2 and Li-β12-borophene. The partial density of states (PDOS) of Li-β12-borophene/H2 is shown in Figure 6. Obvious hybridizations between the Li-s orbit and H-s orbit can be found in 2.0 eV~5.0 eV, which demonstrates a strong interaction between H2 and Li atoms. With an increase in the number of H2 molecules, the peak values of H2 molecules become smaller and further away from the fermi level. This indicates that the interaction between H2 molecules and Li-β12-borophene weakens, which is consistent with the average adsorption energy becoming smaller. Another overlap between the B-p and H-s orbits was found at -10.0 eV~5.0 eV. Upon the adsorption of the second H2 molecule, the H 1s orbit peaks move to the lift, implying an increased stability in the system. This is consistent with the increase in the Eads value after the second H2 addition. With an increase in the number of H2 molecules (an expected in the second H2), the H-s orbits move to the right and the peak values become smaller, which indicates that the interaction between H2 and β12-borophene becomes increasingly weaker. This conclusion is consistent with the decrease in the average adsorption energy ( ). The B-p and Li-s orbits also have hybridization, which implies an interaction between the Li and B atoms. The peaks near and below the fermi surface are mostly contributed by the B-s orbits, which means that the H2 molecules and Li atom have less influence on the β12-borophene. The comparison of the PDOS of a single Li-β12-borophene show that the interaction between Li atom and β12-borophene is weakened due to the adsorption of H2 molecules. The PDOS of two Li-β12-borophene/H2 consistent with this analysis. The bonding strength between atoms can be quantitatively analyzed based on the Mulliken charge population and bond population. Table 2 shows the Mulliken charge population before and after one H 2 moleculebecomes absorbed on the Li-β 12 -borophene. H (1) and H (2) represent the two H atoms of the adsorbed H 2 molecule; while B1, B5 and B6 are three B atoms that transfer the greatest amount of charge in the β 12 -borophene (as shown in Figure 1). The two H atoms have charges of 0.06 e and 0.05 e, respectively. In contrast, the Li atom loses 1.40 e, which occurs mainly in the H and Li atomic orbits. The Li atom transfers charge to the H 2 molecules, resulting in the H 2 molecules carrying more negative charge and Li atom showing positive charge. The interaction between the H 2 molecules and the Li atom is consistent with the conclusion of the PDOS analysis. In addition, the B atoms obtains charge, with this charge transfer mainly occurring in the B-2p orbits and H-s orbits. This is in contrast with the Mulliken charge population of the β 12 -borophene/H 2 , in which the charge transfer mainly occurs in H and B atoms forming a covalent bond of H-B that is not favorable for the desorption of H 2 . Due to the β 12 -borophene being modified by Li atoms, H 2 molecules and B atoms only have small interactions, resulting in the H 2 molecules physically adsorbing on the Li-β 12 -borophene. This is conducive for H 2 desorption and increases the hydrogen storage capacity.

Conclusions
In summary, we performed a study on hydrogen storage properties of pure β 12 -borophene and Li-decorated β 12 -borophene through DFT calculations. It is found that H 2 molecules are mainly adsorbed on pure β 12 -borophene as chemical adsorption with an adsorption energy of −0.536 eV. The H 2 molecules are dissociated into two H atoms, which tend to the bridge of two B site and the H-B bond to form covalent bond. In order to improve the hydrogen storage performance of pure β 12 -borophene and increase the hydrogen storage capacity, we use the Li atom to modify the β 12 -borophene. It is found that a single Li atom adsorbed on the center of Boron ring with the adsorption energy −3.088 eV, the Li-β 12 -borophene can adsorb up to 7 H 2 molecules with the average adsorption energy of −0.210 eV/H 2 . The charge transfer of the Li-β 12 -borophene/H 2 is that H and B atoms lose electron, Li atom get electron. We use two Li atoms to modify β 12 -borophene to increase its hydrogen storage capacity. It is find that the two Li atoms are located at the same position on both sides of the same boron hole. 2Li-β 12 -borophene system can adsorb up to 14 H 2 molecules and the hydrogen storage capacity up to 10.85 wt %. The average adsorption energy is range of −0.381 to −0.220 eV/H 2 , which is necessary for practical application [3,4].