Janus Chromium Dichalcogenide Monolayers with Low Carrier

As expected, the intrinsic dipole moment exists for all Janus CrXY monolayers arising from the broken mirror symmetry, which generates a built-in elec...
0 downloads 0 Views 814KB Size
Subscriber access provided by TULANE UNIVERSITY

C: Surfaces, Interfaces, Porous Materials, and Catalysis

Janus Chromium Dichalcogenides Monolayers with Low Carrier Recombination for Photocatalytic Overall Water-Splitting under Infrared Light Pei Zhao, Yan Liang, Yandong Ma, Baibiao Huang, and Ying Dai J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b11240 • Publication Date (Web): 29 Jan 2019 Downloaded from http://pubs.acs.org on February 5, 2019

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 19 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

Janus Chromium Dichalcogenides Monolayers with Low Carrier Recombination for Photocatalytic Overall Water-Splitting under Infrared Light Pei Zhao, Yan Liang, Yandong Ma*, Baibiao Huang and Ying Dai* School of Physics, State Key Laboratory of Crystal Materials, Shandong University, Shandanan Str. 27, Jinan 250100, People's Republic of China *Corresponding author: [email protected] (Y. M); [email protected] (Y. D)

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 19

ABSTRACT Photocatalytic overall water-splitting is known as one of most promising methods to alleviate energy crisis. Searching for stable and efficient photocatalysts plays a critical role in this process. Here, we propose a novel class of Janus chromium dichalcogenides (CrXY, X/Y = S, Se, Te) monolayers serving as efficient photocatalysts for overall water-splitting under infrared light irradiation. We reveal that these Janus monolayers harbor an intrinsic dipole, which promotes the spatial separation of photo-generated carriers. More significantly, these systems exhibit suitable band gaps as well as band edge positions, enabling preeminent infrared optical absorption and high carrier mobility. Furthermore, the nonradiative recombination of photoinduced charge carriers in CrXY monolayers are evaluated based on time-domain density functional theory. The obtained long-lived excited carriers (~ 2 ns) are even comparable with that in transition-metal dichalcogenides heterostructures, which benefits for the photocatalytic reaction with high efficiency. Our results provide a new guidance for designing brand new photocatalytic systems with broad optical absorption and low carrier recombination.

2

ACS Paragon Plus Environment

Page 3 of 19 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 Searching for clean and renewable energy is an effective way to solve environmental pollution and energy shortage. Hydrogen production utilizing solar energy based on photocatalytic water-splitting technology has attracted widespread attention during the past decades.1,2 In photocatalytic field, exploring robust photocatalysts with high conversion efficiency from solar energy to hydrogen energy is one of the most urgent issue. To this end, excellent optical absorption, spatial separation of photogenerated carriers and energy level are the prerequisites for achieving high photocatalytic efficiency. Nevertheless, conventional three dimensional (3D) photocatlysts, such as TiO2, SrTiO3, CdS and WO3 etc.,3-9 exhibit a wide band gap, which restricts their absorption in the ultraviolet light region, accounting for only about 5% of solar energy. And most of the reactive sites of these 3D photocatalysts cannot be exposed on the surface. Although the carriers can migrate from the interior to the surface of bulk, such migration would increase the possibility of their recombination.10,11 As a result, most of the 3D photocatalysts display a low photocatalytic efficiency. Despite extensive efforts have been devoted to enhancing their performances in photocatalytic reaction (i.e., band gap engineering, surface modification and developing new reactive mechanism, etc.12-14), the photocatalytic water-splitting efficiency still call for further improvement. Since the birth of graphene,15-19 a wide variety of two-dimensional (2D) materials

have been

applied in photocatalytic field, such as, 2D oxides, 2D chalcogenides and 2D g-C3N4 etc..20-25 Compared with 3D photocatalytic materials, 2D photocatalysts hold high carrier mobility, tunable band gap and light absorption, large specific surface area and short migration distance for the generated electrons and holes.21 Recently, Janus MoSSe monolayer was successfully synthesized via chemical vapour deposition (CVD) method.26 Due to the broken structural symmetry, the Janus structure has an intrinsic dipole moment in the out-of-plane direction, thus inducing a built-in electric field. More recently, Janus MoSSe monolayer has been demonstrated to be a promising 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 19

photocatalyst.27-29 However, its large band gap makes it missing the energy of infrared spectrum accounting for even 43% of solar energy. In this case, it is imperative to hunt for novel photocatalysts to fully utilize solar energy and improve the energy conversion efficiency. Inspired by the unique properties of Janus monolayers,30-32 we design a series of novel Janus CrXY monolayers and explore their potential applications in photocatalytic overall water-splitting based on first principles calculations. Their stability is examined by ab initio molecular dynamics (AIMD) and phonon calculations. Interestingly, all these Janus monolayers exhibit suitable band edge positions relative to redox levels of water, and the intrinsic electrostatic potential difference between the top and bottom surfaces could effectively promote spatial separation of photo-generated carriers. The most remarkable point is that Janus monolayers hold excellent optical absorption under the whole solar spectrum including the near-infrared light. Meanwhile, the detailed time-domain ab initio atomistic simulations prove that the electron-hole recombination rate in CrXY is slow, making these materials bear a high potential for applications in photocatalysis.

