Open Access Article

This Open Access Article is licensed under a

Creative Commons Attribution 3.0 Unported Licence

Andrey A.
Kistanov
^{ab},
Yongqing
Cai
*^{b},
Kun
Zhou
*^{a},
Narasimalu
Srikanth
^{c},
Sergey V.
Dmitriev
^{de} and
Yong-Wei
Zhang
*^{b}
^{a}School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798, Singapore. E-mail: kzhou@ntu.edu.sg
^{b}Institute of High Performance Computing, Agency for Science, Technology and Research, Singapore 138632, Singapore. E-mail: caiy@ihpc.a-star.edu.sg; zhangyw@ihpc.a-star.edu.sg
^{c}Energy Research Institute@NTU, Singapore 637141, Singapore
^{d}Institute for Metals Superplasticity Problems, Russian Academy of Sciences, Ufa 450001, Russia
^{e}National Research Tomsk State University, Tomsk 634050, Russia

Received
2nd September 2017
, Accepted 13th December 2017

First published on 13th December 2017

Recently synthesized two-dimensional (2D) boron, borophene, exhibits a novel metallic behavior rooted in the s–p orbital hybridization, distinctively different from other 2D materials such as sulfides/selenides and semi-metallic graphene. This unique feature of borophene implies new routes for charge delocalization and band gap opening. Herein, using first-principles calculations, we explore the routes to localize the carriers and open the band gap of borophene via chemical functionalization, ribbon construction, and defect engineering. The metallicity of borophene is found to be remarkably robust against H- and F-functionalization and the presence of vacancies. Interestingly, a strong odd–even oscillation of the electronic structure with width is revealed for H-functionalized borophene nanoribbons, while an ultra-high work function (∼7.83 eV) is found for the F-functionalized borophene due to its strong charge transfer to the atomic adsorbates.

So far, various atomic models for boron clusters and 2D models have been theoretically proposed.^{13–15} Concerning the structure of borophene, recent experiments reported two dramatically different atomic structures on the Ag substrate: the closely packed structure^{1} and the hole-containing structure.^{2} While experiments claimed the stability of a borophene sheet supported by the Ag substrate, the stability of free-standing borophene is still elusive. First-principles phononic calculations showed that the long-wavelength imaginary modes always exist in isolated triangular sheets and other polymorphs like β_{12} sheets (refer to phonon dispersion in the ESI in ref. 14). Physically or chemically bound species were found to promote the stability of borophene.^{16,17} Recent theoretical studies proposed a viable approach of using substrates^{18,19} or chemical functionalization^{20} to tune the adsorption and energetics of borophene. In principle, proper truncation of the borophene sheet by breaking the lattice periodicity can eliminate these long-wavelength softening phonons, thus stabilizing structures without the need for chemical functionalization.

As the boron atom has three valence electrons, it needs to pair with five additional electrons to satisfy the octet rule. However, in borophene, each boron atom forms bonds with six neighbors, thus favoring the metallic phase according to the band theory. For electronic applications of borophene, an intriguing issue is the localization of its itinerant electrons and ultimately its band gap opening. Approaches for the band gap opening in 2D materials can be categorized into two groups: (i) quantum confinement induced by the construction of finite-sized structures, like ribbons, edges and dots,^{21–23} and (ii) chemical functionalization.^{24,25} However, the effectiveness of both approaches in the borophene band gap opening remains unclear. The formation of ionic bonds in a high-pressure boron phase^{26} suggests a different charge distribution in boron materials in comparison with other 2D materials, especially graphene. The failure of the octet rule in pure borophene implies a new mechanism of the charge localization/delocalization, which is yet to be understood.

In this work, we explore the routes for the charge localization and the band gap opening of borophene via chemical functionalization, defect engineering, and ribbon construction by using first-principles calculations. We show that the metallicity of borophene is remarkably robust against H- and F-functionalization and the presence of vacancies. Interestingly, a strong odd–even oscillation of the electronic structure with the ribbon width is revealed for H-functionalized borophene nanoribbons, and band gap opening occurs only for a specific type of ribbon. Owing to the high density of states (DOS) near the Fermi level, a record-high work function is found for the F-functionalized borophene. These unique features of chemically functionalized borophene sheets and ribbons may indicate many interesting applications.

