Computational Screening of MXene Electrodes for Pseudocapacitive

Dec 5, 2018 - ‡Computational Sciences and Engineering Division and §Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak ...
1 downloads 0 Views 1MB Size
Subscriber access provided by AUSTRALIAN NATIONAL UNIV

C: Energy Conversion and Storage; Energy and Charge Transport

Computational Screening of MXene Electrodes for Pseudocapacitive Energy Storage Cheng Zhan, Weiwei Sun, Paul R.C. Kent, Michael Naguib, Yury Gogotsi, and De-en Jiang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b11608 • Publication Date (Web): 05 Dec 2018 Downloaded from http://pubs.acs.org on December 6, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Computational Screening of MXene Electrodes for Pseudocapacitive Energy Storage Cheng Zhan,1 Weiwei Sun,2 Paul R. C. Kent,2,3 Michael Naguib,4 Yury Gogotsi,5 and De-en Jiang1,* 1Department 2Computational

of Chemistry, University of California, Riverside, California, 92521, United States

Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States

3Center

for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States

4Department

of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States

5Department

of Materials Science and Engineering and A. J. Drexel Nanomaterials Institute, Drexel University, Philadelphia, PA, 19104, United States

*Corresponding author. Email: [email protected]

Abstract: MXenes (two-dimensional transition metal carbides and nitrides) are promising materials for capacitive energy storage due to the large chemical space of existing and potential compositions, but only a few of them have been experimentally explored. In this work, we computationally screen a series of MXene electrodes (Mn+1XnTx: M=Sc, Ti, V, Zr, Nb, Mo; X=C, N; T=O, OH; n=1-3) to simulate their pseudocapacitive performance in the aqueous H2SO4 electrolyte. We find that nitride MXenes exhibit better pseudocapacitive performance than carbide MXenes. Especially, Ti2NTx is predicted to have a high gravimetric capacitance over a wide voltage window, while Zrn+1NnTx MXenes are predicted to possess the best areal capacitive performance. Evaluating the descriptors for the capacitance trends, we find that more positive hydrogen adsorption free energy (weak binding to H) and smaller change of the potential at the point of zero charge after H binding lead to higher capacitance. Our work provides helpful guidance to selectively develop high-performance MXene pseudocapacitors and to further screen MXene electrodes. Keywords: MXene, pseudocapacitance, density functional theory, screening, descriptor

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 15

Introduction MXenes are a large family of two-dimensional (2D) materials including transition metal carbides, nitrides, and carbonitrides.1 The chemical formula of MXene can be expressed as Mn+1XnTx, where M is a transition metal (Sc, Ti, V, Zr, Nb, Mo, etc.), X is C or N, and T represents the surface terminal groups (such as –O, –OH, and –F). MXenes have shown promise in capacitive energy storage,2-4 lithium and sodium ion batteries,5-7 electrocatalysis,8-9 gas sensors,10-11 and many other applications. MXenes’ performance in capacitive energy storage is most promising. For example, Ti3C2Tx has shown high volumetric capacitance and gravimetric capacitance with excellent cyclability.1213

MXene exhibits pseudocapacitance in the aqueous H2SO4 electrolyte due to the reversible

surface redox reaction of hydrogen binding. In aqueous salt or ionic liquid electrolytes, MXene shows mainly the electric double-layer (EDL) capacitance.13-16 In organic electrolytes, MXene can have intercalation pseudocapacitance, especially in the presence of Li+ or Na+.17-19 To achieve high capacitance for MXene electrodes, H2SO4 solution is the most promising electrolyte due to the its contribution surface redox reaction. In Table S1, we listed the reported capacitances of MXene electrodes in comparison with other conventional pseudocapacitor materials, such as RuO2,20-22 poly-pyrrole,23 and benzoquinone.24 One can see that MXene as a pseudocapacitor not only has promising capacitance, but also exhibits much better cyclability. Previous studies have focused on Ti3C2Tx,13-14, 16, 25 in addition to a few other MXenes such as Mo2CTx26 and Mo2TiC2Tx,27 despite the diversity of the large MXene family.1 Because MXene synthesis has many technical challenges, experimentally exploring the capacitive performance of the complete MXene family would be very time-consuming. Thus, exploring the pseudocapacitive performance of MXene from a computational perspective is an effective way to help experimentalists design supercapacitors with improved energy storage ability. In this work, we computationally screen the energy storage of Mn+1XnTx MXenes (Mn+1XnTx: M=Sc, Ti, V, Zr, Nb, Mo; X=C, N; T=O, OH; n=1, 2, 3 ) in H2SO4, based on a recently developed model that incorporates the surface-redox pseudocapacitance in the total capacitive energy storage of MXene in H2SO4.28 Moreover, we will explore descriptors that can correlate the capacitances of MXenes with some easily computable and relevant quantities.

