Next Article in Journal
Multiferroic and Magnetodielectric Effects in Multiferroic Pr2FeAlO6 Double Perovskite
Next Article in Special Issue
Low-Temperature UVO-Sintered ZnO/SnO2 as Robust Cathode Buffer Layer for Ternary Organic Solar Cells
Previous Article in Journal
Topological Resistance-Free One-Way Transport in a Square-Hexagon Lattice Gyromagnetic Photonic Crystal
Previous Article in Special Issue
Ab Initio Study of Structural, Electronic and Magnetic Properties of TM&(B@C60) (TM = V, Cr) Sandwich Clusters and Infinite Molecular Wires
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nanomaterials
Volume 12
Issue 17
10.3390/nano12173012
nanomaterials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

N-Doped CrS2 Monolayer as a Highly-Efficient Catalyst for Oxygen Reduction Reaction: A Computational Study

by
Zengming Qin
1,
Zhongxu Wang
1,
Xiaofeng Li
2,
Qinghai Cai
2,
Fengyu Li
3,* and
Jingxiang Zhao
1,2,*
1
Key Laboratory for Photonic and Electronic Bandgap Materials, Ministry of Education, School of Physics and Electronic Engineering, Harbin Normal University, No. 1, Shida Street, Harbin 150025, China
2
College of Chemistry and Chemical Engineering, Harbin Normal University, Harbin 150025, China
3
School of Physical Science and Technology, Inner Mongolia University, Hohhot 010021, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(17), 3012; https://doi.org/10.3390/nano12173012
Submission received: 18 July 2022 / Revised: 19 August 2022 / Accepted: 29 August 2022 / Published: 30 August 2022
(This article belongs to the Special Issue First-Principles Investigations of Low-Dimensional Nanomaterials)
Browse Figures
Versions Notes

Abstract

:
Searching for low-cost and highly-efficient oxygen reduction reaction (ORR) catalysts is crucial to the large-scale application of fuel cells. Herein, by means of density functional theory (DFT) computations, we proposed a new class of ORR catalysts by doping the CrS2 monolayer with non-metal atoms (X@CrS2, X = B, C, N, O, Si, P, Cl, As, Se, and Br). Our results revealed that most of the X@CrS2 candidates exhibit negative formation energy and large binding energy, thus ensuring their high stability and offering great promise for experimental synthesis. Moreover, based on the computed free energy profiles, we predicted that N@CrS2 exhibits the best ORR catalytic activity among all considered candidates due to its lowest overpotential (0.41 V), which is even lower than that of the state-of-the-art Pt catalyst (0.45 V). Remarkably, the excellent catalytic performance of N@CrS2 for ORR can be ascribed to its optimal binding strength with the oxygenated intermediates, according to the computed linear scaling relationships and volcano plot, which can be well verified by the analysis of the p-band center as well as the charge transfer between oxygenated species and catalysts. Therefore, by carefully modulating the incorporated non-metal dopants, the CrS2 monolayer can be utilized as a promising ORR catalyst, which may offer a new strategy to further develop eligible electrocatalysts in fuel cells.

1. Introduction

The rising energy demand and depletion of fossil fuels have attracted increasing interest in alternative energy sources. This is because the traditional fossil fuels are not renewable and will produce harmful combustion products, such as CO, CO2, NO, and SO2, which thus pose a serious challenge to human health and environmental protection [1,2,3,4,5,6]. In this regard, electrochemical energy storage and conversion technologies, such as fuel cells, were regarded as efficient and clean devices for reducing the global energy shortage [7,8]. However, the efficiency of these energy-related devices was greatly determined by the oxygen reduction reaction (ORR) due to its sluggish kinetics and high overpotential [9,10,11,12,13]. At present, Pt-based materials represent the most promising ORR catalysts [14,15,16,17]. However, their practical application is severely limited by their high price and low abundance. Thus, the search for highly efficient, low-cost alternative electrocatalysts is urgent to boost the ORR.
In recent years, single-layer transition metal dichalcogenides (TMDs) have become a hotspot of theoretical and experimental research owing to their large surface area, stability, and peculiar electronic structure [18,19,20,21,22,23]. Thus, TMD-based materials have been widely employed in nanoelectronics, nanophotonics, the absorber layer in solar cells, anode materials, field-effect transistors, and so on [24,25,26,27,28]. More interestingly, the heteroatom doped TMDs have attracted great interest because they have more active sites, low cost, high stability, and high efficiency as potential electrocatalytic catalysts [29,30,31,32,33,34,35]. For example, Xiong et al. reported that co-doped MoS2 possesses high catalytic activity for oxygen evolution reaction [29]. Moreover, Lv et al. proposed that N-doped MoS2 displayed a high faradaic efficiency and low onset overpotential for CO2 electroreduction to CO [30]. Theoretically, Li et al. demonstrated that B-doping causes the VS2 monolayer to exhibit outstanding catalytic activity towards nitrogen reduction reactions [31]. Moreover, Singh et al. reported that the ORR activity of N-doped WS2 monolayers should originate from the introduction of spin density caused by N doping [32]. In addition, Tian et al. showed that Co doped 1T-TiS2 exhibits significantly enhanced performance toward the ORR [33].
As a representative of 2D TMDs, the chromium disulfide (CrS2) monolayer was synthesized via the chemical vapor deposition (CVD) method in 2019 [21]. Interestingly, due to its unique properties, the CrS2 monolayer has wide applications in spintronic devices [36,37]. For example, Chen et al. suggested that CrS2 has the most diverse electronic and magnetic properties: antiferromagnetic (AFM) metallic, non-magnetic (NM) semiconductor, and ferromagnetic (FM) semiconductor with a Curie temperature of ~1000 K [36]. Moreover, Zhang et al. reported that the magnetic properties of CrC2 monolayer can be effectively tuned by doping metal atoms [37].
Inspired by the interesting properties of the CrS2 monolayer, in this work, we explored the potential of several non-metal-doped CrS2 monolayers (X@CrS2, X = B, C, N, O, Si, P, Cl, As, Se, and Br) as ORR catalysts by performing comprehensive DFT computations. According to the computed formation energies and binding energies of these X@CrS2 systems, we suggested that these doped CrS2 candidates are very likely to be synthesized in experiments and possess extremely high stability. Moreover, the N@CrS2 catalyst was screened out as an eligible ORR catalyst with a rather low overpotential of 0.41 V, which originates from its optimal binding strengths with oxygenated intermediates based on the linear scaling relationships and volcano plot.