Fig. 1 The atomic configuration and the band structure of (a) pristine, (b) H-functionalized, and (c) F-functionalized borophene. The component of the states scales with the radius of the black and cyan circles for H- and F-functionalized cases, respectively. (d) Inward dipole layer built in F-functionalized borophene due to the charge transfer (top panel). Comparison of the work functions of pristine, H-functionalized, and F-functionalized borophene (calculated using the HSE method) with those of the common metals and graphene (middle panel). Schematic plot of the integration of the chemically functionalized borophene for improving the efficiency of the injection and transport of carriers in nanoelectronic devices (bottom panel). The work functions of graphene and other metals are adopted from ref. 34 and 35. |

Fig. 1b and c show the optimized structures and the electronic band structures of pristine, H-functionalized, and F-functionalized borophene sheets with adsorption on both sides. Consistent with previous studies,^{29,30} the pristine borophene shows a metallic behavior (see Fig. 1a) with the Fermi level crossing the electronic levels. Owing to the puckered zigzag structure (side view), the pristine borophene shows an anisotropic electronic structure with half-filled bands along the Γ–X (line-edge) direction but an energy gap along the Y–Γ (zigzag) direction. Such an intriguing electronic property implies orientationally different quantum confinement effects in borophene, which may lead to angle-dependent plasmonic behavior in this ultrathin metallic sheet.

Compared with the pristine borophene, the H-functionalized borophene only shows a small change in its lattice constant. In contrast, an F-functionalized borophene sheet shows significantly deformed B–B bonds, leading to a large distortion of the host borophene lattice. Such a difference can be attributed to the strong subtraction of electrons from the sheet to the anionic adsorbents, in comparison with that of the H- group. Interestingly, the F-functionalization of borophene induces changes in the lattice constant along the zigzag direction, around 1%, while the change along the line-edge direction is only slight. Concerning the electronic properties, surprisingly, all these highly chemically functionalized borophene sheets remain metallic. This is in strong contrast to graphene where hydrogenation and fluorination are well known to lead to band gap openings. Therefore, the charge localization, which is necessary for the band gap opening, is hard to induce via the surface chemical functionalization of borophene. A projected band analysis of the H-functionalized borophene (see Fig. 1b) shows that the states around the Fermi level are predominantly populated with the H states. In contrast, the F states in the F-functionalized borophene are largely distributed below the Fermi level (see Fig. 1c). In addition, upon the functionalization by these atoms, partially occupied levels are formed along the Y–Γ direction, which is otherwise empty in pristine borophene. Therefore, selective atomic functionalization enables the modulation of the anisotropy in the electronic properties of borophene. It is noted that the in-plane doping within the atomic sheet of borophene is able to open the band gap.^{31}

The above-mentioned robust metallicity in these surface functionalized borophene sheets is absent in graphene and transition metal dichalcogenides (TMDs). The finite density of these free carriers at the Fermi level suggests that borophene and its functionalized derivatives are promising for applications as interconnecting and field-emitting materials. Since the work function, which quantifies the ability of electrons to move from the surface of a material to vacuum, is critically important for field emission and rectification of conducting barriers;^{32,33} in the following, we examine the changes in the work functions of these functionalized borophene sheets.

Fig. 1d shows the energetics diagram of the work functions for various functionalized borophene sheets in comparison with other common bulk metals and graphene. From the diagram, the following important features can be identified. Firstly, the work function of a pristine borophene sheet is 5.31 eV (obtained via HSE calculation), which is larger than that of most listed metals, except Pt. Moreover, the work function of pristine borophene is also higher than that of graphene (∼4.5 eV).^{34} This is surprising since a carbon atom has larger electronegativity than a boron atom. The higher work function of borophene could be attributed to the nature of atomic states around the Fermi level. Borophene mainly consists of in-plane s–p hybridized (σ) states, which lie lower than the out-of-plane p_{z} (π) states in the graphene case. Thus, an electron in borophene is harder to knock out than that in graphene. Secondly, the work function of borophene increases slightly to 5.88 eV for the H-functionalized borophene and dramatically to 7.83 eV for the F-functionalized borophene. The underlying origin may arise from the strong dipole layer pointing inward towards the central borophene layer due to its transfer of electrons to the functionalizing atoms (see Fig. 1d). In other 2D materials, like graphene and TMDs, the density of electrons is negligible at the Fermi level, which means that the magnitude of the dipole layer is modest upon the chemical functionalization due to the limited charge transfer. In contrast, a borophene layer has a considerably high density of carriers at the Fermi level because of its intrinsic metallicity, giving rise to a pronounced charge flow and a built-in dipole layer. This great tunability in the work function suggests that the chemically functionalized borophene sheets can be used as a buffer layer for reducing the contact resistance and Schottky barrier at the interface. In addition, the high work function in the F-functionalized borophene sheet is particularly useful for electron collection and hole injection.