ACS Paragon Plus Environment

2

Page 3 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Results and discussion The structure of M2XT2 with the solvation layers is shown in Figure 1 as an example.29 The structure of M3X2T2 and M4X3T2 can be derived from M2X by extending the M and X layers following the ABC stacking order. The termination group T can occupy either the fcc site (Figure 1a,b) or the hcp site (Figure 1c,d). We calculated the relative stability of fcc and hcp configurations for each MXene (Table S2). Since T can be either O or OH, the MXene structures for the redox pair Mn+1XnO2 and Mn+1Xn(OH)2 are selected by the following rules: (1) If both –O and –OH are more stable at the fcc (hcp) site, we use the fcc (hcp) model; (2) if one terminal group is more stable at the fcc site, while the other is more stable at the hcp site, we use the site for which the group has stronger relative stability. This consistent choice of either an fcc or hcp model ensures that we have a most likely surface for reversible hydrogen binding at different coverages, a key step for our prediction of the capacitance.

Figure 1. Model of the M2XT2 MXene with the terminal group T at two different sites: (a) and (b), side and top views of the T groups at the fcc sites (T not directly above M); (c) and (d), side and top views of the T groups at the hcp sites (T directly above M). Color scheme: Red, terminal group T; cyan, transition metal M; gray: carbon or nitrogen X.

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 15

Our first screening criteria are the electronic structure and magnetic property. The results from density functional theory (DFT) calculation are listed in Table S2. To be a high-rate electrode material, metallic MXenes are preferred. Hence, Sc2C(OH)2, Ti2CO2, and Zr2CO2, which are semiconducting (consistent with a previous study30), can be eliminated. In addition, to avoid complication from the magnetic properties, which may be both termination and coverage dependent, we leave out V3C2O2, Sc2NO2, V2NO2, V3N2O2 and V4N3(OH)2 for a future study. Since the pseudocapacitive performance of MXene in the aqueous H2SO4 electrolyte competes with hydrogen evolution reaction (HER) at negative potentials,8 it is necessary to screen MXenes to eliminate those that are too active for HER. The differential hydrogen adsorption free energy is a good descriptor to assess the HER activity.31 We calculated the average hydrogen adsorption free energy (ΔGH, which shows the overall binding trend) and the differential hydrogen adsorption free energy at the optimal H coverage (ΔGHOPT, which reflects the HER activity). From Figure 2, one can see that Sc-based MXenes (Sc2C, Sc3C2, Sc4C3, Sc3N2 and Sc4N3) exhibit relatively stronger binding ability to H atom, with GH close to or less than zero. So we eliminated them for consideration.

Figure 2. Hydrogen adsorption free energy of 2D MXenes: ΔGH , average hydrogen adsorption free energy; ΔGHOPT, the hydrogen adsorption free energy at optimal hydrogen coverage.

ACS Paragon Plus Environment

4

Page 5 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

