Two-dimensional blue-phase CX(X=S,Se)monolayers with high carrier mobility and tunable photocatalytic water splitting capability

Shengi Li,Mengyo Shi,Jihui Yu,Shujun Li,Shuli Lei,*,Lingxu Lin,Jijun Wng,*

a Tianjin Key Laboratory of Structure and Performance for Functional Molecules,Key Laboratory of Inorganic-Organic Hybrid Functional Materials Chemistry,Ministry of Education,College of Chemistry,Tianjin Normal University,Tianjin 300387,China

b Hubei Key Laboratory of Low Dimensional Optoelectronic Materials and Devices,Hubei University of Arts and Science,Xiangyang 441053,China

c ARC Centre of Excellence for Electromaterials Science,Intelligent Polymer Research Institute,Australia Institute for Innovative Materials(AIIM),Innovation Campus,University of Wollongong,Wollongong 2519,Australia

ABSTRACT Photocatalytic water splitting utilizing solar energy is considered as one of the most ideal strategies for solving the energy and environmental issues.Recently,two-dimensional(2D)materials with an intrinsic dipole show great chance to achieve excellent photocatalytic performance.In this work,blue-phase monolayer carbon monochalcogenides(CX,X=S,Se)are constructed and systematically studied as photocatalysts for water splitting by performing first-principles calculations based on density functional theory.After confirming the great dynamical,thermal,and mechanical stability of CX monolayers,we observe that they possess moderate band gaps(2.41 eV for CS and 2.46 eV for CSe)and high carrier mobility(3.23×104 cm2 V-1 s-1 for CS and 4.27×103 cm2 V-1 s-1 for CSe),comparable to those of many recently reported 2D photocatalysts.Moreover,these two monolayer materials are found to have large intrinsic dipole(0.43 D for CS and 0.51 D for CSe),thus the build-in internal electric field can be selfintroduced,which can effectively drive the separation of photongenerated carriers.More importantly,the well-aligned band edge as well as rather pronounced optical absorption in the visible-light and ultraviolet regions further ensure that our proposed CX monolayers can be used as high efficient photocatalysts for water splitting.Additionally,the effects of external strain on the electronic,optical and photocatalytic properties of CX monolayers are also evaluated.These theoretical predictions will stimulate further work to open up the energy-related applications of CX monolayers.

Keywords:First-principle calculations Carbon monochalcogenides Carrier mobility Optical absorption Water splitting

Hydrogen production fromsolar water splitting,which offers a promising technology to solve the global challenge of the depletion of fossil fuels and the associated environmental pollution,has attracted significant interest over the past years[1-4].Since the pioneering work of Fujishima and Honda in 1972[5],numerous inorganic semiconductor materials,including metal oxides,metal sulfides,metal nitrides,and metal phosphides,have been explored for water photolysis[6-12].Unfortunately,most reported inorganic photocatalysts exhibit intrinsic disadvantages such as wide band gap(>3.0 eV),ineffective utilization of solar energy(can only absorb ultraviolet light),and high electron-hole recombination rate induced by the low carrier mobility and long migration distance[13,14].Hence,their photocatalytic performances are still yet to achieve the requirement of practical application.

Compared to the bulk photocatalysts,two-dimensional(2D)photocatalysts possess unique atomic structures and unprecedented properties,thus have been hot topics of extensive experimental and theoretical researches[15-22].For instance,Sun et al.[23]reported that the photocurrent density of singlelayer SnS2reached up to 2.75 mA/cm2at 1.0 V,which was 72 times larger than the bulk counterpart.Cheng et al.[24]found that the hydrogen generation rates of synthesized FePS3sheets were 3 times higher than that of bulk FePS3.Qiao et al.[25]theoretically predicted that PdSeO3monolayer can be easily exfoliated from the bulk phase,and its reduction/oxidation ability was high enough for splitting pure water into hydrogen and oxygen.Most recently,Yang’s group designed a new 2D photocatalyst(F-BNBN-H)with large intrinsic dipole,in which the band gap restriction(1.23 eV)for direct overall water splitting can be broken and thus utilizing the near-infrared light would be achieved[26].Referring to this novel model,some 2D materials with intrinsic electric polarization,such as Janus transition metal dichalcogenides monolayers[27-30],Sc2CO2[31],M2XY(M=Ga,In;X/Y=S,Se,Te)[32-34]and M2X3(M=Al,Ga,In;X=S,Se,Te)[35,36],have been exploited as photocatalyts for water splitting.In virtue of the intrinsic electric field produced by the intrinsic dipole,these 2D materials have not only appropriate band gap and band edge alignments,but also very low recombination possibility of photongenerated carriers,making them be promising candidates as highly efficient photocatalysts.Nevertheless,the amount of such built-in dipoled 2D photocatalysts is still rather rare.It is still necessary to search for other 2D semiconductor materials with good structure stability,appropriate electronic properties,and large intrinsic dipoles simultaneously.

