Monolayered Silicon and Germanium Monopnictide Semiconductors

Jan 29, 2018 - ACS eBooks; C&EN Global Enterprise .... These semiconductors, including the monolayered SiP, SiAs, GeP, and GeAs, possess wide .... Thu...
0 downloads 0 Views 3MB Size
Letter www.acsami.org

Cite This: ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Monolayered Silicon and Germanium Monopnictide Semiconductors: Excellent Stability, High Absorbance, and Strain Engineering of Electronic Properties Ai-Qiang Cheng,†,‡ Zi He,†,‡ Jun Zhao,§,⊥ Hui Zeng,*,‡ and Ru-Shan Chen*,‡ ‡

School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing, Jiangsu 210094, China School of Science, Nanjing University of Posts and Telecommunications, Nanjing, Jiangsu 210023, China ⊥ School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China §

S Supporting Information *

ABSTRACT: The discovery of stable two-dimensional (2D) semiconductors with exotic electronic properties is crucial to the future electronic technologies. Using the first-principles calculations, we predict the monolayered Silicon- and Germanium-monopnictides as a new class of semiconductors owning excellent dynamical and thermal stabilities, prominent anisotropy, and high possibility of experimental exfoliation. These semiconductors, including the monolayered SiP, SiAs, GeP, and GeAs, possess wide bandgaps of 2.08−2.64 eV obtained by hybrid functional calculation. Under small uniaxial strains (−2 to 3%), dramatic modulations of their band structures are observed, and furthermore, all the 2D monolayers (MLs) can be transformed between indirect and direct semiconductors. The monolayered GeAs and SiP exhibits extraordinary optical absorption in the range of visible and ultraviolet (UV) light spectra, respectively. The exfoliation energies of these monolayers are comparable to graphene, implying a strong probability of successful fabrication by mechanical exfoliation. These intriguing properties of the monolayered silicon- and germanium-monopnictides, combined with their highly stable structures, offer tremendous opportunities for electronic and optoelectronic devices working under UV−visible spectrum. KEYWORDS: monolayered monopnictide, strain engineering, optical absorbance, exfoliation energy, indirect-to-direct transition fabricated by using self-flux under high pressure and seeded flux growth.14,15 Experimental characterizations have shown that they crystallize with orthorhombic phase for SiP and monoclinic phase for the others. Moreover, the atoms in each 2D layer are covalently bonded and different layers are stacked by weak van der Waals (vdW) interactions, in analogy to the transition metal dichalcogenides.16 All these IV−V binary compound crystals are semiconductors with bandgaps of 0.52− 1.69 eV.17−20 The anisotropic structures of these monolayered binary compound AB (A = Si, Ge; B = P, As) give rise to highly anisotropic transport properties and high-temperature thermoelectric character.21,22 Previous density functional theory (DFT) calculations revealed that the monoclinic phase is energetically more favorable than the orthorhombic phase,23−25 and the electronic properties of several free-standing MLs have been investigated.25−28 However, the fundamental physical properties related to these MLs are not yet clear, and their potential applications in nanoelectronic and optoelectronic are still lacking. For instance, what about the optical adsorption

T

he rising of graphene has sparked tremendous interests in various two-dimensional (2D) materials consisting of atomically thin sheets,1−3 such as Group-IV elemental analogues,4 transition metal dichalcogenides,5 III−VI layered materials,6 and monolayer MXenes,7 which are promising semiconductor candidates for nanoelectronics and optoelectronics due to their suitable bandgaps. Among them, the monolayered MoS2 and phosphorene, possessing both appreciated fundamental bandgap and moderate intrinsic carrier mobility, have attracted considerable attentions more recently.8−10 Experimental studies have demonstrated intriguing electrical performance in the MoS2- and phosphorene-based nanoelectronic and optoelectronic devices.11 However, fewlayered phosphorene is chemically unstable, which hinder its practical applications in semiconductor device.12 Meanwhile, the presence of various defects in the fabrication of monolayered MoS2 is considered to be inevitable, which significantly deteriorate the carrier mobility and thus impede its device applications.13 In this sense, searching for stable 2D semiconductors with novel electronic properties becomes a challenging problem. More recently, bulk crystals of four IV−V layered materials, i.e., SiP, SiAs, GeP, and GeAs, have been experimentally © 2018 American Chemical Society