Next, we investigate the effects of the atomic vacancies on the electronic properties of borophene (see Fig. 2). We consider both monovacancies (MV) and divacancies (DV) with the loss of one and two boron atoms in the 6 × 5 supercell (56 atoms), respectively. In perfect borophene, each boron atom has a coordination number of six. With the creation of an MV, six peripheral atoms become fivefold coordinated and the defect core has a local symmetry of C_{2v} (see Fig. 2b). For the DV, two different configurations exist: horizontal (C_{2v} symmetry) and tilted (C_{i} symmetry) DV, depending on the relative direction of the deleted boron dimer to the a lattice (see Fig. 2c and d). The horizontal DV is slightly more stable with the energy of 0.07 eV lower than that of the tilted DV. In the DV structures, there are two new-born fourfold boron atoms in the edge in addition to the six fivefold corner atoms. The formation energy (E_{form}), which is defined as E_{form} = E_{perfect} − E_{defect} − nE_{atom}, where E_{perfect}, E_{defect} and E_{atom} are the total energies of perfect and defective borophene, and the single boron atom, respectively, and n is the total number of the removed atoms. The calculated values of E_{form} for MV, and horizontal and tilted DVs are −5.61, −12.58, and −12.51 eV, respectively. The relatively high values of E_{form} suggest that it is difficult to form isolated vacancies in borophene. However, the energy cost for forming an MV in borophene is smaller than that in graphene (−7.57 eV) and that of a boron MV in h-BN (from −7.50 to −10.20 eV).^{36}

It is well known that the electronic properties around the vacancy core may change dramatically owing to the breaking of the lattice periodicity.^{37–39}Fig. 2a shows the DOS in the perfect, MV-, and DV-containing borophene sheets. It is seen that the metallicity of these borophene sheets is robust against the presence of vacancies. Interestingly, the Fermi level significantly shifts upwards for the MV and DV cases compared with perfect borophene (see the arrows in Fig. 2a). The bands of vacancy-containing borophene sheets are non-zero at the Fermi level and are mainly contributed by the p_{z} orbitals of B atoms. Different from the new localized states formed in the band gap-associated vacancies in MoS_{2}^{40} and phosphorene^{41} cases, there are no peaks related to the dangling bond states due to vacancies in borophene. However, by comparing the local DOS (LDOS) of peripheral atoms in the defect core and those of atoms far from the vacancy center (see Fig. 2b–d, bottom panel), we find that there is a significant difference in the LDOS profiles, suggesting that the states are renormalized greatly in the defect center.

To explore the routes to introduce localized states in the intrinsic metallic borophene, we examine one-dimensional (1D) nanostructures of borophene, that is, nanoribbons. Both the line- and zigzag-edge ribbons are considered. Following the normal nomenclature of nanoribbons in graphene, MoS_{2}, and phosphorene, the BNR along the line-edge or zigzag direction is named LE-N_{d} BNR or ZZ-N_{z} BNR according to the number of B–B dimer lines (N_{d}) or zigzag chains (N_{z}) across the ribbon width. Fig. 3a and b show the atomic models and the band structures of pristine LE-9 and LE-10 BNR, selected as the representatives for the odd- and even-width BNRs, respectively. It is seen that both types of BNRs are metallic without any localized states and band gap opening, which is in a sharp contrast with graphene.^{42} The orbital-resolved band structures suggest that the out-of-plane p_{z} and in-plane p_{y} orbitals of B atoms are dominant at the Fermi level, which may account for the quasi-planar structures at the edges (red circles in Fig. 3a and b). The transporting states consisting of p_{x} orbitals (blue lines in Fig. 3a and b) and aligned along the momentum (Γ–X) direction are quantized with an energy gap in the band dispersion. The isosurface plots of electronic densities of the valence band clearly reflect this quantized feature of p_{x} states with regular nodal planes across the width direction of ribbons. The formed quasi-1D p_{x} states are highly delocalized, which may facilitate the stabilization of the LE edge-terminated BNRs. In contrast, pristine ZZ BNRs undergo severe structural distortions and become disordered after structural relaxation (not shown).