2. Materials and Methods

All computations were performed based on spin-polarized density functional theory, as implemented in the Vienna ab initio simulation package (VASP 541) [38,39]. The projector augmented wave (PAW) was adopted to describe the interactions between ions and electrons [40,41]. The generalized gradient approximation (GGA) in the form of Perdew–Burke–Ernzerhof (PBE) was adopted for the exchange-correlation functional [42]. Notably, the PBE functional is, nowadays, the most commonly used functional for solid-state calculations, which remains the best for the solids containing 3d-transition elements [43,44,45]. In particular, our purpose in the present work is to search for the ideal ORR catalysts among various candidates, and we thus mainly focus on the catalytic tendency of these candidates for ORR. Thus, although some more modern and better functionals have been proposed in recent years [46,47], the identified scaling relations of oxygenated species on various electrocatalysts and the derived conclusions (such as the corresponding catalytic activity) will not change. The DFT+D3 method in the Grimme scheme was employed to describe the possible weak interactions of the oxygenated species and the catalysts [48]. An energy cutoff of 500 eV was used for the plane wave (Table S1), and the convergence criteria of force and energy were set to 0.02 eV Å−1 and 10−5 eV, respectively. A k-point of 3 × 3 × 1 was sampled in the Brillouin zones (Table S1), and the vacuum space was set to 15 Å.
To estimate the ORR catalytic performance of these X@CrS2 materials, the change in the Gibbs free-energy change (ΔG) of each elementary reaction step during ORR was computed using the computational hydrogen electrode (CHE) model [15,49]: ΔG = ΔE + ΔZPE − TΔS + ΔGU, in which ΔE is the reaction energy directly obtained from DFT computations, and ΔZPE and ΔS represent the difference of zero-point energy and entropy, respectively, which can be derived from the computations on the vibrational frequencies and the standard thermodynamic data. ΔGU = −eU, where U is the applied potential. Since DFT generally fails to obtain the energy of the O2 molecule, the free energy of O2 (GO2) will be computed via the energies of H2O and H2, namely GO2 = GH2O − 2GH2 + 4.92 eV. Furthermore, the ORR’s catalytic activity was evaluated by computing the corresponding overpotential (η) according to the following equations: η = max{ΔG1, ΔG2, ΔG3, ΔG4}/e + U0, where ΔG represents the free energy change in each elementary reaction, including the formation of OOH*, O*, and OH* and the desorption of OH*, and U0 is the equilibrium potential of ORR, which is equal to 1.23 V. According to the above definition, a smaller η value corresponds to a higher catalytic activity for ORR.

3. Results

3.1. Structures, Stabilities, and Properties of X@CrS2 Catalysts

