Intrinsic Structural, Electrical, Thermal, and Mechanical Properties of

Subscriber access provided by University of Sussex Library

Article

Intrinsic Structural, Electrical, Thermal and Mechanical Properties of the Promising Conductor MoC MXene 2

Xianhu Zha, Jingshuo Yin, Yuhong Zhou, Qing Huang, Kan Luo, Jiajian Lang, Joseph S. Francisco, Jian He, and Shiyu Du J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b04192 • Publication Date (Web): 30 Jun 2016 Downloaded from http://pubs.acs.org on July 1, 2016

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 free 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 accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

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

Page 1 of 26

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

Intrinsic Structural, Electrical, Thermal and Mechanical Properties of the Promising Conductor Mo2C MXene Xian-Hu Zha1, Jingshuo Yin1, Yuhong Zhou1, Qing Huang1, Kan Luo1, Jiajian Lang1, Joseph S. Francisco2, Jian He3, and Shiyu DU1,*

1

Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo

Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo, Zhejiang, 315201, China 2

Departments of Chemistry and Earth and Atmospheric Science, Purdue University,

West Lafayette, IN 47906, USA 3

Center for Translational Medicine, Department of Biotechnology, Dalian Institute of

Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning, 116023, China * Corresponding Authors: [email protected]

1

ACS Paragon Plus Environment


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 26

ABSTRACT Mo2C, the newly synthesized MXene with a large lateral size and superconductivity property, has attracted increasing interest in material science. Employing first-principles density functional calculations, its intrinsic structural, electrical, thermal and mechanical properties are investigated in this work. It is found that this MXene is nonmagnetic with a small molar volume. The electrical conductivity is predicted in the order of 106 Ω−1m−1 , and its value is significantly influenced by doping. For thermal conductivity, both of the electron and phonon contributions are studied. At room temperature, the Mo2C’s thermal conductivity is determined to be 48.4 Wm−1K −1 , which can be further enhanced by increasing temperature and introducing n-type dopants. The specific heat and thermal expansion coefficient are also assessed, their values at room temperature are calculated as 290

Jkg−1K−1 and 2.26 × 10-6 K−1 , respectively. Moreover, the thermal contraction of the MXene is found at low temperatures. Under biaxial strains, the elastic modulus is predicted as 312 ± 10 GPa , and the ideal strength is determined to be 20.8 GPa at a critical strain of 0.086. In view of the small molar volume, super-high electrical conductivity, favorable thermal conductivity, low thermal expansion coefficient and high mechanical strength, the Mo2C MXene generally merits more widespread applications besides superconductors, such as applying to substrates for other layer materials, and candidate materials for batteries and supercapacitors.

2


Page 3 of 26

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


INTRODUCTION

MXenes, a new family of two dimensional materials, have attracted extensive attentions in recent years.1, 2 Since the discovery of Ti3C2Tx in 2011 (Tx referring to surface terminations, which generally are fluorine, hydroxyl and oxygen functional groups),3 over ten MXene members have been successful synthesized: Ti2CTx, Ta4C3Tx, (Ti0.5Nb0.5)2CTx, (V0.5Cr0.5)3C2Tx, Ti3CNTx,4 Nb2CTx, V2CTx,5 Nb4C3Tx,6 Mo2TiC2Tx, Mo2Ti2C3Tx, Cr2TiC2Tx,7 Zr3C2Tx,8 pristine Mo2C,9 Mo2CTx,10,11 Ti4N3Tx.12 Most existing MXenes are prepared by etching corresponding MAX phases and relatives13 with various etchants, such as concentrated hydrofluoric acid,3 ammonium bifluoride (NH4HF2),14 and the combination of common hydrochloric acid and lithium fluoride.15 By means of the etching approaches, various surface groups Tx generally form, and their compositions are related to the choice of etchants and the experimental conditions.16, 17 Since the intrinsic physical properties of the MXenes are found to be significantly dependent on their surface groups,18, 19 variations of surface terminations become an obstacle that hinder MXenes’ practical applications. Recently, Xu et al. successful synthesized a stable and clean Mo2C MXene.9 In contrast to previous approaches, the Mo2C MXene was prepared though the traditional chemical vapor deposition (CVD) method, in which methane was adopted as the carbon source and a copper foil was chosen as the substrate for a molybdenum foil. Using this CVD approach, the lateral size of the synthesized Mo2C MXene was measured to be larger than 100 µ m , and its surface was found to be clean without any functional groups. Additionally, the MXene was determined to present superconductivity under 2.85 K . 3



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 26