Received: November 17, 2017 Accepted: January 29, 2018 Published: January 29, 2018 5133

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces properties of these 2D IV−V binary monolayers (MLs), which is of great importance to their applications in the field of optoelectronic. How are the electronic properties of these MLs modulated by stain engineering ? Furthermore, is it possible to experimentally fabricate these monolayered binary compounds and what is the effective synthesis route ? Motivated by the successful fabrication of few-layered IV−V binary compound in experiments,14,20 in the following part, we perform DFT calculations to systematically study the structures, electronic and optical properties, and stabilities of the monolayered SiP, SiAs, GeP, and GeAs semiconductors. It will be shown that these semiconductors exhibit excellent dynamical and thermal stabilities, and the dramatic modulations of band structures can be obtained by small uniaxial strains. Moreover, the excellent optical absorbance ranging from visible to UV light and the feasibly experimental route to fabricate these 2D MLs from bulk crystal will also be addressed. In this work, we have carried out the first-principles calculations based on DFT within the generalized gradient approximation (GGA) to study the electronic structures. All calculations are implemented in the Vienna ab initio simulation package (VASP).29,30 For structural optimizations, the exchange-correlation potential is described by GGA with Perdew−Burke−Ernzerhof (PBE) version.31 We have adopted the DFT-D3 approach proposed by Grimme et al. to correctly describe the vdW interactions.32 The cutoff energy for planewave basis is 500 eV. For the ionic relaxation, the convergence criteria are 0.01 eV/Å for the remnant force and 1 × 10−6 eV/ atom for the energy. The Brillouin zone (BZ) is sampled by using the Monkhorst−Pack scheme,33 and we adopt the k-grid of 30 × 6 × 1 for geometry optimization and 40 × 8 × 1 for the computation of density of states. To minimize the interlayer interactions, the vacuum spacing of 25 Å was employed for the monolayers. To correct the well-known bandgap underestimation reported by GGA-PBE, we also performed band structure calculations based on screened hybrid functional in the form of Heyd−Scuseria−Ernzerhof (HSE06).34,35 The phonon dispersion spectra are calculated by using density functional perturbation theory utilized in the CASTEP code.36,37 To study the thermal stability, we have carried out ab initio molecular dynamics (AIMD) simulations in the canonical ensemble. The 6 × 1 × 1 ML supercell consisting of 144 atoms is used in simulation, and the temperature is controlled by the Nosé−Hoover algorithm at 1000 and 1500 K.38 To evaluate the optical properties of these monolayered semiconductors, we computated the frequency-dependent dielectric function ϵ(ω) as well as the absorption coefficient α(ω).39 The imaginary part ϵ2(ω) is calculated by using the formula: ϵ2(ω) =

Figure 1. Ball-and-stick model of the monolayered IV−V semiconductors. Top: top view of 2 × 1 monolayered supercell from the Zaxis, and the primitive cell is marked by black rectangle. Bottom: side view of the 2 × 1 supercell from the X-axis. Some bond lengths (A = Si, Ge; B = P, As) are shown and the corresponding values are listed in Table 1.

A(ω) =

(2)