First, we investigated the structures and electronic properties of the pristine CrS2 monolayer with the 2H phase. As shown in Figure 1a, each Cr center is prismatically coordinated by six surrounding S atoms, with the S atoms in the upper layer lying directly above those of the lower layer. Furthermore, the optimized lattice constant for the CrS2 monolayer is 3.04 Å, while the lengths of the formed Cr-S bonds are 2.28 Å. Moreover, the CrS2 monolayer is a direct semiconductor with the band gap of 0.93 eV, in which the conduction band minimum and the valence band maximum located at the K point mainly originate from the contribution of 3d orbitals of the Cr atom (Figure 1b). Notably, these above results on the structure and properties of a pristine CrS2 monolayer are consistent with previous theoretical reports [50], indicating the accuracy of our computational methods to describe the behavior of the CrS2 monolayer.
Based on the optimized CrS2 monolayer, various non-metal atoms were introduced to substitute one of the S atoms within the CrS2 monolayer to construct the doped CrS2 system. After fully geometrical relaxation, the structures of the X@CrS2 monolayer are presented in Figure 2, while their corresponding structural parameters are summarized in Table 1. The results showed that, after the introduction of dopants, the structure of the CrS2 monolayer is changed in various ways, which are highly dependent on the difference in the radius between dopants and S atoms. In detail, for B, C, N, and O dopants with a smaller radius than the S atom, the introduced non-S atoms are concave from the CrS2 plane, whereas the doping of Si, P, Cl, As, Se, and Br induces an unchanged or slightly outward structure due to their comparable or slightly larger radius.
To assess the experimental feasibility of these doped CrS2 monolayers, we computed their formation energies (Ef) based on the following equation: Ef = EX@CrS2 − ECrS2 + μS − μX [51,52], where the EX@CrS2 and ECrS2 are the total energies of doped and pristine CrS2 monolayers, respectively. μX and μS are the chemical potentials of a doping atom and S atom. Notably, μX can be derived from the energy per atom in their bulk X materials. In contrast, μS is dependent on the growth conditions, in which two extreme cases were considered, including Cr-rich and S-rich conditions. According to previous studies [32,53], the μCr and μS should meet the relationship from the thermodynamic equilibrium condition: μCrS2 = μCr + 2μS, where μCrS2 represents the total energy of the pristine CrS2 monolayer per formula unit. Under Cr-rich conditions, the μCr was taken from the bulk Cr. Thus, μS can be obtained by: μS = (μCrS2 − μCr)/2. On the other hand, under S-rich conditions, the μS can be obtained from its bulk S8 state, and then μCr can be determined by: μCr = μCrS2 − 2μS. The computed formation energies of these X@CrS2 materials were listed in Table 1. We can see that the formation energies for B doping are all positive values, indicating that the formation of this kind of doped CrS2 monolayer is difficult. On the contrary, when O, Si, Cl, As, Se, and Br dopants were introduced into the CrS2 monolayer, their Ef values were always negative, implying notable feasibility for their experimental synthesis. Notably, the CrS2 monolayers doped by C, N, and P atoms are more likely to form under Cr-rich conditions due to their negative formation energies. Obviously, by carefully tuning the reaction conditions, these doped CrS2 monolayers can be easily fabricated in experiments, except for B-doping.
To further explore the structural stability of X@CrS2, we computed the binding energies (Ebind) of the introduced dopants on the CrS2 monolayer according to the following equation: Ebind = EX@CrS2 − ECrS2 − EX, where EX@CrS2, ECrS2, and EX represent the total electronic energies of doped CrS2 monolayers, the defective CrS2 substrate, and isolated X atoms in their most stable phase. It can be observed that the Ebind values of these dopants on the CrS2 substrate range from −2.52 eV of Br to −7.16 eV of C, indicative of the strong binding strength, thus ensuring their high stability. The stability of X@CrS2 was further evaluated by using AIMD simulations, where N doping was taken as a representative. As shown from Figure S1, there is still no significant distortion of the geometric structure at 500 K, indicative of its high thermodynamic stability. It should be noted that incorporation of various non-metal atoms into 2D TMDs has been already realized experimentally. For example, Li et al. reported a simple, facile, and effective strategy to fabricate an N-doped MoS2 nanosheet [54], while Xin et al. demonstrated the fabrication of a P-doped MoS2 nanosheet by the one-step hydrothermal method [55]. Therefore, we strongly believe that the as-designed doped CrS2 monolayer holds great promise for synthesis.
Next, we explored the magnetic and electronic properties of these X@CrS2 systems. Our results showed that the pristine CrS2 monolayer is a nonmagnetic material. After introducing these dopants into the CrS2 monolayer, however, different magnetic behaviors can be observed: (1) C-, O-, and Se-doped CrS2 systems still retain their nonmagnetic states due to the absence of unpaired electrons; (2) as for B, N, P, As, or Br-doped systems, the total magnetic moment is close to 1.00 μB, while the Si-doped CrS2 monolayer has a magnetic moment of 2.00 μB. Moreover, we noted that these dopants will loot different amounts of electrons (0.04~0.95 |e−|) from the CrS2 monolayer, except for the Si dopant, which can donate about 0.32 electrons to the CrS2 substrate. As a result, the band gaps of the doped CrS2 monolayer are reduced to different degrees due to the introduction of impurity levels. For example, the N-doped CrS2 system exhibits a smaller band gap of 0.25 eV than that of the pristine one (0.99 eV), suggesting the enhanced electrical conductivity, which may facilitate the rapid charge transfer in electrocatalysis [56]. It is noted that the GGA method usually underestimates band gaps, whereas the hybrid functional method, such as Heyd–Scuseria–Ernzerhof (HSE06), can rectify the band gap [57]. For example, the computed band gaps of the well-established MoS2 monolayer using HSE06 and PBE methods are 2.20 and 1.79 eV [58], respectively. Despite the fact that a more accurate band gap value can be predicted by means of the HSE06 method, the corresponding computational cost is extremely high for the 10 doped CrS2 monolayers, which consist of 25 Cr and 50 S atoms. We will mainly focus on the trend of the band gaps of the CrS2 monolayer after non-metal doping.

3.2. ORR Catalytic Activity

After knowing that these non-metal doped CrS2 candidates exhibit great potential for experimental synthesis, high stability, diverse magnetic moment, and enhanced conductivity, we further explored their catalytic activities towards ORR.
Typically, the O2 molecule can be reduced to H2O along a four-electron (4e−) pathway (Figure 3a): (1) * + O2 (g) + H+ + e− ⭢ OOH*; (2) OOH* + H+ + e− ⭢ O* + H2O; (3) O* + H+ + e− ⭢ OH*; (4) OH* + H+ + e− ⭢ H2O (l) [59,60]. We again took the N@CrS2 monolayer as the representative to compute the ∆G values of the above four elementary steps. As shown in Figure 3b, we found that the formation of OOH* species on the N site of the N@CrS2 catalyst is exothermic, with a ∆G value of −0.84 eV. The length of the formed N–O bond is 1.41 Å, and the O–O bond is elongated to 1.48 Å as compared with that of the free O2 molecule (1.23 Å). Subsequently, the approach of a second hydrogen induces the dissociation of OOH* into (O* + H2O) or the formation of H2O2. Remarkably, the OOH* formation is exothermic by 2.43 eV, which is much larger than that of H2O2 formation (0.16 eV), suggesting that N@CrS2 shows a rather high selectivity towards the 4e− pathway by greatly suppressing the competing 2e− one. Once the O* species is formed, it can be further hydrogenated to form OH* and the second H2O with the ∆G values of −0.83 and −0.82 eV, respectively. According to the CHE model, the final step (i.e., OH* desorption) is identified as the potential-determining step (PDS). Thus, the limiting potential is 0.82 V, which is the smallest applied voltage to make the whole reaction still exergonic, corresponding to the overpotential of 0.41 V.
As for other X@CrS2 candidates, the computed free adsorption energies of oxygenated species and the corresponding ηORR are summarized in Table S2. We found that the computed ηORR increases in the order of As@CrS2 (0.97 V) < C@CrS2 (1.20 V) < O@CrS2 (1.53 V) < Se@CrS2 (1.61 V) < Br@CrS2 (1.74 V) ≈ Cl@CrS2 (1.78 V) < P@CrS2 (2.18 V) < B@CrS2 (2.42 V) < Si@CrS2 (2.82 V), as shown in Figure S2, which are all higher than that of N@CrS2 (0.41 V). In particular, the overpotential of N@VS2 (0.41 V) is even lower than that of the well-established Pt catalyst (0.45 V) [15], implying its excellent ORR catalytic activity. As ORR usually occurs in aqueous solution, we also studied the solvent effect on the ORR activity of the N@CrS2 catalyst using the implicit solvation model in VASPsol with a dielectric constant of 80 [61]. It can be seen from Figure S3 that the PDS locates at the last step, and the computed η value is 0.48 V, which is comparable to that of the value without the solvent effect (0.41 V), implying that the solvent effects left the superior ORR catalytic performance of the N@CrS2 monolayer nearly unchanged. Although implicit solvent models can offer fast and inexpensive starting points for estimating the solvation effects [62,63,64,65], they may not fully describe the interactions of ORR adsorbates with water molecules due to the formation of H-bonds. As an alternative to implicit solvent approaches, explicit solvent models may provide a more comprehensive solution to describe the solvation effects on the ORR catalytic performance, which usually requires sampling thousands of solvent configurations, thus resulting in significant computational expense due to the use of classical molecular dynamic (MD) computations based on force fields. To this end, according to previous studies [66], explicit water layers were employed by placing 10 H2O molecules on the adsorbed oxygenated species (Figure S4). The results showed that the computed η value for ORR on the N@CrS2 catalyst using the explicit models is 0.50 V (Figure S4), which is close to our implicit solvent value (0.48 V). Thus, the implicit solvation model can also provide a reasonable estimation of the solvation energy for ORR intermediates, consistent with previous theoretical studies [67].

