Phase Coexistence and Strain-Induced ... - ACS Publications

Jun 13, 2018 - Hexagonal boron nitride is confirmed as an excellent substrate for supporting the β-BiAs film without perturbing the topological insul...
0 downloads 0 Views 7MB Size
Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

C: Physical Processes in Nanomaterials and Nanostructures

Phase Coexistence and Strain-Induced Topological Insulator in Two-Dimensional BiAs Tamiru Teshome, and Ayan Datta J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b05293 • Publication Date (Web): 13 Jun 2018 Downloaded from http://pubs.acs.org on June 19, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Phase Coexistence and Strain-Induced Topological Insulator in TwoDimensional BiAs Tamiru Teshome and Ayan Datta* Department of Spectroscopy, Indian Association for the Cultivation of Science, 2A, and 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India. Email: [email protected] ABSTRACT: Two-dimensional (2D) binary compounds have been recently reported as promising materials for achieving topological insulator and dissipationless transport devices. Based on first-principle calculations combined with a tight-binding (TB) model, we investigate a new class of 2D BiAs polymorphs which are energetically and dynamically stable. The monolayers of α-, β- and γ-BiAs allotropes show significant direct band-gap, while the δ- and ε-BiAs phases display indirect band gaps. Using strain-engineering along with spin-orbital coupling (SOC), β-BiAs transforms from a normal insulator to nontrivial topological phase. Under tensile strain and SOC, the bands are inverted resulting in a topological phase transition with a sizable band-gap of 0.28 eV. The nontrivial topological state is explicitly confirmed by calculating topological invariant, Z2 = 1 and the characteristic of edge states which are topologically protected in a Dirac cone at the ᴦ-point. Hexagonal boron nitride is confirmed as an excellent substrate for supporting the β-BiAs film without perturbing the topological insulator state. The present results indicate promise for 2D TIs at room temperature. 1.

INTRODUCTION

have not been synthesized till date and would be extremely

Recently, topological insulators (TIs), namely quantum spin

difficult to fabricate devices apart from toxicity of the heavy

Hall (QSH) insulators1,2 have attracted interests as a new state

atoms.9-22 Among these systems, Bi4Br4 and ZrTe5/HfTe5 exist

of quantum materials. These, novel electronic states possess

as three dimensional (3D) layered materials and their

bulk insulating gap and conduct charge and spin in helical

nontrivial band-gaps are practically useful to realize TI at

edge states by time-reversal symmetry. Though first proposed

room temperature.

by Kane and Mele for graphene, a very small opening of band gap ~10-3 meV at the Dirac point remains practically

Several group-V elements in buckled and planar hexagonal

inconsequential due to weak SOC for typical devices.3

structures such as phosphorene,23 arsenene,9 antimonene15-18

Therefore, a major bottleneck in for TIs is to find new

and bismuthene24, 25 exhibit topological phase transition from

materials with sizable band-gap opening at the Dirac-cone due

NI to TIs. At present, Bismuth based systems are an excellent

to controlled SOC.

materials for topological insulator due to their strong SOC. For example; the two well-known binary compounds - Bi2Se3

Recently, several 2D TI materials have been predicted like

and Bi2Te3 are strong TIs.26-27 On the other hand, arsenene in

silicene, germanene and stanene.4 However, so far, only the

α- and β-phase is energetically stable and converts from NI to

5

Bismuthene on a SiC substrate, HgTe/CdTe and InAs/GaSb

TI under suitable strain modification.9, 28-32 Therefore, seeking

quantum wells have been experimentally confirmed at a very

suitable 2D TIs with significantly large band gap is an

low temperature (below 10 K) limited by their small band

important area of research.

gap.6-8 Unfortunately, the majority of these predicted materials

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 11

Figure 1: Equilibrium structures for the polymorphs of Bismuth Arsenic (BiAs) in the top and side views (a) α-BiAs, (b) βBiAs, (c) γ-BiAs, (d) δ-BiAs and (e) ε-BiAs. The shaded regions were shown to be unit-cell of the polymorphs; pink spheres represent Bismuth (Bi) atoms and green Arsenic (As) atoms. Hence, a NI to TI crossover can be in particular interesting in

were performed to verify the stabilities at room temperature.

2D materials system for quantum devices where NI/TI

Wannier

Charge

Center

(WCC)

implemented

in

38

switchability can be tuned under electric or superconducting

WANNIER90 package to identify the topological invariance

fields. To the best of our knowledge, the electronic and

characterized by the Z2 number, which is implemented in Z2

geometric

pack.

phase-diagram

of

a

nontrivial

topological

compound (BiAs) with appreciable band-gap with respect to

3.

its unique polymorphs has not been elucidated till date.

Results and Discussion

In this article, we have investigated 2D BiAs polymorphs and their electronic structures based on first principle

