A new structure of two-dimensional allotropes of group V elements.

The elemental two-dimensional (2D) materials such as graphene, silicene, germanene, and black phosphorus have attracted considerable attention due to their fascinating physical properties. Structurally they possess the honeycomb or distorted honeycomb lattices, which are composed of six-atom rings. Here we find a new structure of 2D allotropes of group V elements composed of eight-atom rings, which we name as the octagonal tiling (OT) structure. First-principles calculations indicate that these allotropes are dynamically stable and are also thermally stable at temperatures up to 600 K. These allotropes are semiconductors with band gaps ranging from 0.3 to 2.0 eV, thus they are potentially useful in near- and mid-infrared optoelectronic devices. OT-Bi is also a 2D topological insulator (TI) with a band gap of 0.33 eV, which is the largest among the reported elemental 2D TIs, and this gap can be increased further by applying compressive strains.

The successful exfoliation of graphene from graphite brings two-dimensional (2D) materials from the theoretical domain to the real world123. Since then, many more 2D materials have been studied. For example, graphene, in a honeycomb structure, has attracted tremendous attention due to its excellent dynamical properties and unique electronic structures3456. Silicene78, germanene7910, stanene1112 and Bi(111) bilayer13141516171819 are analogues of graphene and they have the buckled honeycomb structures. Recent studies7111314 have proposed that these analogues are 2D topological insulators (TIs), also called quantum spin Hall (QSH) insulators, which are characterized by coexistence of gapped bulk states and gapless edge states protected by time reversal symmetry2021. They have attracted much attention due to fundamental interest and potential applications in spintronics and topological quantum computations2021. The nontrivial energy gap of Bi(111) bilayer is about 0.2 eV1415, which is the largest among the reported elemental 2D TIs until now. This is important, because for room temperature utilizations, energy gaps of TIs must be large enough to suppress the effects of thermal perturbations on transport properties. Encouraging signatures have been experimentally observed in this material, such as the existence of edge states1819. Black phosphorus (BP) layers have been investigated for their potentials for optoelectronics and electronics applications. These layers have the puckered honeycomb structure and can be exfoliated from the bulk materials22. They exhibit high-mobility transport anisotropy and linear dichroism2223, and the energy gaps can be tuned from about 0.3 to 1.5 eV by varying the layer thickness24. Recently, BP field-effect transistors have been experimentally fabricated25.

The 2D elemental materials discussed above all have the honeycomb or distorted honeycomb structures, which are composed of six-atom rings. Carbon has been reported to have many other 2D allotropes2627, while the known 2D allotropes for other elements are quite limited. A pseudocubic {012}-oriented allotrope of Bi can be realized on Si(111) surface28, and it was proposed to be in the BP-like puckered-layer structure. In three dimensions, the group V elements exhibit several allotropic forms because their semi-metallic bonding character can be easily shifted to either metallic or covalent side. Different allotropes can have very different physical properties. For example, As and Sb can exist in metallic forms or as semiconductors of different band gaps, depending on their crystal structures29. For two dimensions, it is also of interests to explore new allotropes of the group V elements and investigate their electronic properties.

In this article, we report a new structure of 2D allotropes of the group V elements (P, As, Sb and Bi) composed of eight-atom rings. These allotropes contain unique atomic octagonal tiling (OT) patterns and are shown to be stable by using first-principles calculations. All these OT ultra-thin films are semiconductors with band gaps ranging from 0.3 to 2.0 eV, covering the entire gap range of BP and transition-metal dichalcogenides (TMDCs), thus they are candidates for various near- and mid-infrared optoelectronic applications. In addition, the OT-Bi bilayer is a 2D TI with a band gap of 0.33 eV, which is larger than the gap of the Bi(111) bilayer and is the largest among reported elemental 2D TIs.

Result

The 2D lattice structure model studied in this work is shown in Fig. 1(a), and its unit cell contains eight atoms. The atoms have two different heights, so it can be viewed as a bilayer structure. From top view, the atoms form a beautiful octagonal tiling, which is composed of octagons with alternating long and short edges and smaller squares. In the following, we will show that all the OT-X (X = P, As, Sb and Bi) bilayers are metastable 2D allotropes of group V elements and they exhibit interesting structural and electronic properties.

Figure 1

(a) Atomic structures of the 2D octagonal tiling structure. The two colours indicate atoms in different heights of the bilayer. (b) The calculated phonon spectrum of the OT-Bi bilayer.

