Abstract

The manuscript substantiates the structural and electronic properties of covalent triazine frameworks (CTF)/As van der Waals heterostructures (vdWh) employing the standard first-principles calculation method. The numerical results designate that the CTF/As vdWh has robust crystal structures, a type-II band alignment (BA), and an indirect bandgap of 1.44 eV. The calculated results demonstrate that the strain could lead to interesting indirect-direct semiconductor transitions, while the external electric field could give rise to type-II to type-I BA and semiconductor-metal transitions. The underlined outcomes present the workability of CTF/As vdWhs in unprecedented high-performance optoelectronic equipment.

1. Introduction

Since the successful isolation of black phosphorous [1], the group-V elements of two-dimensional (2D) substances have been largely researched because they have significant physical properties [25]. Monolayer arsenic (As), a standard group-V 2D substance, has been synthesized recently [6]. Recent research suggests that monolayer gray As has a buckled honeycomb crystal structure and is an indirect band-gap semiconductor, contingent on both indirect-direct and semiconducting-metallic transitions under the external electric field or strain [6, 7]. On the other side, organic 2D covalent triazine frameworks (CTF), featured by the crystalline extended organic configurations including covalently bonded building entities, have also attracted recent eloquent consideration [810]. 2D CTF exhibits semiconductor features built by trimerization reaction of carbonitriles and utilized triazine ring (C3N3H3) as the building units. Aside from photocatalysis, 2D CTF has been employed in electronic and optoelectronic properties due to their adjustable band-gap and the correlated structure between their geometric structure and electronic properties [11, 12]. Noted that 2D As and CTF share the same hexagonal crystal structures, and the lattice parameters of CTF are approximately equal to two times those of As, making it applicable to fit CTF/As heterostructures or superlattices in the experiment.

Recently, van der Waals heterostructures (vdWh), where one layered low-dimensional configuration is stacked on the other with the aid of vdW forces, have been built experimentally and theoretically and utilized to establish electronic and optoelectronic devices with novel physical attributes [1321]. For example, theoretical research by Wei et al. [22] suggested that the external strain could constantly tune the electronic structures of the GaS/GaSe vdWh including the indirect band-gap. Our group demonstrated that the germanane/antimonene vdWh endures semiconducting-metallic transitions under external electric fields and transitions from type-II to type-I band alignment (BA) under the in-plane biaxial strain [23]. We also found that the transitions from type-III to type-II and type-II to type-I BA could be conformed to carbon nanotubes/graphene nanoribbons one-dimensional vdWh under axial strains or external electric fields [24]. In particular, vdWh utilizing 2D As or CTF has been also explored. For instance, a theoretical investigation by Su et al. [25] demonstrated that the indirect band-gap of the MoS2/As vdWh can be regulated by varying the interlayer distance. Another theoretical work about As/graphene vdWh demonstrated that altering the interlayer distance can tune the Schottky barriers and contacts of vdWhs [26]. Our group studied the electronic properties of the CTF/GaS vdWh and found that external electric fields and strain can alter the band configurations [27]. All of the above conclusions imply that stacking two 2D materials into vdWhs offers an accessible method to build substitute artificial materials that have enjoyable properties. Therefore, purporting a 2D CTF/As vdWh has been so engaging to see whether they can form a stable 2D CTF/As vdWh and validate their eye-catching electronic properties. The research implements the first-principle calculations to evaluate the structural and electronic properties of the 2D CTF/As vdWhs. The electronic properties and the band-gap of the 2D CTF/As vdWh are adjusted by employing the in-plane biaxial strain and external electric field.

2. Strategies Used for Computations

The optimum crystal structure and electronic properties of CTF/As vdWhs are conducted by using the first-principle calculation through the theory of density functional (DFT) as a part of the package available in the ATK simulation [28]. The generalized gradient approximation (GGA) of the Perdew–Burke–Ernzerhof (PBE) functional is implemented to ascertain the relaxation of the crystal configuration [29, 30]. Afterward, the Kohn–Sham equations are resolved by utilizing the algorithm called the PseudoDojo [31]. 100 Hartree is picked as the value of cutoff energy, and a 15 × 15 × 1 Brillouin area sampling K-mesh is selected. More precise outcomes are reached by utilizing a 21 × 21 × 1 K-mesh to compute the structure of the energy band. The relaxation is obtained by utilizing a 10−4 eV energy measurement and the force of the residual is less than −0.01 eV/Å. An approximately 10 Å thick vacuum layer is selected to evade relations with neighboring layers. vdW relations in a broad range need to be implemented to protect the vdWh, and its relational impacts are denoted by utilizing vdW-DF2 functional [32].

The thermal firmness of the CTF/As vdWhs is validated by running simulations called the AIMD [33], which is the ab initio molecular dynamics and is succeeded by employing the nose-thermostat methodology at a temperature of 300 K for 8 picoseconds (ps) with an increment of 1 femtosecond (fs).

3. Theoretical Results

3.1. The Geometric Structure and Stability of CTF/As vdWhs

