Promising Photocatalysts for Water Splitting in BeN2 and MgN2

of BeN2 and MgN2 monolayers as efficient photocatalysts for water spitting. Our phonon spectra and ab initio molecular dynamics simulations provide el...
1 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Article

Promising Photocatalysts for Water Splitting in BeN2 and MgN2 Monolayers Yining Wei, Yandong Ma, Wei Wei, Mengmeng Li, Baibiao Huang, and Ying Dai J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b01081 • Publication Date (Web): 27 Mar 2018 Downloaded from http://pubs.acs.org on March 27, 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.

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

Promising Photocatalysts for Water Splitting in BeN2 and MgN2 Monolayers Yining Wei, Yandong Ma*, Wei Wei, Mengmeng Li, Baibiao Huang, Ying Dai* School of physics, State Key Laboratory of Crystal Materials, Shandong University, 250100 Jinan, PR China

Corresponding authors: [email protected] (Y.M.); [email protected] (Y.D.)

ABSTRACT Currently, identifying suitable photocatalysts, especially in two-dimensional materials, for photocatalytic water splitting is still challenging. Here, we report the identification of BeN2 and MgN2 monolayers as efficient photocatalysts for water spitting. Our phonon spectra and ab initio molecular dynamics simulations provide eloquent examinations for the dynamical and thermal stabilities of these two monolayers. We find that BeN2 monolayer exhibits a direct band gap of 2.26 eV, while MgN2 monolayer shows an indirect band gap of 2.33 eV. The carrier mobility for BeN2 monolayer is up to 105 cm2V-1s-1 and 104 cm2V-1s-1 for MgN2 monolayer. More importantly, the BeN2 monolayer has appropriate band edge alignment with respect to the water reduction and oxidation levels in water splitting, while MgN2 monolayer satisfy the reduction levels in water splitting only. These results imply that BeN2 monolayer can be a potential photocatalyst for water splitting, while MgN2 monolayer can be a potential hydrogen production material.

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

INTRODUCTION Hydrogen fuel produced from photocatalytic water splitting is one of the promising solutions to address the energy and environmental issues.1,2 To maximize the efficiency of solar energy utilization, i) the material should have a suitable band gap between 1.55 and 3.0 eV, ii) its band edge positions must straddle the reduction and oxidation potentials of water, iii) its conduction band minimum (CBM) should be higher than the hydrogen reduction potential, while the valence band maximum(VBM) should be lower than the water oxidation potential.3-5 Currently, it has been one of the glorious missions of material science to seek for desirable photacatalyst with proper band gap and band alignments. Two-dimensional (2D) materials have aroused the upsurge of research since graphene was first synthesized in 2004.6-15 Many 2D materials have been identified as promising photocatalyst. Compared with traditional 3D photocatalyst, such as TiO2 and Ag-based photocatalyst, 2D materials have its unique superiority: high specific surface areas, good crystallinity, rich options of host–guest species, better charge carrier separation, and abundant surface active sites.16-22 The typical example is C3N4 monolayer, which has been identified as an attractive photacatalyst both theoretically and experimentally.23 Other 2D materials, including GeX (X = S, Se), p-SiC, TiN, and phosphorene also are proven to be suitable photacatalysts in experiment or theory.24-27 Although significant progress has been made in the study of 2D materials for the application in photocatalysis field, new 2D photocatalyst materials are still highly desirable. What’s more, the oxidation reaction of the water splitting is more complicated, and it suffers from a higher overpotential than the reduction reaction.23, 28

So it is of significance to find suitable photocatalyst with lower oxidation potential,

to make sure the water oxidation reaction will happen. Beryllium has been found to be an ideal ligand to stabilize molecules with planar hypercoordinate carbon, such as Cal4Be, CBe5.29, 30 Other stable beryllium-based 2D materials, such as Be2C, Be5C2, BeC and BeN2 monolayer, have been proposed with various properties.31-34 In this work, using first-principles calculations, we identify that BeN2 and MgN2 monolayer are predicted to be stable, and BeN2 can be an ACS Paragon Plus Environment

Page 2 of 19

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

efficient photocatalyst with excellent performances, and MgN2 monolayer is a potential materials to generate hydrogen. Our cohesive energy, phonon spectra calculations and ab initio molecular dynamics (AIMD) simulations provide compelling evidences for energy, thermal and dynamical stabilities of these two monolayers. To firmly demonstrate their photocatalytic performance, we provide the redox capabilities, optical absorption and carrier’s mobility. Also, we investigate the related properties of the double-layers of BeN2 and MgN2.