Interestingly, the structural integrity is largely maintained for both H-functionalized LE (Fig. 3c) and ZZ (Fig. 3d) BNRs, which are created through cutting the 2D H-functionalized borophene, as shown in Fig. 1b. The plots of the band structure show an odd–even oscillation of the band gap with the ribbon width for both H-functionalized LE and ZZ BNRs: the band gap is absent for odd H-functionalized BNRs but present for even H-functionalized BNRs. For instance, the band gap is 0.93 eV for LE-8 H-BNR while zero for LE-9 H-BNR. The band gap is 1.07 eV for ZZ-16 H-BNR but zero for ZZ-15 H-BNR. Therefore, depending on the width of the H-functionalized BNRs, the electronic states can become partially localized and have a band gap opening, as supported by the isosurface plots of the partial charge density of the valence bands in Fig. 3. The orbital-resolved band structure plots (see Fig. 3c) show that the frontier orbitals in LE H-BNRs are still p_{y} and p_{z} orbitals, similar to those in the pristine BNRs (see Fig. 3a and b). For the ZZ H-BNRs, the frontier orbitals mainly consist of p_{x} components, as shown by the blue dispersion lines. By examining the charge density distribution of the LE-8 H-BNRs (Fig. 3c), we find that such band gap opening in the specific type of H-BNR is due to hydrogen-induced Peierls instability of the metallic states in the BNRs. This scenario is evidenced by the strong structural distortion of the BNR lattice and the tilted H–B bonds.

We also investigate the F-functionalized BNRs and find that their electronic properties are insensitive to the ribbon width. Therefore, only one ribbon is selected as a representative for each of LE and ZZ BNRs. Fig. 4a and b show the optimized atomic structure and the band structure of F-functionalized LE-10 and ZZ-16 BNRs, respectively. The orbital-resolved band structures show that F-BNRs remain metallic, the same as pristine BNRs. Different from the graphene case, in which fluorination or ribbon construction can effectively open the band gap, the coexistence of fluorination and ribbon construction is not able to cause a metal–semiconductor transition, implying the robustness of metallicity in borophene. For the F-functionalized ZZ BNR, the in-plane p_{y} orbital of B atoms is dominant at the Fermi level, which is different from the H-functionalized ZZ BNR with p_{x} as the frontier orbital (see Fig. 3d), suggesting that these functional groups can have a selective hybridization of the boron orbitals and alter the orbital population at the Fermi level.

To examine the stability of the considered structures, we calculate the average binding energy (E_{b}) of perfect borophene which is defined as E_{b} = (n_{B}E_{B} − E_{tot})/n, and the E_{b} of H-, F,- and O-functionalized borophene using E_{b} = (n_{B}E_{B} + n_{X}E_{X} − E_{tot})/n, where E_{tot} is the total energy of the functionalized system, E_{B} is the energy of a single boron atom, E_{X} is the energy of a single H or F atom, n_{B} is the total number of B atoms, n_{x} is the total number of H or F atoms, and n is the total number of atoms in the system. The calculated values for the E_{b} for perfect and H- and F-functionalized borophene are 5.86, 4.78, and 5.25 eV, respectively. Clearly, all the three considered structures of the functionalized borophene show better stability than that of the pristine borophene. To verify the stability of BNRs, E_{b} is also calculated and is shown in Table 1. It is seen that the values of E_{b} are positive for both LE-BNRs and ZZ-BNRs, indicating that these nanoribbon derivatives are energetically stable, which is in a good agreement with the recent work.^{43} To confirm the calculation results of E_{b}, we also consider the edge energy (E_{edge}), which is defined as E_{edge} = (E_{BNR} − nE_{B})/L, where E_{BNR} is the total energy of the BNR, L is the length of the ribbon along the periodic direction, E_{B} is the total energy per atom in 2D borophene, and n is the total number of atoms. As shown in Table 1, the calculated values of E_{edge} have the following sequence: 0 < LE BNRs < LE H-BNRs < LE F-BNRs and 0 < ZZ F-BNRs < ZZ H-BNRs, confirming the stability of BNRs.