2. COMPUTATIONAL METHODS First principles calculations are carried out by using projector-augmented wave (PAW) scheme implemented in Vienna ab-initio simulation package (VASP).33,34 For

geometrical optimizations,

we use the generalized gradient approximation (GGA) as formulated by Predew-Burke-Ernzerhof (PBE) to describe the interaction between exchange and correlation interaction between electrons.35,36 Cut-off energy is set as 500 eV for the plane-wave expansion of wave functions, and a vacuum space of about 15 Å perpendicular to 2D plane is included to avoid the interaction between adjacent layers. A 15×15×1 grid for k-point sampling is adopt. The structures are fully relaxed until the residual force and energy difference are, respectively, less than 0.01 eV/Å and 10-5 eV. It is worth mentioning that HSE functional may give wrong electronic structures of CrXY compounds, and for the sake of safety, the electronic properties of CrXY monolayers are checked not only based on PBE and HSE functional but also based on G0W0 method.37-39 In one shot G0W0 calculation, we set 10-8 eV and 150 eV as the 4

ACS Paragon Plus Environment

Page 5 of 19 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

energy convergence and cutoff energy. And the number of empty band is set as 144 to obtain accurate results. Besides, thermal stability is demonstrated with a 5×5×1 supercell at 300 K based on AIMD, and dynamical stability is examined by 6×6×1 supercell based on PHONOPY code.40,41 The nonadiabatic molecular dynamic simulations (NAMD) for electron-hole recombination with decoherence correction in CrXY are carried out using the PYXAID code.42,43 The methodology and justification are presented elsewhere.44,45 A 4 ps MD simulation at 300 K in the microcanonical ensemble is performed with a time step of 1 fs, then the generated trajectories are used to obtain the nonadiabatic (NA) Hamiltonian. 500 initial configurations are randomly sampled from the first 1ps MD trajectory for NAMD simulations.

3. RESULTS AND DISCUSSION 3.1 Structure and stability of the CrXY monolayer The crystal structure of CrXY monolayer is formed by one Cr layer sandwiched by one X (S, Se, Te) layer and one Y layer (X≠Y; S, Se, Te). As shown in Figure 1(a), CrXY monolayer share similar honeycomb structures with that of pristine CrX2 monolayer. To provide more information on the surface of crystal for experiments, the simulated scanning tunneling microscope (STM) is plotted in Figure 1(b). Due to the atomic difference of top and bottom surfaces, the mirror symmetry of the structure is broken, thus leading to a dipole moment along the out-of-plane direction. The optimized lattice constants of pristine CrX2 and Janus CrXY monolayers are listed in Table I. We can see that the lattice constant of CrSSe monolayer is found to be 3.13 Å, approaching to the average value of those of CrS2 and CrSe2. The similar conclusion is also obtained in CrSeTe and CrSTe monolayers. To confirm their stability, AIMD and phonon calculations are performed. Figure 1(d) and S1 display the top and side views of proposed Janus monolayers after heating for 5ps at temperature of 300K. Clearly, there is no bond broken or structural reconstruction in all these monolayers. And the 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 19

vibrations of total energy per unit cell for these monolayers are within the range of 0.2 eV. The phonon dispersions of Janus CrXY monolayers are shown in Figure 1(c) and S2, no imaginary frequencies appear for all proposed monolayers. In addition, we also calculate their formation energies (Ef) based on the following equation: Ef = (Etotal-nECr-nE1-nE2)/3n, where Etotal is the total energy of monolayer structures, ECr, E1 and E2 are the energies of Cr, X and

Figure 1. (a) Top and side views and (b) Simulated STM image for CrSSe monolayer. (c) Phonon dispersion curves of CrSSe monolayer. (d) Side views and total energy per unit cell of CrSSe snapshotted from MD simulations at 300 K during the time of 5ps. (e) The calculated formation energy of CrX2 and CrXY monolayers. The orange arrow indicates the intrinsic electric-field direction. Yellow, orange and light blue balls represent S, Se and Cr atoms, respectively. Y atoms in their stable bulk phases, n is the number of atoms. From Figure 1(e), we find that all the formation energy of these structures are less than the value of the Sb monolayer (~ 0.9 eV/atom), which has been synthesized in experiment.16 These results could serve as convincing evidences to demonstrate the stability of these proposed Janus structures at room temperature. And they are expected to be synthesized in experiment by CVD method. 6

ACS Paragon Plus Environment

Page 7 of 19 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.2 Electronic properties and band alignment As shown in Figure S3, all 2H-CrX2 monolayers are non-magnetic semiconductors with a direct band gap of 1.34, 1.14 and 0.96 eV for CrS2, CrSe2 and CrTe2 monolayers, which is consistent with previous studies.46-48 Due to the mirror symmetry, there is no electrostatic potential difference between top and bottom surfaces. For both CrSSe and CrSeTe monolayers, they are also nonmagnetic semiconductors with a direct band gap of 1.24 and 1.04 eV, respectively, with valence band maximum (VBM) and conduction band minimum (CBM) locating at the K point. While for CrSTe monolayer, it holds an indirect band gap of 0.69 eV, with VBM and CBM locating at  and K points, respectively. In order to obtain the accurate band gaps of Janus CrXY monolayers, the band structures are examined on G0W0 level. From Figure 3, we can see that the band gaps of CrSSe, CrSeTe and CrSTe monolayers increase to 1.57, 1.32 and 1.10 eV, respectively, while the main features of the bands are similar to those based on HSE06 functional.