THEORETICAL METHODS First-principles calculations were performed by using the density functional theory (DFT) calculations in conjunction with the projector augmented wave (PAW) scheme, as implemented in the plane-wave basis code VASP (Vienna ab initio simulation package).35,36 The generalized gradient approximation (GGA) as formulated by Perdew–Burke–Ernzerhof

(PBE)

was

used

for

exchange

and

correlation

contributions.37 Band structures are also examined by the Heyd–Scuseria–Ernzerhof (HSE06) hybrid functional to give more accurate results.38 A kinetic energy cutoff of 630 eV was chosen for the plane-wave expansion of wave functions and the Monkhorst–Pack scheme of k-point sampling was adopted for the integration over the first Brillouin zone.39 A 10 × 10 × 1 grid for k-point sampling was used for geometry optimization, while 15 × 15 × 1 for the static total energy calculations. A vacuum space of 20 Å is applied as periodic boundary conditions along the z direction, normal to the interface, which is large enough to avoid the interactions between the adjacent layers. The convergence criteria of energy and force were set to 10-5 eV and 0.01 eV/Å, respectively. The DFT-D2 semi-empirical dispersion correction approach was employed to correct the van der Waals (vdW) interactions between the layers when considering the double layers.40 Besides, ab initio MD simulations (AIMD) calculations were carried out to examine thermal stability by using 3 × 3 supercells at 800K within each time step of 1 fs.

3. RESULTS AND DISCUSSION 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

As shown in Figure. 1a, the optimized structures of BeN2 and MgN2 monolayers adopt a flat but distorted hexagonal structure. Both monolayers possess a P-62m symmetry, with two metal atoms (Be, Mg) and four N atoms in one primitive cell. The optimized lattice parameters and bond lengths are shown in Table 1. The lattice constants of BeN2 and MgN2 monolayers are 4.54Å and 5.26Å, respectively. The bond length between Be atom and N atom is 1.62 Å, while the bond length between Mg atom and N atom is 2.0 Å for MgN2 monolayer. The bond length between N atoms is ~1.33 Å for both systems. We then investigate the chemical bonding features by analyzing their electron localized functions (ELF). From Figure. 1b, we can see that substantial concentration of electrons are located between Be (Mg) and N atoms, indicating strong covalent bonding between metal atoms and N atoms.

Table 1. The optimized lattice parameters and bond lengths for BeN2 and MgN2 monolayers. Lattice parameter (Å) BeN2 MgN2

4.54 5.26

Be-N bond length(Å) 1.62 2.0

N-N bond length(Å) 1.32 1.33

ACS Paragon Plus Environment

cohesive energy(eV/atom) 4.98 3.89

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

Figure 1. (a) Crystal structures and (b) electron localization functions (ELFs) for BeN2 (left) and MgN2 (right) monolayers.

The cohesive energies for BeN2 and MgN2 monolayers are 4.98eV/atom and 3.89eV/atom, respectively. Both value is larger than that for black phosphorus (3.48eV/atom),34 implying that they both are more stable than black Phosphorus. Since few-layer black phosphorus has been synthesized experimentally, BeN2 and MgN2 monolayers are expected to be synthesized in experiment in the near feature.41 To confirm their dynamical and thermal stabilities, we perform the phonon spectrum and AIMD simulations. The AIMD simulation is carried out after running 5000 steps at the temperature of 800K, there are neither structure reconstruction nor bond broken, suggesting that the structures of BeN2 and MgN2 monolayers are stable even at a temperature of up to 800 K. From the phonon spectrum, we can see that there are no imaginary vibration modes in the entire Brillouin zones for both BeN2 and MgN2 monolayers, which verifies the dynamical stabilities of BeN2 and MgN2 monolayers. The evolution of free energy for BeN2 and MgN2 monolayers during AIMD simulation at 800 K and the phonon spectrum are shown in the Supporting Information (Figure S1). After checking the dynamical and thermal stabilities, we further examine the mechanical properties of BeN2 and MgN2 monolayers. The in-plane elastic stiffness coefficients for 2D materials can be written as: C11 =