2.

calculation, and these new materials labeled as α-, β-, γ -,

Computational Details

The calculations were performed using the Vienna ab initio simulation package (VASP).

33

δ-, and ε-BiAs forms in their stable structures are shown in

The exchange-correlation term

Fig. 1(a-e). These five different phases arise from relative

is described within the generalized gradient approximation

connectives between the atoms in the 2D lattices. There

(GGA) parameterized by the Perdew-Burke-Ernzerhof (PBE)

are four atoms per unit cell in the α-BiAs and γ-BiAs as

functional.34 A vacuum space of 20 Å is set to avoid the

shown in Fig. 1(a, c), while there are eight atoms per unit

interaction between layers caused by the periodic boundary

cell in the δ-BiAs and ε-BiAs phase as shown in Fig. 1(d,

condition. The kinetic-energy cutoff plane-wave expansion is

e). On the other hand, β-BiAs with a hexagonal structure

set to be 500 eV. All the atoms in the unit cell are fully relaxed

contains two atoms per unit cell and the buckling height,

until the force on each atom are less than 0.01 eV/Å and

lattice constant and bond length were found to be 1.56 Å,

-5

convergence threshold was set to be 10 eV of the energy.

2.75 Å and 4.03 Å, respectively as shown in Fig. 1(b). The

Furthermore, Monkhorst-Pack k-point grid of 11 × 11 × 1 was

structural details of all the polymorphs are reported in

35

adopted. The strain modification was described as 𝜀 =

Table 1.

∆𝑎/𝑎! , where 𝑎! is the equilibrium lattice and ∆𝑎 + 𝑎! is the 3.1 Dynamical Stabilities and Electronic properties of

strained modified lattice. SOC effects were included in SCF 36

BiAs polymorphs

electronic structure calculations. The phonon calculations are carried out using the PHONOPY

37

To evaluate the relative stability of 2D-BiAs polymorphs, we

code combined with

calculated the cohesive energies as the following expression:

density functional perturbation theory (DFPT) method in

𝐸!"! = (𝐸!"!!"! − 𝑛𝐸!" − 𝑛𝐸!" )/𝑛

VASP and ab-initio molecular dynamics (MD) simulations

2

ACS Paragon Plus Environment

Page 3 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

where 𝐸!"!!"! , 𝐸!" and 𝐸!" being the total energies of BiAs

The relative cohesive energies of 2D BiAs polymorphs with

allotropes, single Bi atom and single As atom, respectively.

respect to that of α-BiAs are illustrated in Fig. 2(f) and clearly, β-BiAs is the most stable structure. The dynamical stability of the BiAs allotropes was further verified by calculating their vibrational spectra as shown in Fig. 2(a-e). The absences of imaginary frequencies show that all the five BiAs polymorphs are dynamically stable. Additionally, molecular dynamics (MD) simulations were performed on the five BiAs polymorphs for 10 ps with a time step of 1.0 fs at 300 K. In all polymorphs, the T vs time graph remains the unchanged for fluctuations between 290 and 310 K are shown in Fig. S2. The MD simulations indicate that these BiAs polymorphs are also thermally stable at room temperature.

To understand how

significantly the band-gaps depend on the structure of the

Figure 2: Phonon spectrum of (a) α-BiAs, (b) β-BiAs, (c) γ-

polymorph, the band structures are calculated for α-, β-, γ -, δ-,

BiAs, (d) δ-BiAs, (e) ε-BiAs and (f) relative cohesive energy

and ε-BiAs monolayers as shown in Fig 3(a-e) without SOC

with respect to α-BiAs.

and Fig. 3 (f-j) with SOC effect.

Figure 3: Electronic band structures of BiAs polymorphs without and with SOC (a, f) α-BiAs, (b, g) β-BiAs, (c, h) γ-BiAs, (d, i) δBiAs, (e, j) ε- BiAs, respectively.

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 11

Table 1: Equilibrium BiAs polymorphs: Lattice parameters a and b (Å), bond length; d (Å), bond-angle; θ (°), dihedral angle; ɸ (°), buckling height; Δh (Å), g band gap value calculated PBE without and with SOC; (eV), Cohesive energy; Ecoh (eV/atom) and ΔE!"! = 𝐸!"! − 𝐸!"! (α-BiAs); ΔEcoh (eV/atom) relative stability of polymorph with respect to α-BiAs. Model

a

b

d

θ

ɸ

α-BiAs β-BiAs γ-BiAs δ-BiAs ε-BiAs

4.17 4.03 6.91 5.94 6.34

4.81 4.03 4.03 6.14 7.89

2.75 2.79 2.80 2.75 2.74

97.29 92.31 92.99 96.52 98.93

82.03 87.59 87.29 9.76 0.38