Figure 2. Band structures of (a) CrSSe, (b) CrSeTe and (c) CrSTe based on HSE06 (blue line) and G0W0 (red line) functional. EQP is the quasiparticle bandgap and insert picture is the 2D Brillouin zone of the hexagonal structure. The Fermi levels are set as 0 eV. As expected, the intrinsic dipole moment exists for all Janus CrXY monolayers arising from the mirror symmetry broken, which generates a built-in electric field along the out-of-plane direction. As

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 19

displayed in Table I and Figure S4, the dipole moment of CrSSe monolayer is 0.2 D, and a built-in electric field along vertical direction pointing from Se atom to S atom forms, thus inducing electrostatic potential difference of 0.88 eV. Similarly, for CrSeTe and CrSTe, the electrostatic potential difference between the top and bottom surface are 0.98 and 1.82 eV, respectively, which results in a dipole moment of 0.25 and 0.45D. Table I. Optimized lattice constants (a), formation energy Ef, band gaps based on HSE06, band gaps based on G0W0, the quasiparticle and optical band gaps EQP and Eopt, exciton binding energies (Eb), dipole moments (μ) and difference of electrostatic potential (ΔV) between top and bottom surfaces for CrX2 and Janus CrXY monolayers. Materials

a (Å)

Ef(eV/atom)

CrS2 CrSe2 CrTe2 CrSSe CrSeTe CrSTe

3.04 3.21 3.48 3.13 3.35 3.27

-2.13 -1.48 -0.79 -1.79 -1.11 -1.38

EHSE(eV) EGW(eV) EQP(eV) EOpt(eV)

1.34 1.14 0.96 1.24 1.04 0.69

1.57 1.32 1.10

1.57 1.32 1.29

1.15 0.88 1.06

Eb (eV)

μ(debye)

V (eV)

0.42 0.44 0.23

0 0 0 0.20 0.25 0.45

0 0 0 0.88 0.98 1.82

Figure 3. (a) Schematic diagram of absolute positions of band edges with respect to vacuum level according to the band gap center (BGC) theory. Band alignments of (b) CrSSe, (c) CrSeTe and (d) CrSTe monolayers based on G0W0 level. As we know, it is indispensable for fully photocatalytic water-splitting to hold suitable band edge 8

ACS Paragon Plus Environment

Page 9 of 19 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

positions, where CBM is higher than the reduction level of hydrogen and VBM is lower than the oxidation level of oxygen. We initially investigate the band edges of pristine CrX2 monolayers based on HSE06 level. Unfortunately, their band edge positions could not meet the redox potential requirement, reflecting that CrX2 monolayers are not good candidates as being photocatalyts. Then, we examine the band edge positions of Janus CrXY monolayers. According to reaction mechanism proposed by Yang et al.,49 the existence of electrostatic potential difference makes band gap (> 1.23 eV) of photocatalysts no limitation for fully water-splitting. As shown in Figure S5, the proposed Janus monolayers well satisfy the band edge requirements for fully photocatalytic water-splitting. In order to accurately obtain the band edges, GW is a preferred method. However, its well-converged process of band edge is considerately costly. In such situation, approximations that determine the absolute band edge positions but avoid a fully converged GW calculation have been proposed, in which of their validity has yet been verified in TMDs.50-53 As presented in Figure 3(a), the energy of band gap center (EBGC) is obtained by HSE06, Eg is the calculated band gap based on G0W0. Generally, the band edge positions of VBM and CBM can be determined by EVBM = EBGC – 1/2 Eg, ECBM = EBGC +1/2 Eg, which has been verified by previous reports.50,52 As we mentioned above, the CBM and VBM would move to outside under the internal electric field. And the BGC theory tells us that the BGC of one certain material is insensitive to the calculated method. In this case, we redefine the band edge positions (CBM and VBM) as EVBM = EBGC – 1/2 (Eg+ V), ECBM = EBGC +1/2 (Eg + V) for polar systems. As shown in Figure 3(b)-(d), all Janus CrXY monolayers hold suitable band alignments for overall water-splitting, and the energy differences between band edges and redox levels are about 0.5 eV. These results demonstrate that these proposed Janus structures have a good redox ability, indicating that they can be potential photocatalysts for water-splitting. To further explain the redox mechanism for splitting water, we investigate the partial charge density near the VBM and CBM of pristine CrX2 and Janus CrXY monolayers. From Figure S6, we can clearly see that the partial charge density near CBM and VBM of pristine CrX2 monolayers are localized around Cr atoms. While for 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 19