1 ∂ 2U 1 ∂ 2U 1 ∂ 2U = C = C 22 12 2 A0 ∂ε 112 A0 ∂ε 22 A0 ∂ε 11∂ε 22 ,

where A0 is the equilibrium unit cell area, U is the strain energy, ε11 and ε22 are the strain along the x and y directions. ε is set from −0.02 to 0.02 for each direction, with an increment of 0.01. C11 and C22 represent the elastic stiffness along x and y directions, respectively. The coefficients C11, C22 and C12 are calculated using fully relaxed atomic configurations. The in-plane elastic stiffness coefficients can be obtained by fitting the strain energy U to the in-plane strain states (ε11, ε22):

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

1 1 U = C11ε112 + C22ε 222 + C12ε11ε 22 . 2 2 What is more, the Young’s modulus Ex (Ey) is defined as:

Ex = C11 −

C122 C2 E y = C22 − 12 C22 C11 ,

these results related to the C11, C22, C12 and Young’s modulus are listed in Table 2. C11 is almost the same with C22 for both monolayers, resulting from the isotropy geometric structures along the two directions. C11 for BeN2 monolayer is larger than that for MgN2 monolayer, which means the strain can be applied on MgN2 monolayer more easily than BeN2 monolayer. The Young’s modulus of BeN2 and MgN2 are 166 N/m and 92.1 N/m with isotropy, respectively. From the above analyzation we can get that MgN2 monolayer is more flexible than BeN2 monolayer, and both BeN2 and MgN2 monolayers are more flexible than graphene.42, 43

Table 2. In-plane elastic stiffness C11, C22 and C12 and in-plane Young modulus of BeN2 and MgN2 monolayers. The Young modulus of MoS2 monolayer and graphene are also listed for reference. Materials

C11(N/m)

C22(N/m)

C12(N/m)

BeN2 MgN2 graphene MoS2

180.7 112.5

179.6 112.5

51.0 47.9

in-plane Young modulus(N/m) 166.0 92.1 340 123

Next, we investigate the electronic properties of BeN2 and MgN2 monolayers, the corresponding band structures and density of states are plotted in Figure 2. BeN2 is a direct semiconductor with a band gap of 2.26eV, with both the valence band maximum (VBM) and the conduction band minimum (CBM) locating at the Γ point. Different from BeN2, MgN2 presents an indirect-gap band structure, and its VBM located at the K point and CBM lies on the Γ point. From the PDOS and partial charge densities shown in Figure. 2a, we can find that for BeN2 monolayer, the

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

electronic states of the valence band near the Fermi level are mainly dominated by the edge N atoms, while conduction band close to the Fermi level mainly originates from the N center atoms. And for MgN2 monolayer, the VBM is contributed by both N edge atoms and Mg atoms, the CBM originates from the N center atoms. These results agree well with the partial charge density analyze plotted in inset of Figure. 2a and

Figure. 2c. Analysis of the partial DOS, as shown in Figure. 2b and Figure. 2d reveals that the VBM mainly originates from N-pz states for BeN2 monolayer and N-px states for MgN2, while N-pz orbits contribute to the CBM for bot h BeN2 and MgN2 monolayers. We also investigate the double layer structures of BeN2 and MgN2, the geometry and electrical properties are shown in the Supporting Information (Figure S2).

Figure 2. Band structures and projected density of states (PDOS) of (a) BeN2 and (c) MgN2 monolayers, the CBM and VBM of BeN2 and MgN2 monolayers are presented as insets in the PDOS. (b) (d) The orbit contribution PDOS.

Effective electron–hole separation and high charge mobility are two crucial aspects for the application of 2D materials. Herein, we study the effective mass and carrier mobility to identify the carriers transfer properties along specific directions. By fitting

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

parabolic functions to the VBM and CBM of the band structure, we investigate the effective masses of electrons and holes based on the following formula:

d 2 Ek −1 m = ±h ( 2 ) dk , *

2

where k is the wave vector, and Ek is the energy corresponding to the wave vector k. As for mh*, there are two fold degenerate states: the heavier holes correspond to the less dispersive band, the lighter holes are related to the more dispersive band, which means there are two effective masses for holes. We only consider the heavy hole during the

photocatalysis process.44 The carrier mobility in 2D materials is

