Electronic, Magnetic, and Catalytic Properties of ... - ACS Publications

Sep 24, 2017 - However, the overpotential in the average process is far away from zero ..... (64) Pu, Z.; Liu, Q.; Tang, C.; Asiri, A. M.; Sun, X. Ni2...
0 downloads 0 Views 2MB Size
Article Cite This: Chem. Mater. 2017, 29, 8892-8900

pubs.acs.org/cm

Electronic, Magnetic, and Catalytic Properties of Thermodynamically Stable Two-Dimensional Transition-Metal Phosphides Yangfan Shao,†,‡ Xingqiang Shi,*,‡ and Hui Pan*,† †

Institute of Applied Physics and Materials Engineering, University of Macau, Macau SAR, China Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China



S Supporting Information *

ABSTRACT: Two-dimensional (2D) nanomaterials have attracted extensive interest due to their unique properties and versatile applications. In the current work, a new family of 2D nanostructures, 2D transition-metal phosphide monolayer (M2P), is designed to find its multifunctional applications on the basis of density functional theory (DFT) calculations. We show that the 2D M2P monolayers are thermodynamically stable by carrying out molecular dynamic simulations, analyzing their phonon dispersions, and calculating their mechanical properties. We find that all of the stable 2D M2P monolayers are metallic, and their magnetic states strongly depend on the transition metals and geometric phases. We further demonstrate that 2D W2P and Fe2P monolayers in the 2H phase show the best catalytic performance in hydrogen evolution reaction (HER) due to the relatively low overpotentials at high hydrogen coverage. Importantly, oxygen functionalization can efficiently improve the HER activities of the 2D M2P monolayers at low hydrogen coverage. Our systematic study predicts that these new thermodynamically stable 2D transition-metal phosphides could be applicable in nanodevices, spintronics, and electrocatalysis.

1. INTRODUCTION Since the discovery of graphene, two-dimensional (2D) materials have triggered increasing interest because of their unusual physical and chemical properties1−4 and wide applications from nanodevices to catalysis.5−10 For example, 2D transition-metal dichalcogenides (TMDs) showed intriguing electronic and optical properties, leading to versatile applications in logic transistors, chemical sensors, Li-ion batteries, catalysts, etc.11−23 Recently, MXenes, 2D transitionmetal carbides, nitrides, and carbonitrides, have attracted great attention because of their specific metallic and magnetic properties.24−28 The MXene monolayer has a chemical formula of MnXn−1 (n = 2−4),27,29 where M is a transition metal and X is carbon or nitrogen. For example, transition metals in M2X are at the top and bottom, and the carbon or nitrogen layer is in between. MXenes have been fabricated by chemical etching30 and chemical vapor deposition (CVD),31 and have found various applications in energy storage,32−37 hydrogen production,38 and electromagnetic interference shield.39 The CVD is viewed as a very promising and scalable fabrication technique, which has recently made tremendous advances in growth of high-quality and large-scale 2D materials.40,41 MXenes showed high catalytic activity in hydrogen evolution reactions (HERs) in a wide range of hydrogen densities, which is comparable to novel metals.42 Most recently, 2D Tl2O was predicted on the basis of DFT calculations, and the family of metal-shrouded 2D materials has expanded to metal oxides.43 © 2017 American Chemical Society

Recently, M. Yoon et al. reported theoretically that layered transition-metal phosphides (TMPs), Sr2P and Ba2P, are thermodynamically stable and could be applied to 2D electrides.44 However, a systematic study on 2D TMPs and their various applications has not been available. In this work, we investigate the structural stability, and electronic, magnetic, and catalytic properties of a new family of 2D TMP monolayers on the basis of DFT calculations. We consider 24 M2P monolayers, including M = Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W, Mn, Fe, and Ni in both trigonal (T) and hexagonal (H) phases. We find that nine systems are stable by analyzing the dynamic, thermal, and mechanical properties. We show that all of the stable monolayers are metallic, while their magnetic ground states depend on the transition metals and geometrical phases. We further show that the catalytic activity of 2D TMPs depends also on the transition metal and phase, and oxygen functionalization can efficiently improve their HER activities.

2. COMPUTATIONAL METHOD The spin-polarized DFT calculations were performed by using the Vienna ab initio simulation package (VASP).45,46 Projector augmented wave (PAW) potentials were employed to describe the core electrons, and the valence electron orbitals were expanded on the basis of plane Received: September 11, 2017 Revised: September 24, 2017 Published: September 24, 2017 8892

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials waves with a kinetic energy cutoff of 500 eV.47 The Perdew−Burke− Ernzerhof generalized gradient approximation (PBE−GGA)48 exchange and correlation functional was used. For structure relaxation, both lattice constants and atom coordinates were optimized until the force is less than 0.02 eV/Å. The Monkhorst−Pack k-point sampling for the 2D cell was 15 × 15. For density of states (DOS) calculations, Gaussian smearing was 0.05 eV, and Monkhorst−Pack k-point sampling was 41 × 41. The vacuum layer in the vertical direction was 20 Å. The phonon dispersion curves were calculated by using the finite displacement method as implemented in the Phonopy program,49 which is interfaced with the density functional perturbation theory as implemented in VASP. The 3 × 3 supercells were adopted to simulate phonon dispersion curves. The kinetic energy cutoff for phonon calculations was 700 eV.