Janus CrSSe monolayer, the partial charge density near CBM mainly distributes around Se and Cr atoms, and that near VBM mainly resides at S and Cr atoms. This can be attributed to the electrostatic potential difference between two surfaces of CrSSe monolayer. In this case, the photo-generated electrons and holes tend to move towards opposite directions, which leads to the spatial separation of carriers. According to photocatalytic reaction: 4H2O+4e-→2H2+4OH-, 2H2O +4h+→O2+4H+, the hydrogen production mainly occurs in the Se surface, while the oxygen production mainly occurs in S surface, thus being able to improve the photocatalytic efficiency. And similar reactive mechanism is also obtained in CrSeTe and CrSTe monolayers. 3.3 Optical absorption and Photoexcitation Dynamics In view of their photocatalytic applications, it is of great significant to investigate the optical absorption behaviors of CrXY monolayers. It is well known that 2D materials exhibit relatively weak screening, we thus study the optical properties of Janus CrXY monolayers using G0W0-BSE method, which takes the electron-hole Coulomb interactions into consideration. The imaginary part of the dependent complex dielectric function is calculated according to:54

α (ω) = 2ω



ε1 (ω) 2 + ε2 (ω) 2 - ε1 (ω)



1/2

As shown in Figure 4(a), all Janus CrXY monolayers exhibit good optical absorption, which are superior to that of the well-known 2D g-C3N4 and MoSSe photocatalysts.28,55 More remarkably, the initial peaks of the optical absorption for CrSSe, CrSeTe and CrSTe monolayers are 1.15, 0.88 and 1.06 eV, respectively, illustrating that these proposed photocatalysts can be active under near-infrared light region. In addition, we study the exciton binding energy of CrXY monolayers on the basis of G0W0-BSE. The exciton binding energy is the difference between quasiparticle band gaps and optical band gaps, the calculated quasiparticle band gaps and exciton energy are displayed in Table I, which 10

ACS Paragon Plus Environment

Page 11 of 19 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

basically follows the Eb ~ 1/4EQP law proposed in previous study. 56,57 And these values are smaller than those of 2D g-C3N4, BP, and MoS2 etc., which is favorable for the separation of photo-generated electrons and holes. Besides, the carrier mobility is another vital factor determining the photocatalytic efficiency. Therefore, the carrier mobility of Janus CrXY monolayers is calculated using the deformation potential method, which is expressed by the following equation:58

2 D 

2eh3C 2

3K BT m* Ed 2

. Here, e and ℏ are electron charge and reduced Planck constant. Kb, T and Ed represent Boltzmann constant, room temperature (300 K) and deformation potential constant. m  h2 (

d 2 Ek 1 ) is the dk 2

carries effective mass, where Ek and k are total energies and momentum, respectively. C  [ 2 E /  2 ] / S0 is the elastic constant of the 2D system, wherein E and S0 are total energy and

system’s area, respectively. The calculated hole mobility for CrSSe monolayer is about 0.37×103 and 0.31×103 𝑐𝑚2𝑉 ―1𝑠 ―1along zigzag (x) and armchair (y) directions, respectively. While for electrons, the mobilities are 0.53×102 and 0.43×102 𝑐𝑚2𝑉 ―1𝑠 ―1 along x and y directions, respectively. The huge disparity of mobility between electron and hole are commonly beneficial for their separation, further enhancing the photocatalytic efficiency. More detailed descriptions are listed in Table SI.

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 19

Figure 4. (a) The optical absorption coefficients of Janus CrXY monolayers based on G0W0-BSE. (b) Decay of the excited-state population at 300 K. (c) Spectral density obtained from Fourier transforms of the autocorrelation functions for the band gap fluctuations of Janus CrXY monolayers. (d) Averaged absolute nonadiabatic (NA) coupling. The photo-generated carrier lifetime plays a crucial role in affecting photocatalytic efficiency.59 To quantitatively describe the lifetime of photo-generated carriers, the evolutions of the population of the first excited state in these Janus CrXY monolayers are calculated using time-domain density functional theory combined with NAMD, as displayed in Figure 4(b). The electron-hole recombination in CrSSe monolayer needs 71 ps, in the same time scale of TMD monolayers (several hundreds of picoseconds).60,61 Emotionally, the lifetimes of photo-generated carriers in CrSeTe and CrSTe monolayers can reach up to 2.4 and 1.9 ns, respectively, even comparable with those of TMDs heterostructures.62,63 This indicates CrXY prone to possess long-lived photo-generated electrons and holes to be involved in the redox reaction, thus leading to high efficiency in photocatalytic overall water-splitting. 12

ACS Paragon Plus Environment

Page 13 of 19 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

The different recombination rate of these systems can be rationalized by the following two factors. First, spectral density by Fourier transforms of the fluctuations for band gaps characterizes the phonon modes of vibrational motions inducing electron-hole recombination. As depicted in Figure 4(c), the dominant vibrational mode peak of CrXY locate at frequency around 200 cm-1, and the mode intensity vanishes when the frequency is larger than 500 cm-1. Particularly, the vibrational modes for CrSSe monolayer is widespread with highest intensities, while are relatively weak for CrSeTe and CrSTe monolayers. Thus, least intense phonon modes in CrSeTe monolayer coupling to the electronic subsystem participate in the relaxation of excited electrons and create longest recombination time.64,65 Second, the recombination rate is proportional to the NA coupling squared between CBM and VBM,66 smaller NA coupling favors slower nonradiative decay as it determines the hopping probability of hot electrons from one state to another. The NA coupling between two electronic states j and k can be expressed as:67