dominated by acoustic phonon scattering via intra- and inter-valley deformation potential (DP) couplings at room temperature, proposed by Bardeen and Shockley:45, 46

µ=

2eh3C 2

3 K BT m* Ed 2

,

where e, KB, and T represent the electron charge, Boltzmann constant, and temperature, respectively. Ed is the DP constant defined as the shift of the band edges (Eedge) (CBM for electrons and VBM for holes) induced by strain ( ε ): dEedge / d ε . And C =  ∂ 2 E / ∂ε 2  / S 0 is the elastic modulus of 2D system activated by strain, where E is the total energy and S0 is the area of the system. The isotropic strain is defined as ε = ∆c / c0 , where c0 is the lattice constants of the freestanding monolayer. The calculated relative energy and the CBM and VBM positions as a function of the uniaxial strain ε along x and y directions for BeN2 monolayer are shown in the

Supporting Information (Figure S3). It is shown that the linear fittings of relative position of CBM and VBM and the parabolic fitting of total energies fit well with theoretical calculations. The carrier effective mass and mobility of BeN2 and MgN2 along x and y directions are summarized in Table 3. The electron mobility at the CBM along x direction is about 3.24 × 105 cm2V-1s-1, while the holes mobility at the VBM along x direction is about 1.14 × 104 cm2V-1s-1 and 4.48 × 105 cm2V-1s-1 for two

ACS Paragon Plus Environment

Page 8 of 19

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

different degenerate edges. However, for the photocatalysis process, we consider the heavy hole only. The value is higher than most of the 2D materials, including MoS2 (200 cm2V-1s-1), phosphorene (104 cm2V-1s-1) and GeS (2500 cm2V-1s-1).24, 47, 48 For MgN2 monolayer, the mobility along x direction is 1.41 × 104 cm2V-1s-1 for electrons and 172 cm2V-1s-1 for holes. The ultrahigh carriers mobilities suggest fast migration of photo-generated electrons and holes to the surface of the materials to participate in the redox reaction of water splitting. Moreover, the tremendous difference between electron mobility and holes mobility is excellent for the separation of carriers, and could further reduce the probability of the recombination of photo-generated carriers, significantly improving the photocatalytic activity. Furthermore, we investigate the effect of isotropic strains on their electronic structures and band gaps, as shown in

Figure 3. For both BeN2 and MgN2, by exerting isotropic strains ranging from -5% to 5%, there is no change of the position of the CBM and VBM, and the band gap exhibit very little change, which means the electronic properties are not sensitive to the strain.

Table 3. Carrier effective masses (m*) and carrier mobilities of BeN2 and MgN2 monolayers along x and y directions. BeN2

MgN2

Carrier type Electrons (x) Heavy holes (x) Lighter holes(x) Electrons (y) Heavy holes (y) Light holes(y)

m∗/m0 0.27 1.87 0.30 0.26 1.82 0.29

µ(cm2V-1s-1) 3.24 × 105 1.14 × 104 4.48 × 105 3.23 × 105 1.49 × 104 5.88 × 105

Electrons (x) Holes (x) Electrons (y) Holes (y)

0.52 1.78 0.51 1.61

1.41 × 104 1.72 × 102 3.00 × 104 2.62 × 102

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

Figure 3. Band structures of BeN2 monolayer under (a) −5%, (b) -3%, (c) 3% and (d) 5% biaxial strain. And band structures of MgN2 monolayer under (d) −5%, (e) -3%, (f) 3% and (h) 5% biaxial strain. The Fermi level is set to 0 eV.

In addition to the band gap values, another criterion to facilitate the water splitting reaction is that the position of CBM must be higher (more negative) than the hydrogen reduction potential of H+/H2 and the VBM must be lower (more positive) than the water oxidation potential of H2O/O2. We then study the band edges by the HSE functional with respect to the vacuum level. As shown in Figure 4, for BeN2

ACS Paragon Plus Environment

Page 10 of 19

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

monolayer, CBM and VBM satisfy the reduction and oxidation potential of water splitting, respectively. Remarkably, the position of the VBM is -6.284eV, which is much lower than the water oxidation potential of H2O/O2 (-5.67eV). It has been reported that the oxidation reaction has a much higher overpotential, so the lower VBM is essential to the oxidation reaction.28 For MgN2 monolayer, only the CBM position is higher than the reduction level of H+/H2O, which means MgN2 monolayers can sever as potential photoacatalyst to generate hydrogen. The double layers of BeN2 can meet the oxidation potential of water splitting, for MgN2, the reduction potential is satisfied.