Δh

gPBE

gsoc+PBE

2.71 1.56 1.74 2.72 2.11

1.27 1.02 1.01 0.61 0.75

0.58 0.69 0.48 0.45 0.51

Ecoh

ΔEcoh

-1.81 -2.46 -0.88 -0.76 -1.38

0.00 -0.65 0.93 1.10 0.43

Figure 4: Electronic band structures of β-BiAs without (a-d) and with SOC (e-h) by applying different values of tensile strains. The Fermi level set to be 0 eV. We found that α-BiAs, β-BiAs and γ-BiAs are direct band-gap

3.2 Topological Phase Transition in Strained β-BiAs

semiconductors at ᴦ-point, while others δ-BiAs and ε-BiAs are

Strain-engineering is a powerful approach for tuning the

indirect semiconductors regardless of whether SOC effects are

electronic and topological properties of 2D materials.21 When

included or not. For α-BiAs, the VBM and CBM are located at

a biaxial tensile strain is induced, Bi-As is stretched and the

ᴦ-point with the band gap of 1.27 eV and 0.58 eV without and

lattices are relaxed consequently (Fig. S3). We critically

with SOC. Similarly, for β-BiAs, the direct band-gaps are 1.02

evaluate the band structure evolution under biaxial tensile

eV (0.62 eV at HSE06 level, see supp. Inf. Fig. S4) without

strain for β-BiAs as it is the most stable phase as shown in Fig.

and 0.69 eV with SOC. For γ-BiAs without SOC the VBM

4(a-h). It is also interesting to explore if a nontrivial

and CBM are located at ᴦ-point resulting in band-gap = 1.01

topological

eV (0.48 eV with SOC). Clearly, α, β and γ polymorphs of

configuration of β-BiAs under a suitable biaxial tensile strain.

BiAs have reasonable direct band-gaps which make them

The variation of band-gaps under tensile strain with and

prospective candidate for TIs under external perturbations like

without SOC for β-BiAs is shown in Fig. 5(a). We started

strain or electric field.

from the semiconducting phase β-BiAs (ε = 0 %) having a

phase

transition

occurs

for

the

buckled

direct band gap of 1.02 eV at its equilibrium structure (Fig. S4). Till a tensile strain applied, ε < 12 %, it behaves as a

4

ACS Paragon Plus Environment

Page 5 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry !

semiconductors without and with SOC effect as shown in Fig.

𝜀!" = −8.70 𝑒𝑉.40 The optimal values for β-BiAs are

4(a-h). Further increasing the strain, the conduction band

𝑉!!" = −1.162 𝑒𝑉,

progressively shifts downward to the Fermi level. On the other

and 𝑉!!" = −0.672 𝑒𝑉, respectively. With the above set of

hand, the valance band shifts upwards leading to semi-metal at

parameters, the TB bands are in reasonable agreement with the

ε = 12 % without SOC and hence, no band inversion observed

DFT results as shown in Fig. 5(b-d). The p-orbitals near the

as shown in Fig. 5(c). Clearly, at ε > 12 % without SOC, β-

Fermi level split into two groups (𝑝! and 𝑝! ) orbitals

BiAs is a semi-metals. As the strain is increased further a band

𝑉!"# = 1.122 𝑒𝑉,

𝑉!!" = 2.689 𝑒𝑉

and 𝑝! orbital, respectively. Projecting 𝑝! , 𝑝! and 𝑝! bands

inversion occurs at ᴦ-point in presence of SOC (Eg = 0.28 eV).

onto different atomic orbitals, we found that the energy

Therefore, with SOC and ε > 12 %, one observes a phase

spectrum of β-BiAs near the Fermi level is modified mostly to

transition for a normal insulator to topological insulator (Fig.

𝑝! and 𝑝! orbitals.

5(a)). Our predicted nontrivial band topology is in excellent agreement with the results previously reported for Bismuth Oxide film.39 In fact, both Arsenic and Bismuth have strong SOC effects and hence, band inversion can be achieved when a large strain is applied.9, 32To identify the mechanism of the band inversion, we propose a tight-binding model of 𝑝! , 𝑝! and 𝑝! orbitals. The effective Hamiltonian is considered as:

Here,

𝑐!!!

!

𝑡!" ( 𝑐!!! 𝑐! + ℎ. 𝑐. ) !,! ,!,!

!,!

𝜀!! ,

!"

𝜀!! 𝑐!!! 𝑐!! +

𝐻!" =

and

𝑐!!

represent the on-site energy, creation,

and annihilation operators of an electron at the α-orbital of the !"

i-th atom. The parameter 𝑡!" is the nearest-neighbor hopping energy of an electron between α-orbital of i-th atom and βorbital of j-th atom, α, β ∈ (𝑝! , 𝑝! , 𝑝! ), which can be