To comprehend the electronic properties of the CTF/As vdWhs, the geometrical structure configurations of both monolayer CTF and As are needed to be optimized first. The structure of a monolayer As, which contains a lattice in a hexagonal structure and zigzagging alternative configurations based on top and side views, is depicted in Figure 1(a), respectively. The schematic configuration of the monolayer CTF, which has a hexagonal lattice configuration from a top view is described in Figure 1(d). The parameters of the optimal lattice were 3.629 Å for the monolayer As, agreeable with theoretically attained scores of 3.61 reported by Li et al. [34]. The parameters of the optimal lattice were 7.272 Å for the monolayer CTF, agreeable with theoretically and experimentally attained scores of 7.25 and 7.3 Å reported by Jiang et al. [35] and Katekomol et al. [8], respectively. Therefore, we built the CTF/As vdW heterostructure using 2 × 2 × 1 supercell for As and one cell for CTF. The lattice mismatch was defined as follows:where a1 and a2 represent the lattice constants of CTF and As, respectively. The calculated result was only ∼0.2%, indicating that the CTF/As vdW heterostructure could be easily fabricated with very little interfacial stress. The band structures of the monolayers As and CTF are described in Figures 1(b) and 1(e), respectively. For monolayer As, the conduction band minimum (CBM) is put on the Γ-M line near M points, and the valence band maximum (VBM) is put at Γ point, denoting an indirect semiconductor including a 1.69 eV band-gap, analogous to the similar theoretically obtained score by employing the PBE approach [36] (1.6 eV). For monolayer CTF, both the CBM and VBM of CTF were situated at the K point of the 2D hexagonal Brillouin area, showing that the substitute CTF is a direct semiconductor including a 2.59 eV band-gap, analogous to the similar theoretically obtained score by employing the PBE approach [35]. Figure 1(c) depicts the iso-surfaces of charge density indicating the As’s CBM and VBM states that are produced through the As-4p orbitals, while those indicated in Figure 1(f) present the CTF’s CBM and VBM stages that are led by the C/N-2p and N-2s, N/C-2p orbitals, respectively.

To find the most stable structure, cases A through F, the six potential stacking orders, are employed, as depicted in Figure 2 inset. The energy comparison indicates that case A is the most stable layout for the heterostructure. To further examine the energy stability, the formation energy (Ef) was calculated by using the following formula [3739]:where Ef represents the formation energy of the CTF/As vdW heterostructures, and ECTF/As, ECTF, EAs, and n represent the total energy of the CTF/As vdW heterostructures, CTF, As monolayer and the number of atoms in the heterostructure (n = 23), respectively. The calculated values of Ef are −138.55, −138.31, −137.67, −137.33, −137.25, −137.03 meV for case A to case F, respectively. Apparently, all the calculated Ef values were negative, indicating that all the six stacking configurations were stable from an energy perspective. Among them, the case A was the ground structure because it had the lowest formation energy, further demonstrating the ground-state configuration, which was selected as a ground-state configuration of 2D CTF/As vdWh for additional calculations.

A simulation called the AIMD is conducted at 300 K for eight ps to validate the thermodynamic firmness of the heterostructure. In order to include as many atoms as possible in the CTF/As vdWhs systems to acquire the statistical effect for temperature and to obtain the interactions among atoms adequately, the supercell of 4 × 4 × 1 is built for the As monolayers and those of 2 × 2 × 1 for the CTF monolayer and CTF/As vdW heterostructure, respectively. Figure 3 presents the complete energy oscillations versus the time increment. Figure 3 depicts the oscillation versus time increment, which is low based on complete energies (∼0.039 eV/atom) for the approaches used. Also, no distortion in the crystal configuration is compiled because of the CTF/As view after eight ps (see inset of Figure 3). The outcomes depict the thermal firmness of the CTF/As vdWh under real environmental conditions and at room temperature.

Figure 4(a) shows the projected band structures of the CTF/As vdWh. The heterostructure delineates the type-II BA in which the electronic states of CTF and As produce the CBM and VBM states, respectively. Also, the heterostructure contains a 1.44 eV indirect band-gap with the CBM and VBM condensed at K and Γ points. The iso-surfaces of charge density reflecting the C/N-2p and As-4p orbitals that generate the CBM and VBM states, respectively, are shown in Figure 4(b), endorsing the heterostructure with the type-II BA. To obtain a more detailed delineation of the electronic properties of the CTF/As vdWh, the plane-mean charge density disparity, along with the z-direction normal to the heterostructure, is depicted in Figure 4(c), and is represented bywhere ρvdWh (z), ρCTF (z), and ρAs (z) denote the plane-mean charge densities of the CTF/As vdWh, the original CTF, and the original monolayer As, respectively. Also, the red and blue isosurfaces delineate the depletion and accumulation of charge, respectively. It is concluded that the redistribution of the weak charge occurs in the interface region due to the small vdW interaction between the interfaces. Figure 4(d) delineates the plane-mean electrostatic potential along with the z-direction normal to the surface. It is observed that the potential of As is larger than the CTF layer. The potential drop (ΔV) across the bilayer is 1.74 eV (see Figure 4(d)). It implies an electrostatic field across the interface. The underlined electric field could greatly have an impact on the carrier dynamics and lead to an excitation behavior of CTF/As vdWh completely different from that of the isolated CTF layer because it could decouple electrons from holes. Moreover, these bandgaps below are the similar ones underestimated by utilizing the HSE approach. The estimation of the pinpoint inclination for the band properties could be done through a benchmark PBE functional, representing their physical properties. Thus, the electronic properties of the CTF/As vdWhs are obtained by employing the PBE approach while both the external electric field and strain are employed to mitigate the cost of the computation.