BNRs type |
E
_{b}, eV |
E
_{edge}, eV Å^{−1} |
---|---|---|

LE-9 BNRs | 5.89 | 0.062 |

LE-10 BNRs | 5.90 | 0.066 |

LE-8 H-BNRs | 4.85 | 2.170 |

LE-9 H-BNRs | 4.84 | 2.090 |

ZZ-15 H-BNRs | 4.80 | 4.330 |

ZZ-16 H-BNRs | 4.80 | 4.350 |

LE-9 F-BNRs | 5.49 | 2.600 |

LE-10 F-BNRs | 5.51 | 2.670 |

ZZ-15 F-BNRs | 5.52 | 4.040 |

ZZ-16 F-BNRs | 5.52 | 4.190 |

In addition, we also performed ab initio molecular dynamics (AIMD) calculations at 300 K for 8 ps using the Nose–Hoover method to check the stability of the functionalized and vacancy-containing borophene, as well as the pristine and functionalized BNRs. The snapshots of the simulation results are shown in Fig. 5. It is seen that during this long time (in terms of ab initio calculations), all the considered structures are stable. It should be noted that the ZZ H- and F-BNRs (Fig. 5g and i) exhibit a lower stability than the LE H- and F-BNRs.

Currently, it is still a challenge to obtain free-standing borophene experimentally, and thus double-side functionalized borophene. Nevertheless, various attempts are being made to tackle this issue. One possible way to address this challenge is to grow a borophene layer on a substrate with a weaker interaction, and then the attachment and the transfer of the borophene layer may be achieved according to the previous theoretical prediction.^{18} The other way is via surface decoration of supported borophene by creating new bonding states at the surface and edge boron atoms. In this case, dopant atoms may diffuse from the edges to the interior of the borophene–substrate interface, allowing for possible exfoliation and double-side functionalization of borophene.

In this work, we identify possible avenues for the band gap opening and the charge localization through creating vacancies, forming edges, and surface functionalization. These methods are found to be effective in the band gap adjustment in other 2D materials. However, according to our present simulation, there is no band gap opening from vacancies and chemical functionalization of the pristine 2D borophene, signifying its robust metallicity.

To further confirm the absence of a localized state with the surface functionalization and local vacancies, we calculated the electronic localization functions (ELFs) of various borophene derivatives and the results are shown in Fig. 6. The value of the ELF (between 0 and 1) reflects the degree of charge localization in the real space, where 0 represents a free electronic state while 1 represents a perfect localization. The isosurface value of 0.65 is adopted in Fig. 6. In the case of vacancies in borophene (Fig. 6a–c), it can be seen that a significant number of electrons are distributed between the B–B bonds along the [001] direction both around and far from the vacancy, indicating a strong covalent chemical hybridization. There exists a zero ELF along the [110] and [110] directions, indicating highly delocalized electronic states. This is different from other semiconducting 2D materials, like MoS_{2}, where defective states are only localized at the defect cores. A similar situation occurs for the cases of the H- and F-dopants and H-BRNs (Fig. 6d–g), where the ELF is zero at the boron atoms, implying that the electrons around the boron atoms are highly delocalized.

Fig. 6 The ELFs for (a) MV-, (b) horizontal DV-, and (c) tilted DV-containing phosphorene sheets, (d) H- and (e) F-functionalized borophene, and (f) LE-10 and (g) ZZ-16 H-BNRs. |