4. Discussion

Based on the well-accepted Sabatier principle, either too strong or too weak adsorption of reaction intermediates on catalysts can result in poor catalytic activity. This is because adsorption that is too strong will hamper the desorption process, resulting in poisoned catalysts, whereas too weak adsorption will induce insufficient activation of intermediates. Thus, the best catalysts exhibit the optimal adsorption strength with reaction intermediates, which locates at the peak of the volcano plot. Clearly, the ORR’s catalytic activity is intrinsically dependent on the adsorption strength of reaction intermediates with catalysts. To this end, we scaled the adsorption free energies of OOH* (ΔGOOH*), O* (ΔGO*) and OH* (ΔGOH*) on different X@CrS2 systems.
As shown in Figure 4a, obvious linear scaling relationships can be obtained between ΔGOOH* and ΔGOH* by ΔGOOH* = 0.72 ΔGOH* + 3.49 (R2 = 0.97) as well as ΔGO* and ΔGOH* by ΔGO* = 1.10 ΔGOH* + 0.72 (R2 = 0.89). Thus, ΔGOH* can be utilized as an eligible descriptor to describe the catalytic trend of these considered X@CrS2 catalysts. Furthermore, a volcano plot of ORR activity (ηORR) for X@CrS2 with the variation in ΔGOH* is obtained (Figure 4b), in which either a strong Si@CrS2 or a weak binding strength O@CrS2 with OH* species will induce poor ORR activity. On the contrary, the N@CrS2 catalyst displays a moderate binding strength with OH* and exhibits high ORR catalytic activity, making it locate at the peak of the volcano curve, and, thus, it becomes the best ORR catalyst among all the studied systems.
To gain deep insight into the remarkable difference of OH* adsorption on X@CrS2, we turn to exploring the p-band center (εp) model of the non-metal active sites [68], in which Si@CrS2, N@CrS2, and O@CrS2 systems were chosen as the representatives of strong, moderate, and weak OH* adsorption. Notably, according to this model, the position of εp closer to the Fermi level will generally induce a stronger interaction of reaction species with catalysts. Our results demonstrated that the computed εp values of Si, N, and O dopants are −0.98, −2.37, and −3.30 eV, respectively, as shown in Figure 5. The moderate εp value on N@CrS2 suggests its optimal interaction with the oxygenated species, which is responsible for its superior catalytic performance.
In addition, the charge density differences in Si@CrS2, N@CrS2, and O@CrS2 with adsorbed OH* species were computed (Figure 6). Upon OH* adsorption, we found that the charge depletion around the Si, N, and O endures, while the charge accumulation locates at the X–O bonds, indicating the charge transfer from the catalysts to oxygenated intermediates. According to the Bader charges analysis, the charge transfer is 0.71, 0.08, and 0.16 |e−| for OH* adsorption on Si@CrS2, N@CrS2, and O@CrS2, respectively, which is consistent with the binding strengths between them. Thus, the moderate adsorption strength of N@CrS2 endows it with its high ORR catalytic activity.

5. Conclusions

In summary, by performing comprehensive DFT computations, we have systematically investigated the structures, stabilities, and magnetic and electronic properties as well as the ORR catalytic activity of several non-metal doped CrS2 monolayers. Our results demonstrated that most of the doped CrS2 materials generally possess high stability and hold great promise for experimental synthesis. As expected, depending on the kinds of the incorporated dopants, these CrS2 monolayers exhibit different magnetic and electronic properties. Based on the computed free energy changes in all elementary steps during ORR, the N@CrS2 was revealed as a quite promising ORR electrocatalyst due to its lower overpotential (0.41 V) than that of the Pt benchmark (0.45 V). Moreover, obvious scaling linear relationships between oxygenated species can be obtained, which was employed to construct the volcano curve between ORR catalytic activity and OH* binding strength. Understandably, the moderate p-band center and charge transfer from catalyst to oxygenated species render the N@CrS2 catalyst’s optimal interaction with reaction species, thus rationalizing its outstanding ORR catalytic performance. Our findings not only provide a novel strategy for the design of low-cost, highly-efficient ORR electrocatalysts, but also further widen the potential applications of CrS2-based materials.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/nano12173012/s1, Figure S1: Total potential temperature and energy fluctuations after equilibration for N@CrS2 in the AIMD simulation; Figure S2: The computed free energy diagrams for ORR on different X@CrS2 systems; Figure S3: The computed free energy diagrams for ORR on N@CrS2 with solvent effect. Figure S4: The computed free energy diagrams and the corresponding intermediates configurations for ORR on N@CrS2 with the explicit solvent models; Tables S1: The computed the overpotentials (ηORR, V) at various energy cutoff (eV) or k-points; Tables S2: The computed the overpotential (ηORR) and the free adsorption energies (ΔG, eV) for various X@CrS2 materials.