Figure 5: (a) The band gaps of β-BiAs without and with SOC

performed by fitting the DFT results as shown in Fig. 5 (b, c

at the different strain. For ε < 12 % without and with SOC are

and d). According to tight-binding theory, the hopping

normal insulators (Cyan color region), [12, 16] % with SOC

energies can be represented as:

are topological insulators (pink color region) and [12, 16] %

!! 𝑡!! = 𝑉!!"

without SOC are semi-metal materials (yellow region). DFT-

!!

𝑡!" ! = 𝑉!"# × 𝑐𝑜𝑠𝜃

PBE and MLWF fitted band structures at (b) equilibrium state

!!! 𝑡!"

= 𝑉!"# × 𝑐𝑜𝑠𝜑

and (c, d) without and with SOC effect at 12 % strain,

! !!

= 𝑉!!" × 𝑐𝑜𝑠 ! 𝜃 + 𝑉!!" × 𝑠𝑖𝑛! 𝜃

respectively. The TB model data are indicated in red dotted

!! !!

= 𝑉!!" × 𝑐𝑜𝑠 ! 𝜑 + 𝑉!!" × 𝑠𝑖𝑛! 𝜑

𝑡!"! 𝑡!"

!! !!

𝑡!"

line and blue color refers to DFT-PBE (b, c, and d).

= 𝑉!!" − 𝑉!!" ×𝑐𝑜𝑠𝜃 ×𝑐𝑜𝑠𝜑

Topological phase transition in β-BiAs is confirmed by

Where θ and 𝜑 are the angles of the vector pointing from i-th

calculating the topological invariant Z2 index based on the

atom to j-th atom with respect to x and y axis. The on-site

U(2N) non-Abelian Berry connection proposed by Rui Yu and

! = energies of s and p orbitals are set to the values of 𝜀!"

co-worker.41

−17.68 𝑒𝑉,

!

𝜀!" = −8.30 𝑒𝑉,

! 𝜀!" = −15.01 𝑒𝑉,

and

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 11

Figure 6: (a) The calculation of Wilson loop (Wannier Charge Center) along ky for β-BiAs monolayer at 12 % strain yielding Z2 = 1, (b) Dirac edge state of β-BiAs at 12 % biaxial strain and (c) Edge spins polarization.

Figure 7: (a) Side and top view of epitaxial growth β-BiAs of TI on three layers of h-BN substrate, (b) band structure without and (c) with SOC at 12 % strain. The Z2 topological invariant is related by counting the number

utilizing the Green's function method with the tight-binding

of crossing between any arbitrary horizontal reference line and

models from VASP and wanniner90 package based on

evaluation of θ mod 2π, where the odd and even numbers

maximally localized wannier function (MLWFs) to calculate

denoted Z2 = 1 and 0 nontrivial and trivial topological

the edge state. In Fig. 6(b) we show the edge bands perfectly

phase, respectively. Figure 6(a) shows that the WCC

connected the conduction band and valence band with 1D

evolution along ky for the time-reversal plane of β-BiAs at 12

gapless edge states at ᴦ-point. The counter propagating edge

% strain when SOC is considered to yield TI with topological

states exhibit spin up and spin down polarization contribution

invariant, Z2 = 1. Additionally, to confirm the topological

in edge spectral function as shown in Fig. 6(c). We recognize

character of β-BiAs at 12 % strain with SOC, one needs to

that the spin-momenta of those Dirac-type edge states are

also show the existence of helical edge states protected by

locked at the ᴦ-point within the bulk gap thereby confirming

time-reversal symmetry, which is an essential characteristic of

the existence of nontrivial topology leading to absolutely

2D TIs. We have performed the band structure of the ribbon

polarized conductive channels.

6

ACS Paragon Plus Environment

Page 7 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry dynamically under such large tensile strain up to 16 % as

Nevertheless to have practical application, it is highly

presented in Fig. S1(a-c) with the onset of the topological

desirable to experimentally realize a substrate to support 2D

phase transition at 12 % strain. However, for ε = 17.2 % the

films. Here, we choose three layers of hexagonal boron nitride

material shows small imaginary frequency are observed (υ = -

42-45

0.0928 THz) which indicates that the β-BiAs becomes

with a good dielectric constant for the epitaxial growth for β-

unstable beyond such strains. Up to now, graphene have been

BiAs. A lattice mismatch of 0.09 Å and the interlayer

subjected to a large strain by laser shock-induced straining

distance=4.02 Å, results in stable van der Waals (vdW)

process on flexible polydimethylsiloxane (PDMS) substrate,

heterostructure. The binding energy is calculated to be -0.028

and also the indentation method carried out with an atomic

eV per unit cell indicating weak dispersion interactions.