As allotropes of phosphorene,layered group-IV monochalcogenides have been theoretically predicted and experimentally synthesized recently,which show good promise for various applications in these aspects of gas sensors,solar cells,anode materials for ion batteries,and photocatalysts[37-40].In particular,the monolayer materials in the buckled configuration,similar to blue phosphorene,usually generate the intrinsic dipoles due to by the large difference of electronegativity between the top and bottom sides.Monolayer silicon and germanium monochalcogenides(i.e.,SiS,SiSe,GeS and GeSe)have been found to be suitable for overall photocatalytic water splitting[41-44].Inspired by these,we thus wonder whether the monolayer carbon monochalcogenides(CX,X=S,Se)can steadily exist as freestanding sheets and whether they are high efficient photocatalysts for water splitting.Hence,we systematically investigate the structural,electronic,and photocatalytic properties of blue-phase CX(X=S,Se)monolayers by performing first-principles calculations.

First-principles calculations are carried out based on the density functional theory with the projector augmented wave method[45],as implemented in Vienna Ab initio Simulation Package(VASP)[46].The Perdew-Burke-Ernzerhof(PBE)functional[47]is adopted for optimizing the geometric structure,and the Heyd-Scuseria-Ernzerhof(HSE06)functional[48,49]is adopted for giving more accurate electronic and optical properties,in which the 25% PBE exchange part is replaced by the Hartree-Fock exact exchange in the short-range.Namely,the mixing parameter(α)of the HSE06 functional is set to be 0.25.The energy cutoff energy is set to 500 eV,and the Monkhorst-Pack mesh[50]of k-points 21×21×1 points are used to sample the Brillouin zone.A vacuum region of 20.0?along the z direction is applied to avoid the coupling effect of two adjacent image layers.The convergence criterias of total energy and atomic forces are set as 1×10-5eV and 0.02 eV/?,respectively.Dipole correction along the z direction is needed to cancel the errors of electrostatic potential,atomic force,and total energy,caused by periodic boundary condition.

Akin to blue phosphorene,the CX monolayers have a buckled hexagonal structure with the C and S(Se)atoms interlocked and distributed alternately,as shown in Fig.1a.After optimization,the lattice constant of monolayer CS is predicted to be 2.82?with a CS bond length of 1.87?and buckling height of 0.92?,while for monolayer CSe,the lattice parameter,C-Se bond length and buckling height increase to 3.07?,2.06?and 1.04?,respectively.The reason is that the atomic radius of Se atom is larger than that of S atom.To characterize the bonding features in these CX monolayers,the electron localization functions(ELF)are calculated and depicted in Fig.1b and Fig.S1(Supporting information),where the values of 1.0,0.5 and 0.0 represent the fully localized electron,the free electron-gas and very low electron density,respectively.It can be seen that the ELF values around the C and S(Se)sites are different,and a substantial concentration of electrons are localized at the regions between C and S(Se)atoms,revealing the presence of polar covalent bonding in monolayer CS and CSe.The Bader charge analysis further illustrates that 0.23(0.39)electrons are transferred from S(Se)atom to C atom,also confirming the polar chemical bonding.

Fig.1.(a)Geometric structure of CX monolayers from thetop and side views.The hexagonal primitive cell and orthogonal supercell are marked by black and blue lines,respectively.(b)Visualization of the electron localization function for monolayer CS.Calculated phonon dispersions of monolayer(c)CS and(d)CSe.The total energy fluctuations during 10 ps AIMD simulations at room temperature(300 K)for monolayer(e)CS and(f)CSe.Insets in(e)and(f)show the final geometric structures after AIMD simulation.