Figure 4. The locations of VBM and CBM of BeN2 and MgN2 with respect to the vacuum level which is labeled as 0 eV. The positions of the reduction level of H+ to H2 and the oxidation potential of H2O to O2 are indicated by the dashed lines.

For efficient and feasible photocatalysts in practical applications, the efficient utilization of solar energy is of great importance. The optical absorption properties are analyzed based on the imaginary part of the dielectric function. From Figure 5, the 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

absorption edges for both BeN2 monolayer and MgN2 monolayer lie in the visible light region, and the substantial absorption is at ultraviolet area. Since the visible light accounts for more than 50% of solar energy, the optical absorption of BeN2 and MgN2 monolayers would guarantee high efficiency in the utilization of solar energy. Besides, we find that, for both BeN2 and MgN2 double layers the absorption of light have red shift, which is due to the smaller band gap for the double layers.

Figure 5. Imaginary part of the dielectric function of BeN2, MgN2 monolayers and double layers.

CONCLUSION In summary, we theoretically investigate the structural, electronic and optical properties of 2D BeN2 and MgN2 monolayers based on the first-principles calculations. The dynamical stabilities of the two monolayers are confirmed by phonon spectra with no imaginary frequencies. AIMD simulations provide further evidence for their thermal stability, which can withstand temperature even at 800 K ACS Paragon Plus Environment

Page 12 of 19

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

with tiny lattice distortion and bond broken. We reveal that BeN2 and MgN2 monolayers are semiconductor with band gap of 2.26eV and 2.33eV, respectively. Attractively, both monolayers have extremely high carrier mobility: for BeN2 it reaches 105 cm2V-1s-1 and for MgN2, 104 cm2V-1s-1. Remarkably, the band edges of BeN2 monolayers can meet the requirement of the reduction and oxidation levels in water splitting, while MgN2 monolayers can satisfy the reduction levels. Meanwhile, BeN2 and MgN2 monolayers also exhibit good optical absorption in the visible-light wavelengths. Our study shows that the BeN2 monolayer is a promising candidate as a visible-light water splitting photocatalyst, while MgN2 monolayer is a potential hydrogen production material. However, the insight reaction between BeN2, MgN2 with water, and the stability of the two monolayers in aqueous solution are still unknown, which is needed to be solved in the future.

Supporting Information AIMD simulation and Phonon dispersion curves for BeN2 and MgN2 monolayers, optimized configure and band structure for BeN2 and MgN2 double layers, rectangle unit cell, the variation of band edge and energy as a function of the uniaxial strain along x directions for BeN2

ACKNOWLEDGMENTS This work is supported by the National Basic Research Program of China (973 program, under Grant No. 2013CB632401), the National Natural Science foundation of China (under Grant No.21333006, 11404187and 11374190), the Taishan Scholar Program of Shandong Province, the 111 Project B13029, and Qilu Young Scholar Program of Shandong University.

REFERENCES

(1) Iwase, A.; Ng, Y. H.; Ishiguro, Y.; Kudo, A.; Amal, R. Reduced graphene oxide as a solid-state electron mediator in Z-scheme photocatalytic water splitting under

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

visible light. J. Am. Chem. Soc. 2011, 133, 11054-11057. (2) Kudo, A.; Miseki, Y. Heterogeneous photocatalyst materials for water splitting. Chem. Soc. Rev. 2009, 38, 253-278. (3) Gai, Y.; Li, S.-S.; Xia, J.-B.; Wei, S.-H. Design of narrow-gap TiO2: A passivated codoping approach for enhanced photoelectrochemical activity. Phys. Rev. Lett. 2009, 102, 036402. (4) Bai, Y.; Zhang, Q.; Luo, G.; Bu, Y.; Zhu, L.; Fan, L.; Wang, B. GaS0.5Te0.5 monolayer as an efficient water splitting photocatalyst. Phys. Chem. Chem. Phys.

2017, 19, 15394-15402. (5) Guo, Z.; Zhou, J.; Zhu, L.; Sun, Z. MXene: A promising photocatalyst for water splitting. J. Mater. Chem. A 2016, 4, 11446-11452. (6) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric field effect in atomically thin carbon films. Science 2004, 306, 666-669. (7) Avouris, P. Graphene: Electronic and photonic properties and devices. Nano Lett.