nearly single-atomic-layer structure) with a vertical distance of 0.69 Å after geometry optimization (Table S2). 3.1. Dynamic, Thermal, and Mechanical Stability. The structural stability cannot be ensured, even though the structures can be optimized by minimizing the total energy. To further confirm the structural stability of 2D TMPs, we calculate the phonon dispersions of the energetically stable M2P systems, which is an important indication of structural stability. The phonon modes at some K points in the Brillouin zone will be imaginary when a system is dynamically unstable. There are no significant imaginary vibrational frequencies in the calculated phonon dispersion curves of 2D M2P monolayers (M = Ti, Zr, Nb, Ta, Cr, Mo, W, and Fe) in the 2H phase (M2P−2H), indicating that these systems are dynamically stable (Figure S1a,c−g,i,k). For 2D M2P in the 1T phase (M2P−1T), Zr2P−1T, W2P−1T, and Fe2P−1T are found to be dynamically stable from the calculated phonon dispersion (Figure S1b,h,j). However, other systems (Ti2P, V2P, Nb2P, Ta2P, Cr2P, Mo2P, Mn2P, and Ni2P in 1T phases; and Mn2P and Ni2P in 2H phases) are dynamically unstable because of imaginary frequencies in phonon dispersions (Figure S2). The above dynamically stable structures are obtained at T = 0 K. To investigate their thermal stability at higher temperatures, we carry out ab initio molecular dynamics (AIMD) calculations. After the supercells with 108 atoms are annealed at 300 K for 3 ps with a time step of 1 fs, most of the predicted stable structures, except Cr2P−2H and W2P−1T (Figure S3), have no structure reconstruction. We thus conclude that the predicted M2P−2H (M = Ti, Zr, Nb, Ta, Mo, W, and Fe) and M2P−1T (M = Zr and Fe) can be stable at room temperatures. Although the MD simulation confirms the thermodynamic stability of the systems, it is necessary to assess the effect of lattice distortion on structural stability because the supercells are fixed in the MD simulation, which can be confirmed by calculating elastic constants. Mechanical stabilities of these systems are obtained on the basis of Born criteria.50 For a 2D sheet, the four nonzero 2D elastic constants are c11, c22, c12, and c66 using the standard Voigt notation, 1-xx, 2-yy, and 6-xy.51 The hexagonal structure has c11 = c22 because of the symmetry 1 and the additional relation c66 = 2 (c11 − c12). 51 As an indication of a mechanically stable 2D sheet, the elastic constants must satisfy c11c22 − c212 > 0 and c66 > 0. The in-plane planar Young’s and shear moduli are derived from elastic constants as follows:

3. RESULTS AND DISCUSSION In our calculations, we focus on searching new stable 2D transition-metal phosphide monolayers with a sandwich structure (M2P). The M2P monolayer is a three-atom-thick layer in a sequence of M−P−M and has two possible phases (1T and 2H). Both 1T and 2H phases are considered in our study to find the stable structure (Figure 1). The H structure is

Figure 1. Top and side views of 2D M2P in 1T (a) and 2H (b) phases. The purple and blue spheres denote phosphorus and transition metals, respectively.

honeycomb with D3h point-group symmetry, while the T structure is centered honeycomb with C3v symmetry. Here, a new family of 2D M2P monolayers (M = Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W, Mn, Fe, and Ni) is designed. To find their structural stability, we first calculate the total energies of M2P monolayers in both T and H phases, respectively. The lowest total energy can be found by optimizing the supercell size and geometry as well as the atomic coordinates in the supercell. Hf2P monolayers in 1T and 2H phases are unstable because their total energies cannot converge to a local minimum, and will not be considered in the following discussion. The 1T phases of Ti2P, Zr2P, and Mn2P monolayers are more energetically stable than their 2H counterparts because their energies in the 1T phase are lower than those in the 2H phase by 0.210, 0.362, and 0.350 eV per unit cell, respectively. For other 2D M2P structures (M = V, Nb, Ta, Cr, Mo, W, Fe, and Ni), their 2H phases are more energetically stable than the 1T phases (Table S1). The optimized lattice constants of M2P (M = Ti, Zr, V, Nb, Ta, Cr, Mo, W, Mn, Fe, and Ni) are listed in Table S2. We notice that only Fe2P−1T becomes nearly flat (a

Ex2D =

2D 2D 2D 2D c11 c 22 − c12 c 21 2D c 22

,

Ey2D =

2D 2D 2D 2D c11 c 22 − c12 c 21 2D c11

,

2D 2D Gxy = c66 2D 2D 2D 2D where c2D 11 , c12 , c21 , c22 , and c66 are components of the 2D elastic constant. The calculated elastic constants (Table 1) 2D show that the Young’s moduli (E2D x and Ey ), and in-plane 2D planar shear moduli (Gxy ), are positive, indicating that the structures are mechanically stable. We conclude, therefore, that the 2D M2P−1T (M = Zr, W, and Fe) and M2P−2H (M = Ti, Zr, Nb, Ta, Cr, Mo, W, and Fe) are mechanically stable. By calculating the total energies and the dynamical, thermal, and mechanical properties of 2D M2P systems (Figure 2, Figures S1−S4, Table 1, and Table S1), we find: (1) Hf2P in both 1T and 2H phases is unstable; (2) there are 11 M2P monolayers, M2P−1T (M = Zr, W, and Fe) and M2P−2H (M = Ti, Zr, Nb, Ta, Cr, Mo, W, and Fe), that are dynamically and

8893

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials Table 1. In-Plane Planar Elastic Constants c2D 11 , Young’s 2D 2D Moduli E2D x and Ey , and In-Plane Planar Shear Moduli Gxy of 2D M2P−1T (M = Zr, W, and Fe) and M2P−2H (M = Ti, Zr, Nb, Ta, Cr, Mo, W, and Fe) Systems Ti2P−2H Zr2P−1T Zr2P−2H Nb2P−2H Ta2P−2H Cr2P−2H Mo2P−2H W2P−1T W2P−2H Fe2P−1T Fe2P−2H

c2D 11 (N/m)

E2D x (N/m)

E2D y (N/m)

G2D xy (N/m)

51 78 59 129 156 27 114 228 96 48 44

51 55 59 87 112 27 22 48 96 23 41

51 55 59 87 112 27 22 44 21 23 40

27 18 30 28 36 16 6 17 6 7 28

= Zr and Fe) and M2P−2H (M = Ti, Zr, Nb, Ta, Mo, W, and Fe), are stable. Accordingly, their electronic and magnetic properties are investigated. We calculate their spin-polarized density of states (DOS) and band structures to investigate the electronic properties of the stable 2D TMP monolayers (Figures S7−S12). The electronic structures show that all of the stable 2D M2P monolayers are metallic.