Considering the large lateral size, simple and stable configuration, and superconductivity characteristics, the Mo2C MXene has potential applications in the future. For example, a recent theoretical work implied that the Mo2C monolayer is a competitive anode material.20 However, up to date, most intrinsic physical properties of the Mo2C MXene are still unavailable. It is difficult to design functional or structural materials using Mo2C without the important physical properties being obtained. Therefore, to better understand and to utilize this new material, the intrinsic structural, electrical, thermal and mechanical properties of the Mo2C MXene are studied in this work. Moreover, this work provides some critical data for the MXene, which shows promise for a number of practical applications.

COMPUTATIONAL DETAILS First-principles calculations are implemented in the plane-wave VASP code.21 The generalized gradient approximation (GGA) of the Purdue-Burke-Ernzerhof (PBE) approach22 is adopted to describe the exchange-correlation functional on the basis of the projected augmented wave (PAW). Plane-waves with energies up to 500 eV are employed to describe the electronic wavefunctions, in which the C 2s2, 2p2 and Mo 4p6, 5s1, 4d5 electrons are considered as valance states. During optimization, all the structures are relaxed until the forces on each atom are smaller than 1.0 × 10-3 eV/Å, and the criterion for energy convergence is chosen as 1.0 × 10-6 eV/cell. A Γ -centered sampling of 12 × 12 × 1 grid is used in describing the Brillouin zone (BZ). Spin-polarization calculations are performed. To eliminate the interactions of

4


Page 5 of 26

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


neighboring layers, a 30 Å lattice parameter in the c-axis perpendicular to the MXene surface is adopted. To obtain the electrical conductivity and the electron thermal conductivity, we perform the semi-classical Boltzmann transport calculations within the constant relaxation time approximation ( τ = 5.52 ×10−15 s ). The constant relaxation time τ for the Mo2C MXene was determined according to Xu’s work.9 The transport calculations are undertaken by the BoltzTraP code.23 The electronic structure for transport properties are calculated using a very fine k-mesh of 42 × 42 × 1. In the BoltzTraP calculation, the original k-mesh is interpolated onto a mesh twenty times as dense, and the thetrahedron scheme is adopted to obtain the density of states. In regard to the vibrational properties, the phonon spectrum is calculated based on the density functional perturbation theory,24 which is implemented by the combination of the VASP and the Phonopy softwares.25 A 6 × 6 × 1 k-points mesh based on a 4 × 4 × 1 supercell is employed for calculating the dynamical matrix, and a 120 k-points grid is adopted for plotting the phonon dispersion for various directions and the entire BZ. The phonon thermal conductivity is calculated in the framework of the Klemens theory:26-28 κp =

< υ j >4

ρ

∑ ω max, j

ln

here ρ is mass density, which equals to ρ = M / (

ω max, j ω min, j

(1)

3 2 a d) , with M being the mass 2

of the unit cell,27 and a and d denoting the hexagonal lattice constant in the xy-plane and the monolayer thickness in the z-direction, respectively. In order to 5



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 26

obtain the layer thickness, the multilayer Mo2C MXene is optimized, and d is chosen as the monomer thickness, similar to the definition for monolayer MoS2.29 Since the van der Waals correction is crucial in describing layer interactions,30 a damped VDW correction (DFT-D3)31 is adopted in the optimization of multilayer configuration. T denotes temperature; υ j and γ j are the group velocity28 and Grüneisen parameter32 of the j − branch in phonon spectrum; ωmin, j and ωmax, j are the minimum and maximum circular phonon frequency of each branch. Due to the flake length L , the minimum frequency is redefined as ω min, j = (

M ω max, j < υ j > 3 2k B TL < γ j2 >

1

)2

with k B representing the Boltzmann’s constant. The thermal expansion coefficient is investigated based on the Grüneisen approximation:32-34

α=

1 Nk

1 ∂ 2 Es a 02 | ∂a 2 0

∑c

v

( j, k)γ ( j, k)

(2)

j,k

here Nk is the k-point number adopted in plotting the phonon spectrum, which equals to 120 in our calculations; Es is the strain energy; cv ( j, k) is the ( j, k) mode

contribution

c v ( j, k) = k B (

hω j,k k BT

Specific heat33,

c=

to

exp( )

hω j,k

2

[exp( 35, 36

the k BT

hω j,k kBT

heat

capacity,

whose

expression

is

)

, with h denoting the reduced Planck’s constant.

) − 1]2

is proportional to the heat capacity, which is calculated as

1 ∑cv (j,k) . Nk M j,k To predict the elastic modulus and the ideal strength under biaxial strains, the

calculations on stretching strains up to 12.5% with an increment of 0.5% are conducted. To be more precise, a smaller increment of 0.1% in the vicinity of the 6


Page 7 of 26

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