Among the four IV−V bulk crystals, only SiP has the orthorhombic phase with space group Cmc21 (No. 36), while the other three (SiAs, GeP, GeAs) possess the monoclinic C2/ m (No. 12) symmetry. The preferred exfoliation orientation (cleavage planes) for the SiP is [0, 0, 1], whereas the other three crystals is [2, 0,−1]. The primitive cells of the monolayered binary IV−V semiconductors are rectangular, and each primitive cell consist of 12 Group-IV atoms and 12 Group-V atoms. The optimized lattice parameters of each primitive cell are listed in Table 1, which is found to be wellconsistent with the experimental characterizations.17−19,21 For these four monopnictides, every Group-V (denoted as B in Figure 1) atom is coordinated to three Group-IV (denoted as A) atoms, and each A atom is coordinated to three B atoms and one other A atom. There is negligible difference in bond length between the AA1 bond and AA2 bond. However, the bond length is slightly increased from AB1 to AB4. Additionally, the buckling distance δ is also increased as the atomic number is increasing. The thermal and dynamical stabilities of the monolayered semiconductors are evaluated in terms of total energy variation and phonon dispersion, respectively. As shown in Figure 2, the time-dependent evolutions of the total energies are oscillating within a very narrow range, suggesting that these MLs are thermally stable. Hence, the AIMD snapshots have clearly shown that the puckered monoclinic monolayers are well maintained, demonstrating that all 2D MLs are stable at 1000 K. In particular, the SiP and SiAs have superior stabilities and are found to be stable under 1500 K, which is consistent with the experimental measurement.14 The calculated phonon dispersions shown in Figure S1 have clearly demonstrated that there is no imaginary phonon mode, ensuring the dynamical stability of these MLs. Furthermore, the highest frequency obtained for the monolayered SiP is about 400 cm−1, and this value is decreased to about 300 cm−1 for both SiAs and GeP, and 220 cm−1 for the GeAs. This is attributed to the differently elemental IV−V bonds. Thus, these monolayered

4π 2e 2 1 lim Σc , v , k2Wkδ(ϵck − ϵvk − ω)⟨uck + eαq| Ω q → 0 q2 uvk⟩⟨uck + eβq|uvk⟩*

2 ω[ ϵ12(ω) + ϵ22(ω) − ϵ1(ω)]1/2

(1)

where c and v are associated with the conduction and valence band states, respectively. The real part ϵ1(ω) is calculated by using the Kramer-Kroning transformation.40 Furthermore, it is noteworthy to mention that the light incidence along the Z axis, as indicated in Figure 1, and the polarization along X/Y direction. Hence, the frequency-dependent absorption coefficient A(ω) is derived from the dielectric functions: 5134

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces Table 1. Structural Parameters of the Monolayered Binary Compound AB (A = Si, Ge; B = P, As)a

a

MLs

a (Å)

b (Å)

β (deg)

AA1 (Å)

AB1 (Å)

AB4 (Å)

δ (Å)

SiP SiAs GeP GeAs

3.519 3.682 3.648 3.803

20.449 21.246 21.393 22.129

90.00 106.56 100.18 100.38

2.351 2.346 2.479 2.479

2.271 2.394 2.360 2.471

2.300 2.422 2.395 2.504

4.732 4.856 4.982 5.099

The parameters correspond to the representation shown in Figure 1.

Figure 2. Evolution of total energies and snapshots of the 2D monolayered (a) SiP, (b) SiAs, (c) GeP, and (d) GeAs from AIMD simulations at the end of 5 ps.

monopnictide semiconductors exhibit excellent dynamical and thermal stabilities. The electronic structures of the monolayered IV−V semiconductors are shown in Figure 3. All 2D MLs possess suitable bandgaps (Eg), which decrease with the increasing of atomic number. The band structures of each 2D monolayered semiconductor calculated at HSE level are shown in Figure S2.

is generally 0.65 eV larger than the value It is found that EHSE g and EGGA lie in the visible-light reported by PBE. Both EHSE g g region, which is greater than the corresponding value for the germanane (1.53 eV) as well as phosphorene (1.65 eV) and comparable with MoS2 (2.13 eV).16,40−42 As demonstrated in Figure 4a, the valence band maxmium (VBM) is originated