force microscope (AFM).53 Kim et al. obtained under large

Interestingly,

β-BiAs/h-BN

tensile strain of 18.7 % for stretched graphene.54 In an AFM

heterostructure without and with SOC at strain 12 % remains

indentation experiment, MoS2 monolayer could be strained

unperturbed as shown in Fig. 7 (b and c), respectively. There

upto 11%.55 For the purpose of application based on

is no charge transfer between β-BiAs and h-BN as confirmed

electronics and spintronics device technology (spin-controlled

by Bader charge analysis, thus the 𝑝! and 𝑝! orbitals near the

phase transition can be enhanced in spin transport) tunability

fermi level are dominantly contributed by Bi and As atoms

of the spin-splitting by applying large tensile strain has been

and the charge density difference being calculated as shown in

realized in many 2D materials.56-59 In this context, we believe

Fig. S6. Recently, 2D TIs materials such as bismuthene5,

that our prediction of TI in β-BiAs under a tensile strain could

plumbene,46

be realized experimentally.

(h-BN) as shown in Fig. 7(a), which has large band gap

the

band

germanene

structure

and

of

stanene

bilayer,47

SiGe,21

arsenene,32 arsenene oxide48 and Bi monolayer49 have been successfully

grown

on

different

substrates

4.

materials.

In conclusion, we predicted novel of 2D BiAs polymorphs,

Evidently, the β-BiAs/h-BN heterostructure is a robust

which are dynamically and thermally stable. α-, β- and γ-BiAs

topological insulator who’s the band inversion is not perturbed

are direct band-gap semiconductors, while the δ-, and ε-BiAs

by substrate.

phases display indirect band gap. Interestingly, β-BiAs has a nontrivial phase with significantly large band-gap. We expect

Generally, most of the 2D nontrivial topology proposed till

that the excellent TI behavior of β-BiAs with significant band-

date with sizable bulk gaps is composed of heavy metals,

gap can result in QSH devices at room temperature. For the

which have strong SOC interaction, including Pb, Hg, Sn, Sb,

practical application, we also proposed h-BN as excellent

and Bi. Although, 2D TIs has also achieved in some light

candidate substrate for supporting the β-BiAs film without

elements but their bulk gaps are very small induced by weak

perturbing the topological insulator.

SOC, such as, the gaps are ~10-3 meV for graphene,48,49 silicene (1.55-2.9 meV),4 germanene (Eg = 23.9-108 meV),4, 5052

Conclusions

ASSOCIATED CONTENT

arsenene (Eg = 696 meV),9 few-layered black phosphorus (5

Supporting Information

meV),22 BiNH2 (Eg = 0.83 eV)49 and Arsenene Oxide (Eg =

Phonon dispersion, Details of AIMD simulations, variation

232meV).48 In comparison, our first principle calculation

of strain respect to angle, lattice constant and bond-length

showed that 2D β-BiAs, a newly predicted binary compound

for

has a sizable gap up to 0.28 eV, suggesting it as a promising

β-BiAs,

structures

at

equilibrium

state

contribution of orbitals, PBE and HSE06, Partial density of

nontrivial topological insulator at room temperature. As a

states (PDOS) at ε = 12 % CBM and VBM, without and

result of buckling geometry the β-BiAs can be stable

7

Band

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 11

with SOC, respectively, β-BiAs monolayer deposited on h-

Transition in HgTe Quantum Wells, Science 2006, 314,

BN substrate and charge density difference.

1757-1761 9.

Zhang, H. J.; Ma, Y.; Chen, Z. F. Quantum Spin Hall

AUTHOR INFORMATION

Insulators in Strain-Modified Arsenene, Nanoscale 2015,

Corresponding Author

7, 19152-19159. 10. Wang, Z. F.; Su, N. H.; Liu, F. Prediction of a Two-

E-mail: [email protected]. Tel.: +91-33-24734971. Notes

Dimensional Organic Topological Insulator, Nano Lett.

The authors declare no competing financial interest.

2013, 13, 2842-2845. 11. Zhou, L. J.; Kou, L. Z.; Sun, Y.; Felser, C.; Hu, F. M.;

ACKNOWLEDGEMENT

Shan, G. C.; Smith, S. C.; Yan, B. H.; Frauenheim, T.,

We acknowledge The World Academy of Sciences (TWAS) -

New Family of Quantum Spin Hall Insulators in Two-

Indian Association for the Cultivation of Science (IACS) (FR

Dimensional

number: 3240280472) financial support and AD thanks, DST,

Nontrivial Band Gaps, Nano Lett. 2015, 15, 7867-7872.

INSA and BRNS for partial funding.

Transition

Metal

Halide

with

Large