critical strain is undertaken. Due to the vacuum layer, the elastic constants of the MXene is rescaled by h / d , with h being the lattice parameter in c-axis (30 Å in this work). All the structures and electronic density distributions are visualized in the VESTA code.37 All the computational parameters and methods have been well tested. For example, the methods and relevant parameters for calculating the thermal properties have been adopted in our previous works.28, 33 In our benchmark calculations, the predicted thermal conductivities of graphene and ZnSb (4756 and 3.099 Wm −1K −1 at room temperature, respectively) by the current methods have been shown to be in excellent agreement with the measurements.38, 39 Moreover, to verify the semi-classic Boltzmann calculations implemented in the BoltzTraP code, we have calculated the transport coefficient

σ xx

τ , equaled to the electrical conductivity divided by the

scattering time, of the monolayer Sc2C(OH)2. The value at room temperature is predicted to be 1.98 × 1020 Ω −1m −1s −1 at a chemical potential of -0.5 eV , which is consistent with that of the Khazaei’s work.40

RESULTS AND DISCUSSION In order to investigate intrinsic physical properties, the structure and the fundamental electronic properties of the Mo2C MXene are investigated firstly. As shown in Figure 1, the top- and side-views of the MXene are presented in Figure 1(a) and 1(b), respectively. A central carbon layer is sandwiched by two molybdenum monolayers, and a molybdenum atom is on the top-site of the center of three −

neighboring carbon atoms. The space group is determined to be P 3 m1 (No. 164), and its corresponding BZ is given in Figure 1(c). As seen from the figure, the high 7



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 26

symmetry routes ΓΜ and ΓΚ in BZ are correspond to the armchair and zigzag directions in the real-space.28 For the Mo2C MXene, the hexagonal lattice parameter in the xy-plane is optimized to be 3.00 Å. Through the optimization of the multilayer configuration, the monolayer thickness is determined to be 4.64 Å, which is approximately twice that of the bare layer thickness (2.34 Å).41 Evidently, both of the lattice parameter and the layer thickness are smaller than those of functionalized Mo2CT2 (T=O, F, OH) MXenes.19 Based on the small lattice parameter and thin layer thickness, the Mo2C MXene generally possesses a small molar volume. Meanwhile, these structural parameters imply that the molybdenum atoms bond strongly with carbon atoms. To test the bond strength and stability of the Mo2C MXene, the binding energy is also assessed ( E binding = (2E Mo + E C − E Mo2 C ) / 3 , where EMo and EC are the energies of separated molybdenum and carbon atoms, respectively, and E Mo2 C denotes the total energy of Mo2C unit cell), which is determined to be 6.46 eV. This large and positive value partially implicates the stability of the Mo2C MXene. Based on the optimized structure, the electronic density of states (DOS) is studied and presented in Figure 1(d). It shows that the Mo2C MXene is a nonmagnetic and metallic system with a large DOS around the Fermi level. This is markedly different from its fluorine functionalized Mo2CF2 MXene with a small band gap.40 By further analization with the projected DOS, the DOS in the vicinity of the Fermi level is determined to be mainly contributed by the electrons of the molybdenum atoms, especially those valance electrons in the 4d orbitals. The large DOS implies that there exist a high number of migration channels for electrons. For visualization, the electronic wavefunctions with energy in the range of -0.2~0.2 eV which significantly influenced the electrical conductivity are plotted and provided as an inset in Figure 1(d). Wave packets can be seen around the molybdenum atoms, which are mainly related to the 8


Page 9 of 26

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


high degeneracy of the 4d orbitals. The high DOS near Fermi level is a reflection of the superconductivity of the Mo2C MXene.9 It is worth mentioning that we have also considered the influence of spin-orbit coupling (SOC) in the electronic energy bands. It shows that the energy bands calculated with SOC corrections change little from those without them. Therefore, the contribution of SOC may not be significant for the Mo2C MXene system and is not considered in the following sections.

Figure 1. (a) and (b) are the top-view and side-view of the Mo2C MXene. (c) shows the Brillouin zone (BZ) of the 2D hexagonal lattice, in which the high symmetry routes ΓΜ and ΓΚ correspond to the real-space armchair and zigzag directions, respectively. (d) presents the electronic density of states (DOS) of the MXene, in which the projected DOS of each atomic orbital is provided. The inset presents the electronic wavefunctions with energy in the range of -0.2~0.2 eV with the Fermi level is adjusted to zero eV. 9



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 26

Based on the optimized structure and electronic characteristics, the electrical conductivity and the electron thermal conductivity of the Mo2C MXene are further investigated. It is revealed that the two quantities are isotropic in the plane of the monolayer regardless of the temperatures and chemical potentials. To facilitate elaboration, the transport coefficients defined by the electrical conductivity and electron thermal conductivity divided by the scattering time (

σ xx

κ xx ), as τ and τ