Structural optimizations of the octagonal tiling structure were performed for the group V elements P, As, Sb and Bi. The optimized structural parameters such as the lattice constant a, the bilayer thickness d, and the bond lengths l1 and l2 of the long and short edges of the octagons are listed in Table 1. When viewed in terms of the four-atom squares, l1 is the “intra-square” bond length, and l2 is the “inter-square” bond length. These parameters increase as the atomic size increases. Comparing with the reference structures such as the monolayer black P and As, bilayer buckled honeycomb Sb and Bi, the energies of these OT allotropes are slightly higher. Typically, the energy differences are about 70 to 80 meV per atom.

Table 1

The calculated lattice constant a, the bilayer thickness d, and the bond lengths l1 and l2 of the long and short edges of the octagons of the OT allotropes (a, d, l1 and l2 all in units of angstrom).

	a	d	l1	l2
P	6.43	1.27	2.28	2.24
As	7.08	1.42	2.53	2.49
Sb	8.05	1.67	2.91	2.87
Bi	8.47	1.76	3.05	3.02

To check the dynamical stability of these 2D materials, we calculate their phonon spectrums. These OT allotropes have similar phonon dispersions. The phonon spectrum of the OT-Bi bilayer is shown in Fig. 1(b) as an example. Generally speaking, the phonon bands can be separated into four bunches, and quantitatively, lighter materials have higher phonon frequencies. There are no imaginary phonon modes in all cases, indicating that these 2D allotropes are dynamically stable. Additionally, first-principles molecular dynamics simulations were performed for over 9 ps time, corresponding to 6000 simulation steps. All these allotropes maintain their structures at temperatures up to 600 K, showing good thermal stability.

Now we turn to the electronic structures of these 2D allotropes. The calculated energy band structures without and with spin-orbit coupling (SOC) are shown in Fig. 2 for all the 2D OT allotropes. The contributions of the p and p orbitals to the wavefunctions are shown in different colors, and the energy gaps are labeled.

Figure 2

Calculated band structures of OT-P (a), OT-As (b), OT-Sb (c) and OT-Bi (d) without (top panels) and with (bottom panels) SOC. The corresponding global energy gap E is also shown. Red and green symbols denote projections of the p and p states, respectively, and symbol size indicates the contribution weight.

As shown in Fig. 2, all the OT bilayers are insulators. OT-P is an indirect gap insulator. Both its conduction band minimum (CBM) and the valence band maximum (VBM) are located away from Γ point and they are mainly composed of p states. At Γ point, the VBM has contributions from both p and p states, while the CBM mainly comes from p states. OT-As and OT-Sb allotropes are direct gap insulators. Both their VBM and CBM are located at Γ point. The former mainly comes from p states, while the latter mainly comes from p states. Besides, from P to Bi, the energy gap decreases monotonically. This is in agreement with the increasing metallic behaviours from P to Bi. The gaps of OT-P, OT-As and OT-Sb range from about 1 to 2 eV, almost covering the full gap range of TMDCs such as MoS2, MoSe2, WS2 and WSe230, thus they may be alternative candidates for various near-infrared optoelectronic applications. We note that the energy gaps of OT-As and OT-Sb with SOC are smaller than those without SOC. This can be understood as the following. Spin-orbit coupling has the effect to invert the frontier orbits. If the non-SOC energy gap is not too large and the SOC strength is large enough, band inversion usually occurs, which is the typical mechanism of topological insulators2021. The stronger the SOC, the larger the gap opens. Here, As and Sb atoms have weak SOC, and their non-SOC energy gaps are large (all larger than 1.3 eV). Therefore, their SOC is not enough to invert the bands, but only narrowing their gaps. However, the minimum gap (2.11 eV) of OT-P with SOC is slightly larger than the non-SOC gap (1.97 eV), as shown in Fig. 3(a). This is because of the indirect gap nature in OT-P. Considering the direct gap at Γ point, where band inversion usually occurs, the SOC energy gap is 2.38 eV, smaller than the non-SOC gap, which is 2.53 eV.

Figure 3

Electronic structures of the OT-Bi nanoribbon: (a) Energy band dispersion, (b) projected density of states (PDOS) and (c) real space projection of the Dirac point wave functions. Purple symbols in (a) indicate contributions of edge atoms, and symbol size denotes the weight contribution.

When not considering SOC, OT-Bi has similar band structures as OT-P and OT-Sb. It has a direct gap, and the p and p states contribute to the VBM and CBM, respectively, as shown in Fig. 2(d). The non-SOC band gap is 0.59 eV, the smallest among the four OT allotropes due to its strongest metallic properties. When SOC is switched on, things become different. The strong SOC of the heavy Bi atoms inverts the VBM and CBM at Γ, and now the former mainly comes from the p states while the latter mainly from p states. The VBM bands around Γ resemble an M shape, giving the OT-Bi an indirect energy gap of 0.33 eV, very close to the bulk gap of black P (about 0.3 eV)24. Generally, band inversion is a strong indicator of topologically nontrivial band structures, and we will confirm this in the following.