2010, 10, 4285-4294. (8) Pla, J. J.; Tan, K. Y.; Dehollain, J. P.; Lim, W. H.; Morton, J. J.; Zwanenburg, F. A.; Jamieson, D. N.; Dzurak, A. S.; Morello, A. High-fidelity readout and control of a nuclear spin qubit in silicon. Nature 2013, 496, 334-338. (9) Lee, J. H.; Lee, E. K.; Joo, W. J.; Jang, Y.; Kim, B. S.; Lim, J. Y.; Choi, S. H.; Ahn, S. J.; Ahn, J. R.; Park, M. H.; et al. Wafer-scale growth of single-crystal monolayer graphene on reusable hydrogen-terminated germanium. Science 2014, 344, 286−289.

ACS Paragon Plus Environment

Page 14 of 19

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

(10) Radisavljevic, B.; Kis, A. Mobility engineering and a metal-insulator transition in monolayer MoS2. Nat. Mater. 2013, 12, 815-820. (11) Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically thin MoS2: A new direct-gap semiconductor. Phys. Rev. Lett. 2010, 105, 136805. (12) Chhowalla, M.; Liu, Z.; Zhang, H. Two-dimensional transition metal dichalcogenide (TMD) nanosheets. Chem. Soc. Rev. 2015, 44, 2584-2586. (13) Rahman, M. Z.; Kwong, C. W.; Davey, K.; Qiao, S. Z. 2D phosphorene as a water splitting photocatalyst: Fundamentals to applications. Energy Environ. Sci. 2016, 9, 709-728. (14) Naguib, M.; Mashtalir, O.; Carle, J.; Presser, V.; Lu, J.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W. Two-dimensional transition metal carbides. ACS Nano 2012, 6, 1322−1331. (15) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black phosphorus field-effect transistors. Nat. Nanotechnol. 2014, 9, 372-377. (16) Fujishima, A. Electrochemical photolysis of water at a semiconductor electrode. Nature 1972, 238, 37−38. (17) Bard, A. J. Photoelectrochemistry and heterogeneous photo-catalysis at semiconductors. J. Photochem. 1979, 10, 59-75. (18) Bard, A. J. Photoelectrochemistry. Science, 1980, 207, 139-144. (19) Jing, T.; Dai, Y.; Ma, X.; Wei, W.; Huang, B. Electronic structure and photocatalytic water-splitting properties of Ag2ZnSn(S1–xSex)4. J. Phys. Chem. C 2015, 119, 27900-27908.

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

(20) Wang, P.; Huang, B.; Qin, X.; Zhang, X.; Dai, Y.; Wei, J.; Whangbo, M. H. Ag@AgCl: A highly efficient and stable photocatalyst active under visible light. Angew. Chem. Int. Ed. Engl. 2008, 47, 7931-7933. (21) Luo, B.; Liu, G.; Wang, L. Recent advances in 2D materials for photocatalysis. Nanoscale 2016, 8, 6904-6920 (22)

Zhou,

M.;

Lou,

X.

W.; Xie,

Y. Two-dimensional

nanosheets for

photoelectrochemical water splitting: Possibilities and opportunities. Nano Today

2013, 8, 598-618. (23) 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 (24) Lv, X.; Wei, W.; Sun, Q.; Li, F.; Huang, B.; Dai, Y. Two-dimensional germanium monochalcogenides for photocatalytic water splitting with high carrier mobility. Appl. Catal. B 2017, 217, 275-284. (25) Li, X.; Dai, Y.; Li, M.; Wei, W.; Huang, B. Stable Si-based pentagonal monolayers: High carrier mobilities and applications in photocatalytic water splitting. J. Mater. Chem. A 2015, 3, 24055-24063. (26) Sa, B.; Li, Y.-L.; Qi, J.; Ahuja, R.; Sun, Z. Strain engineering for phosphorene: The potential application as a photocatalyst. J. Phys. Chem. C 2014, 118, 26560-26568. (27) Zhou, L.; Zhuo, Z.; Kou, L.; Du, A.; Tretiak, S. Computational dissection of two-dimensional rectangular titanium mononitride TiN: Auxetics and promises for

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