12. Ma, Y. D.; Kou, L. Z.; Li, X.; Dai, Y.; Smith, S. C.; Heine, T. Quantum spin Hall Effect and Topological

References

Phase Transition in Two-Dimensional Square Transition-

1.

Metal Dichalcogenides, Phys. Rev. B 2015, 92, 085427.

Fu, L.; Kane, C. L. Topological Insulators with Inversion Symmetry, Phys. Rev. B 2007, 76, 045302.

2. 3.

13. Nie, S. M.; Song, Z.; Weng, H. M.; Fang, Z. Quantum

Fu, L.; Kane, C. L.; Mele, E. J. Topological Insulators in

Spin Hall effect in Two-Dimensional Transition-Metal

Three Dimensions. Phys. Rev. Lett. 2007, 98, 106803.

Dichalcogenides Haeckelites, Phys. Rev. B 2015, 91,

Kane, C. L.; Mele, E. J. Quantum Spin Hall Effect in

235434.

Graphene. Phys. Rev. Lett. 2005, 95, 226802. 4.

14. Qian, X. F.; Liu, J.W.; Fu, L.; Li, J. Quantum spin Hall

Liu, C.-C.; Feng, W. X.; Yao, Y. G. Quantum Spin Hall

Effect

Effect in Silicene and Two-Dimensional Germanium

Dichalcogenides, Science 2014, 346, 1344-1347.

Hanke, W.; Thomale, R.; Schäfer, J. Claessen, R.

Two-Dimensional Topological Insulators Consisting of

Bismuthene on a SiC substrate: A candidate for a high-

Bilayers of Group III Elements with Bi, Nano Lett. 2014,

temperature quantum spin Hall material, Science 2017,

14, 2505-2508. 16. Ma, Y.; Kou, L.; Dai, Y.; Heine, T. Two-Dimensional

König, M.; Wiedmann, S.; Bröune, C.; Roth, A.;

Topological

Buhmann, H.; Molenkamp, L. W.; Qi, X. L.; Zhang, S. C.

Compounds: M2Te (M = Cu, Ag), Phys. Rev. B 2016, 93,

Quantum spin Hall Insulator State in HgTe Quantum

235451.

Insulators

in

Group-11

Knez, I.; Du, R. R.; Sullivan, G. Evidence for Helical

Topological Phase Diagrams of Bulk and Monolayer TiS2-

Edge Modes in Inverted InAs/GaSb Quantum Wells,

xTex

, Phys. Rev. Lett. 2013, 110, 077202.

18. Ma, Y. Y. D.; Dai, L. Z.; Kou, T. F.; Heine, T. Robust

Bernevig, B. A.; Hughes, T. L.; Zhang, Shou-Cheng,

Two-Dimensional Topological Insulators in Methyl

Quantum Spin Hall Effect and Topological Phase

Functionalized Bismuth, Antimony, and Lead Bilayer Films, Nano Lett. 2015, 15, 1083-1089.

8

Chalcogenide

17. Zhu, Z. Y.; Cheng, Y. C.; Schwingenschlöogl, U.

Phys. Rev. Lett. 2011, 107, 136603. 8.

Metal

C.-H.; Das, T.; Lin, H. B., A., Prediction of Large Gap

Wells, Science 2007, 318, 766-770. 7.

Transition

Reis, F.; Li, G.; Dudy, L.; Bauernfeind, M.; Glass, S.;

357, 287-290. 6.

Two-Dimensional

15. Chuang, F.-C.; Yao, L.-Z.; Huang, Z.-Q.; Liu, Y.-T.; Hsu,

Phys. Rev. Lett. 2011, 107, 076802. 5.

in

ACS Paragon Plus Environment

Page 9 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Semiconducting Properties in Few-Layer Orthorhombic

19. Luo, W.; Xiang, H. J. Room Temperature Quantum Spin

Arsenene, Appl. Phys. Express 2015, 8, 055201.

Hall Insulators with a Buckled Square Lattice, Nano Lett.

30. Wang, Y.; Ding, Y. Unexpected Buckled Structures and

2015, 15, 3230-3235. 20. Jin, K.-H., Jhi, S.-H. Quantum Anomalous Hall and

Tunable Electronic Properties in Arsenic nanosheets:

Quantum Spin-Hall Phases in Flattened Bi and Sb

Insights from First-Principles Calculations, J. Phys.:

Bilayers, Sci. Rep. 2015, 5, 8426.

Condens. Matter 2015, 27, 225304. 31. Zhu, Z.; Guan, J.; Tomanek, D. Strain-Induced Metal

21. Teshome, T.; Datta, A. Topological Insulator in TwoDimensional SiGe Induced by Biaxial Tensile Strain, ACS

Semiconductor Transition in Monolayers and Bilayers of

Omega, 2018, 3(1), 1-7