3.2. The Impact of the Vertical External Electric Field

Both theoretical and experimental results obtained previously designate that the external electric field could importantly have an impact on the electronic properties of the vdWh [4042]. Thus, the electronic structures of the CTF/As vdWhs are validated while a vertical external electric field upright to the CTF/As vdWhs is exerted, while CTF to As has a positive direction. The external electric field is changed in [−1.4, 2.0] V/Å with an increment of 0.2 V/Å. Figure 5 depicts the impact of the external electric field on the electronic properties of CTF/As vdWhs. The CBM is attained from CTF and situated at K points, while the VBM is attained from As and situated at Γ points in [−1.2, 0.0] V/Å, indicating an indirect semiconductor including the type-II BA. Both the CBM and VBM derived from As and CBM are situated at the Γ-M line near the M points, and VBM is situated at Γ points in [0.2, 0.4] V/Å, indicating an indirect semiconductor including the type-I BA. The CBM derived from As was located on the Γ-M line close to the M points, while the VBM acquired from CTF was situated at K points in [0.6, 1.8] V/Å, indicating an indirect semiconductor including the type-II BA. Also, a transition from semiconductor to metal state is attained at −1.4 and 2.0 V/Å. The bandgaps and their BA are given in Figure 6(a). The underlined results suggested that their band structures altered when an external electric field is exerted.

3.3. The Effect of Parallel Strain

Besides, the strain could have an impact on the electronic properties of the vdWhs [40, 43]. Afterward, the impact of in-plane biaxial strain on the electronic structures of the CTF/As vdWhs is examined. The parameters of a crystal lattice are altered to simulate the in-plane biaxial strain on the CTF/As vdWhs computed bywhere a0 and a represent the crystal lattice parameters of vdWhs when both unstrained and strained requirements are under consideration. The [−5% to 5%] interval with a 1% increment was picked for the outside strain. The stability of CTF/As vdWhs under strain was obtained with the strain energy and was calculated using the following equation [44, 45]:where the Es, Estrained, Eunstrained, and n were the strain energy, the energy of pristine heterostructures, the energy of the strained heterostructures, and the number of atoms in the heterostructure (n = 23), respectively. As shown in Figure 6(b), the strain energies are a quadratic function relative to the strains of [−5% to 5%], indicating that the deformation of the CTF/As vdWhs is in the elasticity limit and the geometrical structures are stable. Figure 7 depicts the impacts of the strain on the electronic configurations of the CTF/As vdWhs. The CBM derived from CTF is situated at K points, while the derivation of the VBM based on As is situated at Γ points in [−5%, 1%], indicating an indirect semiconductor including the type-II BA. The CBM derived from CTF is situated at K points, while the derivation of the VBM based on As is situated on the Γ-M line near the M points and the Γ-K line near the K point for 2%, indicating an indirect semiconductor containing the type-II BA. Both the CBM and VBM acquired from As are situated at Γ points in [3%, 5%], indicating a direct semiconductor including the type-I BA. Figure 6(b) presents their bandgaps. When the strain raised in [−5%, 1%], the band-gap raised similarly, while it decreased once the strain is boosted in [2%, 5%]. The outcomes imply that the strain can efficiently adjust the electronic properties of the CTF/As vdWhs.

4. Conclusion

Consequently, the electronic properties of the CTF/As vdWhs were assessed by utilizing the first-principle calculation contingent on both external electric field and strain. Concluded that the CTF/As vdWhs were the indirect semiconductors including a type-II BA changing in [−1.2, 0.0] and [0.6, 1.8] V/Å, while the CTF/As vdWhs were the indirect semiconductors containing a type-I BA altering in [0.2, 0.4] V/Å, and semiconductor-metal transitions occurred at −1.4 and 2.0 V/Å. The CTF/As vdWhs were the indirect semiconductors including a type-II BA changing in [−5%, 2%], and the direct semiconductors including a type-I BA ranging in [3%, 5%]. The outcomes showed that the adjustable band-gap of the CTF/As vdWhs leads to a hopeful perspective for building a substitute nano-electronic and optoelectronic devices.

Data Availability

The data used to support the findings of this study are available upon reasonable request from the corresponding authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China with grant no. 11864011. This work was also supported by the Fundamental Research Program of Shanxi Province (Grant no. 202203021212018) and Foundation of Shanxi Educational Committee (Grant nos. J2021770 and 2021L574). The authors would like to thank the MJEditor (https://www.mjeditor.com) for its linguistic assistance during the preparation of this manuscript.