After our quick screening on electronic structure, magnetic property, and hydrogen binding ability, 24 candidates have been identified. Following our previous work,28 we predicted the capacitive energy storage in the aqueous H2SO4 electrolyte for the selected MXene candidates. Briefly, we computed the free energy of a solvated MXene electrode at different applied potentials and H coverages as well as determined the potentials at the point of zero charge before and after H binding. These values then allowed us to determine the stored charge (Q) as a function of the applied potential (V). From these curves, we then derived the differential and integral capacitances. The details are described in the Methods section. From the simulated gravimetric or specific differential capacitances for carbides (Figure 3a) and nitrides (Figure 3b), one can see that most MXenes exhibit a roughly flat differential capacitance curve across our simulated voltage range, except V2C, Ti2N, and Zr2N. If we use a common voltage window of -0.5 to 0.5 V, we can obtain the specific (Figure 4a,b) and areal (Figure 4c,d) integral capacitances. In the carbide MXenes (Figure 4a,c), V2C exhibit highest specific capacitance of 350 F/g, while Zr4C3 exhibit highest areal capacitance that over 60 μF/cm2. In the nitride MXenes (Figure 4b,d), Ti2N exhibit highest gravimetric capacitance over 450 F/g. In all MXene candidates, Ti2N has the highest gravimetric capacitance due to both Ti’s low atomic weight and the favorable redox chemistry (discussed later), while Zr2N, Zr3N2, and Zr3N2 exhibit the highest and similar areal capacitances. Since experimentally synthesized MXenes have a multilayer structure, a large volumetric capacitance is expected for MXenes that have large areal capacitances. Comparing carbides with nitrides of the same metal and stoichiometry, we found that nitride MXenes have better capacitive performance than carbide MXenes for most cases. Due to the difficulty in the nitride MXene synthesis, so far only Ti4N3 and Ti2N were successfully synthesized.32-33 Thus, there is no available experimental data to directly compare with our predictions. However, there have been recent experimental findings of other nitride pseudocapacitors. For instance, Morel et al. reported the pseudocapacitance of VN34 and Xiao et al. found the promising pseudocapacitive performance of 2D metallic MoN.35 Moreover, Wen and Yang et al. found that the nitrogen doping in Ti3C2Tx can increase the pseudocapacitance.36 These recent progress suggests that nitride MXenes may be available soon for testing our predictions. In terms of structural stability, a recent DFT study analyzed the formation energies of 72 different MXenes and concluded that all the -O and -OH functionalized carbide and nitride MXenes have very negative formation energies, indicating their thermodynamic stability.37

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 15

Figure 3. Predicted differential capacitances of MXene electrodes: (a) specific capacitance of carbide MXenes; (b) specific capacitance of nitride MXenes; (c) areal capacitance of carbide MXenes; (d) areal capacitance of nitride MXenes.

Table 1. Predicted and experimental capacitances (F/g) of MXene in 1M H2SO4. MXene

Theory

Experiment

Ti3C2Tx

230a

235,b 245c

Mo2CTx

150d

190e

V2CTx

350d

380f

aRef.

28; bRef. 14; cRef. 16; dPresent work; eRef. 26; fRef. 38.

The reliability of our model and prediction for the capacitance of Ti3C2Tx in H2SO4 has been discussed previously28 and is summarized in Table 1. In addition, we listed recently measured experimental capacitances for Mo2CTx and V2CTx in H2SO4 for comparison. The reported experimental capacitance of Mo2CTx in H2SO4 at 190 F/g26 also agrees reasonably well with our predicted capacitance of 150 F/g at the same voltage window. A more recent experimental work

ACS Paragon Plus Environment

6

Page 7 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

reported the capacitance of V2C in 1M H2SO4 to be ~380 F/g at the scan rate of 5 mV/s,38 in good agreement with our predicted capacitance of 350 F/g. For other MXenes in Figure 4, no experimental data are available, so our calculations indicate promising MXenes to try.

Figure 4. Predicted integral capacitances of MXene electrodes: (a) specific capacitance of carbide MXenes; (b) specific capacitance of nitride MXenes; (c) areal capacitance of carbide MXenes; (d) areal capacitance of nitride MXenes. The voltage window is from -0.5 V to +0.5 V vs SHE.

To understand what determines the capacitive performance of MXenes, we studied the factors that can be correlated to the charge-storage trend, in order to find a descriptor of pseudocapacitance. As shown in Figure 5, we found that the change or shift of the potential at the point-of-zero-charge (VPZC) from Mn+1XnO2 to Mn+1Xn(OH)2 is closely related to the amount of charge stored: the charge storage (blue dot; right axis) is higher if the VPZC change (gray bar) is smaller. We use ΔVPZC, defined as one half of the VPZC change from Mn+1XnO2 to Mn+1Xn(OH)2, to represent the VPZC shift. The charge stored can be expressed as a function of ΔGH and ΔVPZC as descriptors (see the Methods section). Briefly, ΔGH is computed by hydrogen adsorption simulation via DFT and ΔVPZC is calculated from the Fermi levels of MXene in both reduced and oxidized forms with a continuum solvation model. As shown in Figure S1, the agreement between this descriptor approach and the DFT results (Figure 4) is very good. One can see from the 2D color map (Figure 6) that large ΔGH

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 15