functions of chemical potentials at various temperatures, are provided in Figure 2. As shown in Figure 2(a),

σ xx

τ

is generally independent of the increasing temperatures

without considering the temperature effect on the chemical potential, which is similar to the behaviors in other metallic MXenes.40 Conversely, the chemical potential plays an important role in regulating the electrical conductivity. At zero chemical potential, the

σ xx

τ

value is determined to be 6.51 × 1020 Ω −1m −1s −1 (using τ = 5.52 ×10−15 s ,

the σ xx value is calculated as 3.59 × 106 Ω−1m−1 ), which is substantially consistent with the experiment result ( σ is calculated to be 8.99 × 106 Ω−1m−1 based on a layer thickness of 8.3 nm at 10 K ).9 With a positive chemical potential up to 0.684 eV ,

σ xx

τ

increases with increasing chemical potential, and its maximum value at 0.684

eV is determined to be 10.6 × 1020 Ω −1m −1s −1 (corresponding σ xx is 5.87 × 106 Ω−1m−1 ). Since the positive chemical potential is mainly caused by n-doping, it’s a demonstration that the n-doping is beneficial to increase the MXene’s electrical conductivity. Although a higher chemical potential corresponds to a lower electrical conductivity, it has been beyond the applicability of the current model. In contrast to the circumstance with a positive chemical potential, the negative chemical potential 10


Page 11 of 26

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


can be induced by introducing p-type dopants and increasing temperature.42 The

σ xx

τ

value firstly decreases with the increasing absolute value of the negative

chemical potential, and the minimum value is found to be 5.00 × 1020 Ω −1m −1s −1 at -0.24 eV. Thereafter, the

σ xx

τ value increases. It is interesting to point out that the

lowered Fermi energy by p-doping causes the increase in density of states according to the DOS plot shown in Figure 1(d), however, it can be seen the value of

σ xx

τ

decreases when the chemical potential is slightly lowered. This indicates the electron mobility is lowered when low amount p-doping occurs. In general, the predicted electrical conductivity is much higher than those of Mo2CT2MXenes.40 In regard of

κ xx with different τ values as shown in Figure 2 (b), the function of τ σ κ chemical potentials is similar to that for xx . A remarkable difference is that xx τ τ

the

κ xx

is proportional to the increasing temperature. As a consequence, the metallic Mo2C MXene satisfies the Wiedemann-Franz Law.43 At zero chemical potential, the

κ xx

τ

value is determined to be 47.8 × 1014 Wm −1K −1s −1 at room temperature ( κxx is calculated as 26.4 Wm −1K −1 ). The corresponding value at 500 K is 80.9 × 1014 Wm −1K −1s −1 . At the higher chemical potential of 0.684 eV, the room temperature

κ xx

14 14 −1 −1 −1 τ increases to 77.4 × 10 Wm K s , and the value at 500 K reaches 127 × 10

Wm −1K −1s −1

(corresponding κxx are 42.7 and 70.2 Wm −1K −1 , respectively).

Interestingly, the electron thermal conductivity can be further enhanced by increasing the temperature. Evidently, the Mo2C MXene possesses super-high electrical conductivity and favorable thermal conductivity.

11



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

Figure 2. The transport coefficients (a)

σ xx

Page 12 of 26

κ xx , as functions of the τ and (b) τ

chemical potentials at various temperatures.

To describe the thermal conductivity more precisely, the phonon thermal conductivity of the Mo2C MXene is also investigated. Figure 3(a) and 3(b) show the phonon dispersions along the zigzag ( ΓΜ ) and armchair ( ΓΚ ) directions, respectively. Three acoustic modes [out-of-plane acoustic (ZA), longitudinal acoustic (LA) and transversal acoustic (TA) modes] are denoted in black rectangles, red circles and blue triangles, respectively. The ZA mode shows the lowest phonon frequency. Meanwhile, the phonon dispersions of the TA and LA modes are nearly superposed, especially along the armchair direction. It is worth mentioning that the LA mode is slightly softened at the Μ high-symmetry point along the zigzag direction. Based on the phonon dispersions and a flake length of 5 µ m , the phonon thermal conductivities along the zigzag and armchair directions are calculated and provided in Figure 3(c) 12


Page 13 of 26

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


and 3(d), respectively. From these two graphs, we find that the phonon thermal conductivity is mainly contributed by the LA mode, and its value decreases with increasing temperature in the range of 100-500 K investigated. At room temperature, the phonon thermal conductivity is determined to be 9.72 Wm −1K −1 in the zigzag direction, in which the ZA, TA and LA mode’s contributions are 1.50, 1.93 and 6.29 Wm −1K −1 , respectively. Correspondingly, the phonon thermal conductivity in the