- A. J. Mannix, X. F. Zhou, B. Kiraly, J. D. Wood, D. Alducin, B. D. Myers, X. Liu, B. L. Fisher, U. Santiago, J. R. Guest, M. J. Yacaman, A. Ponce, A. R. Oganov, M. C. Hersam and N. P. Guisinger, Science, 2015, 350, 1513–1516 CrossRef CAS PubMed .
- B. Feng, J. Zhang, Q. Zhong, W. Li, S. Li, H. Li, P. Cheng, S. Meng, L. Chen and K. Wu, Nat. Chem., 2016, 8, 563–568 CrossRef CAS PubMed .
- Z. Zhang, A. J. Mannix, Z. Hu, B. Kiraly, N. P. Guisinger, M. C. Hersam and B. I. Yakobson, Nano Lett., 2016, 16, 6622–6627 CrossRef CAS PubMed .
- D. V. P. Massote, L. Liang, N. Kharche and V. Meunier, Phys. Rev. B, 2016, 94, 195416 CrossRef .
- H. S. Tsai, C. H. Hsiao, Y. P. Lin, C. V. Chen, H. Ouyang and J. H. Liang, Small, 2016, 12, 5251–5255 CrossRef CAS PubMed .
- Z. Zhang, Y. Yang, E. S. Penev and B. I. Yakobson, Adv. Funct. Mater., 2017, 27, 1605059 CrossRef .
- S. Ma, K. Bao, Q. Tao, C. Xu, X. Feng, P. Zhu and T. Cui, Inorg. Chem., 2016, 55, 11140–11146 CrossRef CAS PubMed .
- Z. A. Piazza, H. S. Hu, W. L. Li, Y. F. Zhao, J. Li and L. S. Wang, Nat. Commun., 2014, 5, 3113 Search PubMed .
- S. Jalife, L. Liu, S. Pan, J. L. Cabellos, E. Osorio, C. Lu, T. Heine, K. J. Donald and G. Merino, Nanoscale, 2016, 8, 17639 RSC .
- B. Feng, J. Zhang, R. Y. Liu, T. Limori, C. Lian, H. Li, L. Chen, K. Wu, S. Meng, F. Komori and I. Matsuda, Phys. Rev. B, 2016, 94, 041408 CrossRef .
- Y. Liu, Y. J. Dong, Z. Tang, X. F. Wang, L. Wang, T. Hou, H. Lin and Y. Li, J. Mater. Chem. C, 2016, 4, 6380–6385 RSC .
- H. Wang, Q. Li, Y. Gao, F. Miao, X. F. Zhou and X. G. Wan, New J. Phys., 2016, 18, 073016 CrossRef .
- E. S. Penev, S. Bhowmick, A. Sadrzadeh and B. I. Yakobson, Nano Lett., 2012, 12, 2441–2445 CrossRef CAS PubMed .
- E. S. Penev, A. Kutana and B. I. Yakobson, Nano Lett., 2016, 16, 2522–2526 CrossRef CAS PubMed .
- Y. Yang, Z. Zhang, E. S. Penev and B. I. Yakobson, Nanoscale, 2017, 9, 1805–1810 RSC .
- Y. Jiao, F. Ma, J. Bell, A. Bilic and A. Du, Angew. Chem., Int. Ed., 2016, 128, 10448–10451 CrossRef .
- L. Xu, A. Du and L. Kou, Phys. Chem. Chem. Phys., 2016, 18(39), 27284–27289 RSC .
- Y. Liu, E. S. Penev and B. I. Yakobson, Angew. Chem., Int. Ed., 2013, 52, 3156–3159 CrossRef CAS PubMed .
- X. Sun, X. Liu, J. Yin, J. Yu, Y. Li, Y. Hang, X. Zhou, M. Yu, J. Li, G. Tai and W. Guo, Adv. Funct. Mater., 2017, 27, 1603300 CrossRef .
- Z. Wang, T. Y. Lu, H. Q. Wang, Y. P. Feng and J. C. Zheng, Phys. Chem. Chem. Phys., 2016, 18, 31424–31430 RSC .
- H. Gao, L. Wang, J. Zhao, F. Ding and J. Lu, J. Phys. Chem. C, 2011, 115, 3236–3242 CAS .
- X. Niu, H. Shu, Y. Li and J. Wang, J. Phys. Chem. Lett., 2017, 8, 161–166 CrossRef CAS PubMed .
- A. A. Kistanov, Y. Cai, K. Zhou, S. V. Dmitriev and Y.-W. Zhang, J. Phys. Chem. C, 2016, 120(12), 6876–6884 CAS .
- L. Liu, J. Zhang, H. Gao, L. Wang, X. Jiang and J. Zhao, Comput. Mater. Sci., 2016, 112, 527–546 CrossRef CAS .
- Y. Du, J. Zhuang, H. Liu, X. Xu, S. Eilers, K. Wu, P. Cheng, J. Zhao, X. Pi, K. W. See, G. Peleckis, X. Wang and S. X. Dou, ACS Nano, 2014, 8, 10019–10025 CrossRef CAS PubMed .
- A. R. Oganov, J. Chen, C. Gatti, Y. Ma, C. W. Glass, Z. Liu, T. Yu, O. O. Kurakevych and V. L. Solozhenko, Nature, 2009, 457, 863 CrossRef CAS PubMed .
- G. Kresse and J. Furthmüller, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 11169 CrossRef CAS .
- J. Heyd, G. E. Scuseria and M. Ernzerhof, J. Chem. Phys., 2006, 124, 219906 CrossRef .
- A. Lherbier, A. R. Botello-Méndez and J. C. Charlier, 2D Mater., 2016, 3, 045006 CrossRef .
- H. Shu, F. Li, P. Liang and X. Chen, Nanoscale, 2016, 8, 16284 RSC .
- R. Zhang and J. Yang, 2016, arXiv:1607.08290v1.
- Y. Cai, A. Zhang, Y. P. Feng, C. Zhang, H. F. Teoh and G. W. Ho, J. Chem. Phys., 2009, 131, 224701 CrossRef PubMed .
- Y. Cai, G. Zhang and Y.-W. Zhang, Sci. Rep., 2014, 4, 6677 CrossRef CAS PubMed .
- E. J. H. Lee, K. Balasubramanian, R. T. Weitz, M. Burghard and K. Kern, Nat. Nanotechnol., 2008, 3, 486–490 CrossRef CAS PubMed .
- J. A. Robinson, M. LaBella, M. Zhu, M. Hollander, R. Kasarda, Z. Hughes, K. Trumbull, R. Cavalero and D. Snyder, Appl. Phys. Lett., 2011, 98, 053103 CrossRef .
- Y. Cai, Q. Ke, G. Zhang, B. I. Yakobson and Y.-W. Zhang, J. Am. Chem. Soc., 2016, 138, 10199–10206 CrossRef CAS PubMed .
- Q. Zhou, Q. Chen, Y. Tong and J. Wang, Angew. Chem., Int. Ed., 2016, 128, 11609–11613 CrossRef .
- H. Qiu, T. Xu, Z. Wang, W. Ren, H. Nan, Z. Ni, Q. Chen, S. Yuan, F. Miao, F. Song, G. Long, Y. Shi, L. Sun, J. Wang and X. Wang, Nat. Commun., 2013, 4, 2642 Search PubMed .
- Z. Yu, Y. Pan, Y. Shen, Z. Wang, Z. Y. Ong, T. Xu, R. Xin, L. Pan, B. Wang, L. Sun, J. Wang, G. Zhang, Y.-W. Zhang, Y. Shi and X. Wang, Nat. Commun., 2014, 5, 5290 CrossRef CAS PubMed .
- Y. Cai, H. Zhou, G. Zhang and Y.-W. Zhang, Chem. Mater., 2016, 28(23), 8611–8621 CrossRef CAS .
- A. A. Kistanov, Y. Cai, K. Zhou, S. V. Dmitriev and Y.-W. Zhang, 2D Mater., 2017, 4, 015010 CrossRef .
- W. Wu, W. Guo and X. C. Zeng, Nanoscale, 2013, 5, 9264–9276 RSC .
- F. Meng, X. Chen, S. Sun and J. He, Physica E, 2017, 91, 106–112 CrossRef CAS .

This journal is © The Royal Society of Chemistry 2018 |