Figure 3. (a) Calculated energy difference of the M2P monolayer between ferromagnetic and antiferromagnetic states. Alignment configurations of magnetic moments: ferromagnetic ground state (b) and antiferromagnetic ground state (c).

mechanically stable; (3) Cr2P−2H and W2P−1T are thermally unstable; and (4) nine 2D M2P systems (Ti2P, Nb2P, Ta2P, Mo2P, and W2P in the 2H phase; Zr2P and Fe2P in both 1T and 2H) are thermodynamically stable. The 1T phase of the Zr2P monolayer is more energetically stable than its 2H phase, while the 2H phase of Fe2P is more energetically stable than its 1T counterpart. Our calculations show that the stability of 2D M2P strongly depends on its composition and structure. As reported in literatures, 2D transition-metal dichalcogenides (TMDs) can be stable in 2H or 1T phases. For example, MoS2 is stable in 2H,19 while TiS2 is stable in 1T.52 Thus, we believe that the stability of 2D M2P is determined by the hybridization between d and p electrons from M and P atoms, respectively, and the symmetry. In addition, we calculated the formation energies of the stable 2D M2P monolayer to evaluate the formation possibility in experiments by comparing with stable MxPy compounds that experimentally exist. The negative (positive) formation energy indicates that the process of fabrication is exothermic (endothermic). According to the convex hull, most new 2D M2P monolayers, except Ta2P−2H and W2P−2H, can be synthesized from bulk metal and phosphorus easily (Figure S4). Although the fabrication of Ta2P−2H and W2P−2H is endothermic, the low formation energies indicate the possibility of experimental realization. 3.2. Electronic and Magnetic Properties. On the basis of the analysis of the dynamical, mechanical, and thermal stability, we find that nine 2D M2P monolayers, including M2P−1T (M

The stable M2P monolayers show a rich diversity of magnetic properties. Depending on the transition metal and phase, 2D M2P monolayers can be ferromagnetic (FM), antiferromagnetic (AFM), or nonmagnetic (NM) metals. There are several possible AFM orderings in a supercell with 2 × 2 × 1 unit cells, named as AFM1, AFM2, AFM3, AFM4, and AFM5 (Figure S5). For AFM1, the metal atoms in the same layer have the same magnetic ordering, while the atoms in a different metal layer have opposite magnetic ordering. For AFM2 and AFM3, the two atoms in the same metal layer of the supercell have opposite magnetic directions. For AFM4 and AFM5, the spin of one atom is opposite to those of another three atoms in the same metal layer of a supercell. In the AFM3 and AFM5 configurations, the spins in the different metal layers are antiparallel to each other. In contrast, for AFM2 and AFM4, the spins in the different metal layers are parallel. All of the possible configurations are relaxed to find the stable magnetic ground state. Our calculations show that AFM2, AFM3, AFM4, and AFM5 are destroyed and automatically convert to the FM or AFM1 state after relaxation. The calculated exchange energies (EAFM − EFM) show that the ground magnetic states of Ti2P− 2H, Zr2P−1T, Fe2P−1T, and Fe2P−2H are ferromagnetic

Figure 2. Summary of dynamical, mechanical, and thermal stability analysis, and magnetic properties of 2D M2Ps. M is a transition metal. FM, AFM, and NM refer to ferromagnetic, antiferromagnetic, and nonmagnetic ordering, respectively. 8894

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials (Figure 3a,b) because the exchange energies are positive; Zr2P−2H is antiferromagnetic (Figure 3a,c) because of the negative exchange energy, and the other 2D M2P monolayers (M = Nb, Ta, Mo, and W) are nonmagnetic because of the zero exchange energies. The magnetic moments of P atoms in the ferromagnetic M2P are antiparallel to those of metallic atoms (Table 2 and Figure 3b). For a better understanding of the

Ti2P−2H, Fe2P−1T, and Fe2P−2H monolayers are ferromagnetic, and the d orbitals of metal atoms dominate their magnetic properties. In addition, there is a strong spin polarization around the Fermi level of Ti2P−2H, Fe2P−1T, and Fe2P−2H, and the out-of-plane components Ti_dπ and Fe_dπ make a key contribution (Figure S6c−e). For the Ti2P− 2H monolayer, the spin-up bands are partially occupied, and the spin-down bands are nearly empty, leading to a magnetic moment of 0.77 μB per Ti atom (Figure S6c). For Fe2P−1T and Fe2P−2H, the spin-up bands are fully occupied, while the spin-down bands are partially occupied, resulting in a magnetic moment of 2.13 μB in the 1T phase and 1.49 μB in 2H phase per Fe atom (Figure S6d,e). The magnetic moment is mainly contributed by the Fe_d orbitals. The p electrons of P atoms in ferromagnetic Ti2P−2H, Fe2P−1T, and Fe2P−2H monolayers are also weakly spin-polarized, leading to small magnetic moments of 0.045, 0.032, and 0.081 μB per P atom, respectively, which are antiparallel to those of metal atoms (Figure S6c−e). For magnetic systems, the spin-polarization ratios (SPRs) are calculated to quantify the effect on spin polarization of M2P with ferromagnetic ordering. The spinpolarization ratio (SPR) is defined as SPR(E) = [DOS↑(E) − DOS↓(E)]/[DOS↑(E) + DOS↓(E)], where DOS↑(E) and DOS↓(E) are the DOS for the majority-spin and minority-spin, respectively. There is a strong spin polarization of 50−85% around the Fermi level of Ti2P−2H (Figure 4b). For the Fe2P− 2H structure (Figure S10b), there is also strong spin polarization around the Fermi level, which could be useful to spintronics. 3.3. HER Activity of Pure M2P. Transition-metal phosphides have been reported to be promising catalysts for hydrogen evolution reactions (HERs).58−68 The highly conductive 2D M2P monolayers may lead to high catalytic activity in the HER. To study the catalytic activity of stable 2D TMP monolayers, including M2P−2H (M = Ti, Zr, Nb, Ta, Mo, W, and Fe) and M2P−1T (M = Zr and Fe), we first investigate the adsorption of the hydrogen atom on the monolayers. Three possible sites with high symmetry are

Table 2. Magnetic Momenta (in μB) of M2P−1T (M = Zr and Fe) and M2P−2H (M = Ti, Zr, Cr, and Fe) Structures magnetic moment (μB) M P a

Ti2P−2H

Zr2P−1T

Zr2P−2H

Fe2P−1T

Fe2P−2H

0.77 −0.045

0.20 −0.013

0.32 0

2.13 −0.032

1.49 −0.081

The magnetic moment is mainly from M atoms.

physical origin of the magnetic properties of 2D M2P monolayers, the spin-polarized partial density of states (PDOS) of TM and P atoms are calculated (Figures S6a−e). The TM d orbitals are decomposed to in-plane component dσ and out-of-plane component dπ (Figure S6). The Zr2P monolayer in the 1T phase is ferromagnetic. Zr_dπ dominates the spin polarization around the Fermi level, leading to a moment of 0.20 μB per Zr (Figure S6a). The magnetic moment of the P atom is about 0.013 μB in Zr2P−1T and antiparallel to those of metal atoms, indicating that the carrier-mediated double exchange is the dominant mechanism of ferromagnetism.3,53−55 However, its 2H phase is antiferromagnetic. The PDOS of the Zr2P−2H monolayer show that the electrons of the P atom are spin-unpolarized because of the asymmetric magnetic structures (the two magnetic values are equal in the opposite direction) between two Zr atoms. Both Zr_dσ and Zr_dπ make the main contribution to the spin polarization near the Fermi level and the p−d hybridization between P and Zr atoms also played a significant role (Figure S6b). From the antiparallel magnetic moments between two Zr layers, we see superexchange among Zr atoms via nonmagnetic P atoms.3,56,57

Figure 4. Spin-polarized density of states (DOS) (a), and spin-polarization ratio (SPR; b) of Ti2P−2H. Spin-up (c) and spin-down (d) band structures of Ti2P−2H. 8895

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials

Figure 5. Side and top views of the structure of M2P−2H with hydrogen adsorption at HC (M2P−HC; a) and M2P−2H with hydrogen adsorption at TP (M2P−TP; b) for a 3 × 3 × 1 supercell. (c) Calculated H-adsorption energies at different sites. (d) Calculated Gibbs free energies of W2P− HC and Fe2P−TP as a function of H coverage on pure M2P monolayers. d-ΔG: differential Gibbs free energy; a-ΔG: average Gibbs free energy.

where E(M2P + nH) and E(M2P) are the total energies of the M2P monolayer with and without hydrogen atoms, respectively. n is the number of H atoms adsorbed on single-layer M2P. ΔSH is the difference in entropy and T is room temperature. ΔEZPE is the difference in zero-point energies of hydrogen between the adsorbed and the gas phases. ΔEZPE − TΔSH is about 0.24 eV. The Gibbs free energies are defined as differential d-ΔGH and average a-ΔG H . The individual process describes the production of the hydrogen molecule one by one, which can be assessed by d-ΔGH. The collective process shows that all of the hydrogen atoms on the surface are simultaneously converted to molecules, which can be express by a-ΔGH. To investigate the effect of hydrogen coverage on HER performance, we constructed a 3 × 3 × 1 supercell of pure M2P monolayer for the adsorption of n hydrogen atoms with n changed from 1 to 9. We consider the magnetic properties in the study of HER performance for the magnetic system, including Ti2P−2H, Zr2P−1T, Zr2P−2H, and Fe2P−2H. We see that the ΔGH strongly depends on the H coverage and adsorption site (Figures S14−S17). For the hydrogen adsorption on HC (Figure 5d), the 2D W2P−2H monolayer shows the best HER performance due to the relatively lower overpotentials (absolute value of the calculated Gibbs free energy). Specifically, the d-ΔGH of W2P−2H is 0.051, −0.048, 0.115, and −0.114 eV at n = 3, 6, 8, and 9, respectively, indicating high catalytic activity in the individual process. However, the overpotential in the average process is far away from zero, leading to a poor collective process. For hydrogen adsorbed on the TP site, the calculated Gibbs free energies show that Mo2P−2H and Fe2P−2H are better than other predicted 2D TMPs in HER performance because of the relatively lower overpotentials (Figures S16 and S17). Specifically, the d-ΔGH values of Mo2P−2H and Fe2P−2H are −0.065 and −0.080 eV, approximately equal to 0, at full H coverage (n = 9), indicating good HER activity at high H density in the individual process. Similarly, the average Gibbs

considered for hydrogen atoms adsorbed on the M2P monolayer, including the top of the hexagonal center (HC), top of P atoms (TP), and top of M atoms (TM) (Figure 5a,b, and Figure S13). The adsorption energies (Ead) of different H adsorption sites are calculated to find the stable adsorption configuration, as expressed by Ead = E(M2P + H) − E(M2P) − (1/2)E(H2), where E(M2P + H) and E(M2P) are the total energies of M2P monolayers with and without hydrogen adsorption. The calculated energies show that both HC and TP sites on most of the 2D M2P monolayers are stable to hold hydrogen atoms because of the negative adsorption energy, except Fe2P−1T, for which it is difficult to hold hydrogen atoms. The H atom on the TM site moves away to other sites after relaxation, indicating that the top of the metal is unstable to hold the hydrogen atom. HC is the most stable site for hydrogen adsorption because of the relatively lower adsorption energy (Figure 5c). However, hydrogen can also adsorb on TP because of the small adsorption energy difference between the adsorptions on HC and TP sites. The HER performance of the 2D M2P monolayer can be quantified by the reaction Gibbs free energy of hydrogen adsorption (ΔGH).69,70 In principle, an advanced high-HERactivity catalyst should have a ΔGH near 0 eV, which can be obtained from the following equation: ΔG H = ΔE H + ΔEZPE − T ΔSH

(1)

where ΔEH is the hydrogen chemisorption energy. This can be the differential chemisorption energy, defined as ΔE H = E(M 2P + nH) − E(M 2P + (n − 1)H) −

1 E(H 2) 2 (2)

or the average chemisorption energy, defined as ⎛ ⎞ n ΔE H = ⎜E(M 2P + nH) − E(M 2P) − E(H 2)⎟ /n ⎝ ⎠ 2

(3) 8896

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials

Figure 6. Side and top views of the structure of M2PO2−HC (a) and M2PO2-TP (b) for a 3 × 3 × 1 supercell. Calculated differential Gibbs free energies of oxidized M2P monolayers as a function of H coverage: individual (c) and collective (d) processes.

the HER activity of a few 2D TMPs, including Ti2P−2H, Ta2P−2H, and Nb2P−2H, at low H coverage. However, the HER performances of Mo2P−2H and Fe2P−2H are greatly reduced by oxidation because of large overpotentials (Figures S21−S24). In addition, we investigate the effect of oxygen coverage on HER performance. We find that low O coverage leads to reduction of HER performance, and high/full O coverage shows better catalytic ability (Figure S26).

free energies of Mo2P−2H and Fe2P−2H for hydrogen at the TP site are far away from 0 eV, indicating that it is more difficult for the collective processes to take place than the individual process. By comparing the Gibbs free energies at different adsorption sites, we see that Mo2P−2H shows generally better HER performance when hydrogen atoms are at HC sites than when at TP sites, except n = 9, where the situation is reversed. However, the HER performance of Fe2P− 2H is much better for the hydrogen atom at TP sites than for the one at HC sites. Furthermore, the effect of defects in Fe2P− 2H and W2P−2H on HER performance is investigated. We find that their HER performances are slightly reduced by introducing defects (Figure S18). 3.4. HER Activity of Oxidized 2D M2P. It is well-known that MXenes can be easily oxidized because the out-layers are transition metal.71−73 To give a comprehensive understanding on catalytic properties of the 2D TMPs, we investigate the HER performance of oxidized 2D TMPs (M2PO2), where oxygen atoms can adsorb on HC or TP sites (Figure 6a,b) because of the strong binding energy. A 2 × 2 × 1 supercell of the M2PO2 monolayer is used to investigate the effect of hydrogen coverage on the HER performance of the oxidized TMPs. From calculated differential Gibbs free energies, the ΔGH of the M2PO2 monolayer strongly depend also on the H coverage and ΔGH increase with the H coverage in both individual and collective processes. For hydrogen adsorption at the HC site, the d-ΔGH of Nb2PO2−2H is −0.046 eV at the lowest H coverage (n = 1), indicating high HER performance in the individual process at low hydrogen coverage. For hydrogen adsorption at the TP site, the d-ΔGH of Ti2PO2−2H for n = 1, Ta2PO2−2H for n = 1, and Nb2PO2−2H for n = 2 are 0.077, −0.079, and 0.08 eV, respectively, indicating good HER activity at low H density in the individual process. Their HER activities are comparable with those of MXenes and MX2 monolayers when considering the H-coverage density (Table S3).16,42,74−76 Therefore, the oxygen functionalization can efficiently improve

4. CONCLUSIONS In summary, a new family of 2D transition-metal phosphides is designed on the basis of DFT calculations. We predict that these 2D M2P monolayers (M = Ti, Zr, Nb, Ta, Mo, W, and Fe) are thermodynamically stable by the analysis of phonon dispersions, molecular dynamic simulations, and investigation of mechanical properties. We find that all of the stable 2D M2P monolayers are metallic, and the magnetic states strongly depend on the transition metal and phase. We further demonstrate that pure W2P−2H and Fe2P−2H monolayers show the best catalytic performance in the hydrogen evolution reaction at high hydrogen coverage because of the relatively lower overpotentials, and oxygen functionalization can efficiently improve the HER activities of Nb2PO2, Ti2PO2, and Ta2PO2 at low H coverage. It is expected that these 2D M2P monolayers can be applied in spintronics and electrocatalysis.



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.chemmater.7b03832. Optimized lattice constants, phonon dispersion curves, snapshots of atomic configurations at the end of AIMD simulations, spin-polarized density of states, and Gibbs free energies (PDF) 8897

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials



(13) Zhu, J.; Schwingenschlögl, U. P and Si Functionalized Mxenes for Metal-Ion Battery Applications. 2D Mater. 2017, 4, 025073. (14) Jing, Y.; Zhou, Z.; Cabrera, C. R.; Chen, Z. Graphene, Inorganic Graphene Analogs and Their Composites for Lithium Ion Batteries. J. Mater. Chem. A 2014, 2, 12104−12122. (15) Deng, J.; Li, H.; Xiao, J.; et al. Triggering the electrocatalytic hydrogen evolution activity of the inert two-dimensional MoS2 surface viasingle-atom metal doping. Energy Environ. Sci. 2015, 8, 1594−1601. (16) Pan, H. Metal Dichalcogenides Monolayers: Novel Catalysts for Electrochemical Hydrogen Production. Sci. Rep. 2014, 4, 5348. (17) Pan, H. Tension-Enhanced Hydrogen Evolution Reaction on Vanadium Disulfide Monolayer. Nanoscale Res. Lett. 2016, 11, 1−6. (18) Qu, Y.; Pan, H.; Kwok, C. T.; Wang, Z. A First-Principles Study on the Hydrogen Evolution Reaction of VS2 Nanoribbons. Phys. Chem. Chem. Phys. 2015, 17, 24820−24825. (19) Splendiani, A.; Sun, L.; Zhang, Y.; Li, T.; Kim, J.; Chim, C.-Y.; Galli, G.; Wang, F. Emerging Photoluminescence in Monolayer MoS2. Nano Lett. 2010, 10, 1271−1275. (20) Voiry, D.; Yang, J.; Chhowalla, M. Recent Strategies for Improving the Catalytic Activity of 2D TMD Nanosheets toward the Hydrogen Evolution Reaction. Adv. Mater. 2016, 28, 6197−6206. (21) Zeng, Z.; Tan, C.; Huang, X.; Bao, S.; Zhang, H. Growth of Noble Metal Nanoparticles on Single-Layer TiS2 and TaS2 Nanosheets for Hydrogen Evolution Reaction. Energy Environ. Sci. 2014, 7, 797− 803. (22) Zhou, J.; Qin, J.; Guo, L.; Zhao, N.; Shi, C.; Liu, E.; He, F.; Ma, L.; Li, J.; He, C. Scalable Synthesis of High-Quality Transition Metal Dichalcogenide Nanosheets and Their Application as Sodium-Ion Battery Anodes. J. Mater. Chem. A 2016, 4, 17370−17380. (23) Ataca, A.; Ş a hin, H.; Ciraci, S. Stable, Single-Layer MX2Transition-Metal Oxides and Dichalcogenides in a HoneycombLike Structure. J. Phys. Chem. C 2012, 116, 8983−8999. (24) Boota, M.; Pasini, M.; Galeotti, F.; Porzio, W.; Zhao, M.; Halim, J.; Gogotsi, Y. Interaction of Polar and Nonpolar Polyfluorenes with Layers of Two-Dimensional Titanium Carbide (Mxene): Intercalation and Pseudocapacitance. Chem. Mater. 2017, 29, 2731−2738. (25) Hu, J.; Xu, B.; Ouyang, C.; Yang, S.; Yao, Y. Investigations on V2C and V2CX2 (X= F, OH) Monolayer as a Promising Anode Material for Li Ion Batteries from First-Principles Calculations. J. Phys. Chem. C 2014, 118, 24274−24281. (26) Naguib, M.; Mochalin, V.; Barsoum, M. W.; Gogotsi, Y. 25th Anniversary Article: Mxenes: A New Family of Two-Dimensional Materials. Adv. Mater. 2014, 26, 992−1005. (27) Wang, X.; Shen, X.; Gao, Y.; Wang, Z.; Yu, R.; Chen, L. AtomicScale Recognition of Surface Structure and Intercalation Mechanism of Ti3C2X. J. Am. Chem. Soc. 2015, 137, 2715−2721. (28) Lei, J.; Zhang, X.; Zhou, Z. Recent advances in MXene: Preparation, properties, and applications. Front. Phys. 2015, 10, 276− 286. (29) Liu, Y.; Xiao, H.; Goddard, W. A. Schottky-Barrier-Free Contacts with Two-Dimensional Semiconductors by Surface-Engineered Mxenes. J. Am. Chem. Soc. 2016, 138, 15853−15856. (30) Naguib, M.; Kurtoglu, M.; Presser, V.; Lu, J.; Niu, J.; Heon, M.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W. Two-Dimensional Nanocrystals Produced by Exfoliation of Ti3AlC2. Adv. Mater. 2011, 23, 4248−4253. (31) Xu, C.; Wang, L.; Liu, Z.; Chen, L.; Guo, J.; Kang, N.; Ma, X.; Cheng, H.; Ren, W. Large-Area High-Quality 2D Ultrathin Mo2C Superconducting Crystals. Nat. Mater. 2015, 14, 1135−1141. (32) Anasori, B.; Lukatskaya, M.; Gogotsi, Y. 2D Metal Carbides and Nitrides (Mxenes) for Energy Storage. Nat. Rev. Mater. 2017, 2, 16098. (33) Mashtalir, O.; Lukatskaya, M. R.; Zhao, M.; Barsoum, M. W.; Gogotsi, Y. Amine-Assisted Delamination of Nb2C Mxene for Li-Ion Energy Storage Devices. Adv. Mater. 2015, 27, 3501−3506. (34) Pan, H. Electronic Properties and Lithium Storage Capacities of Two-Dimensional Transition-Metal Nitride Monolayers. J. Mater. Chem. A 2015, 3, 21486−21493. (35) Tang, Q.; Zhou, Z.; Shen, P. Are MXenes Promising Anode Materials for Li Ion Batteries? Computational Studies on Electronic

AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. ORCID

Hui Pan: 0000-0002-6515-4970 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The work is supported by the Science and Technology Development Fund from Macau SAR (FDCT-068/2014/A2, FDCT-132/2014/A3, and FDCT-110/2014/SB), the MultiYear Research Grant (MYRG2017-00027-FST and MYRG2015-00017-FST) from the Research & Development Office at University of Macau, the NSF of China (Grants 11474145, and 11334003), the Nanshan Key Lab on Nonvolatile Memory Grant (KC2015ZDYF0003A), and the special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase) under Grant U1501501. The DFT calculations were partially performed at High Performance Computing Cluster (HPCC) of Information and Communication Technology Office (ICTO) at University of Macau.