photocatalysis. Nano Lett. 2017, 17, 4466-4472. (28) Wirth, J.; Neumann, R.; Antonietti, M.; Saalfrank, P. Adsorption and photocatalytic splitting of water on graphitic carbon nitride: A combined first principles and semiempirical study. Phys. Chem. Chem. Phys. 2014, 16, 15917-15926. (29) Jimenez-Halla, J. O.; Wu, Y. B.; Wang, Z. X.; Islas, R.; Heine, T.; Merino, G. CAl4Be and CAl3Be2: Global minima with a planar pentacoordinate carbon atom. Chem. Commun. 2010, 46, 8776-8778. (30) Guo, J. C.; Ren, G. M.; Miao, C. Q.; Tian, W. J.; Wu, Y. B.; Wang, X. CBe5Hnn-4 (n = 2-5): Hydrogen-stabilized CBe5 pentagons containing planar or quasi-planar pentacoordinate carbons. J. Phys. Chem. A 2015, 119, 13101-13106. (31) Wang, Y.; Li, F.; Li, Y.; Chen, Z. Semi-metallic Be5C2 monolayer global minimum with quasi-planar pentacoordinate carbons and negative Poisson's ratio. Nat. Commun. 2016, 7, 11488. (32) Li, Y.; Liao, Y.; Chen, Z. Be2C monolayer with quasi-planar hexacoordinate carbons: A global minimum structure. Angew. Chem. Int. Ed. Engl. 2014, 53, 7248-7252. (33) Liu, C. S.; Zhu, H. H.; Ye, X. J.; Yan, X. H. Prediction of a new BeC monolayer with perfectly planar tetracoordinate carbons. Nanoscale 2017, 9, 5854-5858. (34) Zhang, C.; Sun, Q. A honeycomb BeN2 sheet with a desirable direct band gap and high carrier mobility. J. Phys. Chem. Lett. 2016, 7, 2664-2670. (35) 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,

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

15. (36) Kresse, G.; Hafner, J. Norm-conserving and ultrasoft pseudopotentials for first-row and transition elements. J. Phys.: Condens. Matter 1994, 6, 8245−8257. (37) Venkataraman, L.; Klare, J. E.; Nuckolls, C.; Hybertsen, M. S.; Steigerwald, M. L. Dependence of single-molecule junction conductance on molecular conformation. Nature 2006, 442, 904−907. (38) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid functionals based on a screened coulomb potential. J. Chem. Phys. 2003, 118, 8207-8215. (39) Monkhorst, H. J.; Pack, J.D. Special points for brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188-5192. (40) 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. (41) Castellanos-Gomez, A.; Vicarelli, L.; Prada, E.; Island, J. O.; Narasimha-Acharya, K. L.; Blanter, S. I.; Groenendijk, D. J.; Buscema, M.; Steele, G. A.; Alvarez, J. V.; Zandbergen, H. W.; Palacios, J. J.; van der Zant, H. S. J. Isolation and characterization of few-layer black phosphorus. 2D Mater. 2014, 1, 025001. (42) Lee, C.; Wei, X.D.; Kysar, J. W.; Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science, 2008, 321, 385-388. (43) Yue, Q.; Kang, J.; Shao, Z.; Zhang, X.; Chang, S.; Wang, G.; Qin, S.; Li, J. Mechanical and electronic properties of monolayer MoS2 under elastic strain. Phy. Lett. A 2012, 376, 1166-1170.

ACS Paragon Plus Environment

Page 18 of 19

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

(44) Ma, X.; Dai, Y.; Guo, M.; Huang, B. The role of effective mass of carrier in the photocatalytic behavior of silver halide-based Ag@AgX (X=Cl, Br, I): A theoretical study. ChemPhysChem 2012, 13, 2304−2309. (45) Kaasbjerg, K.; Jauho, A. P.; Thygesen, K. S. Acoustic phonon-limited mobility in two-dimensional MoS2: Deformation potential and piezoelectric scattering from first principles. Phys. Rev. B 2013, 87, 235312−235326. (46) 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. (47) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-layer MoS2 transistors. Nat. Nanotechnol. 2011, 6, 147-150. (48) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tomanek, D.; Ye, P. D. Phosphorene: An unexplored 2D semiconductor with a high hole mobility. ACS Nano,

2014, 8, 4033-4041.

TOC Graphic

ACS Paragon Plus Environment