and low ΔVpzc lead to higher charge storage per unit of formula. Placing the MXene candidates onto this map, one can see that Zr-based nitride MXenes (Zr2N, Zr3N2 and Zr4N3) show the highest charge storage per unit of formula, consistent with the areal capacitance trends in Figure 4.

Figure 5. Shift in the potential at the point-of-zero-charge (VPZC) from Mn+1XnO2 to the reduced Mn+1Xn(OH)2 form (bars; left axis) and charge storage per formula unit (blue dot; right axis) of 24 MXene electrode candidates.

Figure 6. Charge storage per formula unit vs the shift in the potential at the point-of-zero-charge (ΔVPZC) and hydrogen adsorption free energy (ΔGH). Voltage window is from -1V to 1V vs SHE. Values (in (|e| per Mn+1XnTx) for the contour lines are given at right.

ACS Paragon Plus Environment

8

Page 9 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The correlation of ΔVpzc with the redox charge storage of MXenes can be understood from the Fermi-level shift from Mn+1XnO2 to Mn+1Xn(OH)2, as ΔVpzc is proportional to the Fermi-level shift. For MXenes with higher capacitances, it takes smaller shifts in the Fermi level to have the complete two-electron transfer; in other words, the more metallic the MXenes are, the higher the capacitances. Further test of this insight in more MXenes such as bimetallic, doped, and defected ones is warranted. To understand why most nitride MXenes have higher capacitances than their carbide counterparts, we compare their Fermi-level shifts from Mn+1XnO2 to Mn+1Xn(OH)2 in Figure S2. We find that changing X from C to N for the same Mn+1XnO2 MXene lowers ΔVpzc, thereby making the system more metallic. This explains the nitrides’ higher capacitances. Summary and conclusions In summary, we have conducted a computational screening of the pseudocapacitive Mn+1XnTx MXene electrodes in the H2SO4 electrolyte. 24 MXene electrode candidates were selected by considering electronic structure, magnetism, and hydrogen adsorption energy calculated via DFT. Simulated pseudocapacitive charging curves showed that most of the candidates exhibited better capacitive performance than Ti3C2Tx, a prototypical MXene comprehensively reported before. Our screening indicated that the nitrides tends to have better capacitive performance than carbides: Zr2N, Zr3N2, and Zr4N3 were predicted to have highest charge storage per unit formula and Ti2N to have the highest specific capacitance. Moreover, the hydrogen adsorption energy (ΔGH) and the shift in the potential at PZC (ΔVPZC) were found to be the key quantities that dominate MXenes’ charge storage. A promising MXene pseudocapacitor should have large ΔGH and low ΔVPZC. This rule can be applied to pseudocapacitive MXene electrode design and high-throughput screening in the future. Methods Geometry optimization and electronic structure calculations with DFT were performed using the Quantum Espresso package.39 The Perdew–Burke-Ernzerhof functional (GGA-PBE) of generalized-gradient approximation was used for electron exchange and correlation.40 Ultrasoft pseudopotentials were used to describe the nuclei-electron interaction.41 A plane wave basis with the energy cutoff of 40 Ry and a k-mesh of 7 × 7 × 1 (for sampling the Brillouin zone) were used for structure optimizations, and a higher cutoff of 60 Ry and a k-mesh of 21 × 21 × 1 were used

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 15

for single-point calculations of optimized geometries. Spin polarization was applied to test the magnetic property of MXenes. To assess the hydrogen binding ability of MXene, we calculated the average and differential hydrogen adsorption free energy. The average hydrogen adsorption free energy is defined by: 1

∆𝐺𝐻 = 2(𝐺(𝑀𝑛 + 1𝑋𝑛(𝑂𝐻)2) ―𝐺(𝑀𝑛 + 1𝑋𝑛𝑂2) ―𝐺(𝐻2))

(1)

The free energy difference between Mn+1Xn(OH)2 and Mn+1XnO2 was obtain the following: 𝐺(𝑀𝑛 + 1𝑋𝑛(𝑂𝐻)2) ―𝐺(𝑀𝑛 + 1𝑋𝑛𝑂2) = 𝐸(𝑀𝑛 + 1𝑋𝑛(𝑂𝐻)2) ―𝐸(𝑀𝑛 + 1𝑋𝑛𝑂2) +∆𝑍𝑃𝐸 (2) The free energy of a hydrogen molecule is defined as: 7

𝐺(𝐻2) = 𝐸(𝐻2) +𝑍𝑃𝐸(𝐻2) + 2𝑘𝐵𝑇 ― 𝑇𝑆𝐻2

(3)

All quantities in equation (3) can be obtained from DFT and standard thermodynamic database. To assess the HER activity of MXene, we also computed the differential hydrogen adsorption free energy at coverages of 0.25 ML, 0.5 ML, 0.75 ML and 1.00 ML. When the calculated differential hydrogen adsorption free energy is closest to 0 eV, the corresponding coverage is treated as the optimal coverage. For the energy and PZC calculation of a solvated electrode, we performed DFT calculation with an implicit solvation model implemented in the JDFTx code.42-43 Identical parameters and settings were used for the electronic DFT part as used in the single-point Quantum Espresso calculations. For the solvation part, we used the charge-asymmetric nonlocally determined localelectric (CANDLE) model, which has shown excellent reliability in predicting the solvation energy and PZC of metal electrodes.44 To simulate the ionic strength of 1M H2SO4, our implicit solvation model was based on a 3M monovalent aqueous solution. The pseudocapacitive charging simulation was carried out by following our previously proposed method;28 namely, the free energy function of any intermediate state (partially oxidized/reduced MXene layer) at a fixed electrode potential 𝜑 can be expressed as: 𝐺(𝑥,𝜑) = 𝐸(𝑥) +𝑥𝐸𝑍𝑃𝐸 +𝑄(𝑉(𝑥,𝜑))𝜑 + 𝐸𝐸𝐷𝐿(𝑉(𝑥,𝜑)) +(1 ― 𝑥)𝜇𝐻 +

,

(4)

where E(x) is the total electronic energy of a solvated electrode with the H coverage of x in zero surface charge, while EZPE is the zero-point energy (ZPE) difference between x=0 and x=1 states. We use the ZPE of M2XTx to represent the ZPE of M3X2Tx and M4X3Tx if they share the same “M” and “X”. The ZPE contributions calculated by DFT for different M2XTx systems were listed in Table S3. 𝑄𝜑 is the electrical work to move the charge Q (net charge on the electrode) from

ACS Paragon Plus Environment

10

Page 11 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

zero potential (in the bulk electrolyte) to the electrode with the potential 𝜑. EEDL is the energy electric double layer contributed by the electrode charge Q. The last term is the chemical potential of a solvated proton in the electrolyte. In the third and fourth terms, V(x, 𝜑) is the relative potential with respect to the potential at the PZC at coverage x and the electrode potential 𝜑, given by: 𝑉(𝑥,𝜑) = 𝜑 ― 𝜑𝑃𝑍𝐶(𝑥)

(5)

Since our previous study has proved that the E(x) and 𝜑𝑃𝑍𝐶(𝑥) can be approximately described by linear equations, in this work, we only performed DFT calculation on Mn+1XnO2 (x=0) and Mn+1Xn(OH)2 (x=1). 𝐺(𝑥,𝜑) and 𝜑𝑃𝑍𝐶(𝑥) values for x between 0 and 1 were derived from linear interpolation. More technical details of our method can be found in our previous work.28 One can also obtain 𝜑𝑃𝑍𝐶(𝑥) from: 𝜑𝑃𝑍𝐶(𝑥) = 𝐸0 +∆𝑉𝑃𝑍𝐶 ―𝑥 ∗ [(𝐸0 + ∆𝑉𝑝𝑧𝑐) ― (𝐸0 ― ∆𝑉𝑃𝑍𝐶)] ,

(6)

where E0 is the standard redox potential calculated from ∆𝐺𝐻 = ―𝐹𝐸0 and ∆𝑉𝑃𝑍𝐶 is one half of the 𝜑𝑃𝑍𝐶 change from Mn+1XnO2 to Mn+1Xn(OH)2. ΔGH and ΔVPZC are descriptors that can be directly obtained from DFT. With the estimated values of 𝜑𝑃𝑍𝐶(𝑥), we can then use Eqs. (4) and (5) to derive the charge stored in the MXene electrode at an applied voltage. Acknowledgements This research is sponsored by the Fluid Interface Reactions, Structures, and Transport (FIRST) Center, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC0205CH11231. Supporting Information Reported performance of conventional pseudocapacitors; lattice parameters and phase stability/property for MXenes examined; zero-point-energy (ZPE) differences between M2XO2 and M2X(OH)2 calculated by DFT; differential hydrogen adsorption free energy of MXene electrode candidates at optimal hydrogen coverage; comparison of charge storage predicted by the

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 15

descriptor approach and DFT; definition of electron chemical potential (Fermi energy) in our DFT calculation and its relation to electrode potential. References 1. Anasori, B.; Lukatskaya, M. R. Gogotsi, Y. 2D Metal Carbides and Nitrides (MXenes) for Energy Storage. Nat. Rev. Mater. 2017, 2, 16098. 2. Naguib, M.; Kurtoglu, M.; Presser, V.; Lu, J.; Niu, J. J.; Heon, M.; Hultman, L.; Gogotsi, Y. Barsoum, M. W. Two-Dimensional Nanocrystals Produced by Exfoliation of Ti3AlC2. Adv. Mater. 2011, 23, 4248-4253. 3. Zhao, M. Q.; Ren, C. E.; Ling, Z.; Lukatskaya, M. R.; Zhang, C. F.; Van Aken, K. L.; Barsoum, M. W. Gogotsi, Y. Flexible Mxene/Carbon Nanotube Composite Paper with High Volumetric Capacitance. Adv. Mater. 2015, 27, 339-345. 4. Ling, Z.; Ren, C. E.; Zhao, M. Q.; Yang, J.; Giammarco, J. M.; Qiu, J. S.; Barsoum, M. W. Gogotsi, Y. Flexible and Conductive Mxene Films and Nanocomposites with High Capacitance. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 16676-16681. 5. Naguib, M.; Halim, J.; Lu, J.; Cook, K. M.; Hultman, L.; Gogotsi, Y. Barsoum, M. W. New TwoDimensional Niobium and Vanadium Carbides as Promising Materials for Li-Ion Batteries. J. Am. Chem. Soc. 2013, 135, 15966-15969. 6. Naguib, M.; Come, J.; Dyatkin, B.; Presser, V.; Taberna, P. L.; Simon, P.; Barsoum, M. W. Gogotsi, Y. Mxene: A Promising Transition Metal Carbide Anode for Lithium-Ion Batteries. Electrochem. Commun. 2012, 16, 61-64. 7. Er, D. Q.; Li, J. W.; Naguib, M.; Gogotsi, Y. Shenoy, V. B. Ti3C2 MXene as a High Capacity Electrode Material for Metal (Li, Na, K, Ca) Ion Batteries. ACS Appl. Mater. Interfaces 2014, 6, 1117311179. 8. Seh, Z. W.; Fredrickson, K. D.; Anasori, B.; Kibsgaard, J.; Strickler, A. L.; Lukatskaya, M. R.; Gogotsi, Y.; Jaramillo, T. F. Vojvodic, A. Two-Dimensional Molybdenum Carbide (MXene) as an Efficient Electrocatalyst for Hydrogen Evolution. ACS Energy Lett. 2016, 1, 589-594. 9. Ma, T. Y.; Cao, J. L.; Jaroniec, M. Qiao, S. Z. Interacting Carbon Nitride and Titanium Carbide Nanosheets for High-Performance Oxygen Evolution. Angew. Chem., Int. Ed. 2016, 55, 1138-1142. 10. Yu, X. F.; Li, Y. C.; Cheng, J. B.; Liu, Z. B.; Li, Q. Z.; Li, W. Z.; Yang, X. Xiao, B. Monolayer Ti2co2: A Promising Candidate for NH3 Sensor or Capturer with High Sensitivity and Selectivity. ACS Appl. Mater. Interfaces 2015, 7, 13707-13713. 11. Chen, J.; Chen, K.; Tong, D. Y.; Huang, Y. J.; Zhang, J. W.; Xue, J. M.; Huang, Q. Chen, T. CO2 and Temperature Dual Responsive "Smart'' MXene Phases. Chem. Commun. 2015, 51, 314-317. 12. Lukatskaya, M. R.; Mashtalir, O.; Ren, C. E.; Dall'Agnese, Y.; Rozier, P.; Taberna, P. L.; Naguib, M.; Simon, P.; Barsoum, M. W. Gogotsi, Y. Cation Intercalation and High Volumetric Capacitance of TwoDimensional Titanium Carbide. Science 2013, 341, 1502-1505. 13. Lukatskaya, M. R.; Kota, S.; Lin, Z.; Zhao, M.; Shpigel, N.; Levi, M. D.; Halim, J.; Taberna, P. L.; Barsoum, M. W.; Simon, P., et al. Ultra-High-Rate Pseudocapacitive Energy Storage in Two-Dimensional Transition Metal Carbides. Nat. Energy 2017, 2, 17105. 14. Lukatskaya, M. R.; Bak, S. M.; Yu, X. Q.; Yang, X. Q.; Barsoum, M. W. Gogotsi, Y. Probing the Mechanism of High Capacitance in 2D Titanium Carbide Using in Situ X-Ray Absorption Spectroscopy. Adv. Energy Mater. 2015, 5, 1500589. 15. Wang, X.; Kajiyama, S.; Iinuma, H.; Hosono, E.; Oro, S.; Moriguchi, I.; Okubo, M. Yamada, A. Pseudocapacitance of MXene Nanosheets for High-Power Sodium-Ion Hybrid Capacitors. Nat. Commun. 2015, 6, 6544. 16. Ghidiu, M.; Lukatskaya, M. R.; Zhao, M. Q.; Gogotsi, Y. Barsoum, M. W. Conductive TwoDimensional Titanium Carbide 'Clay' with High Volumetric Capacitance. Nature 2014, 516, 78. 17. Mashtalir, O.; Naguib, M.; Mochalin, V. N.; Dall'Agnese, Y.; Heon, M.; Barsoum, M. W. Gogotsi, Y. Intercalation and Delamination of Layered Carbides and Carbonitrides. Nat. Commun. 2013, 4, 1716.

ACS Paragon Plus Environment

12

Page 13 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

18. Mashtalir, O.; Lukatskaya, M. R.; Zhao, M. Q.; Barsoum, M. W. Gogotsi, Y. Amine-Assisted Delamination of Nb2C Mxene for Li-Ion Energy Storage Devices. Adv. Mater. 2015, 27, 3501-3506. 19. Naguib, M.; Unocic, R. R.; Armstrong, B. L. Nanda, J. Large-Scale Delamination of Multi-Layers Transition Metal Carbides and Carbonitrides "MXenes". Dalton Trans. 2015, 44, 9353-9358. 20. Bi, R.-R.; Wu, X.-L.; Cao, F.-F.; Jiang, L.-Y.; Guo, Y.-G. Wan, L.-J. Highly Dispersed RuO2 Nanoparticles on Carbon Nanotubes: Facile Synthesis and Enhanced Supercapacitance Performance. J. Phys. Chem. C 2010, 114, 2448-2451. 21. Amir, F.; Pham, V. Dickerson, J. Facile Synthesis of Ultra-Small Ruthenium Oxide Nanoparticles Anchored on Reduced Graphene Oxide Nanosheets for High-Performance Supercapacitors. RSC Adv. 2015, 5, 67638-67645. 22. Kuratani, K.; Kiyobayashi, T. Kuriyama, N. Influence of the Mesoporous Structure on Capacitance of the Ruo2 Electrode. J. Power Sources 2009, 189, 1284-1291. 23. Sharma, R.; Rastogi, A. Desu, S. Pulse Polymerized Polypyrrole Electrodes for High Energy Density Electrochemical Supercapacitor. Electrochem. Commun. 2008, 10, 268-272. 24. Anjos, D. M.; McDonough, J. K.; Perre, E.; Brown, G. M.; Overbury, S. H.; Gogotsi, Y. Presser, V. Pseudocapacitance and Performance Stability of Quinone-Coated Carbon Onions. Nano Energy 2013, 2, 702-712. 25. Hu, M. M.; Li, Z. J.; Hu, T.; Zhu, S. H.; Zhang, C. Wang, X. H. High-Capacitance Mechanism for Ti3C2Tx Mxene by in Situ Electrochemical Raman Spectroscopy Investigation. ACS Nano 2016, 10, 1134411350. 26. Halim, J.; Kota, S.; Lukatskaya, M. R.; Naguib, M.; Zhao, M. Q.; Moon, E. J.; Pitock, J.; Nanda, J.; May, S. J.; Gogotsi, Y., et al. Synthesis and Characterization of 2D Molybdenum Carbide (MXene). Adv. Funct. Mater. 2016, 26, 3118-3127. 27. Anasori, B.; Xie, Y.; Beidaghi, M.; Lu, J.; Hosler, B. C.; Hultman, L.; Kent, P. R. C.; Gogotsi, Y. Barsoum, M. W. Two-Dimensional, Ordered, Double Transition Metals Carbides (MXenes). ACS Nano 2015, 9, 9507-9516. 28. Zhan, C.; Naguib, M.; Lukatskaya, M.; Kent, P. R. C.; Gogotsi, Y. Jiang, D. E. Understanding the MXene Pseudocapacitance. J Phys. Chem. Lett. 2018, 9, 1223-1228. 29. Khazaei, M.; Arai, M.; Sasaki, T.; Chung, C. Y.; Venkataramanan, N. S.; Estili, M.; Sakka, Y. Kawazoe, Y. Novel Electronic and Magnetic Properties of Two-Dimensional Transition Metal Carbides and Nitrides. Adv. Funct. Mater. 2013, 23, 2185-2192. 30. Liu, Y. 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. 31. Greeley, J.; Jaramillo, T. F.; Bonde, J.; Chorkendorff, I. B. Norskov, J. K. Computational HighThroughput Screening of Electrocatalytic Materials for Hydrogen Evolution. Nat. Mater. 2006, 5, 909-913. 32. Urbankowski, P.; Anasori, B.; Makaryan, T.; Er, D.; Kota, S.; Walsh, P. L.; Zhao, M.; Shenoy, V. B.; Barsoum, M. W. Gogotsi, Y. Synthesis of Two-Dimensional Titanium Nitride Ti4N3 (MXene). Nanoscale 2016, 8, 11385-11391. 33. Soundiraraju, B. George, B. K. Two-Dimensional Titanium Nitride (Ti2N) Mxene: Synthesis, Characterization, and Potential Application as Surface-Enhanced Raman Scattering Substrate. ACS Nano 2017, 11, 8892-8900. 34. Morel, A.; Borjon-Piron, Y.; Porto, R. L.; Brousse, T. Bélanger, D. Suitable Conditions for the Use of Vanadium Nitride as an Electrode for Electrochemical Capacitor. J. Electrochem. Soc. 2016, 163, A1077-A1082. 35. Xiao, X.; Yu, H.; Jin, H.; Wu, M.; Fang, Y.; Sun, J.; Hu, Z.; Li, T.; Wu, J. Huang, L. Salt-Templated Synthesis of 2D Metallic MoN and Other Nitrides. ACS Nano 2017, 11, 2180-2186. 36. Wen, Y.; Rufford, T. E.; Chen, X.; Li, N.; Lyu, M.; Dai, L. Wang, L. Nitrogen-Doped Ti3C2Tx Mxene Electrodes for High-Performance Supercapacitors. Nano Energy 2017, 38, 368-376. 37. Pandey, M. Thygesen, K. S. Two-Dimensional MXenes as Catalysts for Electrochemical Hydrogen Evolution: A Computational Screening Study. J. Phys. Chem. C 2017, 121, 13593-13598.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 15

38. Shan, Q.; Mu, X.; Alhabeb, M.; Shuck, C. E.; Pang, D.; Zhao, X.; Chu, X.-F.; Wei, Y.; Du, F. Chen, G. Two-Dimensional Vanadium Carbide (V2C) Mxene as Electrode for Supercapacitors with Aqueous Electrolytes. Electrochem. Commun. 2018, 96, 103-107. 39. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I., et al. Quantum Espresso: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter 2009, 21, 395502. 40. Perdew, J. P.; Burke, K. Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 41. Garrity, K. F.; Bennett, J. W.; Rabe, K. M. Vanderbilt, D. Pseudopotentials for High-Throughput Dft Calculations. Comput. Mater. Sci. 2014, 81, 446-452. 42. Sundararaman, R.; Letchworth-Weaver, K.; Schwarz, K. A.; Gunceler, D.; Ozhabes, Y. Arias, T. A. Jdftx: Software for Joint Density-Functional Theory. SoftwareX 2017, 6, 278-284. 43. Letchworth-Weaver, K. Arias, T. A. Joint Density Functional Theory of the Electrode-Electrolyte Interface: Application to Fixed Electrode Potentials, Interfacial Capacitances, and Potentials of Zero Charge. Phys. Rev. B 2012, 86, 075140. 44. Sundararaman, R. Goddard, W. A. The Charge-Asymmetric Nonlocally Determined Local-Electric (Candle) Solvation Model. J. Chem. Phys. 2015, 142, 064107.

ACS Paragon Plus Environment

14

Page 15 of 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC graphic:

ACS Paragon Plus Environment

15