Author Contributions

Conceptualization, writing—original draft preparation, Z.Q.; validation, Z.W.; writing—review and editing, formal analysis, X.L.; data curation, Q.C.; writing—review and editing, F.L.; Supervision, project administration, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Funds (NSF) for Distinguished Young Scholar of Heilongjiang Province (No. JC2018004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Topsøe, N.-Y. Mechanism of the Selective Catalytic Reduction of Nitric Oxide by Ammonia Elucidated by in Situ On-Line Fourier Transform Infrared Spectroscopy. Science 1994, 265, 1217–1219. [Google Scholar] [CrossRef]
  2. Gattuso, J.-P.; Magnan, A.; Billé, R.; Cheung, W.W.L.; Howes, E.L.; Joos, F.; Allemand, D.; Bopp, L.; Cooley, S.R.; Eakin, C.M.; et al. Contrasting futures for ocean and society from different anthropogenic CO2 emissions scenarios. Science 2015, 349, aac4722. [Google Scholar] [CrossRef]
  3. Meng, D.; Zhan, W.; Guo, Y.; Guo, Y.; Wang, L.; Lu, G. A Highly Effective Catalyst of Sm-MnOx for the NH3-SCR of NOx at Low Temperature: Promotional Role of Sm and Its Catalytic Performance. ACS Catal. 2015, 5, 5973–5983. [Google Scholar] [CrossRef]
  4. Chu, S.; Cui, Y.; Liu, N. The path towards sustainable energy. Nat. Mater. 2017, 16, 16–22. [Google Scholar] [CrossRef]
  5. Obama, B. The irreversible momentum of clean energy. Science 2017, 355, 126–129. [Google Scholar] [CrossRef]
  6. Wang, L.; Chen, W.; Zhang, D.; Du, Y.; Amal, R.; Qiao, S.; Wu, J.; Yin, Z. Surface strategies for catalytic CO2 reduction: From two-dimensional materials to nanoclusters to single atoms. Chem. Soc. Rev. 2019, 48, 5310–5349. [Google Scholar] [CrossRef]
  7. Dresselhaus, M.S.; Thomas, I.L. Alternative energy technologies. Nature 2001, 414, 332–337. [Google Scholar] [CrossRef]
  8. Panwar, N.L.; Kaushik, S.C.; Kothari, S. Role of renewable energy sources in environmental protection: A review. Renew. Sustain. Energy Rev. 2011, 15, 1513–1524. [Google Scholar] [CrossRef]
  9. Liu, H.; Logan, B.E. Electricity Generation Using an Air-Cathode Single Chamber Microbial Fuel Cell in the Presence and Absence of a Proton Exchange Membrane. Environ. Sci. Technol. 2004, 38, 4040–4046. [Google Scholar] [CrossRef]
  10. Winter, M.; Brodd, R.J. What Are Batteries, Fuel Cells, and Supercapacitors? Chem. Rev. 2004, 104, 4245–4270. [Google Scholar] [CrossRef] [Green Version]
  11. Li, Y.; Dai, H. Recent advances in zinc–air batteries. Chem. Soc. Rev. 2014, 43, 5257–5275. [Google Scholar] [CrossRef]
  12. Wen, X.; Zhang, Q.; Guan, J. Applications of metal–organic framework-derived materials in fuel cells and metal-air batteries. Coord. Chem. Rev. 2020, 409, 213214. [Google Scholar] [CrossRef]
  13. Huang, Z.-F.; Wang, J.; Peng, Y.; Jung, C.-Y.; Fisher, A.; Wang, X. Design of Efficient Bifunctional Oxygen Reduction/Evolution Electrocatalyst: Recent Advances and Perspectives. Adv. Energy Mater. 2017, 7, 1700544. [Google Scholar] [CrossRef]
  14. Marković, N.M.; Schmidt, T.J.; Stamenković, V.; Ross, P.N. Oxygen Reduction Reaction on Pt and Pt Bimetallic Surfaces: A Selective Review. Fuel Cells 2001, 1, 105–116. [Google Scholar] [CrossRef]
  15. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J.R.; Bligaard, T.; Jónsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  16. Lim, B.; Jiang, M.; Camargo, P.H.C.; Cho, E.C.; Tao, J.; Lu, X.; Zhu, Y.; Xia, Y. Pd-Pt Bimetallic Nanodendrites with High Activity for Oxygen Reduction. Science 2009, 324, 1302–1305. [Google Scholar] [CrossRef]
  17. Nie, Y.; Li, L.; Wei, Z. Recent advancements in Pt and Pt-free catalysts for oxygen reduction reaction. Chem. Soc. Rev. 2015, 44, 2168–2201. [Google Scholar] [CrossRef]
  18. Zhang, H.; Liu, L.-M.; Lau, W.-M. Dimension-dependent phase transition and magnetic properties of VS2. J. Mater. Chem. A 2013, 1, 10821–10828. [Google Scholar] [CrossRef]
  19. Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive photodetectors based on monolayer MoS2. Nat. Nanotechnol. 2013, 8, 497–501. [Google Scholar] [CrossRef]
  20. Kou, L.; Du, A.; Chen, C.; Frauenheim, T. Strain engineering of selective chemical adsorption on monolayer MoS2. Nanoscale 2014, 6, 5156–5161. [Google Scholar] [CrossRef] [Green Version]
  21. Habib, M.R.; Wang, S.; Wang, W.; Xiao, H.; Obaidulla, S.M.; Gayen, A.; Khan, Y.; Chen, H.; Xu, M. Electronic properties of polymorphic two-dimensional layered chromium disulphide. Nanoscale 2019, 11, 20123–20132. [Google Scholar] [CrossRef]
  22. Zhao, R.; Wang, T.; An, Y.; Dai, X.; Xia, C. VS2 nanosheet as a promising candidate of recycle and reuse NO2 gas sensor and capturer: A DFT study. J. Phys. Condens. Matter 2021, 33, 165501. [Google Scholar] [CrossRef]
  23. Tursun, M.; Wu, C. Electrocatalytic Reduction of N2 to NH3 Over Defective 1T′-WX2 (X=S, Se, Te) Monolayers. ChemSusChem 2022, 15, e202200191. [Google Scholar] [CrossRef]
  24. Qin, Z.; Wang, Z.; Zhao, J. Computational screening of single-atom catalysts supported by VS2 monolayers for electrocatalytic oxygen reduction/evolution reactions. Nanoscale 2022, 14, 6902–6911. [Google Scholar] [CrossRef]
  25. Qin, Z.; Zhao, J. 1 T-MoSe2 monolayer supported single Pd atom as a highly-efficient bifunctional catalyst for ORR/OER. J. Colloid Interface Sci. 2022, 605, 155–162. [Google Scholar] [CrossRef]
  26. Ling, F.; Xia, W.; Li, L.; Zhou, X.; Luo, X.; Bu, Q.; Huang, J.; Liu, X.; Kang, W.; Zhou, M. Single Transition Metal Atom Bound to the Unconventional Phase of the MoS2 Monolayer for Catalytic Oxygen Reduction Reaction: A First-Principles Study. ACS Appl. Mater. Interfaces 2021, 13, 17412–17419. [Google Scholar] [CrossRef]
  27. Meza, E.; Diaz, R.E.; Li, C.W. Solution-Phase Activation and Functionalization of Colloidal WS2 Nanosheets with Ni Single Atoms. ACS Nano 2020, 14, 2238–2247. [Google Scholar] [CrossRef]
  28. Zhang, J.; Xu, X.; Yang, L.; Cheng, D.; Cao, D. Single-Atom Ru Doping Induced Phase Transition of MoS2 and S Vacancy for Hydrogen Evolution Reaction. Small Methods 2019, 3, 1900653. [Google Scholar] [CrossRef]
  29. Xiong, Q.; Zhang, X.; Wang, H.; Liu, G.; Wang, G.; Zhang, H.; Zhao, H. One-step synthesis of cobalt-doped MoS2 nanosheets as bifunctional electrocatalysts for overall water splitting under both acidic and alkaline conditions. Chem. Commun. 2018, 54, 3859–3862. [Google Scholar] [CrossRef]
  30. Lv, K.; Suo, W.; Shao, M.; Zhu, Y.; Wang, X.; Feng, J.; Fang, M.; Zhu, Y. Nitrogen doped MoS2 and nitrogen doped carbon dots composite catalyst for electroreduction CO2 to CO with high Faradaic efficiency. Nano Energy 2019, 63, 103834. [Google Scholar] [CrossRef]
  31. Li, Q.; Guo, Y.; Tian, Y.; Liu, W.; Chu, K. Activating VS2 basal planes for enhanced NRR electrocatalysis: The synergistic role of S-vacancies and B dopants. J. Mater. Chem. A 2020, 8, 16195–16202. [Google Scholar] [CrossRef]
  32. Singh, Y.; Back, S.; Jung, Y. Activating Transition Metal Dichalcogenides by Substitutional Nitrogen-Doping for Potential ORR Electrocatalysts. ChemElectroChem 2018, 5, 4029–4035. [Google Scholar] [CrossRef]
  33. Tian, S.; Tang, Q. Activating transition metal dichalcogenide monolayers as efficient electrocatalysts for the oxygen reduction reaction via single atom doping. J. Mater. Chem. C 2021, 9, 6040–6050. [Google Scholar] [CrossRef]
  34. Saraf, D.; Chakraborty, S.; Kshirsagar, A.; Ahuja, R. In pursuit of bifunctional catalytic activity in PdS2 pseudo-monolayer through reaction coordinate mapping. Nano Energy 2018, 49, 283–289. [Google Scholar] [CrossRef]
  35. Zhou, W.; Jiang, Z.; Chen, M.; Li, Z.; Luo, X.; Guo, M.; Yang, Y.; Yu, T.; Yuan, C.; Wang, S. Directly anchoring non-noble metal single atoms on 1T-TMDs with tip structure for efficient hydrogen evolution. Chem. Eng. J. 2022, 428, 131210. [Google Scholar] [CrossRef]
  36. Chen, K.; Deng, J.; Yan, Y.; Shi, Q.; Chang, T.; Ding, X.; Sun, J.; Yang, S.; Liu, J.Z. Diverse electronic and magnetic properties of CrS2 enabling strain-controlled 2D lateral heterostructure spintronic devices. NPJ Comput. Mater. 2021, 7, 79. [Google Scholar] [CrossRef]
  37. Zhang, J.; Zheng, H.; Han, R.; Du, X.; Yan, Y. Tuning magnetic properties of CrS2 monolayer by doping transition metal and alkaline-earth atoms. J. Alloys Compd. 2015, 647, 75–81. [Google Scholar] [CrossRef]
  38. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558. [Google Scholar] [CrossRef] [PubMed]
  39. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef]
  40. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. [Google Scholar] [CrossRef] [Green Version]
  41. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. [Google Scholar] [CrossRef]
  42. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  43. Tran, F.; Laskowski, R.; Blaha, P.; Schwarz, K. Performance on molecules, surfaces, and solids of the Wu-Cohen GGA exchange-correlation energy functional. Phys. Rev. B 2007, 75, 115131. [Google Scholar] [CrossRef]
  44. Ropo, M.; Kokko, K.; Vitos, L. Assessing the Perdew-Burke-Ernzerhof exchange-correlation density functional revised for metallic bulk and surface systems. Phys. Rev. B 2008, 77, 195445. [Google Scholar] [CrossRef]
  45. Haas, P.; Tran, F.; Blaha, P. Calculation of the lattice constant of solids with semilocal functionals. Phys. Rev. B 2009, 79, 085104. [Google Scholar] [CrossRef]
  46. Haas, P.; Tran, F.; Blaha, P.; Schwarz, K.; Laskowski, R. Insight into the performance of GGA functionals for solid-state calculations. Phys. Rev. B 2009, 80, 195109. [Google Scholar] [CrossRef]
  47. Pedroza, L.S.; da Silva, A.J.R.; Capelle, K. Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids. Phys. Rev. B 2009, 79, 201106. [Google Scholar] [CrossRef]
  48. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. Phys. Rev. B 2006, 27, 1787. [Google Scholar] [CrossRef]
  49. Peterson, A.A.; Abild-Pedersen, F.; Studt, F.; Rossmeisl, J.; Nørskov, J.K. How copper catalyzes the electroreduction of carbon dioxide into hydrocarbon fuels. Energy Environ. Sci. 2010, 3, 1311–1315. [Google Scholar] [CrossRef]
  50. Tian, X.-H.; Zhang, J.-M. First-Principles Prediction of the Structural, Electronic, and Magnetic Properties of Nonmetal Atoms Doped Single-Layer CrS2. Phys. Status Solidi B 2019, 256, 1900149. [Google Scholar] [CrossRef]
  51. Dolui, K.; Rungger, I.; Das Pemmaraju, C.; Sanvito, S. Possible doping strategies for MoS2 monolayers: An ab initio study. Phys. Rev. B 2013, 88, 075420. [Google Scholar] [CrossRef]
  52. Zhou, J.; Lin, J.; Sims, H.; Jiang, C.; Cong, C.; Brehm, J.A.; Zhang, Z.; Niu, L.; Chen, Y.; Zhou, Y.; et al. Synthesis of Co-Doped MoS2 Monolayers with Enhanced Valley Splitting. Adv. Mater. 2020, 32, 1906536. [Google Scholar] [CrossRef]
  53. Zhang, H.; Tian, Y.; Zhao, J.; Cai, Q.; Chen, Z. Small Dopants Make Big Differences: Enhanced Electrocatalytic Performance of MoS2 Monolayer for Oxygen Reduction Reaction (ORR) by N– and P–Doping. Electrochim. Acta 2017, 225, 543–550. [Google Scholar] [CrossRef]
  54. Li, R.; Yang, L.; Xiong, T.; Wu, Y.; Cao, L.; Yuan, D.; Zhou, W. Nitrogen doped MoS2 nanosheets synthesized via a low-temperature process as electrocatalysts with enhanced activity for hydrogen evolution reaction. J. Power Sources 2017, 356, 133–139. [Google Scholar] [CrossRef]
  55. Xin, X.; Song, Y.; Guo, S.; Zhang, Y.; Wang, B.; Wang, Y.; Li, X. One-step synthesis of P-doped MoS2 for efficient photocatalytic hydrogen production. J. Alloys Compd. 2020, 829, 154635. [Google Scholar] [CrossRef]
  56. Liu, M.; Zhang, R.; Chen, W. Graphene-Supported Nanoelectrocatalysts for Fuel Cells: Synthesis, Properties, and Applications. Chem. Rev. 2014, 114, 5117–5160. [Google Scholar] [CrossRef]
  57. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207–8215. [Google Scholar] [CrossRef]
  58. Wang, Z.; Zhao, J.; Cai, Q.; Li, F. Computational screening for high-activity MoS2 monolayer-based catalysts for the oxygen reduction reaction via substitutional doping with transition metal. J. Mater. Chem. A 2017, 5, 9842–9851. [Google Scholar] [CrossRef]
  59. Li, S.; Miao, P.; Zhang, Y.; Wu, J.; Zhang, B.; Du, Y.; Han, X.; Sun, J.; Xu, P. Recent Advances in Plasmonic Nanostructures for Enhanced Photocatalysis and Electrocatalysis. Adv. Mater. 2021, 33, 2000086. [Google Scholar] [CrossRef]
  60. Tong, M.; Sun, F.; Xie, Y.; Wang, Y.; Yang, Y.; Tian, C.; Wang, L.; Fu, H. Operando Cooperated Catalytic Mechanism of Atomically Dispersed Cu−N4 and Zn−N4 for Promoting Oxygen Reduction Reaction. Angew. Chem. Int. Ed. 2021, 60, 14005–14012. [Google Scholar] [CrossRef]
  61. Mathew, K.; Sundararaman, R.; Letchworth-Weaver, K.; Arias, T.; Hennig, R.G. Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways. J. Chem. Phys. 2014, 140, 084106. [Google Scholar] [CrossRef] [PubMed]
  62. Xu, H.; Cheng, D.; Cao, D.; Zeng, X.C. A universal principle for a rational design of single-atom electrocatalysts. Nat. Catal. 2018, 1, 339–348. [Google Scholar] [CrossRef]
  63. Chen, Y.; Cheng, T.; Goddard Iii, W.A. Atomistic Explanation of the Dramatically Improved Oxygen Reduction Reaction of Jagged Platinum Nanowires, 50 Times Better than Pt. J. Am. Chem. Soc. 2020, 142, 8625–8632. [Google Scholar] [CrossRef]
  64. Mitchell, C.G.; Rachel, M.A.; Derek, J.H.; Yasemin, B.; Donald, F.R.; Steven, A.P.; John, A.K. Doped Amorphous Ti Oxides To Deoptimize Oxygen Reduction Reaction Catalysis. J. Phys. Chem. C 2017, 121, 16825–16830. [Google Scholar]
  65. Li, F.; Ai, H.; Shen, S.; Geng, J.; Lo, K.H.; Pan, H. Two-Dimensional Dirac Nodal Line Carbon Nitride to Anchor Single-Atom Catalyst for Oxygen Reduction Reaction. ChemSusChem 2022, 15, e202102537. [Google Scholar] [CrossRef]
  66. Svane, K.L.; Hansen, H.A.; Vegge, T. A comparison of single and double Co sites incorporated in N-doped graphene for the oxygen reduction reaction. J. Catal. 2021, 393, 230–237. [Google Scholar] [CrossRef]
  67. Duan, Z.; Henkelman, G. Identification of Active Sites of Pure and Nitrogen-Doped Carbon Materials for Oxygen Reduction Reaction Using Constant-Potential Calculations. J. Phys. Chem. C 2020, 124, 12016–12023. [Google Scholar] [CrossRef]
  68. Zhang, X.; Shi, S.; Gu, T.; Li, L.; Yu, S. The catalytic activity and mechanism of oxygen reduction reaction on P-doped MoS2. Phys. Chem. Chem. Phys. 2018, 20, 18184–18191. [Google Scholar] [CrossRef]