d jk   j |

 j |  R H | k   | k  R, k   j t

so NA coupling is mainly determined by energy difference term  k   j , electron-phonon (e-ph) coupling term

 j |  R H | k



and the nuclear velocity term R . As shown in Figure 4(d), the

strength of NA coupling is in the order of CrSSe > CrSTe > CrSeTe, mainly because the atomic mass is in the order of Te> Se> S, and the e-ph coupling is in the sequence of CrSSe > CrSTe > CrSeTe. In comparison, less active vibrations and small NA coupling cause longest lifetime of photogenerated carriers in CrSeTe.

4. CONCLUSIONS In summary, we propose a series of novel Janus CrXY structures, and demonstrate that all these CrXY monolayers are promising photocatalysts under near-infrared light. The thermal and dynamic stabilities are confirmed. These Janus structurs are all semiconductors with a band gap of 1.57, 1.32 and 1.10 eV for CrSSe, CrSeTe and CrSTe monolayers, respectively. Due to the broken symmetry, 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

Page 14 of 19

Janus CrXY monolayers display an intrinsic electric field. This, makes the reactions of hydrogen and oxygen production happening on different surfaces, thus improving the photocatalytic efficiency. Moreover, all these Janus monolayers hold suitable band edges, which ensures their redox abilities in fully water-splitting. The optical absorption spectrum demonstrates that these Janus monolayers can be active even under near-infrared light, which is of great significance for increasing the utilization of solar energy. Besides, the electron-hole recombination of Janus CrXY monolayers is evaluated as well. All these systems hold long–lived photogenerated carriers, especially for CrSeTe and CrSTe monolayers. Given all this, the proposed Janus monolayers can be promising candidates with high efficiency applied in photocatalytic overall water-splitting.

CONFLICTS OF INTEREST There are no conflicts to declare.

SUPPORTING INFORMATION The carrier mobility of Janus CrXY monolayers, MD simulation and phonon dispersion for CrSeTe and CrSTe monolayers, the band structures of CrX2 monolayer based on HSE level, the electrostatic potential and band alignments of CrX2 and CrXY monolayers based on HSE level, partial charge density of CrX2 and CrXY monolayers.

ACKNOWLEDGEMENT This work is supported by the National Natural Science foundation of China (No. 11804190 and 21333006), Qilu Young Scholar Program of Shandong University, Taishan Scholar Program of Shandong Province, and 111 Project (No. B13029).

REFERENCES (1) Dincer, I. Renewable Energy and Sustainable Development: A Crucial Review. Renew. Sustain. Ener. Rev. 2000, 4, 157-175. (2) Christoforidis, K. C.; Fornasiero, P. Photocatalytic Hydrogen Production: A Rift into the Future Energy Supply. Chem. Cat. Chem. 2017, 9, 1523-1544. (3) Schneider, J.; Matsuoka, M.; Takeuchi, M.; Zhang, J.; Horiuchi, Y.; Anpo, M.; Bahnemann, D. W. Understanding TiO2 Photocatalysis: Mechanisms and Materials. Chem. Rev. 2014, 114, 99199986. 14

ACS Paragon Plus Environment

Page 15 of 19 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