We now discuss the topological properties of the band structures in these 2D allotropes. To do this, the topological invariant, i.e., the Z2 number is calculated using the parity method31. With this method, the Z2 invariant can be determined through computing the parities of wave functions at time-reversal invariant momenta (TRIMs) Γ in the Brillouin zone as the following31,

Here, ξ = ±1 is the parity eigenvalues of wave functions at four TRIMs, N the number of occupied bands, and ν (mod 2) the Z2 number. If Z2 = 1, the system is nontrivial, otherwise it is trivial.

The calculated Z2 numbers for the four OT allotropes are listed in Table 2. Obviously, Z2 = 1 for the OT-Bi allotrope, thus it is a 2D TI, while Z2 = 0 for OT-P, OT-As and OT-Sb, indicating that they are trivial insulators. This is in agreement with the band structures shown in Fig. 2, where band inversion only occurs in OT-Bi. The following discussions on OT-Bi are all based on SOC band structure calculations.

Table 2

Calculated parities (δ) and Z2 invariants of the four OT allotropes.

	P	As	Sb	Bi
δ(Γ)	+	+	+	−
3δ(M)	+	+	+	+
ν	0	0	0	1

As shown in the bottom panel of Fig. 2(d), OT-Bi has an indirect gap, which usually occurs in most reported 2D TIs111532. The indirect gap here is 0.33 eV (direct gap at Γ is 0.40 eV), which is large enough for room temperature observations and utilizations. As far as we know, this gap is the largest among the reported elemental 2D TIs. For elemental TIs without strains or chemical decorations, the largest gap before the present work occurs in the buckled honeycomb Bi(111) bilayer, which is about 0.2 eV1415. Gaps of other 2D TIs range from a few meV to 0.1 eV, for examples, 1.55 meV in silicene, 23.9 meV in germanene7, 0.1 eV in tin films11, and 0.1 eV in ZrTe5 and HfTe5 single-layers32. The rather large gap here is due to the strong spin-orbit coupling of the heavy Bi atom. Besides, the interesting octagonal tiling structure here may provide proper Bi-Bi distances suitable for the strong SOC.

2D TIs are bulk insulating, with Dirac-type gapless states at the edges2021. To confirm this we construct a nanoribbon from the OT-Bi bilayer and calculate its electronic structures. The band structures of the ribbon are shown in Fig. 3(a), where gapless edge states clearly appear. The Fermi velocity at the Dirac point is 6.18 × 105 m/s, which is larger than the 4.4 × 105 m/s in stanene but it is similar to the 6.8 × 105 m/s in fluorinated stanene11. Figure 3(b) shows the projected density of states (PDOS) on both the edge and the bulk atoms. As shown in the figure, edge atoms contribute almost all the DOS around the Fermi level, indicating that the edge is metallic. On the other hand, bulk atoms have negligible contributions near the Fermi level, consistent with the insulating bulk states. Real space projection of wave functions of Dirac states at Γ point (Dirac point) in Fig. 3(c) shows that the edge states here have a large penetration depth of about 3 nm. This is larger than that in the buckled honeycomb Bi(111) bilayer, where the depth is less than 2 nm18. Further analysis indicates that spins and momenta at the Dirac states are always perpendicularly locked. These features of the ultrathin OT-Bi film are consistent with their topologically nontrivial properties2021.

Nontrivial TIs are robust, i.e., they are able to keep their nontrivial properties against nonmagnetic perturbations, such as strain. Strain is also a widely-used method to improve or manipulate the band gaps of 2D TIs71132. In this work, we apply the following external in-plane biaxial strain

where a0 is the equilibrium lattice constant, and a is the strained lattice constant. So ε > 0 (<0) means tensile (compressive) strains. Figure 4(a) shows the band gap of the OT-Bi bilayer as the function of the applied strain ε. The energy gap remains open for strains between ε = −9% and 6%, indicating that no topological phase transition occurs within this range.

Figure 4

Evolutions of (a) energy gaps and (b–f) band structures of the OT-Bi ultra-thin film under different external strains.