REFERENCES

(1) Geim, A. K.; Novoselov, K. S. The Rise of Graphene. Nat. Mater. 2007, 6, 183−191. (2) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Katsnelson, M. I.; Grigorieva, I. V.; Dubonos, S. V.; Firsov, A. TwoDimensional Gas of Massless Dirac Fermions in Graphene. Nature 2005, 438, 197−200. (3) Pan, H. Magnetic and Electronic Evolutions of Hydrogenated VTe2 Monolayer under Tension. Sci. Rep. 2014, 4, 7524. (4) Ghidiu, M.; Kota, S.; Halim, J.; Sherwood, A. W.; Nedfors, N.; Rosen, J.; Mochalin, V. N.; Barsoum, M. W. Alkylammonium Cation Intercalation into Ti3C2 (MXene): Effects on Properties and IonExchange Capacity Estimation. Chem. Mater. 2017, 29, 1099−1106. (5) Seok, J.; Lee, J.; Cho, S.; et al. Active Hydrogen Evolution through Lattice Distortion in Metallic MoTe2. 2D Mater. 2017, 4, 025061. (6) Tang, Q.; Jiang, D. Stabilization and Band-Gap Tuning of the 1TMoS2 Monolayer by Covalent Functionalization. Chem. Mater. 2015, 27, 3743−3748. (7) Zhou, M.; Zhang, A.; Dai, Z.; Feng, Y.; Zhang, C. StrainEnhanced Stabilization and Catalytic Activity of Metal Nanoclusters on Graphene. J. Phys. Chem. C 2010, 114, 16541−16546. (8) Rakhi, R. B.; Ahmed, B.; Hedhili, M. N.; Anjum, D. H.; Alshareef, H. N. Effect of Postetch Annealing Gas Composition on the Structural and Electrochemical Properties of Ti2CTx MXene Electrodes for Supercapacitor Applications. Chem. Mater. 2015, 27, 5314−5323. (9) Zhou, M.; Liu, Z.; Ming, W.; Wang, Z.; Liu, F. sd2 Graphene: Kagome Band in a Hexagonal Lattice. Phys. Rev. Lett. 2014, 113, 236802. (10) Lin, Z.; McCreary, A.; Briggs, N.; et al. 2D Materials Advances: From Large Scale Synthesis and Controlled Heterostructures to Improved Characterization Techniques, Defects and Applications. 2D Mater. 2016, 3, 042001. (11) Fang, Y.; Lv, Y.; Gong, F.; Elzatahry, A.; Zheng, G.; Zhao, D. Synthesis of 2D-Mesoporous-Carbon/MoS2 Heterostructures with Well-Defined Interfaces for High-Performance Lithium-Ion Batteries. Adv. Mater. 2016, 28, 9385−9390. (12) Tang, Q.; Jiang, D. Mechanism of hydrogen evolution reaction on 1T-MoS2 from first principles. ACS Catal. 2016, 6, 4953−4961. 8898

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials

(58) Callejas, J. F.; Read, C. G.; Roske, C. W.; Lewis, N. S.; Schaak, R. E. Synthesis, Characterization, and Properties of Metal Phosphide Catalysts for the Hydrogen-Evolution Reaction. Chem. Mater. 2016, 28, 6017−6044. (59) Hu, G.; Tang, Q.; Jiang, D. CoP for hydrogen evolution: implications from hydrogen adsorption. Phys. Chem. Chem. Phys. 2016, 18, 23864−23871. (60) Jiang, P.; Liu, Q.; Liang, Y.; Tian, J.; Asiri, A. M.; Sun, X. A CostEffective 3D Hydrogen Evolution Cathode with High Catalytic Activity: Fep Nanowire Array as the Active Phase. Angew. Chem., Int. Ed. 2014, 53, 12855−12859. (61) Kibsgaard, J.; Jaramillo, T. F. Molybdenum Phosphosulfide: An Active, Acid-Stable, Earth-Abundant Catalyst for the Hydrogen Evolution Reaction. Angew. Chem., Int. Ed. 2014, 53, 14433−14437. (62) McEnaney, J. M; Chance Crompton, J.; Callejas, J. F.; Popczun, E. J.; Read, C. G.; Lewis, N. S.; Schaak, R. E. Electrocatalytic Hydrogen Evolution Using Amorphous Tungsten Phosphide Nanoparticles. Chem. Commun. 2014, 50, 11026−11028. (63) Popczun, E. J.; McKone, J. R.; Read, C. G.; Biacchi, A. J.; Wiltrout, A. M.; Lewis, N. S.; Schaak, R. E. Nanostructured Nickel Phosphide as an Electrocatalyst for the Hydrogen Evolution Reaction. J. Am. Chem. Soc. 2013, 135, 9267−9270. (64) Pu, Z.; Liu, Q.; Tang, C.; Asiri, A. M.; Sun, X. Ni2P Nanoparticle Films Supported on a Ti Plate as an Efficient Hydrogen Evolution Cathode. Nanoscale 2014, 6, 11031−11034. (65) Tian, J.; Liu, Q.; Cheng, N.; Asiri, A. M.; Sun, X. Self-Supported Cu3P Nanowire Arrays as an Integrated High-Performance ThreeDimensional Cathode for Generating Hydrogen from Water. Angew. Chem., Int. Ed. 2014, 53, 9577−9581. (66) Tian, L.; Yan, X.; Chen, X. Electrochemical Activity of Iron Phosphide Nanoparticles in Hydrogen Evolution Reaction. ACS Catal. 2016, 6, 5441−5448. (67) Xu, Y.; Wu, R.; Zhang, J.; Shi, Y.; Zhang, B. Anion-Exchange Synthesis of Nanoporous FeP Nanosheets as Electrocatalysts for Hydrogen Evolution Reaction. Chem. Commun. 2013, 49, 6656−6658. (68) Yang, H.; Zhang, Y.; Hu, F.; Wang, Q. Urchin-Like Cop Nanocrystals as Hydrogen Evolution Reaction and Oxygen Reduction Reaction Dual-Electrocatalyst with Superior Stability. Nano Lett. 2015, 15, 7616−7620. (69) Nørskov, J. K.; Bligaard, T.; Logadottir, A.; Kitchin, J. R.; Chen, J. G.; Pandelov, S.; Stimming, U. Trends in the Exchange Current for Hydrogen Evolution. J. Electrochem. Soc. 2005, 152, 23−26. (70) Voiry, D.; Salehi, M.; Silva, R.; Fujita, T.; Chen, M. W.; Asefa, T.; Shenoy, V. B.; Eda, G.; Chhowalla, M. Conducting MoS2 nanosheets as catalysts for hydrogen evolution reaction. Nano Lett. 2013, 13, 6222−6227. (71) Hu, M.; Li, Z.; Hu, T.; Zhu, S.; Zhang, C.; Wang, X. HighCapacitance Mechanism for Ti3C2Tx Mxene by in Situ Electrochemical Raman Spectroscopy Investigation. ACS Nano 2016, 10, 11344− 11350. (72) Naguib, M.; Mashtalir, O.; Lukatskaya, M. R.; Dyatkin, B.; Zhang, C.; Presser, V.; Gogotsi, Y.; Barsoum, M. W. One-Step Synthesis of Nanocrystalline Transition Metal Oxides on Thin Sheets of Disordered Graphitic Carbon by Oxidation of MXenes. Chem. Commun. 2014, 50, 7420−7423. (73) Wang, H.; Naguib, M.; Page, K.; Wesolowski, D.; Gogotsi, Y. Resolving the Structure of Ti3C2Tx MXenes through Multilevel Structural Modeling of the Atomic Pair Distribution Function. Chem. Mater. 2016, 28, 349−359. (74) Gao, G.; O’Mullane, A. P.; Du, A. 2D MXenes: A new family of Promising Catalysts for the Hydrogen Evolution Reaction. ACS Catal. 2017, 7, 494−500. (75) Chou, S.; Sai, N.; Lu, P.; Coker, E.; Liu, S.; Artyushkova, K.; Luk, T.; Kaehr, B.; Brinker, C. Understanding Catalysis in a Multiphasic Two-dimensional Transition Metal Dichalcogenide. Nat. Commun. 2015, 6, 8311. (76) Voiry, D.; Yamaguchi, H.; Li, J.; Silva, R.; Alves, D.; Fujita, T.; Chen, M.; Asefa, T.; Shenoy, V.; Eda, G.; Chhowalla, M. Enhanced