(4) Yang, H.; Li, J.; Yu, L.; Huang, B.; Ma, Y.; Dai, Y. A theoretical Study on the Electronic Properties of In-plane CdS/ZnSe Heterostructures: Type-II Band Alignment for Water Splitting. J. Mater. Chem. A 2018, 6, 4161-4166. (5) Zhang, L. J.; Li, S.; Liu, B. K.; Wang, D. J.; Xie, T. F. Highly Efficient CdS/WO3 Photocatalysts: Z-Scheme Photocatalytic Mechanism for Their Enhanced Photocatalytic H2 Evolution under Visible Light. ACS Catal. 2014, 4, 3724-3729. (6) Fujishima, A.; Honda, K. Electrochemical Photolysis of Water at a Semiconductor Electrode. Nature 1972, 238, 37-38. (7) Zhang, H.; Zhu, Y. F. Significant Visible Photocatalytic and Antiphotocorrosion Performance of CdS Photocatalysts after Monolayer Polyaniline Hybridization. J. Phys. Chem. C 2010, 114, 5822-5826. (8) Wang, D.; Liu, Z. P.; Yang, W. M. Proton-Promoted Electron Transfer in Photocatalysis: Key Step for Photocatalytic Hydrogen Evolution on Metal/Titania Composites. ACS Catal. 2017, 7, 2744-2752. (9) Wang, D.; Liu, Z. P.; Yang, W. M. Revealing the Size Effect of Platinum Cocatalyst for Photocatalytic Hydrogen Evolution on TiO2 Support: A DFT Study. ACS Catal. 2018, 8, 72707278. (10)Hisatomi, T.; Kubota, J.; Domen, K. Recent Advances in Semiconductors for Photocatalytic and Photoelectrochemical Water Splitting. Chem. Soc. Rev. 2014, 43, 7520-7535. (11)Ma, X. C.; Dai, Y.; Yu, L.; Huang, B. B. Energy Transfer in Plasmonic Photocatalytic Composites. Light: Sci. & Appl. 2016, 5, e16017. (12)Chen, S.; Qi, Y.; Hisatomi, T.; Ding, Q.; Asai, T.; Li, Z.; Ma, S. S.; Zhang, F.; Domen, K.; Li, C. Efficient Visible-Light-Driven Z-Scheme Overall Water Splitting Using a MgTa2O(6x)N(y)/TaON Heterostructure Photocatalyst for H2 Evolution. Angew. Chem. Int. Ed. 2015, 54, 8498-8501. (13)Khan, M. M.; Ansari, S. A.; Pradhan, D.; Ansari, M. O.; Han, D. H.; Lee, J.; Cho, M. H. Band Gap Engineered TiO2 Nanoparticles for Visible Light Induced Photoelectrochemical and Photocatalytic Studies. J. Mater. Chem. A 2014, 2, 637-644. (14)Nolan, M.; Iwaszuk, A.; Lucid, A. K.; Carey, J. J.; Fronzi, M. Design of Novel Visible Light Active Photocatalyst Materials: Surface Modified TiO2. Adv. Mater. 2016, 28, 5425-5446. (15)Nair, P. R.; Blake, P.; Grigorenko, A. N.; Novoselov, K. S.; Booth, T. J.; Stauber, T.; Peres, M. R.; Geim, A. K. Fine Structure Constant Defines Visual Transparency of Graphene. Science, 2008, 320, 1308. (16)Ji, J.; Song, X.; Liu, J.; Yan, Z.; Huo, C.; Zhang, S.; Su, M.; Liao, L.; Wang, W.; Ni, Z.; Hao, Y.; Zeng, H. Two-Dimensional Antimonene Single Crystals Grown by van der Waals Epitaxy. Nat. Commun. 2016, 7, 13352. (17)Kou, L.; Ma, Y.; Sun, Z.; Heine, T.; Chen, C. Two-Dimensional Topological Insulators: Progress and Prospects. J. Phys. Chem. Lett. 2017, 8, 1905-1919. (18)Xu, W. P.; Jin, Y. J.; Zhwng, B. B.; Xu, H. Extremely High Mobilities in Two-Dimensional Group-VA Binary Compounds with Large Conversion Efficiency for Solar Cells. J. Phys. Chem. C 2017, 139, 11694-11697. (19)Kou, L. Z.; Ma, Y. D.; Yan, B. H.; Tan, X.; Chen, C. F.; Smith, S. C. Encapsulated Silicene: A Robust Large-Gap Topological Insulator. ACS Appl. Mater. Interfaces 2015, 7, 19226-19233. (20)Singh, A. K.; Mathew, K.; Zhuang, H. L.; Hennig, R. G. Computational Screening of 2D Materials for Photocatalysis. J. Phys. Chem. Lett. 2015, 6, 1087-1098. 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 19