As shown in Fig. 4(a), the gap value is very sensitive to strain. We can understand this by analyzing the evolution of the four frontier bands, C1, C2 and V1, V2 as indicated in Fig. 4(b–f). With compressive strains, V1 and V2 bands change very slowly, C1 and C2 bands shift significantly. The C1 band shifts upward, while the C2 band shifts downward, so the energy gaps increase at the initial stage of compression. At about ε = −5%, the C1 and C2 bands touch each other at Γ point, where the direct and indirect gaps reach their maximum values of 0.53 and 0.48 eV, respectively. Under further compressive strains, the C2 band becomes the lowest conduction band and shifts down, and the energy gaps decrease with increasing compressive strain. Under tensile strain, the C1 and C2 bands change very slowly, while the V1 and V2 bands shift upward. So the energy gaps decrease with tensile strain. The V2 band moves faster than the V1 band, and they touch at Γ point when ε = 1.5%. Under further tensile strain, the two bands exchange order around Γ, now the direct gap at Γ point is between the C1 and V2 bands. So the direct gap decreases at a faster rate with tensile strain beyond 1.5%. The indirect gap is still determined by the C1 and V1 bands, because the band exchange only occurs in a small area around the Γ point of the Brillouin zone. In this strain process, the energy gap never reaches zero, so the Z2 number of the system does not change. Strain is an often used method for manipulating the gaps of 2D TIs. For example, it has been reported that a 20% tensile strain can increase the energy gap of the single-layer ZrTe5 from 0.1 eV to more than 0.2 eV32, and a compressive strain of −6% can change the gap of silicene from 1.55 meV to 2.90 meV7.

In addition, the OT-Bi quadruple layer is also calculated to be a 2D TI with an energy gap of 0.15 eV (the direct gap is 0.37 eV). Films with even more layers become metallic. It is possible to modify their electronic properties by methods such as surface adsorptions16, which is of interest for future studies.

Discussion

Different crystal structures of the same element can possess different physical properties. For example, in contrast to the OT-P bilayer studied here, the black phosphorus bilayer in puckered honeycomb structure is a direct gap insulator23. Its VBM and CBM are both located at Γ. The OT-P bilayer has a larger band gap than the BP bilayer, which was calculated to be about 1.51 eV. For Bi, both the OT-Bi bilayer and the Bi(111) bilayer are indirect gap insulators14. However, the energy gap of the OT-Bi bilayer is also larger than the gap of the Bi(111) bilayer in distorted honeycomb structure.

These OT allotropes have the square lattices. For experimental synthesis on a substrate, candidate substrates are zincblende (ZB) CdTe, InSb, and HgTe (for OT-P), ZB-BN, BAs (for OT-As), Ge, ZB-AlAs, GaAs (for OT-Sb), and ZB-CdSe, GaSb, InAs, ZnTe (for OT-Bi). All these substrates can provide square lattices with lattice mismatch less than 2% for the corresponding 2D octagonal tiling allotropes. The experimental fabrications of these OT allotropes, the band structure engineering by methods such as surface adsorption, and properties of thicker films are all of great interests, and further theoretical and experimental explorations along these directions are needed.

In conclusion, we predict a new structure of 2D allotropes of group V elements in octagonal tiling structures. They are insulators with energy gaps ranging from about 0.3 to 2.0 eV, covering almost the whole gap range of BP and TMDCs, so they may be promising candidates for various near and mid-infrared optoelectronic applications, such as photodetectors and modulators. In particular, the octagonal tiling Bi allotrope is predicted to be a 2D topological insulator, and its energy gap is larger than any previously reported elemental 2D TIs. We hope that the present findings will motivate further experimental and theoretical efforts.

Methods

This work is performed using first-principles method within the framework of density functional theory of generalized gradient approximations (DFT-GGA)3334, as implemented in the VASP codes35. We have also performed calculations using the local-density approximation (LDA), and the optimized in-plane lattice parameters and corresponding band gaps when using LDA and GGA are very similar. The spin-orbit coupling as implemented in VASP is based on relativistic spin-density functional theory. Calculations including spin-orbit coupling are performed in the non-collinear mode. A good description of the spin-orbit coupling as implemented in VASP can be found in the paper by Blonski et al.36. Phonon spectrum is obtained by the calculation of the interatomic force constants (IFCs) in VASP, with the help of the phonopy package37. The phonon calculations are performed using both the 4 × 4 × 1 and the 6 × 6 × 1 supercells, and the obtained phonon spectra are practically identical. We have checked that the IFCs are short range, which guarantees that the Fourier interpolated phonon frequencies are highly accurate38. A vacuum of more than 20 Å is added between the layers in the supercell calculations of the 2D bulks or 1D ribbons. Both the in-plane lattice constants and the atom positions of 2D bulks are relaxed until the force on each atom is less than 0.01 eV/Å. A 17 × 17 k-point mesh is used for all the relaxation and static calculations.

Additional Information

How to cite this article: Li, P. and Luo, W. A new structure of two-dimensional allotropes of group V elements. Sci. Rep. 6, 25423; doi: 10.1038/srep25423 (2016).