Figure 3. Calculated electronic band structures from GGA-PBE for the monolayered (a) SiP, (b) SiAs, (c) GeP, and (d) GeAs. The dashed line indicates the Fermi level.

Figure 4. Partial density of states (PDOS) for the monolayered (a) SiP and (b) GeP. The Fermi level is set to zero. 5135

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces Table 2. Electronic Structure Parameters of the IV−V Binary Monopnictide MLsa MLs

type

EGGA (eV) g

EHSE (eV) g

Φ (eV)

mxe (m0)

mye (m0)

mxh (m0)

myh (m0)

SiP SiAs GeP GeAs

direct direct indirect direct

1.899 1.694 1.610 1.423

2.641 2.353 2.309 2.077

4.783 4.496 4.752 4.457

0.111 0.050 0.053 0.118

0.062 0.003 0.004 0.004

−0.351 −0.137 −0.175 −0.152

−0.050 −0.002 −0.004 −0.002

(eV) and EHSE (eV) are the bandgap calculated by GGA and HSE, respectively; Φ is the work function. me and mh are the effective masses of g electron and hole, respectively, along the x and y directions, and m0 is the rest mass of electron. a GGA Eg

Figure 5. Electronic bandgaps and total energies of the monolayered binary compound AB (A = Si, Ge; B = P, As) varied with respect to the strain along y axis. (a−d) Monolayered (a) SiP, (b) SiAs, (c) GeP, and (d) GeAs.

experimentally by using strain engineering.47 We thus studied the influence of a uniaxial strain (εy) along the y axis on the electronic structures of monolayered binary compounds, as presented in Figure 5. The uniaxial strain εy is defined by

from the hybridization of p-orbitals of A and B atoms. In contrast, the s- and p-orbital of A and B atoms almost have equal contributions to the conducting band minimum (CBM). Thus, except for the GeP, the other three MLs are direct bandgap semiconductors. It is noted that the CBM of GeP is located at 4/5 along Γ → X, in striking contrast to the SiP case. It is attributed to the fact that CBM of the former is governed by the s-orbital of the Ge atom, as confirmed by its PDOS result shown in Figure 4b. To estimate the transport properties of these MLs, we have calculated the effective mass (m*) of charge carriers by using the expression m* = ℏ2

(

∂ 2E(k) ∂k 2

εy =

where b0 and b are the unstrained and strained

lattice constants along the y direction, respectively. The staindriven direct (indirect)-to-indirect (direct) bandgap transition is observed for all 2D MLs. More importantly, such bandgap transition is achievable for experimental realization because the critical tensile strain is only 1, 3, and 1% for the monolayered SiP, SiAs, GeAs, respectively. Hence, 2% compressive strain allows indirect-to-diect bandgap transition for the monolayered GeP, which is crucial for its optoelectronic applications. It is interesting to find that the direct bandgap of SiAs could be effectively modulated from 1.3−1.9 eV by the uniaxial strain εy, which is similar to the cases in SmSe and GeSe.48,49 Hence, the bandgaps of the monolayered SiP and SiAs are linearly decreased with respect to the compressive strain. The uniaxial stain, which can be imposed to the 2D nanomaterials by different substrates in practice, allow us to tailoring the band structures of these MLs for potential applications in electronics.50 The photon energy absorption spectra of these MLs are shown in Figure 6. For the SiP, SiAs, and GeP monolayers, these semiconductors exhibit significant absorbance starting at approximately 2.0 eV, which is explicitly presented in the zoomed inset. In contrast, the optical absorbance for the GeAs

−1

)

b − b0 , b0

, and the calculated results