Gray Arsenic: A Computational Study, Phys. Rev. B 2015, 91, 161404.

22. Liu, Q. H.; Zhang, X.W.; Abdalla, L. B.; Fazzio, A. Zunger, A. Switching a Normal Insulator into a

32. Wang, Y. P.; Zhang, C. W.; Ji, W. X.; Zhang, R. W.; Li, P.;

Topological Insulator via Electric Field with Application

Wang, P. J.; Ren, M. J.; Chen, X. L.; Yuan, M. Tunable

to Phosphorene, Nano Lett. 2015, 15, 1222-1228.

Quantum Spin Hall Effect via Strain in Two-Dimensional Arsenene Monolayer, J. Phys. D: Appl. Phys., 2016, 49,

23. Cheng, L.; Liu, H. J.; Tan, X. J.; Zhang, J.; Wei, J.; Lv, H.

055305.

Y.; Shi, J.; Tang, X. F. Thermoelectric Properties of a

33. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized

Monolayer Bismuth, J. Phys. Chem. C 2014, 118, 904-

Gradient Approximation Made Simple, Phys. Rev. Lett.

910. 24. Freitas, R. R. Q.; Rivelino, R..; Mota, F. de Brito, Castilho,

1996, 77, 3865.

C. M. C. de, Kakanakova-Georgieva, A.; Gueorguiev, G.

34. Markov, G.; Payne, M. C. Periodic Boundary Conditions

K. Topological Insulating Phases in Two-Dimensional

in Ab Initio Calculations, Phys. Rev. B 1995, 51, 4014-

Bismuth-Containing

4022.

Single

Layers

Preserved

by

Hydrogenation, J. Phys. Chem. C, 2015, 119, 23599-23606.

35. Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-

25. Zhang, H.; Liu, C. X.; Qi, X. L.; Dai, X.; Fang, Z.; Zhang,

Zone Integrations, Phys. Rev. B 1976, 13, 5188. 36. Qi, X. L.; Zhang, S. C. Topological Insulators, and

S. C. Topological Insulators in Bi2Se3, Bi2Te3, and

Superconductors, Rev. Mod. Phys. 2011, 83, 1057.

Sb2Te3 with A Single Dirac Cone on the Surface, Nat. Phys., 2009, 5, 438-442.

37. Togo, A.; Tanaka, I. First Principles Phonon Calculations

26. Shi, H.; Parker, D.; Du, M.-H.; Singh, D. Connecting

in Materials Science, Scr. Mater. 2015, 108, 1-5

Thermoelectric Performance, and Topological-Insulator

38. Mostofi, A. A.; Yates, J. R.; Lee, Y. S.; Souza, I.;

Behavior: Bi2Te3 and Bi2Te2Se from First Principles, J.

Vanderbilt, D.; Marzari, N. Wannier90: A tool for

Phys. Rev. Appl. 2015, 3, 014004.

obtaining

Wannier

Functions,

Comput. Phys. Commun. 2008, 178, 685.

27. Crowley, J. M.; Tahir-Kheli, J.; Goddard III, W. A.; Accurate Ab Initio Quantum Mechanics Simulations of

39. Miao, L.; Yao, M. Y.; Ming, W.; Zhu, F.; Han, C. Q.;

Bi2Se3 and Bi2Te3 Topological Insulator Surfaces, J.

Wang, Z. F.; Guan, D. D.; Gao, C. L.; Liu, C.; Liu, F.;

Phys. Chem. Lett., 2015, 6 (19), 3792-3796.

Qian, D.; Jia, J. F. Evolution of the Electronic Structure in

28. Kamal, C.; Ezawa, M. Arsenene: 2D Buckled and

Ultrathin Bi(111) Films, Phys. Rev. B 2015, 91, 205414

Puckered Honeycomb Arsenic Systems, Phys. Rev. B

40. Vogl, P.; Hjamarson, H. P.; Dow, J. D. A Semi-Empirical

2015, 91, 085423.

Tight-Binding Theory of the Electronic Structure of

29. Zhang, Z.; Y.; Xie, J. F.; Yang, D. Z.; Wang, Y. H.;

Semiconductors, J. Phys. Chem. Solids 1983, 44, 365-

Si, M. S.; Xue, D. S. Manifestation of Unexpected

378.

9

Maximally-Localized

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

41. Yui, R., Qi, X. L., Bernevig, A., Fang, Z., Dai, X. Equivalent Expression of Ζ2 Topological Invariant for Band

Insulators

using

the

Nnon-Abelian

Berry

42. Liu, Z.; Song, L.; Zhao, S.; Huang, J.; Ma, L.; Zhang, J.; J.;

Ajayan,

P.

M.

Direct