Before moving on to other properties,we first focus on the structural stability of the CX monolayers,including their dynamical,thermal,and mechanical stability.To assess their dynamical stabilities,the phonon dispersion curves are calculated based on the frozen-phonon method as implemented in the PHONOPY package[51].The 5×5×1 supercell and 4×4×1Γcentered kpoint sampling are used for the phonon calculations.As shown in Figs.1c and d,we find that the low-frequency optical and acoustic modes separate apart from each other around theΓpoint,and the frequencies of all phonon modes in the entire Brillouin zone are positive.The highest frequencies of monolayer CS and CSe reach up to 680 cm-1and 580 cm-1,respectively,which are comparable to those of silicene(550 cm-1)[52]and blue phosephorene(520 cm-1)[53].These results demonstrate that our proposed two monolayer materials are dynamically stable.Furthermore,we examine their thermal stability by performing ab initio moleculardynamics(AIMD)simulations at 300 K using a 5×5×1 supercell.As plotted in Figs.1e and f,the total energies have small fluctuations and the initial geometric structures are well maintained after 10 ps simulation with the NVT ensemble,suggesting that monolayer CS and CSe can be stable at the room-temperature.In addition,we also investigate the effect of the O2molecule on their structural stabilities.As shown in Fig.S2(Supporting information),we can see that the distances from the O2molecule to the CS and CSe monolayers are 2.67 and 2.52?,respectively,and the corresponding binding energy(Ea)are 151 meV and 223 meV.This indicates that O2molecule is bonded to the monolayers mainly by vdW force.As a result,our proposed CX monolayers can be stable in the presence of O2.

Additionally,we systematically examine the elastic properties of the CX monolayers to analyze the mechanical stability.Based on the strain-energy method,the calculated elastic constants of monolayer CS(CSe)are C11=C22=132.64(95.59)N/m,C12=14.69(16.20)N/m and C66=58.97(39.69)N/m.Clearly,the Born-Huang criteria of mechanical stability for 2D hexagonal materials(e.g.,C11C22-C122>0,C66>0)is nicely satisfied,reflecting the mechanical robustness of monolayer CS and CSe.Given these elastic constants,we also calculate the in-plane Young’s modulus(Y)and Poisson’s ratio(υ),which are defined as Y=(C112-C122)/C11and υ=C12/C11,respectively.The in-plane Young’s modulus(Y)and Poisson's ratio(υ)of monolayer CS(CSe)are found to be 131.01(92.84)N/m and 0.11(0.17),respectively.These results are lower than graphene[54]but comparable to MoS2[55],indicate that the bonding between atoms in these two monolayers are strong and would avail to avoid the curling.

Encouraged by the excellent stability of the CX monolayers,we then turn to discuss their electronic properties by computing the band structures as well as the density of states(DOSs).As shown in Fig.2,both of CS and CSe monolayers present indirect band gap semiconductor characteristics.The valence band maximum(VBM)of CS monolayer lies at the point between K andΓpoints and their conduction band minimum(CBM)is situated at the M point.As for CSe monolayer,the VBM is still located on the K→Γpath while the CBM is located between theΓand M points.Although a direct gap semiconductor is preferred,the materials with indirect band gaps can also be used as promising photocatalysts for water splitting because they advantageously reduce the recombination possibility of photogenerated electron-hole pairs.On basis of the HSE06 calculation,the predicted band gap values are 2.41 eV for CS and 2.46 eV for CSe.Noted that these values are lower than the band gap of metal-free graphitic carbon nitride(g-C3N4),which has been experimentally identified as a visible-light-responsive photocatalyst[56,57].From the computed partial DOSs depicted in the Fig.2,we can see that the VBMs of monolayer CS and CSe are mainly contributed by X 3p,C 2 s and C 2p states,while their CBMs are primarily originated from the hybridization of X 3p and C 2p states.

To further explore the charge transport properties for potential energy conversion application,the carrier mobility of CX monolayers are estimated according to the deformation potential(DP)theory proposed by Bardeen and Shockley[58,59].The carrier mobility for 2D materials is described as

Fig.2.HSE06 calculated electronic band structures and projected density of states of monolayer(a)CS and(b)CSe.Here,the Fermi energy levels are scaled to zero and shown as red dashed lines.

Table 1 Carrier effective mass(m*),2D elastic modulus(C2D),deformation potential(Ed),and carrier mobility(μ)of CX monolayers.

Now,we explore the possibility of using CX monolayers as photocatalysts for water splitting.Owing to the electronegativity difference between C and X atoms,the C atoms in the top layer are negatively charged whereas the X atoms in the bottom layer are positively charged.Thus,an intrinsic dipole(P)is generated,which are found to be 0.43 D for CS and 0.51 D for CSe.As we known,these relatively large dipole can introduce internal electric field(Eeff)in monolayer CX,which are along with the perpendicular plane direction and point from X side to C side.Moreover,as described by the Yang group,the energy levels of the monolayer CS and CSe would bend along the direction of the built-in Eeff,accompanied by the presence of the difference(ΔΦ)between the vacuum energy levels on the two respective sides[26].As shown in Figs.3a and b,the top and bottom vacuum regions are different by 2.34 and 2.38 eV for monolayer CS and CSe,respectively,which are as large as some widely known 2D materials with intrinsic electricpolarization[27-36],such as MXene,germanium monochalcogenide and Janus metal chalcogenides.