armchair direction is calculated as 16.2 Wm −1K −1 . The lower thermal conductivity in the zigzag direction can be ascribed to the softening of the LA mode at the M point, which introduces a smaller group velocity and a higher Grüneisen parameter. In general, a long flake of a 2D material may reduce boundary scattering of phonons, leading to high thermal conductivity. With a larger flake length of 100 µ m approaching to the experimental result,9 the phonon thermal conductivities increase to 12.9 and 22.0 Wm −1K −1 in the zigzag and armchair directions, respectively. Considering both of the electron and phonon contributions, the thermal conductivity at room temperature is predicted as 48.4 Wm −1K −1 in the armchair direction, which increases to 64.7 Wm −1K −1 by means of a heavy n-doping (with chemical potential equal to 0.68 eV ). Moreover, a higher temperature will further enhance the thermal conductivity. For example, the value is determined to be 92.2 Wm −1K −1 at 500 K in view of the n-doping, which is proximately close to that of pure iron.44 As to the thermal conductivity in the zigzag direction, the value is approximately 9 Wm −1K −1 lower at room temperature.

13



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 26

Figure 3. (a) and (b) are the phonon dispersions in the zigzag and armchair directions, respectively. (c) and (d) show the phonon thermal conductivities along these two directions.

Since the Mo2C MXene shows excellent electrical and favorable thermal conductivities, it is capable as a conductive material at various temperatures. Herein, the thermal expansion coefficient and the specific heat become two crucial factors in determining its practical applications. Based on the phonon dispersion in the entire BZ as shown in Figure 4(a), these two quantities of the Mo2C MXene are also investigated. The specific heat and thermal expansion coefficient as functions of the temperatures are given in Figure 4(b) and 4(c), respectively. As seen from the figures, the specific heat increases with increasing temperature in the entire range of 14


Page 15 of 26

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


temperatures investigated, and its room temperature value is determined to be 290

Jkg−1K−1 . At 800 K , the highest temperature investigated, the specific heat reaches 351 Jkg−1K−1 . With respect to the thermal expansion behavior, the Mo2C MXene is distinct from those of functionalized MXene and many bulk systems.45-47 This MXene contracts with the increasing temperature when temperature is low, and it expands afterwards. The negative thermal expansion coefficient with the largest absolute value is determined to be -3.27 × 10-6 K−1 at 25 K , and the maximum thermal contraction is found at 125 K where the thermal expansion coefficient reaches zero. After that, the value of the thermal expansion coefficient is positive and increases with increasing temperature. At 800 K , the value is determined to be 3.65 × 10-6 K−1 . Evidently, the magnitude of the thermal expansion coefficient is much smaller than those of most metals, such as copper and nickel.46 In order to gain a deeper understanding of the thermal expansion behavior, the Grüneisen parameter γ

32, 34

shown each mode contribution in the entire BZ are investigated and presented in Figure 4(d). This figure shows that the absolute values of γ for acoustic modes are much larger than those of optical modes, which implies that the acoustic modes undertake the main contributions. The ZA mode shows large negative γ near the BZ center, and positive values at the BZ boundary. Thus, this mode both contributes to the thermal contraction and expansion of the Mo2C MXene. In contrast to the ZA mode, the TA mode mainly contributes to the thermal contraction, and the LA mode generally causes the thermal expansion. Based on the predicted low thermal expansion coefficient of the Mo2C MXene, a large thermal mismatch is expected to 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

Page 16 of 26

appear between the Mo2C MXene and the pure molybdenum,47 which probably is the main reason that we can obtain the clean Mo2C MXene free from Mo nanoparticles by a rapid cooling approach.9 Moreover, the relatively low thermal expansion coefficient can be beneficial in the Mo2C’s practical applications because the structure is robust to temperature variations.

Figure 4. (a) The phonon dispersion in the entire Bollouin zone (BZ); (b) the specific heat as a function of temperature; (c) the thermal expansion coefficient as a function of temperature; (d) the Grüneisen parameter in the entire BZ.

Since biaxial strains can inevitably emerge due to existence of substrates when using layer materials, it is interesting to investigate the mechanical properties of the Mo2C MXene under the biaxial strains. Figure 5 shows the relationship between the stress and the biaxial strain. By fitting the stress strain function based on linear 16


Page 17 of 26

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


regression up to 1.50% (as shown in the inset), we obtained a biaxial elastic modulus

Y2D = 312 ± 10GPa , which is slightly higher than that of monolayer MoS2.29,

48