are listed in Table 2. As shown in Figure 3, the CBM significantly disperse along x direction, indicating small effective mass mxe in the four MLs. The value of mxe is smaller than that of monolayered phosphorene (0.17 m0)43 and antimonene.44 The effective mass along y direction mye is much lighter than the corresponding mxe value, implying strong anisotropy in electron transport. For the hole carrier, the mxh value is between 0.137 and 0.351 m0, which are also smaller than that of most 2D materials.45,46 Hence, the hole carrier along the y direction is lowered to 0.002 m0. Overall, these monolayered IV−V semiconductors are anticipated to have excellent and highly anisotropic transport properties, leading to a very promising perspective for future nanoelectronics. Recent investigations show that dramatic modulation of electronic properties in 2D semiconductors has been realized 5136

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces

2.08−2.64 eV according to accurate hybrid HSE06 functional. The semiconductors of the IV−V MLs exhibit high anisotropy in the effective mass of the charge carriers, suggesting strongly anisotropic electrical transport properties. These 2D monolayers can transformed between indirect and direct semiconductors and their bandgaps can be effectively tuned by applying an achievable uniaxial strain. The exfoliation energy of the IV−V monopnictide monolayers is found be be energetically favorable, therefore these semiconductor MLs are expected to be experimentally fabricated by using mechanical cleavage in the near future. The monolayered SiP, SiAs, GeP, and GeAs semiconductors can enrich the family of 2D nanomaterials and could have promising potentials for novel electronic and optoelectronic devices.

Figure 6. Calculated absorbance of these monolayered semiconductors as a function of photon energy. The magnified inset shows the light absorption at the region of solar spectrum.