To assess whether monolayer CX possess strong reduction/oxidation ability,we further calculate their accurate VBM/CBM locations by considering the vacuum level corrections,and then align them relative to the redox potential levels for H+/H2(-4.44 eV at pH 0)and O2/H2O(-5.67 at pH 0),which are illustrated in Figs.3c and d.It should be pointed that the photogenerated electrons and holes can aggregate separately on disparate sides of monolayer CX because of the built-in electric field,meaning that the hydrogen reduction reaction and water oxidation reaction happen on the X and C sides,respectively.On the X atomic layer side,the CBM potentials of CX monolayers(-3.47 eV for CS and-3.12 eV for CSe)lie above hydrogen reduction potential,showing that H+can be reduced to H2by the photogenerated electrons.On the C atomic layer side,the VBMs(-8.21 eV for CS and-7.96 eV for CSe)lie below water oxidation potential,showing that H2O can be oxidized to O2by the photogenerated holes.Following a previous study,the energy difference between the CBM and the hydrogen reduction potential is defined as the potential of photogenerated electrons for hydrogen evolution reaction(Ue),while the energy difference between the VBM and the water oxidation potential is calculated as the potential of photogenerated holes for oxygen evolution reaction(Uh).Hence,the Ueand Uhare calculated to be 0.97/2.54 eV for CS and 1.32/2.29 eV for CSe,suggesting that both of them are capable to work as highly efficient photocatalysts for water splitting without using sacrificial reagents or cocatalysts.

Fig.3.Electrostatic potentials of monolayer(a)CS and(b)CSe.Band alignments with respect to the redox potentials of water for monolayer(c)CS and(d)CSe.Eeff represents the internal effective electric field generated by the intrinsic dipole.

In practical applications for energy storage and conversion,the external strain would inevitably occur and usually play an important role in modulating the properties of 2D photocatalyst materials.As such,we further study the effects of in-plane biaxial strains in the range from-5%to 5%on the electronic structures and photocatalytic properties of monolayer CS and CSe.The biaxial strain is defined aswhere a0and a are the lattice parameters of nonstrained and strained systems,respectively.The positive value indicates tensile strain,while the negative value denotes compressive strain.Fig.4 gives the HSE06 calculated band structures of monolayer CX under compression(-1%,-3% and-5%)and tension(1%,3%and 5%).Comparing with the strain-free cases,we find that the band structure topologies of monolayer CS and CSe at different values of strains are similar except that the CBM of CSe shifts toΓpoint with the applied tensile strain larger than 3%.More importantly,the band gaps of monolayer CX can be significantly influenced by the biaxial strain,decreasing with increasing tensile strain and increasing with an increase of compressive strain.The band gaps of monolayer CS/CSe can be reduced to 2.23/1.77 eV under 5% tensile strain and enlarged to 2.71/2.77 eV when the biaxial compressive strain reaches-5%.Such controllable band gap properties of monolayer CX,easily achieved via biaxial strain engineering,would be beneficial for extending the visible-light adsorption.

Fig.4.Band structures of monolayer(a)CS and(b)CSe at-5%,-3%,-1%,1%,3%and 5% biaxial strain.The Fermi levels(red dashed lines)are set to zero.

Meanwhile,the dipole moment P and electrostatic potential difference of monolayer CS and CSe are expected to be effectively modulated by applying external biaxial strain.To explore this issue,the variations of the dipole moment and electrostatic potential difference as a function of the biaxial strain are studied.As shown in Fig.5a,the dipole moment of monolayer CS/CSe slightly decrease with the biaxial strain changing from-5% to 5%,for example,0.46/0.54 D,0.45/0.53 D,0.41/0.50 D and 0.39/0.49 D at-5%,-3%,3% and 5% strain.The applied strain has similar influences on the electrostatic potential difference of monolayer CX,which can be engineered from 2.48/2.54 eV to 2.12/2.16 eV for the CS and CSe monolayers,respectively.These values under the compressive and tensile strain states are still sufficiently large to induce strong intrinsic electric field,and thus remarkable charge separation can also occur in these strained CX systems.

Fig.5.Strain effects on(a)electrostatic potential differences and dipole moments,(b)band edge positions of monolayer CS and CSe.