Incidentally, the large elastic modulus can be attributed to the strong interactions between the molybdenum and carbon atoms. With increasing strains, the stress increases up to the critical strain of 0.086, and then the MXene shows a creep deformation. The ideal strength at the critical strain is determined to be 20.8 GPa . Interestingly, the predicted ideal strength also approaches to that of the monolayer MoS2 (23.8 GPa ), though the critical strain is smaller than that of MoS2(0.195).29 These data imply that the Mo2C MXene is a virtually strong elastic 2D material which is robust to applied strains.

Figure 5. The calculated stress versus biaxial strain for the Mo2C MXene. The inset shows the linear regression of the initial stress strain curve for calculating the elastic modulus. 17



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 18 of 26

From this work, the structure, electrical, thermal and mechanical properties of the Mo2C are comprehensively studied. It is found that the MXene presents a small molar volume, super-high electrical conductivity, favorable thermal conductivity, low thermal expansion coefficient and strong mechanical strength. These properties suggest that this MXene may have broader applications. For example, the Mo2C MXene can be used as conductive material based on its super-high electrical conductivity and favorable thermal conductivity. Moreover, it may be used as a substrate for other layer systems49 due to its strong conductivity and the structure robustness to temperature variation and strains. Considering the strong capacity on adsorption of ions,50 small molar volume and super-high conductivity, the Mo2C MXene can also act as a good candidate material for batteries and supercapacitors.

CONCLUSIONS In summary, we have investigated various intrinsic physical properties of the promising Mo2C MXene. The MXene is found to present a small molar volume, super-high electrical conductivity and favorable thermal conductivity. Moreover, it is robust to temperature variations and strains. Explicitly, the lattice parameter in the xy-plane is determined to be 3.00 Å, and the layer thickness is predicted as 4.64Å. The electrical conductivity is of order of 106 Ω−1m−1 . The thermal conductivity increases with the increasing temperature; the intrinsic value at room temperature is predicted to be 48.4 Wm −1K −1 , which increases to 64.7 Wm −1K −1 by introducing a heavy n-doping. The maximum thermal expansion coefficient is determined below 3.65 × 10-6 K−1 in the entire range of temperatures up to 800 K . The biaxial elastic modulus is predicted as 312 ± 10 GPa by a linear fitting, and the ideal strength is 18


Page 19 of 26

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


found to be 20.8 GPa at the critical strain of 0.086. Based on these predicted properties, the Mo2C MXene merits promise for a number of potential applications.

ACKNOWLEDGEMENTS The authors acknowledge the support of the Division of Functional Materials and Nanodevices, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, the National Natural Science of Foundations of China (Grant Nos. 51372046, 51479037 and 91226202), the Ningbo Municipal Natural Science Foundation (Nos. 2014A610006 and 2016A610272), ITaP at Purdue University for computing resources and the key technology of nuclear energy, 2014, CAS Interdisciplinary Innovation Team.

19



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

Page 20 of 26

REFERENCES (1) Naguib, M.; Mochalin, V. N.; Barsoum, M. W.; Gogotsi, Y., 25th Anniversary Article: MXenes: A New Family of Two-Dimensional Materials. Adv. Mater. 2014, 26, 992-1005. (2) Lei, J.-C.; Zhang, X.; Zhou, Z., Recent Advances in MXene: Preparation, Properties, and Applications. Front. Phys-Beijing 2015, 10, 276-286. (3) 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. (4) Naguib, M.; Mashtalir, O.; Carle, J.; Presser, V.; Lu, J.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W., Two-Dimensional Transition Metal Carbides. ACS Nano 2012, 6, 1322-1331. (5) Naguib, M.; Halim, J.; Lu, J.; Cook, K. M.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W., New Two-Dimensional Niobium and Vanadium Carbides as Promising Materials for Li-Ion Batteries. J. Am. Chem. Soc. 2013, 135, 15966-15969. (6) Ghidiu, M.; Naguib, M.; Shi, C.; Mashtalir, O.; Pan, L. M.; Zhang, B.; Yang, J.; Gogotsi, Y.; Billinge, S. J. L.; Barsoum, M. W., Synthesis and Characterization of Two-Dimensional Nb4C3 (MXene). Chem. Commun. 2014, 50, 9517-9520. (7) Anasori, B. X., Y. Beidaghi, M. Lu, J. Hosler, B. C. Hultman, L. Kent, P. R. C. Gogotsi, Y. and Barsoum, M. W., Two-Dimensional, Ordered, Double Transition Metals Carbides (MXenes). ACS Nano 2015, 9, 9507. (8) Zhou, J.; Zha, X.-H.; Chen, F.; Ye, Q.; Eklund, P.; Du, S.; Huang, Q., A New Two-Dimensional Zirconium Carbide MXene by Selective Etching of Al3C3 from Nanolaminated Zr3Al3C5. Angew. Chem. Int. Edit. 2016, 55, 5008-5013. (9) Xu, C.; Wang, L. B.; Liu, Z. B.; Chen, L.; Guo, J. K.; Kang, N.; Ma, X. L.; Cheng, 20