Properties and Li Storage Capability of Ti3C2 and Ti3C2X2 (X = F, OH) Monolayer. J. Am. Chem. Soc. 2012, 134, 16909−16916. (36) Xie, Y.; Naguib, M.; Mochalin, V.; Barsoum, M. W.; Gogotsi, Y.; Yu, X.; Nam, K. W.; Yang, K. Q.; Kolesnikov, A.; Kent, P. Role of Surface Structure on Li-Ion Energy Storage Capacity of TwoDimensional Transition-Metal Carbides. J. Am. Chem. Soc. 2014, 136, 6385−6394. (37) Zhang, C.; Kim, S.; Ghidiu, M.; Zhao, M.; Barsoum, M. W.; Nicolosi, V.; Gogotsi, Y. Layered Orthorhombic Nb2O5@Nb4C3Tx and TiO2@Ti3C2Tx Hierarchical Composites for High Performance Li-Ion Batteries. Adv. Funct. Mater. 2016, 26, 4143−4151. (38) Ling, C.; Shi, L.; Ouyang, Y.; Chen, Q.; Wang, J. Transition Metal-Promoted V2CO2 (MXenes): A New and Highly Active Catalyst for Hydrogen Evolution Reaction. Adv. Sci. 2016, 3, 1600180. (39) Shahzad, F.; Alhabeb, M.; Hatter, C. B.; Anasori, B.; Man Hong, S.; Koo, C. M.; Gogotsi, Y. Electromagnetic Interference Shielding with 2D Transition Metal Carbides (Mxenes). Science 2016, 353, 1137. (40) Samorì, P.; Palermo, V.; Feng, X. Chemical Approaches to 2D Materials. Adv. Mater. 2016, 28, 6027−6029. (41) Karvonen, L.; Säynätjoki, A.; Huttunen, M.; et al. Rapid visualization of grain boundaries in monolayer MoS2 by multiphoton microscopy. Nat. Commun. 2017, 8, 15714. (42) Pan, H. Ultra-High Electrochemical Catalytic Activity of MXenes. Sci. Rep. 2016, 6, 32531. (43) Ma, Y.; Kuc, A.; Heine, T. Single-Layer Tl2O: A Metal-Shrouded 2D Semiconductor with High Electronic Mobility. J. Am. Chem. Soc. 2017, 139, 11694. (44) Ming, W.; Yoon, M.; Du, M.-H.; Lee, K.; Kim, S. W. FirstPrinciples Prediction of Thermodynamically Stable Two-Dimensional Electrides. J. Am. Chem. Soc. 2016, 138, 15336−15344. (45) Hafner, J. Ab-Initio Simulations of Materials Using Vasp: Density-Functional Theory and Beyond. J. Comput. Chem. 2008, 29, 2044−2078. (46) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169. (47) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758. (48) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. (49) Togo, A.; Tanaka, I. First Principles Phonon Calculations in Materials Science. Scr. Mater. 2015, 108, 1−5. (50) Ding, Y.; Wang, Y. Density Functional Theory Study of the Silicene-Like SiX and XSi3 (X = B, C, N, Al, P) Honeycomb Lattices: The Various Buckled Structures and Versatile Electronic Properties. J. Phys. Chem. C 2013, 117, 18266−18278. (51) Andrew, R. C.; Mapasha, R. E.; Ukpong, A. M.; Chetty, N. Mechanical Properties of Graphene and Boronitrene. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 125428. (52) Qu, Y.; Pan, H.; Kwok, C. T. Hydrogenation-controlled phase transition on two-dimensional transition metal dichalcogenides and their unique physical and catalytic properties. Sci. Rep. 2016, 6, 34186. (53) Akai, H. Ferromagnetism and Its Stability in the Diluted Magnetic Semiconductor (in, Mn)As. Phys. Rev. Lett. 1998, 81, 3002− 3005. (54) Behan, A. J.; Mokhtari, A.; Blythe, H. J.; Score, D.; Xu, X. H.; Neal, J. R.; Fox, A. M.; Gehring, G. A. Two Magnetic Regimes in Doped Zno Corresponding to a Dilute Magnetic Semiconductor and a Dilute Magnetic Insulator. Phys. Rev. Lett. 2008, 100, 047206. (55) Pan, H.; Feng, Y. P.; Wu, Q. Y.; Huang, Z. G.; Lin, J. Magnetic Properties of Carbon Doped CdS: A First-Principles and Monte Carlo Study. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 125211. (56) Anderson, P. W. Antiferromagnetism. Theory of Superexchange Interaction. Phys. Rev. 1950, 79, 350−356. (57) Goodenough, J. B. Theory of the Role of Covalence in the Perovskite-Type Manganites [La, M(II)]MnO3. Phys. Rev. 1955, 100, 564−573. 8899

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900

Article

Chemistry of Materials Catalytic Activity in Strained Chemically Exfoliated WS2 Nanosheets for Hydrogen Evolution. Nat. Mater. 2013, 12, 850.

8900

DOI: 10.1021/acs.chemmater.7b03832 Chem. Mater. 2017, 29, 8892−8900