To explore the strain effect on the reduction/oxidation ability of monolayer CS and CSe,their band edge alignments as a function of the biaxial strain are further calculated and shown in Fig.5b.It should be noted that the redox potentials of water splitting depend on the pH value of the solution,which increase with pH by pH×0.059 eV[63].Thus,the band edge alignments of the CS and CSe monolayers as well as the reduction/oxidation ability will change in different solution environments.Consequently,the water redox potential at the acidic(pH 0)and neutral environment(pH 7)are computed and plotted in Fig.5b by light gray and green horizontal dashed lines,respectively.We can see that applying compressive strain makes the Ueand Uhdistinctly increase,while they gradually decrease when introducing external tensile strain.More fascinatingly,the CBM energies of all strained monolayer CX systems are higher than the reduction potential of H+/H2and the corresponding VBM energies are lower than the oxidation potential of O2/H2O.The band edge positions are all situated at the favorable positions regardless of acidic or neutral environment,satisfying the requirements of photocatalyst for overall water splitting.

It is well-known that the light-harvesting performance is also an essential ingredient in determining the photocatalytic efficiency of water splitting.So that,we finally simulate the optical absorption spectrums of monolayer CX systems based on the HSE06 functional.According to the frequency-dependent dielectric function ε（ω）=ε1（ω）+iε2（ω）,the optical absorption coefficient αabsas a function of photon energy can be calculated by using the following equation[64],

Fig.6.Optical absorption spectrums of monolayer(a)CS and(b)CSe within the strain range of-5% to 5%.

where ε1（ω）and ε2（ω）are the real and imaginary parts of dielectric function,respectively.Fig.6 shows the obtained results for monolayer CS and CSe within the strain range of-5%～5%.Here,only the absorption spectrums that the polarization directions are parallel to the planes of monolayer CX are considered.For the CS and CSe monolayers at strain-free states,the optical absorptions can be activated by the incoming light irradiation with photon energies higher than about 2.40 eV,which are coincident with the values of the band gaps obtained in pervious calculations.Moreover,they exhibit strong light absorption in both the ultraviolet and visible-light region,and the absorption intensities are as large as about 105cm-1.As for monolayer CX systems under compressive strains,blue-shifted absorption edges in their absorption spectrums are observed,which are due to their larger band gaps.However,when the biaxial tensile strains are applied at monolayer CS and CSe,the absorption curves are considerably extended toward the visible-light region because of the narrowing of band gaps.These above results indicate that our proposed CX monolayers have outstanding absorption behaviours and can absorb as much visible-light as possible by tensile strain engineering,suggesting that they are promising metal-free photocatalyst for visible-light-driven water splitting.

For the practical application of photocatalytic water splitting,of course require the semiconductor material to be hydrophilic,which is another significant parameter characterizing its performance.To investigate how strong the interactions are between the H2O molecule and our proposed CX monolayers,several possible adsorption configurations are considered and geometrically optimized.The most energetically favorable configurations for H2O molecule adsorbed on CX monolayers are shown in Fig.S5(Supporting information).The distances from the H2O molecule to the CS and CSe monolayers are 2.29 and 2.09?,respectively,and the corresponding binding energy(Ea)are 191 meV and 235 meV.The relatively large Easuggests a strong interaction between the H2O molecule and the monolayers,ensuring better photocatalytic water-splitting activity.

In conclusion,by performing first-principles calculations,we systematically investigate the structural,electronic and photocatalytic properties of blue-phase CX(X=S,Se)monolayers.From the calculated phonon spectrums and elastic constants as well as AIMD simulations,we first verify that the proposed CX monolayers are dynamically,thermally,and mechanical stable.Then,CS and CSe monolayers are found to be indirect semiconductors with band gap of 2.41 and 2.46 eV,respectively.The calculated hole mobility along the b direction can be as high 3.23×104cm2V-1s-1for CS and 4.27×103cm2V-1s-1for CSe,which are comparable to those of many recently reported 2D photocatalysts.Furthermore,internal electrical fields are induced in these predicted monolayer materials due to the large intrinsic dipole,reducing the possibility of electron-hole recombination.More fascinatingly,CS and CSe monolayers not only possess favorable band edge positions with respect to water redox potentials and but also have strong capacity of optical absorption with the intensity up to 105cm-1.Finally,we find that applying external strain can further adjust the band gaps,band alignments,and optical absorption behaviours of CX monolayers.These theoretical findings suggest that our proposed CS and CSe monolayers have great potential in the field of photocatalytic water splitting.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China(No.21503149),by the Program for Innovative Research Team in University of Tianjin(No.TD13-5074),by the Project of Hubei University of Arts and Science(Nos.xk2020043,xk2020044,2020kypyfy015).

Appendix A.Supplementary data

Supplementary material related to this articlecanbe found, in the online version,at doi:https://doi.org/10.1016/j.cclet.2020.09.056.