Page 21 of 26

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


H. M.; Ren, W. C., Large-Area High-Quality 2D Ultrathin Mo2C Superconducting Crystals. Nat. Mater. 2015, 14, 1135. (10) Meshkian, R.; Näslund, L.-Å.; Halim, J.; Lu, J.; Barsoum, M. W.; Rosen, J., Synthesis of Two-Dimensional Molybdenum Carbide, Mo2C, from the Gallium Based Atomic Laminate Mo2Ga2C. Scripta Mater. 2015, 108, 147-150. (11) Halim, J.; Kota, S.; Lukatskaya, M. R.; Naguib, M.; Zhao, M.-Q.; Moon, E. J.; Pitock, J.; Nanda, J.; May, S. J.; Gogotsi, Y.; Barsoum, M. W., Synthesis and Characterization of 2D Molybdenum Carbide (MXene). Adv. Funct. Mater. 2016, 26, 3118-3127. (12) 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. (13) Barsoum, M. W.; Radovic, M., Elastic and Mechanical Properties of the MAX Phases. Annu. Rev. Mater. Res. 2011, 41, 195-227. (14) Halim, J., Lukatskaya, M. R.; Cook, K. M.; Lu, J.; Smith, C. R.; Näslund, L.-Å.; May, S. J.; Hultman, L.; Gogotsi, Y.; Eklund, P. et al., Transparent Conductive Two-Dimensional Titanium Carbide Epitaxial Thin Films. Chem. Mater. 2014, 26, 2374-2381. (15) Ghidiu, M.; Lukatskaya, M. R.; Zhao, M.-Q.; Gogotsi, Y.; Barsoum, M. W., Conductive Two-Dimensional Titanium Carbide 'Clay' with High Volumetric Capacitance. Nature 2014, 516, 78-81. (16) Hope, M. A.; Forse, A. C.; Griffith, K. J.; Lukatskaya, M. R.; Ghidiu, M.; Gogotsi, Y.; Grey, C. P., NMR Reveals the Surface Functionalisation of Ti3C2 MXene. Phys. Chem. Chem. Phys. 2016, 18, 5099-5102. (17) Xie, Y.; Naguib, M.; Mochalin, V. N.; Barsoum, M. W.; Gogotsi, Y.; Yu, X.; Nam, 21



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 22 of 26

K.-W.; Yang, X.-Q.; Kolesnikov, A. I.; Kent, P. R. C., Role of Surface Structure on Li-Ion Energy Storage Capacity of Two-Dimensional Transition-Metal Carbides. J. Am. Chem. Soc. 2014, 136, 6385-6394. (18) 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. (19) Zha, X.-H.; Luo, K.; Li, Q.; Huang, Q.; He, J.; Wen, X.; Du, S., Role of the Surface Effect on the Structural,Electronic and Mechanical Properties of the Carbide MXenes. EPL- Europhys. Lett. 2015, 111, 26007. (20) Cakir, D.; Sevik, C.; Gulseren, O.; Peeters, F. M., Mo2C as a High Capacity Anode Material: A First-Principles Study. J. Mater. Chem. A 2016, 4, 6029-6035. (21) 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-11186. (22) Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (23) Madsen, G. K. H.; Singh, D. J., Boltztrap. A Code for Calculating Band-Structure Dependent Quantities. Comput. Phys. Commun. 2006, 175, 67-71. (24) Gonze, X.; Lee, C., Dynamical Matrices, Born Effective Charges, Dielectric Permittivity Tensors, and Interatomic Force Constants from Density-Functional Perturbation Theory. Phys. Rev. B 1997, 55, 10355-10368. (25) Togo, A.; Oba, F.; Tanaka, I., First-Principles Calculations of the Ferroelastic Transition between Rutile-Type and CaCl2-Type SiO2 at High Pressures. Phys. Rev. B 2008, 78, 134106. (26) Klemens, P. G.; Pedraza, D. F., Thermal Conductivity of Graphite in the Basal 22


Page 23 of 26

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