Calculations. Phys. Rev. B: Condens. Matter, 2007, 75, 041401. 51. Gmitra; M., Konschuh; S., Ertler; C., Ambrosch-Draxl; C., Fabian; J. Band-Structure Topologies of Graphene:

Connection, Phys. Rev. B 2011, 84, 075119. Lou,

Page 10 of 11

Growth

of

Graphene/Hexagonal Boron Nitride Stacked Layers. Nano

Spin-Orbit Coupling Effects from First Principles. Phys. Rev. B: Condens. Matter, 2009, 80, 2354. 52. Liu, J. W.; Qian, X. F.; Fu, L. Crystal Field Effect Induced

Lett. 2011, 11, 2032-2037. 43. Ö zcelik, V. O.; Cahangirov, S.; Ciraci, S. Epitaxial Growth; ̧ Mechanisms of Graphene and Effects of

Topological

Crystalline

Insulators

in

Monolayer IV-VI semiconductors. Nano Lett. 2015, 15, 2657-2661. 53. Li, J.; Chung, T.-F.; Chen, Y. P.; Cheng, G. J.;

Substrates. Phys. Rev. B 2012, 85, 235456. I.;

Nanoscale Strainability of Graphene by Laser Shock-

Lichtenstein, A. I. Adhesion and Electronic Structure of

Induced Three Dimensional Shaping, Nano Lett.,

Graphene on Hexagonal Boron Nitride Substrates. Phys.

2012, 12, 4577-4583.

44. Sachs,

B.;

Wehling,

T.

O.;

Katsnelson,

M.

54. Lee, C.; Wei, X.; Kysar, J. W.; Hone, J. Measurement

Rev. B 2011, 84, 195414. 45. Teshome, T.; Datta, A. Two-Dimensional Graphene-Gold

of Elastic Properties and Intrinsic Strength of

Interfaces Serve as Robust Templates for Dielectric

Monolayer Graphene, Science, 2008, 321, 385-388.

Capacitors. ACS Appl. Mater. Interfaces 2017, 9, 34213-

55. Bertolazzi, S.; Brivio, J.; Kis, A. Stretching and Breaking of Ultrathin MoS2, ACS Nano, 2011, 5,

34220. 46. Zhao, H.; Zhang, C.; Ji, W.; Zhang, R.; Li, S.; Yan, S.

9703-9709.

Unexpected Giant-Gap Quantum Spin Hall Insulator in

56. Kim, K. S.; Zhao, Y.; Jang, H.; Lee, S. Y.; Kim, J.

Chemically Decorated Plumbene Monolayer, Sci. Rep.

M.; Kim, K. S.; Ahn, J.-H.; Kim, P.; Choi, J.-Y.;

2016, 6, 20152.

Hong, B. H. Large-Scale Pattern Growth of Graphene

47. Huang, C.; Zhou, J.; Wu, H.; Deng, K.; Jena, P.; Kan, E. Quantum Phase transition in Germanene and Stanene Bilayer: From Normal metal to Topological Insulator. J.

Films for Stretchable Transparent Electrodes, Nature, 2009, 457, 706-710. 57. Ma, Y.; Dai, Y.; Wei, W.; Li, X.; Huang, B. Emergence of Electric Polarity in BiTeX (X = Br and

Phys. Chem. Lett. 2016, 7, 1919-1925. 48. Wang, Y.; Ji, W.; Zhang, C.; Li, P.; Zhang, S.; Wang, P.;

I) Monolayers and the Giant Rashba Spin Splitting,

Li, S.; Yan, S. Two-Dimensional Arsenene oxide: A

Phys. Chem. Chem. Phys. 2014, 16(33), 17603-17606.

Realistic Large-Gap Quantum Spin Hall Insulator, Appl.

58. Ouyang, B.; Song, J. Strain Engineering of Magnetic States

Phys. Lett. 2017, 110, 213101. 49. Li, S.; Ji, W.; Hu, S.; Zhang, C.; Yan, S. Effect of Amidogen Functionalization on Quantum Spin Hall

Vacancy-Decorated

Hexagonal

Boron

Nitride, Appl. Phys. Lett., 2013, 103, 102401-10240. 59. Mirό,

P.;

Ghorbani-Asl,

M.;

Materials

beyond

Heine, MoS2:

T.

Two Noble-

Effect in Bi/Sb(111) Films, ACS Appl. Mater. Interfaces

Dimensional

2017, 9, 41443-41453.

Transition-Metal Dichalcogenides, Angew. Chem. Int. Ed. 2014, 53, 3015-3018

50. Yao; Y., Ye; F., Qi; X. L., Zhang; S. C., Fang; Z. Spin-Orbit

Gap

of

Graphene:

First-Principles

10

of

ACS Paragon Plus Environment

Page 11 of 11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC Graphic

11

ACS Paragon Plus Environment