is about 1.5 eV. This is because significant absorbance takes place as ℏω > Eg. Therefore, the monolayered GeAs has prior visible light absorption to the other three MLs, especially in the energy range of 1.5−3.6 eV. In general, the calculated absorbance spectra are in accordance with their corresponding bandgaps shown in Figure 3. In addition, the monolayered SiP exhibit strong ultraviolet absorption, and its absorbance peak is located at 5.0 eV. Thus, the vdW heterostructure consisting of SiP and GeAs could be promising 2D semiconductor candidate for photovoltaics as well as optoelectronics.51 To investigate the possibility of exfoliating these 2D MLs from the surface of its layered crystal, we have examined the exfoliation process and calculated the exfoliation energy varied with respect to the separation, as exhibited in Figure 7. The

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.7b17560. Calculated phonon dispersion spectra of the monolayered semiconductors; electronic band structures obtained from HSE for the monolayered semiconductors (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. ORCID

Jun Zhao: 0000-0001-7118-4992 Hui Zeng: 0000-0002-7657-6714 Author Contributions †

A.-Q.C. and Z.H. contributed equally to this work.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We are grateful to Dr. N. Xu for electronic structure simulations as well as inspirational discussions. We thank L. Zhang, Y. Sun, and Q.-Z. Tian for analyzing the simulation data. This work is financially supported by National Natural Science Foundation of China (Grants 11404037, 61701232, and 11304022).

Figure 7. Calculated exfoliation energies as a function of separation distance d, where d0 is the vdW distance between adjacent layers in bulk crystal.



exfoliation energy is only 0.26 J/m2 for the SiP and increased to 0.37 J/m2 for the GeAs. It is worthy to mention that the exfoliation energy of graphite is 0.32 J/m2, which is higher than that of Silicon monopnictides (i.e., SiP and SiAs) and comparable with Germanium monopnictides (i.e., GeP and GeAs).52,53 Therefore, the monolayered Silicon and Germanium monopnictides could be experimentally fabricated from its bulk crystal by mechanical cleavage. Moreover, compared to those of MoS2 and phosphorene, the exfoliation energy of the IV−V monopnictide monolayers is fairly favorable, implying their great potentials in constructing 2D nanoelectronic devices.54 In summary, we have studied the structures, electronic properties, stability, and possible exfoliation of the monolayered SiP, SiAs, GeP, and GeAs semiconductors. These monolayered monopnictide semiconductors are found to have excellent dynamical and thermal stabilities, and their bandgaps are 1.42− 1.90 eV calculated by GGA-PBE functional and are increased to

REFERENCES

(1) Miro, P.; Audiffred, M.; Heine, T. An Atlas of Two-Dimensional Materials. Chem. Soc. Rev. 2014, 43, 6537−6554. (2) Tan, C.; Cao, X.; Wu, X.-J.; He, Q.; Yang, J.; Zhang, X.; Chen, J.; Zhao, W.; Han, S.; Nam, G.-H.; Sindoro, M.; Zhang, H. Recent Advances in Ultrathin Two-Dimensional Nanomaterials. Chem. Rev. 2017, 117, 6225−6331. (3) Mannix, A. J.; Kiraly, B.; Hersam, M. C.; Guisinger, N. P. Synthesis and chemistry of elemental 2D materials. Nat. Rev. Chem. 2017, 1, 14. (4) Balendhran, S.; Walia, S.; Nili, H.; Sriram, S.; Bhaskaran, M. Elemental Analogues of Graphene: Silicene, Germanene, Stanene, and Phosphorene. Small 2015, 11, 640−652. (5) Choi, W.; Choudhary, N.; Han, G. H.; Park, J.; Akinwande, D.; Lee, Y. H. Recent development of two-dimensional transition metal dichalcogenides and their applications. Mater. Today 2017, 20, 116− 130. (6) Xu, K.; Yin, L.; Huang, Y.; Shifa, T. A.; Chu, J.; Wang, F.; Cheng, R.; Wang, Z.; He, J. Synthesis, properties and applications of 2D

5137

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces layered MIIIXVI (M = Ga, In; X = S, Se, Te) materials. Nanoscale 2016, 8, 16802−16818. (7) Gao, G.; Ding, G.; Li, J.; Yao, K.; Wu, M.; Qian, M. Monolayer MXenes: promising half-metals and spin gapless semiconductors. Nanoscale 2016, 8, 8986−8994. (8) Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically Thin MoS2: A New Direct-Gap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805. (9) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black phosphorus field-effect transistors. Nat. Nanotechnol. 2014, 9, 372−377. (10) Fan, Z.-Q.; Jiang, X.-W.; Luo, J.-W.; Jiao, L.-Y.; Huang, R.; Li, S.S.; Wang, L.-W. In-plane Schottky-barrier field-effect transistors based on 1T/2H heterojunctions of transition-metal dichalcogenides. Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 96, 165402. (11) Mak, K. F.; Shan, J. Photonics and optoelectronics of 2D semiconductor transition metal dichalcogenides. Nat. Photonics 2016, 10, 216−226. (12) Carvalho, A.; Wang, M.; Zhu, X.; Rodin, A. S.; Su, H.; Castro Neto, A. H. Phosphorene: from theory to applications. Nat. Rev. Mater. 2016, 1, 16061. (13) Chhowalla, M.; Jena, D.; Zhang, H. Two-dimensional semiconductors for transistors. Nat. Rev. Mater. 2016, 1, 16052. (14) Barreteau, C.; Michon, B.; Besnard, C.; Giannini, E. Highpressure melt growth and transport properties of SiP, SiAs, GeP, and GeAs 2D layered semiconductors. J. Cryst. Growth 2016, 443, 75−80. (15) Li, C.; Wang, S.; Zhang, X.; Jia, N.; Yu, T.; Zhu, M.; Liu, D.; Tao, X. Controllable seeded flux growth and optoelectronic properties of bulk o-SiP crystals. CrystEngComm 2017, 19, 6986−6991. (16) Manzeli, S.; Ovchinnikov, D.; Pasquier, D.; Yazyev, O. V.; Kis, A. 2D transition metal dichalcogenides. Nat. Rev. Mater. 2017, 2, 17033. (17) Beck, C. G.; Stickler, R. Crystallography of SiP and SiAs Single Crystals and of SiP Precipitates in Si. J. Appl. Phys. 1966, 37, 4683− 4687. (18) Chu, T. L.; Kelm, R. W., Jr.; Chu, S. S. C. Crystal Growth of Silicon Arsenide. J. Appl. Phys. 1971, 42, 1169−1173. (19) Kunioka, A.; Ho, K. K.; Sakai, Y. Optical properties of SiAs single crystals. J. Appl. Phys. 1973, 44, 1895−1896. (20) Lee, K.; Synnestvedt, S.; Bellard, M.; Kovnir, K. GeP and (Ge1−xSnx)(P1-yGey) (x ≈ 0.12, y ≈ 0.05): Synthesis, structure, and properties of two-dimensional layered tetrel phosphides. J. Solid State Chem. 2015, 224, 62−70. (21) Rau, J. W.; Kannewurf, C. R. Optical Absorption, Reflectivity, and Electrical Conductivity in GeAs and GeAs2. Phys. Rev. B 1971, 3, 2581−2587. (22) Lee, K.; Kamali, S.; Ericsson, T.; Bellard, M.; Kovnir, K. GeAs: Highly Anisotropic van der Waals Thermoelectric Material. Chem. Mater. 2016, 28, 2776−2785. (23) Kutzner, J.; Kortus, J.; Pätzold, O.; Wunderwald, U.; Irmer, G. Phonons in SiAs: Raman scattering study and DFT calculations. J. Raman Spectrosc. 2011, 42, 2132−2136. (24) Huang, B.; Zhuang, H.; Yoon, M.; Sumpter, B. G.; Wei, S.-H. Highly stable two-dimensional silicon phosphides: Different stoichiometries and exotic electronic properties. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 121401. (25) Shojaei, F.; Kang, H. S. Electronic structure of the germanium phosphide monolayer and Li-diffusion in its bilayer. Phys. Chem. Chem. Phys. 2016, 18, 32458−32465. (26) Wu, P.; Huang, M. Stability, bonding, and electronic properties of silicon and germanium arsenides. Phys. Status Solidi B 2016, 253, 862−867. (27) Zhang, S.; Guo, S.; Huang, Y.; Zhu, Z.; Cai, B.; Xie, M.; Zhou, W.; Zeng, H. Two-dimensional SiP: an unexplored direct band-gap semiconductor. 2D Mater. 2017, 4, 015030. (28) Zhou, L.; Guo, Y.; Zhao, J. GeAs and SiAs monolayers: Novel 2D semiconductors with suitable band structures. Phys. E 2018, 95, 149−153.

(29) Kresse, G.; Furthmuller, J. Efficient Iterative Schemes for ab initio Total-Energy Calculations using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (30) Kresse, G.; Furthmuller, J. Efficiency of ab initio total energy calculations for metals and semiconductors using a planewave basis set. Comput. Mater. Sci. 1996, 6, 15−50. (31) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (32) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (33) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188−5192. (34) Heyd, J.; Scuseria, G. E. Efficient hybrid density functional calculations in solids: Assessment of the Heyd-CScuseria-CErnzerhof screened Coulomb hybrid functional. J. Chem. Phys. 2004, 121, 1187− 1192. (35) Heyd, J.; Peralta, J. E.; Scuseria, G. E.; Martin, R. L. Energy band gaps and lattice parameters evaluated with the Heyd-ScuseriaErnzerhof screened hybrid functional. J. Chem. Phys. 2005, 123, 174101. (36) Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 2001, 73, 515−562. (37) Segall, M. D.; Lindan, P. J. D.; Probert, M. J.; Pickard, C. J.; Hasnip, P. J.; Clark, S. J.; Payne, M. C. First-Principles Simulation: Ideas, Illustrations and the CASTEP code. J. Phys.: Condens. Matter 2002, 14, 2717−2774. (38) Nosé. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511−519. (39) Gajdoš, M.; Hummer, K.; Kresse, G.; Furthmüller, J.; Bechstedt, F. Linear Optical Properties in the Projector-Augmented Wave Methodology. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 73, 045112. (40) Zhao, J.; Zeng, H. Two-dimensional Germanane and Germanane Ribbons: Density Functional Calculation of Structural, Electronic, Optical and Transport Properties and the Role of Defects. RSC Adv. 2016, 6, 28298−28307. (41) Bianco, E.; Butler, S.; Jiang, S.; Restrepo, O. D.; Windl, W.; Goldberger, J. E. Stability and Exfoliation of Germanane: A Germanium Graphane Analogue. ACS Nano 2013, 7, 4414−4421. (42) Kou, L.; Chen, C.; Smith, S. C. Phosphorene: Fabrication, Properties, and Applications. J. Phys. Chem. Lett. 2015, 6, 2794−2805. (43) Qiao, J.; Kong, X.; Hu, Z.-X.; Yang, F.; Ji, W. High-mobility transport anisotropy and linear dichroism in few-layer black phosphorus. Nat. Commun. 2014, 5, 4475. (44) Wang, Y.; Huang, P.; Ye, M.; Quhe, R.; Pan, Y.; Zhang, H.; Zhong, H.; Shi, J.; Lu, J. Many-body Effect, Carrier Mobility, and Device Performance of Hexagonal Arsenene and Antimonene. Chem. Mater. 2017, 29, 2191−2201. (45) Bhimanapati, G. R.; et al. Recent Advances in Two-Dimensional Materials beyond Graphene. ACS Nano 2015, 9, 11509−11539. (46) Guo, Y.; Zhou, S.; Zhang, J.; Bai, Y.; Zhao, J. Atomic structures and electronic properties of phosphorene grain boundaries. 2D Mater. 2016, 3, 025008. (47) Ahn, G. H.; Amani, M.; Rasool, H.; Lien, D.-H.; Mastandrea, J. P.; Ager, J. W., III; Dubey, M.; Chrzan, D. C.; Minor, A. M.; Javey, A. Strain-engineered growth of two-dimensional materials. Nat. Commun. 2017, 8, 608. (48) Kuroda, M. A.; Jiang, Z.; Povolotskyi, M.; Klimeck, G.; Newns, D. M.; Martyna, G. J. Anisotropic strain in SmSe and SmTe: Implications for electronic transport. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 245124. (49) Fan, Z.-Q.; Jiang, X.-W.; Wei, Z.; Luo, J.-W.; Li, S.-S. Tunable Electronic Structures of GeSe Nanosheets and Nanoribbons. J. Phys. Chem. C 2017, 121, 14373−14379. 5138

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139

Letter

ACS Applied Materials & Interfaces (50) Sun, Y.; Wang, R.; Liu, K. Substrate induced changes in atomically thin 2-dimensional semiconductors: Fundamentals, engineering, and applications. Appl. Phys. Rev. 2017, 4, 011301. (51) Novoselov, K. S.; Mishchenko, A.; Carvalho, A.; Castro Neto, A. H. 2D materials and van der Waals heterostructures. Science 2016, 353, aac9439. (52) Zacharia, R.; Ulbricht, H.; Hertel, T. Interlayer cohesive energy of graphite from thermal desorption of polyaromatic hydrocarbons. Phys. Rev. B: Condens. Matter Mater. Phys. 2004, 69, 155406. (53) Ziambaras, E.; Kleis, J.; Schröder, E.; Hyldgaard, P. Potassium intercalation in graphite: A van der Waals density-functional study. Phys. Rev. B: Condens. Matter Mater. Phys. 2007, 76, 155425. (54) Ajayan, P.; Kim, P.; Banerjee, K. Two-dimensional van der Waals materials. Phys. Today 2016, 69, 38−44.

5139

DOI: 10.1021/acsami.7b17560 ACS Appl. Mater. Interfaces 2018, 10, 5133−5139