Plane. Carbon 1994, 32, 735-741. (27) Klemens, P. G., Theory of Thermal Conduction in Thin Ceramic Films. Int. J. Thermophys. 2001, 22, 265-275. (28) Zha, X.-H.; Zhou, J.; Zhou, Y.; Huang, Q.; He, J.; Francisco, J. S.; Luo, K.; Du, S., Promising Electron Mobility and High Thermal Conductivity in Sc2CT2(T = F, OH) MXenes. Nanoscale 2016, 8, 6110-6117. (29) Li, T., Ideal Strength and Phonon Instability in Single-Layer MoS2. Phys. Rev. B 2012, 85, 235407. (30) Rydberg, H.; Dion, M.; Jacobson, N.; Schröder, E.; Hyldgaard, P.; Simak, S. I.; Langreth, D. C.; Lundqvist, B. I., Van der Waals Density Functional for Layered Structures. Phys. Rev. Lett. 2003, 91, 126402. (31) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H., A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (32) Mounet, N.; Marzari, N., First-Principles Determination of the Structural, Vibrational and Thermodynamic Properties of Diamond, Graphite, and Derivatives. Phys. Rev. B 2005, 71, 205214. (33) Zha, X.-H.; Zhang, R.-Q.; Lin, Z., Point Defect Weakened Thermal Contraction in Monolayer Graphene. J. Chem. Phys. 2014, 141, 064705. (34) Zha, X.-H.; Zhang, R.-Q.; Lin, Z., Tuning Thermal Expansions of Zinc Oxide Sheets by Varying the Layer Thickness. EPL-Europhys. Lett. 2014, 107, 26007. (35) Upadhyay Kahaly, M.; Waghmare, U. V., Size Dependence of Thermal Properties of Armchair Carbon Nanotubes: A First-Principles Study. Appl. Phys. Lett. 2007, 91, 023112. (36) Zha, X.-H.; Li, S.; Zhang, R.-Q.; Lin, Z., Remarkable Thermal Contraction in 23



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 24 of 26

Small Size Single-Walled Boron Nanotubes. Commun. Comput. Phys. 2014, 16, 201-212. (37) Momma, K.; Izumi, F., VESTA 3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272-1276. (38) Balandin, A. A.; Ghosh, S.; Bao, W. Z.; Calizo, I.; Teweldebrhan, D.; Miao, F.; Lau, C. N., Superior Thermal Conductivity of Single-Layer Graphene. Nano Lett. 2008, 8, 902-907. (39) Eklöf, D.; Fischer, A.; Wu, Y.; Scheidt, E. W.; Scherer, W.; Häussermann, U., Transport Properties of the II–V Semiconductor ZnSb. J. Mater. Chem. A 2013, 1, 1407-1414. (40) Khazaei, M.; Arai, M.; Sasaki, T.; Estili, M.; Sakka, Y., Two-Dimensional Molybdenum Carbides: Potential Thermoelectric Materials of the MXene Family. Phys. Chem. Chem. Phys. 2014, 16, 7841. (41) Kurtoglu, M.; Naguib, M.; Gogotsi, Y.; Barsoum, M. W., First Principles Study of Two-Dimensional Early Transition Metal Carbides. MRS Commun. 2012, 2, 133-137. (42) Zhang, Z.; Xie, Y.; Peng, Q.; Chen, Y., A Theoretical Prediction of Super High-Performance

Thermoelectric

Materials

Based

on

MoS2/WS2

Hybrid

Nanoribbons. Sci. Rep.-UK 2016, 6, 21639. (43) Kim, W.; Carbotte, J. P., Transport Properties in the d-Density Wave State: Wiedemann–Franz Law. Phys. Rev. B 2002, 66, 033104. (44) Lide D. R., CRC Handbook of Chemistry and Physics (84th Ed.), CRC Press, Boca Raton, Florida, 2003. (45) Zha, X.-H.; Huang, Q.; He, J.; He, H.; Zhai, J.; Francisco, J. S.; Du, S., The Thermal and Electrical Properties of the Promising Semiconductor MXene Hf2CO2. 24


Page 25 of 26

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


Sci. Rep.-UK 2016, 6, 27971. (46) Nix, F. C.; MacNair, D., The Thermal Expansion of Pure Metals: Copper, Gold, Aluminum, Nickel, and Iron. Phys. Rev. 1941, 60, 597-605. (47) Nix, F. C.; MacNair, D., The Thermal Expansion of Pure Metals. II: Molybdenum, Palladium, Silver, Tantalum, Tungsten, Platinum, and Lead. Phys. Rev. 1942, 61, 74-78. (48) Bertolazzi, S.; Brivio, J.; Kis, A., Stretching and Breaking of Ultrathin MoS2. ACS Nano 2011, 5, 9703-9709. (49) Gan, L.-Y.; Zhao, Y.-J.; Huang, D.; Schwingenschlögl, U., First-Principles Analysis

of

MoS2/Ti2C

and

MoS2/Ti2CY2

(Y=F

and

OH)

All-2D

Semiconductor/Metal Contacts. Phys. Rev. B 2013, 87, 245307. (50) 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 Two-Dimensional Titanium Carbide. Science 2013, 341, 1502-1505.

25



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

Table of contents (TOC) Image

26


Intrinsic Structural, Electrical, Thermal, and Mechanical Properties of

Recommend Documents