Nanomaterials 12 03012 g001 550
Figure 1. (a) The obtained stable configurations and (b) computed band structure for the pristine CrS2 monolayer.
Figure 1. (a) The obtained stable configurations and (b) computed band structure for the pristine CrS2 monolayer.
Nanomaterials 12 03012 g001
Nanomaterials 12 03012 g002 550
Figure 2. The obtained stable configurations for X@CrS2 monolayers (N@CrS2 was chosen as a representative).
Figure 2. The obtained stable configurations for X@CrS2 monolayers (N@CrS2 was chosen as a representative).
Nanomaterials 12 03012 g002
Nanomaterials 12 03012 g003 550
Figure 3. (a) The involved reaction pathway and (b) the computed free energy diagrams for ORR pathways on the N@CrS2 monolayer at zero and applied voltages.
Figure 3. (a) The involved reaction pathway and (b) the computed free energy diagrams for ORR pathways on the N@CrS2 monolayer at zero and applied voltages.
Nanomaterials 12 03012 g003
Nanomaterials 12 03012 g004 550
Figure 4. (a) Scaling relations between the free adsorption energies of intermediates (ΔGO* vs. ΔGOH* and ΔGOOH* vs. ΔGOH*) and (b) ηORR vs. ΔGOH* on different X@CrS2 systems.
Figure 4. (a) Scaling relations between the free adsorption energies of intermediates (ΔGO* vs. ΔGOH* and ΔGOOH* vs. ΔGOH*) and (b) ηORR vs. ΔGOH* on different X@CrS2 systems.
Nanomaterials 12 03012 g004
Nanomaterials 12 03012 g005 550
Figure 5. The computed partial density of states for Si@CrS2, N@CrS2, and O@CrS2. The black and pink vertical dotted lines represent the Fermi energy level and p-band center (εp), respectively.
Figure 5. The computed partial density of states for Si@CrS2, N@CrS2, and O@CrS2. The black and pink vertical dotted lines represent the Fermi energy level and p-band center (εp), respectively.
Nanomaterials 12 03012 g005
Nanomaterials 12 03012 g006 550
Figure 6. The charge difference density of OH* adsorption on (a) Si@CrS2N, (b) N@CrS2, and (c) O@CrS2 with the isosurface of 0.003 e/Å3.Cyane and yellow bubbles represent positive and negative charges, respectively.
Figure 6. The charge difference density of OH* adsorption on (a) Si@CrS2N, (b) N@CrS2, and (c) O@CrS2 with the isosurface of 0.003 e/Å3.Cyane and yellow bubbles represent positive and negative charges, respectively.
Nanomaterials 12 03012 g006
Table
Table 1. The computed bond’s length of X-Cr (dX-Cr, Å), formation energies (Ef, eV), binding energies (Ebind, eV), magnetic moment (μB), charge transfer (Q, e−), and band gaps (Egap, eV) for various X@CrS2 monolayers.
Table 1. The computed bond’s length of X-Cr (dX-Cr, Å), formation energies (Ef, eV), binding energies (Ebind, eV), magnetic moment (μB), charge transfer (Q, e−), and band gaps (Egap, eV) for various X@CrS2 monolayers.
dX-CrEfEbindμBQEgap
S-richCr-rich
pristine2.28///0.00/0.93
B1.972.060.94−5.621.000.090.20
C1.920.76−0.36−7.160.010.770.90
N1.870.16−0.96−5.751.000.940.25
O1.93−2.70−3.83−6.960.010.950.96
Si2.47−3.40−4.53−3.932.000.320.27
P2.350.15−0.97−4.151.000.270.95
Cl2.41−2.48−3.61−3.141.000.470.18
As2.49−0.02−1.14−3.700.960.040.00
Se2.42−1.17−2.30−4.760.010.280.98
Br2.56−0.46−1.58−2.521.000.310.19
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qin, Z.; Wang, Z.; Li, X.; Cai, Q.; Li, F.; Zhao, J. N-Doped CrS2 Monolayer as a Highly-Efficient Catalyst for Oxygen Reduction Reaction: A Computational Study. Nanomaterials 2022, 12, 3012. https://doi.org/10.3390/nano12173012

AMA Style

Qin Z, Wang Z, Li X, Cai Q, Li F, Zhao J. N-Doped CrS2 Monolayer as a Highly-Efficient Catalyst for Oxygen Reduction Reaction: A Computational Study. Nanomaterials. 2022; 12(17):3012. https://doi.org/10.3390/nano12173012

Chicago/Turabian Style

Qin, Zengming, Zhongxu Wang, Xiaofeng Li, Qinghai Cai, Fengyu Li, and Jingxiang Zhao. 2022. "N-Doped CrS2 Monolayer as a Highly-Efficient Catalyst for Oxygen Reduction Reaction: A Computational Study" Nanomaterials 12, no. 17: 3012. https://doi.org/10.3390/nano12173012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop