Li Xiao-Hong1, Su Xiang-Ying1, Zhang Rui-Zhou1. 1. College of Physics and Engineering, Henan University of Science and Technology Luoyang 471003 China lorna639@126.com.
Two-dimensional (2D) materials have been an active research area since the discovery of graphene.[1] Graphene is a significant 2D material[2] with zero-band gap semi-metallic properties, which hampers its application in highly integrated electronic components. We can tune the band gap of graphene, but the available range of the band gap is not enough to satisfy the requirement of the application of nanoelectronic, optoelectronic devices. Therefore, it is necessary to design or search for new layered materials with intrinsic bandgaps, high thermal stability and conductivity, and high carrier mobility.MXene is a new group of low-dimensional materials[3-5] and has attracted intensive attention in the field of supercapacitors,[6,7] lead adsorption,[8] catalysis[9] and so on. MXene is obtained by selectively etching layers of A elements from its host, layered ternary MAX materials (MAX), where n = 1, 2, 3; M is a transition metal, and X is C or N.[10] The frequently used etchant is hydrofluoric acid (HF), HF-containing etchants and so on. F, O, and/or OH functional groups are often generated on the external surface of the exfoliated layers. So functionalized MXenes are denoted as MXT, in which T stands for the terminating F, O, and/or OH functional groups, x is the termination number.[11] These functional groups are very important in the applications of energy storage, catalysis etc.[12-14]Defects have been investigated in 2D materials and can influence the performance of materials and devices.[15,16] It is important to understand the layered structure and point defects, which are crucial for exploration of the physiochemical properties of MXenes. The properties of 2D materials can be engineered by defects. For the low-dimensional materials, defects have a much stronger impact, when compared to their bulk counterparts. And defect engineering has attracted extensive attention[17] and is often used to modify the electronic structures of many layered structures.[18] Zhang et al.[19] theoretically investigated the carbon-vacancy ordering in Nb4AlC3−, and reported the importance roles of carbon vacancies in the structure stability and order–disorder phase transformation. Carbon vacancies also have important influences on the superconducting properties of materials.[20-22]Vacancies are universally present in the MAX phases.[23,24] So it is believed that MXenes inherit their hosts' defect. Ti2CO2 is the thinnest O-functionalized Ti-based MXene, and has many promising applications, such as catalyst, optical devices, and gas sensor.[25-27] Wang et al.[28] investigated the stabilities and electronic properties of vacancy-doped Ti2CO2. Their results indicated that the formation energies of C-vacancy are relatively small. Hu et al.[29] reported the influence of C-vacancy on the structural stability, electronic properties of a 2 × 2 × 1 Ti2CT2 (T = O, F, and OH) supercell. They thought that C-vacancy in MXene is much easier to form when compared with graphene and other two-dimensional material. Bandyopadhyay et al.[30] investigated the structural and magnetoelectronic behavior of a 4 × 4 × 1 Ti2CO2 supercell with single vacancy and double vacancies. They thought that the defect formation energies are greatly dependent on local chemical bonding. Sang et al.[31] reported the atomic defects in monolayer Ti3C2T MXene. They reported that the formation energy of Ti-vacancy on Ti3C2O (7.74 eV), and the defects can't strongly influence the metallic conductivity, but can influence the surface morphology and termination groups. To the best of our knowledge, limited researches report the impact of different vacancy defect on the properties of Ti2CO2 monolayer.Work function (WF) is an important parameter for electrode materials. A material with higher WF can be as an anode, and a material with low WF can be treated as a cathode. Research[32] indicates that the materials with higher work function can reduce the Schottky barrier, which can cause increased contact resistances and limit the performance of devices. Therefore, the modulation of WF is important to improve the device performance. Up to now, the report on the modulation of WF of Ti2CO2 monolayer is unavailable.Here, we investigated the effect of carbon vacancy (VC), Ti vacancy (VTi), and oxygen vacancy (VO) on structure, electronic properties, and work function of Ti2CO2 monolayer. The comparisons of related properties between pristine and vacancy Ti2CO2 monolayer and the modulation of vacancy defect on WF are performed in our study, which are not discussed in other people's works.[29-31] Our investigation will provide an alternative method to engineer the physical and chemical properties of MXenes.
Computational details
All density functional theory (DFT) calculations are carried out using the Dmol3 code.[33] The Perdew–Burke–Ernzerhof (PBE) functional[34] generalized gradient approximation (GGA)[35] is used. The used basis set is double numerical plus polarization (DNP), which corresponds to a double-ζ quality basis set with p-type polarization functions added to hydrogen, and d-type polarization functions added to heavier atoms. DNP basis set is comparable with the Gaussian 6-31G (d, p) basis set and exhibits a better accuracy.[36] A 3 × 3 × 1 Ti2CO2 supercell is constructed with one carbon, titanium, or oxygen removal to highlight the influence of VC, VTi, and VO, respectively. The three models of Ti2CO2-VC, Ti2CO2-VTi, and Ti2CO2-VO monolayers are presented in Fig. 1. In order to have a comparison, the structure of pristine Ti2CO2 monolayer (PTM) is also presented in Fig. 1.
Fig. 1
The atomic structures of pristine Ti2CO2 (a), Ti2CO2-VTi (b), Ti2CO2-VO (c), and Ti2CO2-VC (d) from top and side view.
In order to simulate the isolated monolayer, a large vacuum space of 25 Å was set in the direction normal to MXene layers to avoid any unnecessary interaction between the MXene layers. In our calculation, the Brillouin zone (BZ) integration was sampled by using a 5 × 5 × 1 k-mesh according to Monkhorst–Pack method.[37] The convergence in energy, force and displacement were 1 × 10−6 Ha, 0.001 Ha Å−1, and 0.005 Å, respectively. A denser 9 × 9 × 1 mesh is used to calculate the electronic density of states (DOS).
Results and discussion
The structure of pristine and defect-Ti2CO2 monolayer
In order to determine what defects are likely to form, the vacancy formation energies Eform for Ti, C, or O defect are obtained by the following equation:where Eform is the vacancy formation energy, Evac is the total energy of Ti2CO2 monolayer with mono-vacancy, Uatom is the total energy of an isolated atom, Epristine is the total energy of pristine Ti2CO2 monolayer (PTM). The calculated O-vacancy formation energy is 6.59 eV, lower than Ti-vacancy formation energy (9.67 eV), and higher than C-vacancy formation energy (5.31 eV). So C-vacancy is energetically more favorable than Ti-vacancy and O-vacancy, which is consistent with the conclusion of other researches.[29,30,38] Ti-vacancy and O-vacancy are not easy to occur when compared with C-vacancy, but they can be induced in non-equilibrium processes such as e-beam lithography,[39] which is used to decorate the vacancy defects in the MoS2 monolayer.[40] Ti-vacancy deformation is more energy consuming when compared with O-vacancy and C-vacancy, which is consistent with the other researchers' results.[29,31] Sang et al.[31] reported the formation energy of Ti-vacancy on Ti3C2O (7.74 eV), which indicates that Ti-vacancy is more difficult to occur on Ti2CO2 monolayer (9.67 eV) than Ti3C2O.The calculated crystal parameter of PTM is 9.0483 Å, while the crystal parameters of Ti2CO2 monolayers with Ti, C, or O defect are 9.0737, 9.0719, and 9.0636 Å, respectively. So the crystal parameters of Ti2CO2-VC, Ti2CO2-VTi, and Ti2CO2-VO monolayers change slightly when compared with that of PTM. The volume of PTM is 1665.59 Å3. Compared with the volume of PTM, the volume of Ti2CO2-VC monolayer increases 1.99%. The creation of C vacancy makes the Ti atoms near C vacancy away from the vacancy by the remaining strengthened pd bonds, which results in the volume increase of Ti2CO2-VC. Compared with the volume of PTM, the volume of Ti2CO2-VTi increases by 2.01%. The reason is that the creation of Ti vacancy makes the formation of unsaturated electrons on nearest C/O atoms, which result in the stronger interactions between C/O atoms and their nearest Ti atoms. The volume of Ti2CO2-VO increases by 1.90% when compared to that of PTM. The creation of O vacancy on the surface of Ti2CO2 monolayer makes the charge redistribution, and results in the formation of the stronger Ti–O and Ti–C bonds near O vacancy.
The electronic properties of pristine and defect-Ti2CO2 monolayers
Fig. 2 presents the structures of pristine and defect Ti2CO2 monolayer along the symmetry directions Γ–M–K–Γ. The valence band maximum (VBM) and conduction band minimum (CBM) are labelled in Fig. 2. PTM is an M–Γ indirect semiconductor with the band gap of 0.27 eV, which is close to the results of Xie et al. (0.24 eV by PBE)[41] and Zha et al. (0.261 eV by PBE).[42] The introduction of vacancy defect results in the transition from semiconductor to metal.
Fig. 2
The band structure of pristine and vacancy defect Ti2CO2 monolayer.
From Fig. 2(b), the introduction of Ti-vacancy makes the Fermi energy level shift downward compared with that of PTM, and the valence band passes through the Fermi level, which results in the metallic character of Ti2CO2-VTi. Ti defect level lies in the valence band region of PTM. From Fig. 2(c), the introduction of O vacancy defect makes the Fermi energy level shift upward. The defect energy level (band 1 in Fig. 2(c)) appears because of the removal of one oxygen atom, which passes through the Fermi energy level and results in the metallization of Ti2CO2-VO. This indicates that the O vacancy has greater effect on the photocatalyst properties of Ti2CO2 monolayer. In addition, it is noted that the defect energy level fluctuates slightly and changes gently, which indicates that the electrons in this energy level have larger effective mass, stronger localization and electron binding. From Fig. 2(d), the introduction of C vacancy makes the Fermi level shift upward compared with that of PTM. The defect energy level (band 1 in Fig. 2(d)) is introduced below the Fermi level. So the vacancy defect can effectively engineer the band structure of Ti2CO2 monolayer.Fig. 3 presents the total density of state (DOS) of pristine and vacancy defect Ti2CO2 monolayer. It is noted that the introduction of vacancy defect imposes the significant impact on the electronic structure of Ti2CO2 monolayer. The improved electronic conductivities of Ti2CO2-VTi, Ti2CO2-VO, and Ti2CO2-VC monolayers are observed due to the significantly increased DOS at the Fermi level (EF). So the introduction of vacancy defect can improve the electronic conductivities of functionalized MXenes.
Fig. 3
The total density of states of pristine and vacancy defect Ti2CO2 monolayer.
Fig. 4 presents the partial density of state (PDOS) of Ti, O, and C atoms of vacancy defect Ti2CO2 monolayer in order to illustrate the origin of electronic conductivity. The PDOS of PTM is also plotted in order to have a comparison. For PTM, the valence band is mainly from the contribution of Ti-d, C-p and O-p states, while the conduction band is mainly from the contribution of Ti-d states. Similar to PTM, Ti-d, O-p, and C-p states provide the mainly contribution of valence band for Ti2CO2-VTi monolayer, while C-p states mainly contribute to the conduction band. For Ti2CO2-VO and Ti2CO2-VC monolayers, Ti-d, O-p and C-p states provide the mainly contribution to valence band and conduction band. Combined Fig. 2 and 4, for Ti2CO2-VTi monolayer, the band at Γ point is mainly localized on C-p state, the introduction of Ti vacancy results in the unsaturated C and O atoms. For Ti2CO2-VO monolayer, the introduction of O vacancy results in three unsaturated Ti-3d states, which have the main contribution to the band at Γ point. The introduction of C vacancy results in the six unsaturated Ti atoms for Ti2CO2-VC monolayer, so Ti-d state have the main contribution to the band at Γ point and VBM from Fig. 4(d).
Fig. 4
The partial density of states of pristine and vacancy defect Ti2CO2 monolayer.
It is noted that the introduction of C and O vacancy makes Ti-d state split into several peaks in the energy range of 0 to 2 eV because of the structural modification. While the introduction of Ti vacancy makes Ti-d state weaker and C-p state stronger. Combined with Fig. 3 and 4, the valley of the total DOS at −3.5 eV in Fig. 3 is mainly from Ti-d state for Ti2CO2-VTi monolayer, Ti-d, O-p and C-p states for Ti2CO2-VO and Ti2CO2-VC monolayers. In the energy range of −3 to 0 eV, there exist strong interactions between Ti-d and C-p states, Ti-d and O-p states for all studied Ti2CO2 monolayers.
The analysis of work function and atomic charge
Work function (WF) is an important physical parameter of materials and plays the significant role in the light emitting diodes and the field effect transistors. The work function of pristine Ti2CO2 monolayer is 6.39 eV, while the work functions of Ti2CO2-VC, Ti2CO2-VTi, and Ti2CO2-VO monolayers are 7.27, 6.69, and 6.12 eV, respectively. This indicates that the work function can be effectively engineered by vacancy defect. According to the report of Ma et al.,[32] Ti2CO2-VC and Ti2CO2-VTi monolayers with higher WFs can reduce the Schottky barrier and decrease the field emission performance, while Ti2CO2-VO monolayer with lower WF can improve the field emission performance. In the other hand, Ti2CO2-VC and Ti2CO2-VTi monolayers can be considered as anode because of the higher WFs, while Ti2CO2-VO monolayer can be considered as cathode due to the lower WF (Fig. 5).
Fig. 5
The work function of pristine and vacancy defect Ti2CO2 monolayer.
The analysis of Mulliken charge is performed in order to further investigate the effect of vacancy defect on the charge distribution. Fig. 6 presents the electrostatic potential of Ti2CO2-VTi, Ti2CO2-VO, and Ti2CO2-VC monolayers. For PTM, the charge of Ti atom in supercell is 1.202e−, while the charges of the adjacent C and O atoms are −0.980 and −0.712e−, respectively. For Ti2CO2-VTi monolayer, the removal of Ti atom results in the charge localization on the defect site, adjacent unsaturated O and C atoms. The negative charge accumulation is exhibited around Ti vacancy defect. The charge aggregation is exhibited around O atoms, while the loss of negative charge is exhibited around C atoms. The charges of the adjacent C and O atoms are −0.704 and −0.584e−, respectively. For Ti2CO2-VO monolayer, the positive charge accumulation is exhibited around O vacancy defect. The charge of the adjacent Ti atom is 1.187e−, while the charge of the adjacent C atom has little change because of the longer distance from O vacancy. From Fig. 6(c), there is a clear charge accumulation around C vacancy defect for Ti2CO2-VC monolayer. The charge of the adjacent Ti atom is 1.189e−, while the charge of adjacent O atom increases and is −0.708 e−. The introduction of C vacancy weakens the adjacent Ti–O bonds.
Fig. 6
The electrostatic potentials of Ti2CO2-VTi, Ti2CO2-VO, and Ti2CO2-VC monolayers.
Conclusion
In this paper, the structure, electronic properties of pristine and vacancy defect Ti2CO2 monolayer have been investigated by the first principles calculation. The analysis of formation energies indicates that C-vacancy is energetically more favorable than Ti-vacancy and O-vacancy. PTM is a semiconductor. The introduction of vacancy defect results in the metallic character of Ti2CO2 monolayer and improves the electronic conductivity of Ti2CO2 monolayer. The work function can be effectively engineered by vacancy defect. O vacancy defect can improve the field emission performance, while Ti or C vacancy defect can decrease the field emission performance for Ti2CO2 monolayer. The analysis of Mulliken charge is further performed.
Authors: K S Novoselov; D Jiang; F Schedin; T J Booth; V V Khotkevich; S V Morozov; A K Geim Journal: Proc Natl Acad Sci U S A Date: 2005-07-18 Impact factor: 11.205
Authors: Romain Bourrellier; Sophie Meuret; Anna Tararan; Odile Stéphan; Mathieu Kociak; Luiz H G Tizei; Alberto Zobelli Journal: Nano Lett Date: 2016-06-16 Impact factor: 11.189
Authors: Yu Xie; Michael Naguib; Vadym N Mochalin; Michel W Barsoum; Yury Gogotsi; Xiqian Yu; Kyung-Wan Nam; Xiao-Qing Yang; Alexander I Kolesnikov; Paul R C Kent Journal: J Am Chem Soc Date: 2014-04-21 Impact factor: 15.419