(21)Su, T.; Shao, Q.; Qin, Z.; Guo, Z.; Wu, Z. Role of Interfaces in Two-Dimensional Photocatalyst for Water Splitting. ACS Catal. 2018, 8, 2253-2276. (22)Wang, X.; Maeda, K.; Thomas, A.; Takanabe, K.; Xin, G.; Carlsson, J. M.; Domen, K.; Antonietti, M. A Metal-Free Polymeric Photocatalyst for Hydrogen Production from Water under Visible Light. Nat. Mater. 2009, 8, 76-80. (23)Zhao, P.; Ma, Y. D.; Lv, X. S.; Li, M. M.; Huang, B. B.; Dai, Two-dimensional group In2-VI3: promising photocatalysts for overall water splitting under infrared light spectrum. Nano Energy 2018, 51, 533-538. (24)Gao, G.; Jiao, Y.; Waclawik, E. R.; Du, A. Single Atom (Pd/Pt) Supported on Graphitic Carbon Nitride as an Efficient Photocatalyst for Visible-Light Reduction of Carbon Dioxide. J. Am. Chem. Soc. 2016, 138, 6292-6297. (25)Zhang, X.; Zhao, X.; Wu, D.; Jing, Y.; Zhou, Z. MnPSe3 Monolayer: A Promising 2D VisibleLight Photohydrolytic Catalyst with High Carrier Mobility. Adv. Sci. 2016, 3, 1600062. (26)Lu, A. Y.; Zhu, H.; Xiao, J.; Chu, C. P.; Han, Y.; Chiu, M. H.; Cheng, C. C.; Yang, C. W.; Wei, K. H.; Yang, Y.; Wang, Y.; Sokaras, D.; Nordlund, D.; Yang, P.; Muller, D. A.; Chou, M. Y.; Zhang, X.; Li, L. J. Janus Monolayers of Transition Metal Dichalcogenides. Nat. Nanotech. 2017, 12, 744-749. (27)Ma, X.; Wu, X.; Wang, H.; Wang, Y. A Janus MoSSe Monolayer: A Potential Wide SolarSpectrum Water-Splitting Photocatalyst with A Low Carrier Recombination Rate. J. Mater. Chem. A 2018, 6, 2295-2301. (28)Guan, Z.; Ni, S.; Hu, S. Tunable Electronic and Optical Properties of Monolayer and Multilayer Janus MoSSe as A Photocatalyst for Solar Water Splitting: A First-Principles Study. J. Phys. Chem. C 2018, 122, 6209-6216. (29)Ji, Y.; Yang, M.; Lin, H.; Hou, T.; Wang, L.; Li, Y.; Lee, S. T. Janus Structures of Transition Metal Dichalcogenides as the Heterojunction Photocatalysts for Water Splitting. J. Phys. Chem. C 2018, 122, 3123-3129. (30)Li, F.; Wei, W.; Zhao, P.; Huang, B.; Dai, Y. Electronic and Optical Properties of Pristine and Vertical and Lateral Heterostructures of Janus MoSSe and WSSe. J. Phys. Chem. Lett. 2017, 8, 5959-5965. (31)Tang, X.; Li, S. S.; Ma, Y. D.; Du, A. J.; Liao, T.; Gu, Y. T.; Kou, L. Z. Distorted Janus Transition Metal Dichalcogenides: Stable Two-Dimensional Materials with Sizable Band Gap and Ultrahigh Carrier Mobility. J. Phys. Chem. C 2018, 122, 19153-19160. (32)Jin, H.; Wang, T.; Gong, Z.; Long, C.; Dai, Y. Prediction of Extremely Long Exciton Lifetime in Janus-MoSTe Monolayer. Nanoscale 2018, 10, 19310-19315. (33)Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. (34)Kresse, G. Efficient Interative Schemes for Ab Initio Total-Energy Calculations Using A PlaneWave Basis Set. Phys. Rev. B 1996, 54, 11169. (35)Perdew, J. P.; Wang, Y. Accurate and Simple Analytic Representation of the Electron-Gas Correlation Energy. Phys. Rev. B 1992, 45, 13244. (36)Perdew, J. P.; Bruke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. (37)Paier, J.; Marsman, M.; Hummer, K.; Kresse, G.; Gerber, I. C.; Angyan, J. G. Screened Hhybrid Density Functionals Applied to Solids. J. Chem. Phys. 2006, 124, 154709. (38)Hybertsen, M. S.; Louie, S. G. Electron Correlation in Semiconductors and Insulators: Band Gaps and Quasiparticle Energies. Phys. Rev. B 1986, 34, 5390-5413. 16

ACS Paragon Plus Environment

Page 17 of 19 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

(39)Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. Wannier90: A Tool for Obtaining Maximally-Localized Wannier Functions. Comput. Phys. Commun. 2008, 178, 685-699. (40)Barnett, R. N.; Landman, U. Born-Oppenheimer Molecular-Dynamics Simulations of Finite Systems: Structure and Dynamics of (H2O)2. Phys. Rev. B 1993, 48, 2081-2097. (41)Togo, A.; Tanaka, I. First Principles Phonon Calculations in Materials Science. Scripta Mater. 2015, 108, 1-5. (42)Alexey V. A.; Prezdo, O. V. Advanced Capabilities of the PYXAID Program: Integration Schemes Decoherence Effect, Multiexcitonic States and Field-Matter Interaction. J. Chem. Theory Comput. 2014, 10, 789-804. (43)Alexey V. A.; Prezdo, O. V. The PYXAID Program for Non-Adiabatic Molecular Dynamics in Condensed Matter Systems. J. Chem. Theory Comput. 2013, 9, 4959-4972. (44)He, J. L; Fang, W. H; Long, R; Prezhdo, O. V. Increased Lattice Stiffness Suppresses Nonradiactive Charge Recombination in MAPbI3 Doped with Larger Cations: Time Domain Ab Initio Analysis. ACS Energy Lett. 2018, 3, 2070-2076. (45)Long, R.; Prezhdo, O. V. Ab Initio Nanodiabatic Molecular Dynamics of Ultrafast Electron Injection from a PbSe Quantum Dot into the TiO2 Surface. J. Am. Chem. Soc. 2011, 133, 1924019249. (46)Wang, C.; Zhou, X. Y.; Pan, Y. H.; Qiao, J. S.; Kong, X. H.; Kaun, C. C.; Ji, W. Layer and Doping Tunable Ferromagnetic Order in Two-Dimensional CrS2 Layers. Phys. Rev. B 2018, 97, 245409. (47)Guo, H. Y.; Lu, N.; Wang, L.; Wu, X. J.; Zeng, X. C. Tuning Electronic and Magnetic Properties of Early Transition-Metal Dichalcogenides via Tensile Strain. J. Phys. Chem. C 2014, 118, 72427249. (48)Zhuang, H. L.; Johannes, M. D.; Blonsky, M. N.; Hennig, R. G. Computational Prediction and Characterization of Single-Layer CrS2. Appl. Phys. Lett. 2014, 104, 022116. (49)Li, X.; Li, Z.; Yang, J. Proposed Photosynthesis Method for Producing Hydrogen from Dissociated Water Molecules Using Incident Near-Infrared Light. Phys. Rev. Lett. 2014, 112, 018301. (50)Jiang, H. Electronic Band Structures of Molybdenum and Tungsten Dichalcogenides by the GW Approach. J. Phys. Chem. C 2012, 116, 7664-7671. (51)Liang, Y.; Huang, S.; Soklaski, R.; Yang, L. Quasiparticle Band-Edge Energy and Band Offsets of Monolayer of Molybdenum and Tungsten Chalcogenides. Appl. Phys. Lett. 2013, 103, 042106. (52)Pan, Y.; Wang, Y.; Wang, L.; Zhong, H.; Quhe, R.; Ni, Z.; Ye, M.; Mei, W. N.; Shi, J.; Guo, W.; Yang, J.; Lu, J. Graphdiyne-Metal Contacts and Graphdiyne Transistors. Nanoscale 2015, 7, 2116-2127. (53)Toroker, M. C.; Kanan, D. K.; Alidoust, N.; Isseroff, L. Y.; Liao, P.; Carter, E. A. First Principles Scheme to Evaluate Band Edge Positions in Potential Transition Metal Oxide Photocatalysts and Photoelectrodes. Phys. Chem. Chem. Phys. 2011, 13, 16644-16654. (54)Lalitha, S.; Karazhanov, S. Z.; Ravindran, P.; Senthilarasu, S.; Sathyamoorthy, R.; Janabergenov, J. Electronic Structure, Structural and Optical Properties of Thermally Evaporated CdTe Thin Films. Phys. B: Condensed Matter 2007, 387, 227-238. (55)Wei, W.; Jacob, T. Strong Excitonic Effects in the Optical Properties of Graphitic Carbon Nitride g-C3N4 from First Principles. Phys. Rev. B 2013, 87, 085202. (56)Choi, J. H.; Cui, P.; Lan, H. P.; Zhang, Z. Y. Linear Scaling of the Exciton Binding Energy versus 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 19

