Flexible and Anisotropic Properties of Monolayer MX2 (M = Tc and Re

Oct 2, 2017 - Mechanical, electronic, and dielectric properties of VIIB-VIA 2D TMDs have been performed through the Vienna ab initio Simulation Packag...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV NEW ORLEANS

Article

Flexible and Anisotropic Properties of Monolayer MX (M = Tc and Re; X = S, Se) 2

Qiyi Zhao, Yaohui Guo, Yixuan Zhou, Xiang Xu, Zhaoyu Ren, Jin Tao Bai, and Xin Long Xu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b07939 • Publication Date (Web): 02 Oct 2017 Downloaded from http://pubs.acs.org on October 2, 2017

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 31

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

Flexible and Anisotropic Properties of Monolayer MX2 (M = Tc and Re; X = S, Se) Qiyi Zhao, Yaohui Guo, Yixuan Zhou, Xiang Xu, Zhaoyu Ren, Jintao Bai, Xinlong Xu*

Shaanxi Joint Lab of Graphene, State Key Lab Incubation Base of Photoelectric Technology and Functional Materials, International Collaborative Center on Photoelectric Technology and Nano Functional Materials, Institute of Photonics & Photon-Technology, Northwest University, Xi'an 710069, China

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 31

ABSTRACT: Flexible and anisotropic response of layered transition metal dichalcogenides MX2 (M = Tc and Re; X = S, Se) are important for wearable and polarized optoelectronics. Herein, the elastic, electronic, and optical dielectric properties of these two-dimensional (2D) materials have been investigated by density functional theory (DFT) with different van der Waals correction and Heyd-Scuseria-Ernzerhof hybrid functional. The Young’s modulus of these materials is low, which indicates that they are favorable for the flexible optoelectronic devices. The band gaps fall in between 1.70 eV and 2.12 eV with the d states of transition metal atoms playing an important role in conduction and valence bands. In addition, the appearance of band nesting implies that there are strong light-matter interactions in these materials, indicating they are suitable for photovoltaic and photocatalytic applications. Unlike the traditional 2D materials such as MoS2, the optical dielectric properties manifest highly in-plane anisotropic in the infrared and visible region, which is suitable for on-chip polarization manipulation with these materials. This work promotes the understanding of flexible and anisotropic response of these materials and their potential applications in new types of optoelectronic devices.

2

ACS Paragon Plus Environment

Page 3 of 31

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

1. INTRODUCTION

Flexible devices are important as wearable and consumer optoelectronics, and will afford humankind more flexibility in daily life. As such, displays1, thin-film transistors2, touch sensors3, photovoltaic devices4-5, and energy-storage devices6 have recently caughtmuch attention in the flexible format. Two-dimensional (2D) materials, especially 2D transition metal dichalcogenides (TMDs), have been paid wide attentions for their excellent performance in transparent and flexible electronics due to their high crystallinity in atomic thickness and their potential for different types of electronics and optoelectronics. To achieve flexibility, the materials must comply with bending to some degree without losing their function. The hunt for flexible materials suggests that Young’s modulus and shear modulus of them should be relatively small, while bulk modulus of them should be large. It has been reported that large scale molybdenum disulfide (MoS2) film is more suitable for transparent and flexible electronic compared to conventional amorphous silicon or organic films due to its high optical transparency, high electrical mobility and high on/off ratio7. Transparent and high flexible MoS2/hexagonal boron nitride/graphene heterostructures have been utilized in display logic circuits with high mobility and low power consumption8. The flexible and transparent thin film transistors and photodetector have also been achieved based on large-area tungsten selenide (WSe2) film9-10, which display high current on/off ratio. Anisotropic response is a physical property of crystal orientation dependence with different properties along different axes. This usually happens in one dimensional materials as the 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 31

electrons or excitons have different mobility along the longitude and axial axes. Taking single-walled carbon nanotubes as an example, their anisotropic dynamic dielectric response can be used to control THz polarization11, which show great application in THz polarized devices12. 2D materials usually do not show the anisotropic response. However, there are some 2D materials which exhibit anisotropic properties due to their low symmetry. Black phosphorus (BP), as the most stable allotrope of phosphorus, exhibits strong intrinsic in-plane anisotropy with different Hall mobility for holes along the light (x) and heavy (y) effective mass directions. BP thin films also exhibit large and anisotropic in-plane optical conductivity and on/off current ratio13. The calculated electron and hole effective masses show MX3 (M = Ti, Zr; X = S, Se, Te) have anisotropic conductive properties and monolayer TiS3 exhibits high anisotropy in the light absorption14. Unlike general TMDs such as MoS2 and zirconium sulfide (ZrS2), VIIB group TcX2 and ReX2 (X = S, Se) as a new family of 2D TMDs demonstrate an additional valence electron. This additional valence electron inTcX2 and ReX2 (X = S, Se) results in special electronic structures. Using ReS2 as an example, there is no charge couple between layers15, because the charge is confined by the additional valence electron in layers. As a consequence, ReS2 belongs to direct bandgap semiconductor no matter whether in bulk or in monolayer structures15-16. Hence, the special electronic properties and symmetry splitting caused by distorted phase make ReS2 a promising candidate for electron band engineering. The band gaps of TcX2 and ReX2 (X = S, Se) locate in the range of visible-infrared region and they have great applications in photovoltaics17-19. 4

ACS Paragon Plus Environment

Page 5 of 31

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

There are some theoretical calculations of VIIB-VIA 2D TMDs20-24, which mainly emphasize the investigation on the electronic properties, optical absorption, and Raman response. The optical absorption spectra and Raman spectra of monolayer ReX2 (X = S, Se) show in-plane anisotropic response. However, early works are short of an in-depth study on mechanical and optical dielectric properties of TcX2 and ReX2 (X = S, Se). On one hand, flexible electronics requires information of the Young’s modulus, shear modulus, and bulk modulus of the 2D materials. On the other hand, the dielectric functions of materials provide a bridge to connect macroscopic properties and microscopic electronic structures and have a profound effect on the optoelectronic applications of materials. In this paper, the mechanical properties and dielectric functions of VIIB-VIA 2D TMDs are investigated by hybrid density functional. Their elasticity modulus suggests they belong to flexible materials which could be used for flexible devices. Their band gaps mainly depend on d states of transition metal atoms. The band nesting phenomenon occurs in their energy bands manifest the strong light-matter interactions. Highly in-plane anisotropy is observed due to their special electronic structures and crystal structures, through observing the optical dielectric properties of materials.

2. COMPUTATIONAL METHOD

Mechanical, electronic, and dielectric properties of VIIB-VIA 2D TMDs have been performed through Vienna ab initio Simulation Package (VASP)25-26. The Green-Wannier (GW) version of 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 31

the Perdew, Burke, and Ernzerh (PBE) parameterization has been used for the generalized gradient approximation (GGA)27-28. Projector augmented wave (PAW) pseudo-potentials29 act as exchange correlation and electron-ion interaction for all calculations. To obtain reasonable structures, three different kinds of van der Waals correction terms including van der Waals density functional (vdW-DF)30-32, van der Waals density functional dispersion correction (vdW-D3)33, and van der Waals Tkatchenko-Scheffler method with iterative Hirshfeld partitioning (vdW-TS/HI)34-37 are used to describe non-bonding interactions respectively. Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional38 is utilized for the calculation of optical dielectric functions within the framework of PAW method and band gap correction without changing the band character. Frequency dependent dielectric matrix contains the imaginary part of dielectric functions can be obtained after the electronic ground state is determined by hybridization calculations38-39. Afterwards, the real part of dielectric functions could be obtained through Kramers-Kronig relations40. The cut off energy of the plane-wave basis set is 650 eV. The energy relaxation is taken as 1.0 × 10-5 eV and the Hellmann-Feynman force between each atom is set to less than 0.01 eV/Å. In the k-points sampling routine, 5×5×5 and 5×5×1 Monkhorst−Pack grid41 are adopted for bulk and monolayer structures, respectively. In order to eliminate the spurious interaction between neighboring layers, a vacuum layer with a thickness of 3 nm is adopted for the calculation of monolayer structures. All the parameters meet the stability and accuracy criterions.

6

ACS Paragon Plus Environment

Page 7 of 31

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

3. RESULTS AND DISCUSSION

3.1 Crystal Structure and Stability. It is known that the majority of TMDs belong to a hexagonal system with either D6h point-group symmetry or D3d point-group symmetry, while the VIIB-VIA 2D TMDs are triclinic structures with P 1 point-group symmetry. Their stable phases belong to distorted octahedral (T) phase. The chalcogen layers are buckled and the metal atoms form zigzag chains along the direction of lattice vector16. In exception to that, the triclinic unit cells of TcS2, TcSe2, and ReSe2 contain one sandwich which is constituted by four formula units. However, the unit cell of ReS2 contains two sandwiches which are related by a center of symmetry42. This suggests ReS2 has a unit cell with eight formula units. The VIIB atoms locate in a distorted octahedral coordination formed by VIA atoms42-43. Figure 1 illustrates the structures of MX2 and their flexibility.

7

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 31

Figure 1. (a) Side view of MX2. (b) Top view of MX2. (c) MX2 on flexible substrate. The green represent chalcogenide atoms and gray represent transition metal atoms.

The most reasonable crystal structures compared with experimental data are used for the study on mechanical, electronic, and optical dielectric properties. The parameters of relaxed structures are listed in Table S1 in supporting information. The relaxed structures are different in different van der Waals correction terms. From Table S1, it could be found that the relaxed structures of TcS2 from vdW-D3 agree well with the experimental data, while the relaxed structures of ReS2 from vdW-TS/HI are in good agreement with the experimental data. 3.2 Mechanical Properties. The stability of structures under the influence of external 8

ACS Paragon Plus Environment

Page 9 of 31

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

distorting force has a deep influence on the flexible optoelectronics. Therefore, it is important to calculate the elastic constants Cij, which could build the bridge between mechanical and dynamical behavior of structures. In addition, there is a lot of other information, such as the Young’s modulus, Poisson’s ratio etc, could be obtained from the elastic constants. Hence, elastic constants play an important role in determining the stability and stiffness of the materials44-47.

Table 1. Calculated Elastic Constants for VIIB-VIA TMDs; Unit of each Value is GPa

TcS2

C1

C1

C1

C1

C1

C1

C2

C2

C2

C2

C2

C3

C3

1

2

3

4

5

6

2

3

4

5

6

3

4

25

53

36

-6.

12

4.

15

41

7.

58

-0.

56

-8.

29.

-9.

77.

-0.

27.

9

9

4

7

8

7

6

8

1.

14

-3.

11.

-1

65

-0.

22.

5

7

4

8

2

2

8

-5

24

0.3

-0.

-6.

90.

-2

-1.

8

0

7

2

4

5

2

-1.

-2.

-6.

76

-2

-1.

1

4

6

1

6

1 TcS

20

e2

3

ReS

20

2

9

ReS

17

e2

7

4 50

28

-1.

2.8

6 51

29

-5.

0.3

9 42

20

0 43

6.

1 41

9 12

-0.

37

3

2.

-6.

-1

3

3

3

-5.

-1.

-4

19

7.

1.

-7.

-1

8

4

8

7

3

4

3

2

33

11

C35

C3

C44

6

C4

C46

5

C5

C5

C6

5

6

6

39

-1

47

0 29

-1

36

0 19

0.

30

8 17

0.

25

7

As triclinic structures, VIIB-VIA 2D TMDs have 21 independent elastic constants. Table 1 shows the elastic constants of MX2. The eigenvalues of elastic constants matrixes are all positive and meet the stability criterion48 so that the structures would be stable under mechanical distortion. The elastic constants C11, C22, and C33 could indicate the resistance to linear compression from of a-, b-, and c-direction (Figure 1), respectively. Moreover, they also reflect sound propagation along these directions. Hence, it could be found that all VIIB-VIA 2D TMDs have large resistance to linear compression along a-and b-direction, which are larger than that of c-direction. In other words, the materials are stiffer for strain along the a- and b-direction than

9

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 10 of 31

that along the c-direction. This is the reason there exists strong ionic bonds in a- and b-direction, while only weak interaction in c-direction. Besides, the C11 of TcS2 and TcSe2 is bigger than C22 of TcS2 and TcSe2, respectively. While C11 of ReS2 and ReSe2 is lower than C22 of ReS2 and ReSe2, respectively. These suggest that TcX2 (X = S, Se) are easy to slip along the b-direction, while ReX2 (X = S, Se) are easy to slip along the a-direction. The elastic constant C44 plays the most important role in governing the indentation hardness of a material. This parameter reflects the ability of resisting the monoclinic shear distortion in (1 0 0) plane, and the parameter C66 presents the resistance to shear in the direction. It could be found that the C44 of TcS2 and ReS2 is higher than that of TcSe2 and ReSe2, respectively. Besides, the C66 of TcS2, TcSe2, ReS2 and ReSe2 decreased in turn. The resistance of VIIB-VIA 2D TMDs to linear compression along basis vector of lattices is similar to that of MoS245. While the resistance of VIIB-VIA 2D TMDs to shear in direction is lower than that of MoS2, it suggests that VIIB-VIA 2D TMDs are more likely to produce elastic deformation than MoS2. It is known that Young’s modulus (Y), bulk modulus (B), shear modulus (G) and Poisson’s ratio ( ν ), which are regarded as the most concern elastic properties for applications, especially for the investigation on the hardness of polycrystalline materials. The computational process is presented in supporting information. The values of Y, B, G and ν of VIIB-VIA 2D TMDs are listed in Table 2.

10

ACS Paragon Plus Environment

Page 11 of 31

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

Table 2. Calculated Bulk Modulus (B in GPa), Shear Modulus (G in GPa), Young’s Modulus (Y in GPa) and Poisson’s Ratio ( ν ) for VIIB-VIA TMDs from Single Crystal Elastic Constants ν B G Y B/G TcS2 50.591 39.309 93.804 0.193 1.287 TcSe2 50.355 30.510 76.151 0.248 1.650 ReS2 50.736 36.974 89.243 0.207 1.372 33.667 28.524 66.727 0.170 1.180 ReSe2 Young’s modulus is utilized for the measurement of the stiffness of materials, which is defined as the ratio of stress and strain. Compared with other 2D materials, the Young’s modulus of VIIB-VIA 2D TMDs is far lower than that of MoS2 (199.525 GPa)45, but similar to that of BP (70.3 GPa) 49 and bismuth sulfide (Bi2S3) (106.9GPa)50. It has been reported that MoS2 could be used in transparent and flexible electronic devices7-8. Therefore, the VIIB-VIA 2D TMDs, which are more flexible than MoS2, also have potential applications in flexible optoelectronics. Poisson’s ratio is used to reflect the degree of directionality of the covalent bonds. For the covalent materials, the Poisson’s ratio is relatively small ( ν = 0.1) and the typical value of ν for ionic materials is 0.2551. The Poisson’s ratio of VIIB-VIA 2D TMDs is about 0.2, which is lower than those of BP and Bi2S349-50. It means that ionic contribution plays a major role in inter-atomic bonding for VIIB-VIA 2D TMDs. The bulk modulus (B) as the resistance to fracture and shear modulus (G) represents the resistance to plastic deformation. Pugh used the ratio of B and G to measure the ductility52. A material is brittle if B/G is less than 1.75. Contrarily, the material is ductile. 3.3 Electronic Properties. In order to describe the optical dielectric properties accurately, the accurate description of electronic properties of structures should first be obtained. The electronic 11

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 12 of 31

band structures of monolayer VIIB-VIA 2D TMDs are calculated in Brillouin zones (BZ) along the lines connecting high-symmetry points Γ, X, S, X1, Y, X2 , Y1, Γ in the k-space, whose coordination are shown in Figure S1 in supporting information. Figure 2 shows the band structures of monolayer VIIB-VIA 2D TMDs calculated by PBE method.

Figure 2. Band structures of (a) monolayer TcS2, (b) monolayer TcSe2, (c) monolayer ReS2, (d)

monolayer ReSe2. The Fermi level is indicated as the dash line at E=0.0 eV. 12

ACS Paragon Plus Environment

Page 13 of 31

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

It can be found that the monolayer TcS2, TcSe2, and ReSe2 are all indirect gap semiconductors. The CBM of them are between Γ and X of high symmetry k-points and their VBM are between X1 and Y of high symmetry k-points. The ReS2 exhibits a direct gap at Γ. These results are consistent with previous works22-24. Then HSE06 hybrid exchange-correlation (XC) functional is adopted for the correction of band gaps. Table 3 shows the band gaps of monolayer VIIB-VIA TMDs calculated by PBE, PBE with spin-orbit coupling (SOC) and HSE06. Table 3. Calculated Band Gaps in eV for Monolayer VIIB-VIA TMDs PBE PBE-SOC TcS2 1.34 1.28 TcSe2 1.22 1.02 1.46 1.36 ReS2 ReSe2 1.28 1.18

HSE06 1.87 1.70 2.12 1.78

As shown in Table 3, the band gaps of monolayer materials decrease a little bit when SOC is included. The reason is that spin-orbit splitting occur in valence band. While their band gaps increase significantly when the band gaps are corrected by HSE06 hybrid functional. These characteristics are similar to other TMDs such as MoS2, ZrS2 and so on53-55. Besides, the values of band gaps we obtained by HSE06 hybrid functional are similar to previous references24, 56. On account of the fact that band gap defines the threshold of photon absorption, which is related to the electron transition from the valence band to the conduction band. However, the largest absorption usually does not appear at band-gap edge. Therefore, in order to learn about the significant optical response, band nesting of materials is investigated. In 2D materials, the band nesting means the divergence of the joint density of states (JDOS), which implies high optical conductivity57.

According

to

the

definition,

band

nesting

would

appear

once

13

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

∇k ( Ec − E v )

Page 14 of 31

1 ( 2π a ) ( 2π a is the modulus of the reciprocal lattice vector). The lines

connecting high-symmetry points of all structures are searched for exploring the scope where the band nesting occurs. Figure 3 illustrates the modulus of gradient of difference between conduction band and valence band for VIIB-VIA TMDs. The lowest unoccupied band, the highest occupied band, and the second highest occupied band of materials are symbolized by Ec1, Ev1, and Ev2, respectively.

14

ACS Paragon Plus Environment

Page 15 of 31

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

Figure 3. Modulus of gradient of difference for monolayer (a) TcS2, (b) TcSe2, (c) ReS2 and (d) ReSe2. The scopes where the band nesting occurs are summarized in Table S2 in supporting information. It could be found that the band nesting mainly occurs between X and Y1 of the high symmetry k-points. Because of the difference of materials, scopes are slightly different. It means there exist singularities of JDOS in the scopes that embody strong light-matter interactions and high optical conductivity in monolayer structures. Besides, the total and partial densities of states of monolayer VIIB-VIA TMDs are calculated and the results are displayed in Figure 4 and Figure 5 for the understanding of the contribution of different atoms to the band structures of materials.

15

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 16 of 31

Figure 4. (a, b) DOS of TcS2 and TcSe2; (c, d) DOS of S atoms and Se atoms of monolayer TcS2 and TcSe2; (e, f) DOS of Tc atoms of monolayer TcS2 and TcSe2. It could be found that the highest occupied valance bands and the lowest unoccupied conduction bands are mainly dominated by Tc d states for monolayer TcS2 and TcSe2 in Figure 4.

16

ACS Paragon Plus Environment

Page 17 of 31

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

Figure 5. (a, b) DOS of ReS2 and ReSe2; (c, d) DOS of S atoms and Se atoms of monolayer ReS2 and ReSe2; (e, f) DOS of Tc atoms of monolayer ReS2 and ReSe2. In Figure 5, the obvious difference for monolayer ReS2 and ReSe2 is that the highest occupied valance bands are dominated by p states of chalcogenide atoms and d states of Re atoms. The lowest unoccupied conduction bands are dominated by d states of transition metal atoms, which are similar to TcS2 and TcSe2. Compared with Figure 4, the electronic properties of VIIB-VIA 17

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 18 of 31

TMDs could provide information for the manipulating electron band.

3.4 Optical Dielectric Properties. Optical dielectric properties of monolayer VIIB-VIA TMDs are relatively less reported in literatures. However, they are important for the investigation on other optical properties of materials such as absorption spectrum, photoconductivity and so on. Figure 6 illustrates the imaginary part of dielectric functions of monolayer VIIB-VIA TMDs. Because of the in-plane anisotropy, the imaginary part of dielectric functions are indicated by

ε 2 ⊥ XY , ε 2 ⊥ X and ε 2 ⊥ Y , repectively. They indicate imaginary part of dielectric function perpendicular to the plane (the plane defined by the X-axis and Y-axis and as shown in Figure S1 in the supporting information), X-axis, and Y-axis respectively. The peaks of dielectric function are related to the Van Hove singularity of the high symmetry k-points and high symmetry lines, which are listed in Table S3 in supporting information.

18

ACS Paragon Plus Environment

Page 19 of 31

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

Figure 6. Imaginary part of optical dielectric functions of monolayer (a) TcS2, (b) TcSe2, (c) ReS2 and (d) ReSe2. Compared with the imaginary part of optical dielectric functions of monolayer MoS258, imaginary part of optical dielectric functions of monolayer VIIB-VIA TMDs shows obvious in-plane anisotropy after 2.5 eV. What is noteworthy is that TcS2 and TcSe2 tend to be in-plane isotropic in the relatively high energy range (> 5 eV), while ReS2 and ReSe2 always show 19

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 20 of 31

in-plane anisotropic. Table S3 in the supporting information shows the transition of each peak modes for MX2. In monolayer MoS2, the band transition is from valence bands of S p states to conduction bands of Mo d states. While for TcS2 and TcSe2, the transition is from valence bands of Tc d states to conduction bands of Tc d states. However, for ReS2 and ReSe2, the transition is from valence bands of Re d states and p states of chalcogenide atoms to conduction bands of Re d states. Figure 6 also shows that VIIB-VIA TMDs have excellent absorption properties in visible range. Hence, they could be regarded as promising candidates for the third generation solar cell. Through Kramers-Kronig relations40, the real part of dielectric functions could be obtained and the results are displayed in Figure 7.

20

ACS Paragon Plus Environment

Page 21 of 31

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

Figure 7. Real part of optical dielectric functions of monolayer (a) TcS2, (b) TcSe2, (c) ReS2 and (d) ReSe2. The real part of optical dielectric functions of monolayer VIIB-VIA TMDs shows relatively highly in-plane isotropic except the energy range between 2 eV and 3 eV. The complex dielectric constant can be expressed as ε ⊥ (ω) = ε1⊥ (ω ) + i ε2⊥ (ω ) and ε (ω) = ε1 (ω) + i ε2 (ω) , from which absorption, conductivity, transmission, and reflection could be obtained. Compared with 21

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 22 of 31

the bulk materials, the dielectric functions difference is bigger than those of lithium niobate and kalium niobate with the photon energy below 4 eV59. The anisotropy suggests they could be used for the on-chip wave plate, polarizer, and other polarization sensitive devices.

4. CONCLUSION

In summary, we study the structural, electronic, mechanical and optical dielectric properties of the VIIB-VIA 2D TMDs utilizing the ab initio approach. As one kind of new 2D materials that contains in-plane anisotropy, VIIB-VIA 2D TMDs is much more flexible than most studied common TMDs such as MoS2, which have been applied in flexible devices. Therefore, they could be regarded as ideal candidates for flexible electronics. In addition, the band gaps and band nesting of VIIB-VIA monolayer TMDs indicate that they could be the ideal materials for third photovoltaic applications and photo-catalysis. Besides, through the analysis of total and partial densities of states, it could provide guidance for tuning electronic properties of materials. Finally, the optical dielectric functions are obtained. Their highly in-plane anisotropic in the infrared and visible region suggests they could be utilized for on-chip polarization manipulation. The results pave the way for the flexible and anisotropic optoelectronic devices with VIIB-VIA 2D TMDs.

AUTHOR INFORMATION Corresponding Author * Corresponding author: [email protected]. 22

ACS Paragon Plus Environment

Page 23 of 31

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

Supporting Information Lattice constants of relaxed structures, computational formulas of elasticity modulus, Brillouin zones with high symmetry k-points, band nesting regions, transitions and peak positions in dielectric functions spectra

ACKNOWLEDGMENTS This work was supported by National Natural Science Foundation of China (No. 11374240, 11774288, 61605160), Natural Science Foundation of Shaanxi Province (2017KCT-01, 2016JQ1010), Natural Science Research Plan of Shaanxi Education Department (16JK1781), Young Talent Plan from Institute of Science and Technology of University in Shaanxi Province (20160114).

REFERENCES

(1) Choi, G.-M.; Jin, J.; Shin, D.; Kim, Y. H.; Ko, J.-H.; Im, H.-G.; Jang, J.; Jang, D.; Bae, B.-S. Flexible Hard Coating: Glass-Like Wear Resistant, yet Plastic-Like Compliant, Transparent Protective Coating for Foldable Displays. Adv. Mater. 2017, 29, 1700205. (2) Petti, L.; Faber, H.; Münzenrieder, N.; Cantarella, G.; Patsalas, P. A.; Tröster, G.; Anthopoulos, T. D. Low-Temperature Spray-Deposited Indium Oxide for Flexible Thin-Film Transistors and Integrated Circuits. Appl. Phys. Lett. 2015, 106, 092105. (3) Zha, J.-W.; Huang, W.; Wang, S.-J.; Zhang, D.-L.; Li, R. K. Y.; Dang, Z.-M. Difunctional Graphene–Fe3O4 Hybrid Nanosheet/Polydimethylsiloxane Nanocomposites with High Positive Piezoresistive and Superparamagnetism Properties as Flexible Touch Sensors. Adv. 23

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 24 of 31

Mater. Interfaces 2016, 3, 1500418. (4) Leung, S.-F.; Gu, L.; Zhang, Q.; Tsui, K.-H.; Shieh, J.-M.; Shen, C.-H.; Hsiao, T.-H.; Hsu, C.-H.; Lu, L.; Li, D. Roll-to-Roll Fabrication of Large Scale and Regular Arrays of Three-Dimensional Nanospikes for High Efficiency and Flexible Photovoltaics. Sci. Rep.

2014, 4, 4243. (5) Garner, S.; Glaesemann, S.; Li, X. Ultra-Slim Flexible Glass for Roll-to-Roll Electronic Device Fabrication. Appl. Phys. A 2014, 116, 403-407. (6) Yuan, L.; Yao, B.; Hu, B.; Huo, K.; Chen, W.; Zhou, J. Polypyrrole-Coated Paper for Flexible Solid-State Energy Storage. Energy Environ. Sci. 2013, 6, 470-476. (7) Park, J.; Choudhary, N.; Smith, J.; Lee, G.; Kim, M.; Choi, W. Thickness Modulated MoS2 Grown by Chemical Vapor Deposition for Transparent and Flexible Electronic Devices. Appl.

Phys. Lett. 2015, 106, 012104. (8) Lee, G.-H.; Yu, Y.-J.; Cui, X.; Petrone, N.; Lee, C.-H.; Choi, M. S.; Lee, D.-Y.; Lee, C.; Yoo, W. J.; Watanabe, K.; et al. Flexible and Transparent MoS2 Field-Effect Transistors on Hexagonal Boron Nitride-Graphene Heterostructures. ACS Nano 2013, 7, 7931-7936. (9) Zheng, Z.; Zhang T.; Yao, J,; Zhang, Y.; Xu, J.; Yang, G. Flexible, Transparent and Ultra-Broadband Photodetector Based on Large-Area WSe2 Film for Wearable Devices.

Nanotechnology 2016, 27, 225501. (10) Das, S.; Gulotty, R.; Sumant, A. V.; Roelofs, A. All Two-Dimensional, Flexible, Transparent, and Thinnest Thin Film Transistor. Nano Lett. 2014, 14, 2861-2866. (11) Xu, X. L.; Parkinson, P.; Chuang, K. C.; Johnston, M. B.; Nicholas, R. J.; Herz, L. M. 24

ACS Paragon Plus Environment

Page 25 of 31

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

Dynamic Terahertz Polarization in Single-Walled Carbon Nanotubes. Phys. Rev. B: Condens.

Matter Mater. Physic. 2010, 82, 085441. (12) Xu, X.; Yao, Z.; Jin, Y. Texture and Light-Induced Anisotropic Terahertz Properties of Free-Standing Single-Walled Carbon Nanotube Films with Random Networks. Mater. Chem.

Phys. 2015, 162, 743-747. (13) Xia, F.; Wang, H.; Jia, Y. Rediscovering Black Phosphorus as an Anisotropic Layered Material for Optoelectronics and Electronics. Nat. Commun. 2014, 5, 4458. (14) Jin, Y.; Li, X.; Yang, J. Single Layer of MX3 (M= Ti, Zr; X= S, Se, Te): A New Platform for Nano-Electronics and Optics. Phys. Chem. Chem. Phys. 2015, 17, 18665-18669. (15) Tongay, S.; Sahin, H.; Ko, C.; Luce, A.; Fan, W.; Liu, K.; Zhou, J.; Huang, Y.-S.; Ho, C.-H.; Yan, J.; et al. Monolayer Behaviour in Bulk ReS2 Due to Electronic and Vibrational Decoupling. Nat. Commun. 2014, 5, 3252. (16) Fujita, T.; Ito, Y.; Tan, Y.; Yamaguchi, H.; Hojo, D.; Hirata, A.; Voiry, D.; Chhowalla, M.; Chen, M. Chemically Exfoliated ReS2 Nanosheets. Nanoscale 2014, 6, 12458-12462. (17) Corbet, C. M.; McClellan, C.; Rai, A.; Sonde, S. S.; Tutuc, E.; Banerjee, S. K. Field Effect Transistors with Current Saturation and Voltage Gain in Ultrathin ReS2. ACS Nano 2015, 9, 363-370. (18) Liu, E.; Fu, Y.; Wang, Y.; Feng, Y.; Liu, H.; Wan, X.; Zhou, W.; Wang, B.; Shao, L.; Ho, C.-H.; et al. Integrated Digital Inverters Based on Two-Dimensional Anisotropic ReS2 Field-Effect Transistors. Nat. Commun. 2015, 6, 6991. (19) Wolverson, D.; Hart, L. S. Lattice Dynamics of the Rhenium and Technetium 25

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 26 of 31

Dichalcogenides. Nanoscale Res. Lett. 2016, 11, 250. (20) Abdulsalam, M.; Joubert, D. Structural, Electronic and Optical Properties of TcX2 (X = S, Se, Te) from First Principles Calculations. Comput. Mater. Sci. 2016, 115, 177-183. (21) Yu, Z. G.; Cai, Y.; Zhang, Y.-W. Robust Direct Bandgap Characteristics of One- and Two-Dimensional ReS2. Sci. Rep. 2015, 5, 13783. (22) Zhong, H.-X.; Gao, S.; Shi, J.-J.; Yang, L. Quasiparticle Band Gaps, Excitonic Effects, and Anisotropic Optical Properties of the Monolayer Distorted 1T Diamond-Chain Structures ReS2 and ReSe2. Phys. Rev. B: Condens. Matter Mater. Physic. 2015, 92, 115438. (23) Wolverson, D.; Crampin, S.; Kazemi, A. S.; Ilie, A.; Bending, S. J. Raman Spectra of Monolayer, Few-Layer, and Bulk ReSe2: An Anisotropic Layered Semiconductor. ACS Nano

2014, 8, 11154-11164. (24) Jiao, Y.; Zhou, L.; Ma, F.; Gao, G.; Kou, L.; Bell, J.; Sanvito, S.; Du, A. Predicting Single-Layer Technetium Dichalcogenides (TcX2, X = S, Se) with Promising Applications in Photovoltaics and Photocatalysis. ACS Appl. Mater. Interfaces 2016, 8, 5385-5392. (25) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15-50. (26) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Physic.

1996, 54, 11169-11186. (27) Ernzerhof,

M.;

Scuseria,

G.

E.

Assessment

of

the

Perdew–Burke–Ernzerhof

Exchange-Correlation Functional. J. Chem. Phys. 1999, 110, 5029-5036. 26

ACS Paragon Plus Environment

Page 27 of 31

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

(28) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At. Mol. Opt. Phys. 1988, 38, 3098-3100. (29) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater.

Physic. 1994, 50, 17953-17979. (30) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (31) Klimeš, J.; Bowler, D. R.; Michaelides, A. Van Der Waals Density Functionals Applied to Solids. Phys. Rev. B: Condens. Matter Mater. Physic. 2011, 83, 195131. (32) Román-Pérez, G.; Soler, J. M. Efficient Implementation of a Van Der Waals Density Functional: Application to Double-Wall Carbon Nanotubes. Phys. Rev. Lett. 2009, 103, 096102. (33) 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. (34) Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data. Phys. Rev. Lett. 2009, 102, 073005. (35) Kerber, T.; Sierka, M.; Sauer, J. Application of Semiempirical Long-Range Dispersion Corrections to Periodic Systems in Density Functional Theory. J. Comput. Chem. 2008, 29, 2088-2097. (36) Bučko, T.; Lebègue, S.; Hafner, J.; Ángyán, J. G. Improved Density Dependent Correction 27

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 28 of 31

for the Description of London Dispersion Forces. J. Chem. Theory Comput. 2013, 9, 4293-4299. (37) Bučko, T.; Lebègue, S.; Ángyán, J. G.; Hafner, J. Extending the Applicability of the Tkatchenko-Scheffler Dispersion Correction Via Iterative Hirshfeld Partitioning. J. Chem.

Phys. 2014, 141, 034114. (38) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207-8215. (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.

Physic. 2006, 73, 045112. (40) Lucarini, V.; Saarinen, J. J.; Peiponen, K.-E.; Vartiainen, E. M. Kramers-Kronig Relations in

Optical Materials Research; Springer Science & Business Media, 2005; Vol. 110. (41) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Physical

Review B 1976, 13, 5188-5192. (42) Lamfers, H. J.; Meetsma, A.; Wiegers, G. A.; de Boer, J. L. The Crystal Structure of Some Rhenium and Technetium Dichalcogenides. J. Alloys Compd. 1996, 241, 34-39. (43) Wildervanck, J. C.; Jellinek, F. The Dichalcogenides of Technetium and Rhenium. J.

Less-Common Met. 1971, 24, 73-81. (44) Deligoz, E.; Ozisik, H.; Colakoglu, K.; Surucu, G.; Ciftci, Y. O. Mechanical and Phonon Properties of the Superhard LuB2, LuB4, and LuB12 Compounds. J. Alloys Compd. 2011, 509, 1711-1715. 28

ACS Paragon Plus Environment

Page 29 of 31

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

(45) Wei, L.; Jun-fang, C.; Qinyu, H.; Teng, W. Electronic and Elastic Properties of MoS2. Phys.

B 2010, 405, 2498-2502. (46) Ravindran, P.; Fast, L.; Korzhavyi, P. A.; Johansson, B.; Wills, J.; Eriksson, O. Density Functional Theory for Calculation of Elastic Properties of Orthorhombic Crystals: Application to TiSi2. J. Appl. Phys. 1998, 84, 4891-4904. (47) Zhao, Q.; Guo, Y.; Si, K.; Ren, Z.; Bai, J.; Xu, X. Elastic, Electronic, and Dielectric Properties of Bulk and Monolayer ZrS2, ZrSe2, HfS2, HfSe2 from Van Der Waals Density‐ Functional Theory. Phys. Status Solidi B 2017, 1700033. (48) Mouhat, F.; Coudert, F.-X. Necessary and Sufficient Elastic Stability Conditions in Various Crystal Systems. Phys. Rev. B: Condens. Matter Mater. Physic. 2014, 90, 224104. (49) Appalakondaiah, S.; Vaitheeswaran, G.; Lebègue, S.; Christensen, N. E.; Svane, A. Effect of Van Der Waals Interactions on the Structural and Elastic Properties of Black Phosphorus.

Phys. Rev. B: Condens. Matter Mater. Physic. 2012, 86, 035105. (50) Koc, H.; Ozisik, H.; Deligöz, E.; Mamedov, A. M.; Ozbay, E. Mechanical, Electronic, and Optical Properties of Bi2S3 and Bi2Se3 Compounds: First Principle Investigations. J. Mol.

Model. 2014, 20, 2180. (51) Bannikov, V.; Shein, I.; Ivanovskii, A. Electronic Structure, Chemical Bonding and Elastic Properties of the First Thorium‐Containing Nitride Perovskite TaThN3. Phys. Status Solidi

RRL 2007, 1, 89-91. (52) Pugh, S., Xcii. Relations between the Elastic Moduli and the Plastic Properties of Polycrystalline Pure Metals. The London, Edinburgh, and Dublin Philosophical Magazine 29

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 30 of 31

and Journal of Science 1954, 45, 823-843. (53) Kang, J.; Tongay, S.; Zhou, J.; Li, J.; Wu, J. Band Offsets and Heterostructures of Two-Dimensional Semiconductors. Appl. Phys. Lett.2013, 102, 012111. (54) Li, Y.; Kang, J.; Li, J. Indirect-to-Direct Band Gap Transition of the ZrS2 Monolayer by Strain: First-Principles Calculations. RSC Adv. 2014, 4, 7396-7401. (55) Kang, J.; Sahin, H.; Peeters, F. M. Mechanical Properties of Monolayer Sulphides: A Comparative Study between MoS2, HfS2 and TiS3. Phys. Chem. Chem. Phys. 2015, 17, 27742-27749. (56) Liu, H.; Xu, B.; Liu, J. M.; Yin, J.; Miao, F.; Duan, C.-G.; Wan, X. G. Highly Efficient and Ultrastable Visible-Light Photocatalytic Water Splitting over Res2. Phys. Chem. Chem. Phys.

2016, 18, 14222-14227. (57) Carvalho, A.; Ribeiro, R.; Neto, A. C. Band Nesting and the Optical Response of Two-Dimensional Semiconducting Transition Metal Dichalcogenides. Phys. Rev. B:

Condens. Matter Mater. Physic. 2013, 88, 115205. (58) Kumar, A.; Ahluwalia, P. K. A First Principle Comparative Study of Electronic and Optical Properties of 1H – MoS2 and 2H – MoS2. Mater. Chem. Phys. 2012, 135, 755-761. (59) Wiesendanger, E.; Güntherodt, G. Optical Anisotropy of LiNbO3 and KNbO3 in the Interband Transition Region. Solid State Commun.1974, 14, 303-306.

30

ACS Paragon Plus Environment

Page 31 of 31

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

31

ACS Paragon Plus Environment