the Band Gap of Two-Dimensional Materials. Phys. Rev. Lett. 2015, 115, 066403. (57)Jiang, Z.; Liu, Z.; Li, Y.; Duan, W. Scaling Universality between Band Gap and Exciton Binding Energy of Two-Dimensional Semiconductors. Phys. Rev. Lett. 2017, 118, 266401. (58)Xi, J.; Long, M.; Tang, L.; Wang, D.; Shuai, Z. First-principles Prediction of Charge Mobility in Carbon and Organic Nanomaterials. Nanoscale 2012, 4, 4348-4369. (59)Wei, Y.; Li, L.; Fang, W.; Long, R.; Prezhdo, O. V. Weak Donor-Acceptor Interaction and Interface Polarization Define Photoexcitation Dynamics in the MoS2/TiO2 Composite: TimeDomain Ab Initio Simulation. Nano Lett. 2017, 17, 4038-4046. (60)Korn, T.; Heydrich, S.; Hirmer, M; Schmutzler, J.; Schuller, C. Low-Temperature Photocarrier Dynamics in Monolayer MoS2. Appl. Phys. Lett. 2011, 99, 102109. (61)Sun, D. Z.; Rao, Y.; Reider, G. A.; Chen, G.; You, Y. M.; Brezin, L.; Harutyunyan, A. R.; Heinz, T. F. Observed of Rapid Exciton-Exciton Annihilation in Monolayer Molybdnum Disulfide. Nano Lett. 2014, 14, 5625-5629. (62)Li, L.; Long, R.; Prezhdo, O. V. Charge Separation and Recombination in Two-Dimensional MoS2/WS2: Time-Domain Ab Initio Modeling. Chem. Mater. 2016, 29, 2466-2473. (63)Liang, Y.; Li, J.; Jin, H.; Huang, B.; Dai, Y. Photoexcitation Dynamics in Janus-MoSSe/WSe2 Heterobilayers: Ab Initio Time-Domain Study. J. Phys. Chem. Lett. 2018, 9, 2797-2802. (64)Meggiolaro, D.; Mosconi, E.; De Angelis, F. Mechanism of Reversible Trap Passivation by Molecular Oxygen in Lead-Halide Perovskites. ACS Energy Lett. 2017, 2, 2794-2798. (65)Quarti, C.; Grancini, G.; Mosconi, E.; Bruno, P.; Ball, J.; Lee, M.; Snaith, H.; Petrozza, A.; De Angelis, F. The Raman Spectrum of the CH3NH3PbI3 Hybrid Perovskite: Interplay of Theory and Experiment. J. Phys. Chem. Lett. 2014, 5, 279-284. (66)Li, W.; Long. R.; Hou, Z.; Tang, J.; Prezhdo, O. V. Influence of Encapsulated Water on Luminescence Energy, Line Width and Lifetime of Carbon Nanotubes: Time Domain Ab Initio Analysis. J. Phys. Chem. Lett. 2018, 9, 4006-4013. (67)Reeves, K. G.; Schleife, A.; Correa, A. A.; Kanai, Y. Role of Surface Termination on Hot Electron Relaxation in Silicon Quantum Dots: A First-Principles Dynamics Simulation Study. Nano Lett. 2015, 15, 6429-6433.

TOC Graphic

18

ACS Paragon Plus Environment

Page 19 of 19 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

19

ACS Paragon Plus Environment